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 |
if(2.eq.1 )then |
252 |
print *,' In step fizhi corr, new fizhi fields ',bi,' dt= ',dt |
print *,' In step fizhi corr, new fizhi fields ',bi,' dt= ',dt |
253 |
do L = 1,Nrphys |
do L = 1,Nrphys |
254 |
do j = jdim1,jdim2 |
do j = jdim1,jdim2 |
280 |
tempij(i,j) = sphy(i,j,L,bi,bj) |
tempij(i,j) = sphy(i,j,L,bi,bj) |
281 |
enddo |
enddo |
282 |
enddo |
enddo |
283 |
c print *,' sphy at level ',l,' ',tempij |
print *,' sphy at level ',l,' ',tempij |
284 |
enddo |
enddo |
285 |
endif |
endif |
286 |
|
|
287 |
call qcheck (im1,im2,jm1,jm2,Nrphys,idim1,idim2,jdim1,jdim2, |
call qcheck (im1,im2,jm1,jm2,Nrphys,Nsx,Nsy,idim1,idim2, |
288 |
. dpphys(1,1,1,bi,bj),pephy(1,1,1,bi,bj),sphy(1,1,1,bi,bj)) |
. jdim1,jdim2,bi,bj,dpphys,sphy) |
289 |
|
|
290 |
|
if(2.eq.1 )then |
291 |
|
print *,' In step fizhi corr after qcheck ',bi |
292 |
|
do L = 1,Nrphys |
293 |
|
do j = jdim1,jdim2 |
294 |
|
do i = idim1,idim2 |
295 |
|
tempij(i,j) = sphy(i,j,L,bi,bj) |
296 |
|
enddo |
297 |
|
enddo |
298 |
|
print *,' sphy after qcheck at level ',l,' ',tempij |
299 |
|
enddo |
300 |
|
endif |
301 |
|
|
302 |
|
|
303 |
enddo |
enddo |
304 |
enddo |
enddo |
305 |
|
|