54 minidx,maxidx,segid,icoff,ijcoff,inoff,ijnoff, &
55 ds,segm,cccofa,cccofb )
68 INTEGER :: segid,icoff,ijcoff,inoff,ijnoff
69 REAL(RFREAL) :: ds(minidx:maxidx)
70 REAL(RFREAL),
POINTER :: segm(:,:),cccofa(:,:),cccofb(:,:)
73 INTEGER ::
i,
j,
k, ijkc, ijkn, ijkn1
76 CHARACTER(CHRLEN) :: rcsidentstring
77 REAL(RFREAL) :: rwidth1, rwidth2, fac1, fac2, fac3
81 rcsidentstring =
'$RCSfile: TURB_floLesGenCoCCUtil.F90,v $'
84 'TURB_floLesGenCoCCUtil.F90' )
92 ijkn = indijk(
i ,
j ,
k ,inoff,ijnoff)
93 ijkn1 = indijk(
i ,
j ,
k+1 ,inoff,ijnoff)
94 ds(
i) = 0.5_rfreal*(segm(segid,ijkn)+segm(segid,ijkn1))
101 ijkc = indijk(
i ,
j ,
k ,icoff,ijcoff)
102 rwidth1= 1._rfreal/(ds(
i-1)+ds(
i))
103 rwidth2= 1._rfreal/(ds(
i+1)+ds(
i))
105 fac2 = ds(
i-1)*rwidth1+ds(
i+1)*rwidth2+2._rfreal
107 cccofa(1,ijkc) = 0.25_rfreal*fac1
108 cccofa(2,ijkc) = 0.25_rfreal*fac2
109 cccofa(3,ijkc) = 0.25_rfreal*fac3
111 rwidth1= 1._rfreal/(ds(
i-1)+2._rfreal*ds(
i)+ds(
i+1))
112 fac1 = (ds(
i-1)+ds(
i))*rwidth1
113 fac3 = (ds(
i+1)+ds(
i))*rwidth1
114 cccofb(1,ijkc) = 0.5_rfreal*fac1
115 cccofb(2,ijkc) = 0.5_rfreal
116 cccofb(3,ijkc) = 0.5_rfreal*fac3
121 ELSEIF (ijk==dirj)
THEN
125 ijkn = indijk(
i ,
j ,
k ,inoff,ijnoff)
126 ijkn1 = indijk(
i+1 ,
j ,
k ,inoff,ijnoff)
127 ds(
j) = 0.5_rfreal*(segm(segid,ijkn)+segm(segid,ijkn1))
134 ijkc = indijk(
i ,
j ,
k ,icoff,ijcoff)
135 rwidth1= 1._rfreal/(ds(
j-1)+ds(
j))
136 rwidth2= 1._rfreal/(ds(
j+1)+ds(
j))
138 fac2 = ds(
j-1)*rwidth1+ds(
j+1)*rwidth2+2._rfreal
140 cccofa(1,ijkc) = 0.25_rfreal*fac1
141 cccofa(2,ijkc) = 0.25_rfreal*fac2
142 cccofa(3,ijkc) = 0.25_rfreal*fac3
144 rwidth1= 1._rfreal/(ds(
j-1)+2._rfreal*ds(
j)+ds(
j+1))
145 fac1 = (ds(
j-1)+ds(
j))*rwidth1
146 fac3 = (ds(
j+1)+ds(
j))*rwidth1
147 cccofb(1,ijkc) = 0.5_rfreal*fac1
148 cccofb(2,ijkc) = 0.5_rfreal
149 cccofb(3,ijkc) = 0.5_rfreal*fac3
154 ELSEIF (ijk==dirk)
THEN
158 ijkn = indijk(
i ,
j ,
k ,inoff,ijnoff)
159 ijkn1 = indijk(
i ,
j+1 ,
k ,inoff,ijnoff)
160 ds(
k) = 0.5_rfreal*(segm(segid,ijkn)+segm(segid,ijkn1))
164 ds(kend+
k) = ds(kend)
167 ijkc = indijk(
i ,
j ,
k ,icoff,ijcoff)
168 rwidth1= 1._rfreal/(ds(
k-1)+ds(
k))
169 rwidth2= 1._rfreal/(ds(
k+1)+ds(
k))
171 fac2 = ds(
k-1)*rwidth1+ds(
k+1)*rwidth2+2._rfreal
173 cccofa(1,ijkc) = 0.25_rfreal*fac1
174 cccofa(2,ijkc) = 0.25_rfreal*fac2
175 cccofa(3,ijkc) = 0.25_rfreal*fac3
177 rwidth1= 1._rfreal/(ds(
k-1)+2._rfreal*ds(
k)+ds(
k+1))
178 fac1 = (ds(
k-1)+ds(
k))*rwidth1
179 fac3 = (ds(
k+1)+ds(
k))*rwidth1
180 cccofb(1,ijkc) = 0.5_rfreal*fac1
181 cccofb(2,ijkc) = 0.5_rfreal
182 cccofb(3,ijkc) = 0.5_rfreal*fac3
197 minidx,maxidx,segid,icoff,ijcoff,inoff,ijnoff, &
198 ds,segm,cccofa,cccofb )
206 #include "Indexing.h"
211 INTEGER :: segid,icoff,ijcoff,inoff,ijnoff
212 REAL(RFREAL) :: ds(minidx:maxidx)
213 REAL(RFREAL),
POINTER :: segm(:,:),cccofa(:,:),cccofb(:,:)
216 INTEGER ::
i,
j,
k, ijkc, ijkn, ijkn1
219 CHARACTER(CHRLEN) :: rcsidentstring
220 REAL(RFREAL) :: rwidth1, rwidth2, fac1, fac2, fac3, fac4, fac5
224 rcsidentstring =
'$RCSfile: TURB_floLesGenCoCCUtil.F90,v $'
227 'TURB_floLesGenCoCCUtil.F90' )
235 ijkn = indijk(
i ,
j ,
k ,inoff,ijnoff)
236 ijkn1 = indijk(
i ,
j ,
k+1 ,inoff,ijnoff)
237 ds(
i) = 0.5_rfreal*(segm(segid,ijkn)+segm(segid,ijkn1))
244 ijkc = indijk(
i ,
j ,
k ,icoff,ijcoff)
245 rwidth1= 1._rfreal/(ds(
i-1)+2._rfreal*ds(
i)+ds(
i+1))
246 fac1 = (ds(
i-1)+ds(
i))*rwidth1
247 fac3 = (ds(
i+1)+ds(
i))*rwidth1
248 cccofa(1,ijkc) = 0.5_rfreal*fac1
249 cccofa(2,ijkc) = 0.5_rfreal
250 cccofa(3,ijkc) = 0.5_rfreal*fac3
252 rwidth1= 1._rfreal/ &
253 (ds(
i-2)+ds(
i+2)+2._rfreal*(ds(
i-1)+ds(
i)+ds(
i+1)))
254 fac1 = (ds(
i-2)+ds(
i-1))*rwidth1
255 fac2 = (ds(
i-2)+2._rfreal*ds(
i-1)+ds(
i))*rwidth1
256 fac3 = (ds(
i-1)+2._rfreal*ds(
i)+ds(
i+1))*rwidth1
257 fac4 = (ds(
i)+2._rfreal*ds(
i+1)+ds(
i+2))*rwidth1
258 fac5 = (ds(
i+1)+ds(
i+2))*rwidth1
259 cccofb(1,ijkc) = 0.5_rfreal*fac1
260 cccofb(2,ijkc) = 0.5_rfreal*fac2
261 cccofb(3,ijkc) = 0.5_rfreal*fac3
262 cccofb(4,ijkc) = 0.5_rfreal*fac4
263 cccofb(5,ijkc) = 0.5_rfreal*fac5
268 ELSEIF (ijk==dirj)
THEN
272 ijkn = indijk(
i ,
j ,
k ,inoff,ijnoff)
273 ijkn1 = indijk(
i+1 ,
j ,
k ,inoff,ijnoff)
274 ds(
j) = 0.5_rfreal*(segm(segid,ijkn)+segm(segid,ijkn1))
281 ijkc = indijk(
i ,
j ,
k ,icoff,ijcoff)
282 rwidth1= 1._rfreal/(ds(
j-1)+2._rfreal*ds(
j)+ds(
j+1))
283 fac1 = (ds(
j-1)+ds(
j))*rwidth1
284 fac3 = (ds(
j+1)+ds(
j))*rwidth1
285 cccofa(1,ijkc) = 0.5_rfreal*fac1
286 cccofa(2,ijkc) = 0.5_rfreal
287 cccofa(3,ijkc) = 0.5_rfreal*fac3
289 rwidth1= 1._rfreal/ &
290 (ds(
j-2)+ds(
j+2)+2._rfreal*(ds(
j-1)+ds(
j)+ds(
j+1)))
291 fac1 = (ds(
j-2)+ds(
j-1))*rwidth1
292 fac2 = (ds(
j-2)+2._rfreal*ds(
j-1)+ds(
j))*rwidth1
293 fac3 = (ds(
j-1)+2._rfreal*ds(
j)+ds(
j+1))*rwidth1
294 fac4 = (ds(
j)+2._rfreal*ds(
j+1)+ds(
j+2))*rwidth1
295 fac5 = (ds(
j+1)+ds(
j+2))*rwidth1
296 cccofb(1,ijkc) = 0.5_rfreal*fac1
297 cccofb(2,ijkc) = 0.5_rfreal*fac2
298 cccofb(3,ijkc) = 0.5_rfreal*fac3
299 cccofb(4,ijkc) = 0.5_rfreal*fac4
300 cccofb(5,ijkc) = 0.5_rfreal*fac5
305 ELSEIF (ijk==dirk)
THEN
309 ijkn = indijk(
i ,
j ,
k ,inoff,ijnoff)
310 ijkn1 = indijk(
i ,
j+1 ,
k ,inoff,ijnoff)
311 ds(
k) = 0.5_rfreal*(segm(segid,ijkn)+segm(segid,ijkn1))
315 ds(kend+
k) = ds(kend)
318 ijkc = indijk(
i ,
j ,
k ,icoff,ijcoff)
319 rwidth1= 1._rfreal/(ds(
k-1)+2._rfreal*ds(
k)+ds(
k+1))
320 fac1 = (ds(
k-1)+ds(
k))*rwidth1
321 fac3 = (ds(
k+1)+ds(
k))*rwidth1
322 cccofa(1,ijkc) = 0.5_rfreal*fac1
323 cccofa(2,ijkc) = 0.5_rfreal
324 cccofa(3,ijkc) = 0.5_rfreal*fac3
326 rwidth1= 1._rfreal/ &
327 (ds(
k-2)+ds(
k+2)+2._rfreal*(ds(
k-1)+ds(
k)+ds(
k+1)))
328 fac1 = (ds(
k-2)+ds(
k-1))*rwidth1
329 fac2 = (ds(
k-2)+2._rfreal*ds(
k-1)+ds(
k))*rwidth1
330 fac3 = (ds(
k-1)+2._rfreal*ds(
k)+ds(
k+1))*rwidth1
331 fac4 = (ds(
k)+2._rfreal*ds(
k+1)+ds(
k+2))*rwidth1
332 fac5 = (ds(
k+1)+ds(
k+2))*rwidth1
333 cccofb(1,ijkc) = 0.5_rfreal*fac1
334 cccofb(2,ijkc) = 0.5_rfreal*fac2
335 cccofb(3,ijkc) = 0.5_rfreal*fac3
336 cccofb(4,ijkc) = 0.5_rfreal*fac4
337 cccofb(5,ijkc) = 0.5_rfreal*fac5
**********************************************************************Rocstar Simulation Suite Illinois Rocstar LLC All rights reserved ****Illinois Rocstar LLC IL **www illinoisrocstar com **sales illinoisrocstar com WITHOUT WARRANTY OF ANY **EXPRESS OR INCLUDING BUT NOT LIMITED TO THE WARRANTIES **OF FITNESS FOR A PARTICULAR PURPOSE AND **NONINFRINGEMENT IN NO EVENT SHALL THE CONTRIBUTORS OR **COPYRIGHT HOLDERS BE LIABLE FOR ANY DAMAGES OR OTHER WHETHER IN AN ACTION OF TORT OR **Arising OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE **USE OR OTHER DEALINGS WITH THE SOFTWARE **********************************************************************INTERFACE SUBROUTINE ibeg
subroutine registerfunction(global, funName, fileName)
**********************************************************************Rocstar Simulation Suite Illinois Rocstar LLC All rights reserved ****Illinois Rocstar LLC IL **www illinoisrocstar com **sales illinoisrocstar com WITHOUT WARRANTY OF ANY **EXPRESS OR INCLUDING BUT NOT LIMITED TO THE WARRANTIES **OF FITNESS FOR A PARTICULAR PURPOSE AND **NONINFRINGEMENT IN NO EVENT SHALL THE CONTRIBUTORS OR **COPYRIGHT HOLDERS BE LIABLE FOR ANY DAMAGES OR OTHER WHETHER IN AN ACTION OF TORT OR **Arising OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE **USE OR OTHER DEALINGS WITH THE SOFTWARE **********************************************************************INTERFACE SUBROUTINE knode iend
subroutine turb_flolesgencocclo(global, ijk, ibeg, iend, jbeg, jend, kbeg, kend, minIdx, maxIdx, segId, iCOff, ijCOff, iNOff, ijNOff, ds, segm, ccCofA, ccCofB)
**********************************************************************Rocstar Simulation Suite Illinois Rocstar LLC All rights reserved ****Illinois Rocstar LLC IL **www illinoisrocstar com **sales illinoisrocstar com WITHOUT WARRANTY OF ANY **EXPRESS OR INCLUDING BUT NOT LIMITED TO THE WARRANTIES **OF FITNESS FOR A PARTICULAR PURPOSE AND **NONINFRINGEMENT IN NO EVENT SHALL THE CONTRIBUTORS OR **COPYRIGHT HOLDERS BE LIABLE FOR ANY DAMAGES OR OTHER WHETHER IN AN ACTION OF TORT OR **Arising OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE **USE OR OTHER DEALINGS WITH THE SOFTWARE **********************************************************************INTERFACE SUBROUTINE knode jend
**********************************************************************Rocstar Simulation Suite Illinois Rocstar LLC All rights reserved ****Illinois Rocstar LLC IL **www illinoisrocstar com **sales illinoisrocstar com WITHOUT WARRANTY OF ANY **EXPRESS OR INCLUDING BUT NOT LIMITED TO THE WARRANTIES **OF FITNESS FOR A PARTICULAR PURPOSE AND **NONINFRINGEMENT IN NO EVENT SHALL THE CONTRIBUTORS OR **COPYRIGHT HOLDERS BE LIABLE FOR ANY DAMAGES OR OTHER WHETHER IN AN ACTION OF TORT OR **Arising OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE **USE OR OTHER DEALINGS WITH THE SOFTWARE **********************************************************************INTERFACE SUBROUTINE knode jbeg
subroutine turb_flolesgencocchi(global, ijk, ibeg, iend, jbeg, jend, kbeg, kend, minIdx, maxIdx, segId, iCOff, ijCOff, iNOff, ijNOff, ds, segm, ccCofA, ccCofB)
**********************************************************************Rocstar Simulation Suite Illinois Rocstar LLC All rights reserved ****Illinois Rocstar LLC IL **www illinoisrocstar com **sales illinoisrocstar com WITHOUT WARRANTY OF ANY **EXPRESS OR INCLUDING BUT NOT LIMITED TO THE WARRANTIES **OF FITNESS FOR A PARTICULAR PURPOSE AND **NONINFRINGEMENT IN NO EVENT SHALL THE CONTRIBUTORS OR **COPYRIGHT HOLDERS BE LIABLE FOR ANY DAMAGES OR OTHER WHETHER IN AN ACTION OF TORT OR **Arising OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE **USE OR OTHER DEALINGS WITH THE SOFTWARE **********************************************************************INTERFACE SUBROUTINE knode kbeg
subroutine deregisterfunction(global)