191 |
C ========== |
C ========== |
192 |
IF ( alarm('turb') ) THEN |
IF ( alarm('turb') ) THEN |
193 |
|
|
|
if(2.eq.1)then |
|
|
print *,' In fizhi driver, before turb - bi = ',bi |
|
|
do j = 1,jm |
|
|
do i = 1,im |
|
|
tempij(i,j) = radswt(i,j) |
|
|
enddo |
|
|
enddo |
|
|
print *,' radswt ',tempij |
|
|
do j = 1,jm |
|
|
do i = 1,im |
|
|
tempij(i,j) = radswg(i,j) |
|
|
enddo |
|
|
enddo |
|
|
print *,' radswg ',tempij |
|
|
c do j = 1,jm |
|
|
c do i = 1,im |
|
|
c tempij(i,j) = st4(i,j) |
|
|
c enddo |
|
|
c enddo |
|
|
c print *,' st4 ',tempij |
|
|
c do j = 1,jm |
|
|
c do i = 1,im |
|
|
c tempij(i,j) = dst4(i,j) |
|
|
c enddo |
|
|
c enddo |
|
|
c print *,' dst4 ',tempij |
|
|
c do j = 1,jm |
|
|
c do i = 1,im |
|
|
c tempij(i,j) = radlwg(i,j) |
|
|
c enddo |
|
|
c enddo |
|
|
c print *,' radlwg ',tempij |
|
|
c do j = 1,jm |
|
|
c do i = 1,im |
|
|
c tempij(i,j) = fdifpar(i,j) |
|
|
c enddo |
|
|
c enddo |
|
|
c print *,' fdifpar ',tempij |
|
|
c do j = 1,jm |
|
|
c do i = 1,im |
|
|
c tempij(i,j) = fdirpar(i,j) |
|
|
c enddo |
|
|
c enddo |
|
|
c print *,' fdirpar ',tempij |
|
|
c do j = 1,jm |
|
|
c do i = 1,im |
|
|
c tempij(i,j) = tgz(i,j) |
|
|
c enddo |
|
|
c enddo |
|
|
c print *,' tgz ',tempij |
|
|
c do j = 1,jm |
|
|
c do i = 1,im |
|
|
c tempij(i,j) = tg0(i,j) |
|
|
c enddo |
|
|
c enddo |
|
|
c print *,' tg0 ',tempij |
|
|
endif |
|
|
|
|
194 |
call turbio (im,jm,lm,istrip,nymd,nhms,bi,bj,ndturb, |
call turbio (im,jm,lm,istrip,nymd,nhms,bi,bj,ndturb, |
195 |
. ptop,p,u,v,t,q,ntracer,ptracer,pl,ple,dpres,pkht,pkl, |
. ptop,p,u,v,t,q,ntracer,ptracer,pl,ple,dpres,pkht,pkl, |
196 |
. ctmt,xxmt,yymt,zetamt,xlmt,khmt,tke, |
. ctmt,xxmt,yymt,zetamt,xlmt,khmt,tke, |