49 |
integer im1, im2, jm1, jm2, idim1, idim2, jdim1, jdim2 |
integer im1, im2, jm1, jm2, idim1, idim2, jdim1, jdim2 |
50 |
_RL dt |
_RL dt |
51 |
|
|
52 |
|
_RL tempij(sNx,sNy) |
53 |
|
|
54 |
im1 = 1-OLx |
im1 = 1-OLx |
55 |
im2 = sNx+OLx |
im2 = sNx+OLx |
56 |
jm1 = 1-OLy |
jm1 = 1-OLy |
248 |
. Nrphys,Nsx,Nsy,1,sNx,1,sNy,bi,bj, |
. Nrphys,Nsx,Nsy,1,sNx,1,sNy,bi,bj, |
249 |
. uphytemp,vphytemp,thphytemp,sphytemp) |
. uphytemp,vphytemp,thphytemp,sphytemp) |
250 |
|
|
251 |
|
if(1.eq.1 )then |
252 |
|
print *,' In step fizhi corr, new fizhi fields ',bi,' dt= ',dt |
253 |
|
do L = 1,Nrphys |
254 |
|
do j = jdim1,jdim2 |
255 |
|
do i = idim1,idim2 |
256 |
|
tempij(i,j) = uphy(i,j,L,bi,bj) |
257 |
|
enddo |
258 |
|
enddo |
259 |
|
c print *,' uphy at level ',l,' ',tempij |
260 |
|
enddo |
261 |
|
do L = 1,Nrphys |
262 |
|
do j = jdim1,jdim2 |
263 |
|
do i = idim1,idim2 |
264 |
|
tempij(i,j) = vphy(i,j,L,bi,bj) |
265 |
|
enddo |
266 |
|
enddo |
267 |
|
c print *,' vphy at level ',l,' ',tempij |
268 |
|
enddo |
269 |
|
do L = 1,Nrphys |
270 |
|
do j = jdim1,jdim2 |
271 |
|
do i = idim1,idim2 |
272 |
|
tempij(i,j) = thphy(i,j,L,bi,bj) |
273 |
|
enddo |
274 |
|
enddo |
275 |
|
print *,' thphy at level ',l,' ',tempij |
276 |
|
enddo |
277 |
|
do L = 1,Nrphys |
278 |
|
do j = jdim1,jdim2 |
279 |
|
do i = idim1,idim2 |
280 |
|
tempij(i,j) = sphy(i,j,L,bi,bj) |
281 |
|
enddo |
282 |
|
enddo |
283 |
|
c print *,' sphy at level ',l,' ',tempij |
284 |
|
enddo |
285 |
|
endif |
286 |
|
|
287 |
|
|
288 |
enddo |
enddo |
289 |
enddo |
enddo |
290 |
|
|