124 |
logical alarm |
logical alarm |
125 |
external alarm |
external alarm |
126 |
|
|
127 |
save lwdt,lwdtclr,swdt,swdtclr,turbu,turbv,turbt,turbq |
c save lwdt,lwdtclr,swdt,swdtclr,turbu,turbv,turbt,turbq |
128 |
save moistu,moistv,moistt,moistq |
c save moistu,moistv,moistt,moistq |
129 |
save radswg,swgclr,fdirpar,fdifpar,osr,osrclr,tg0,radlwg |
c save radswg,swgclr,fdirpar,fdifpar,osr,osrclr,tg0,radlwg |
130 |
save st4,dst4,dlwdtg,rainlsp,raincon,snowfall,iras |
c save st4,dst4,dlwdtg,rainlsp,raincon,snowfall,iras |
131 |
save nlwcld,cldtot_lw,cldras_lw,cldlsp_lw,nlwlz,lwlz |
c save nlwcld,cldtot_lw,cldras_lw,cldlsp_lw,nlwlz,lwlz |
132 |
save nswcld,cldtot_sw,cldras_sw,cldlsp_sw,nswlz,swlz |
c save nswcld,cldtot_sw,cldras_sw,cldlsp_sw,nswlz,swlz |
133 |
save imstturbsw,imstturblw,qliqavesw,qliqavelw,fccavesw,fccavelw |
c save imstturbsw,imstturblw,qliqavesw,qliqavelw,fccavesw,fccavelw |
134 |
save qq |
c save qq |
135 |
save pl,ple,dpres,pkle,pkl |
c save pl,ple,dpres,pkle,pkl |
136 |
|
|
137 |
|
common /saver/ lwdt,lwdtclr,swdt,swdtclr,turbu,turbv,turbt,turbq |
138 |
|
common /saver/ moistu,moistv,moistt,moistq |
139 |
|
common /saver/ radswg,swgclr,fdirpar,fdifpar,osr,osrclr,tg0,radlwg |
140 |
|
common /saver/ st4,dst4,dlwdtg,rainlsp,raincon,snowfall,iras |
141 |
|
common /saver/ nlwcld,cldtot_lw,cldras_lw,cldlsp_lw,nlwlz,lwlz |
142 |
|
common /saver/ nswcld,cldtot_sw,cldras_sw,cldlsp_sw,nswlz,swlz |
143 |
|
common /saver/ imstturbsw,imstturblw,qliqavesw,qliqavelw,fccavesw |
144 |
|
common /saver/ fccavelw |
145 |
|
common /saver/ qq |
146 |
|
common /saver/ pl,ple,dpres,pkle,pkl |
147 |
|
|
148 |
C*********************************************************************** |
C*********************************************************************** |
149 |
C Unshadow input arrays (and make 'fizhi theta' from true theta) |
C Unshadow input arrays (and make 'fizhi theta' from true theta) |
150 |
C*********************************************************************** |
C*********************************************************************** |
249 |
call pkappa (im2,jm2,Nrphys,ple(1,1,1,bi,bj),pkle(1,1,1,bi,bj), |
call pkappa (im2,jm2,Nrphys,ple(1,1,1,bi,bj),pkle(1,1,1,bi,bj), |
250 |
. pkl(1,1,1,bi,bj)) |
. pkl(1,1,1,bi,bj)) |
251 |
|
|
252 |
if(bi.eq.1)then |
if(bi.eq.1 .and. 1.eq.1 )then |
253 |
print *,' In do fizhi, before fizhi driver - bi = ',bi |
print *,' In do fizhi, before fizhi driver - bi = ',bi |
254 |
do L = 1,Nrphys+1 |
do L = 1,Nrphys+1 |
255 |
do j = jm1,jm2 |
do j = jm1,jm2 |
377 |
. qliqavelw(1,1,1,bi,bj),fccavesw(1,1,1,bi,bj), |
. qliqavelw(1,1,1,bi,bj),fccavesw(1,1,1,bi,bj), |
378 |
. fccavelw(1,1,1,bi,bj),qq(1,1,1,bi,bj)) |
. fccavelw(1,1,1,bi,bj),qq(1,1,1,bi,bj)) |
379 |
|
|
|
endif |
|
|
|
|
380 |
if(bi.eq.1)then |
if(bi.eq.1)then |
381 |
print *,' In do fizhi, after fizhi driver - bi = ',bi |
print *,' In do fizhi, after fizhi driver - bi = ',bi |
382 |
do L = 1,Nrphys |
do L = 1,Nrphys |
513 |
enddo |
enddo |
514 |
enddo |
enddo |
515 |
|
|
516 |
|
endif |
517 |
|
|
518 |
call fizhi_step_diag(myid,p,uphy,vphy,thphy,sphy,qq,pkl,dpres, |
call fizhi_step_diag(myid,p,uphy,vphy,thphy,sphy,qq,pkl,dpres, |
519 |
. radswt,radswg,swgclr,osr,osrclr,st4,dst4,tgz,tg0,radlwg,lwgclr, |
. radswt,radswg,swgclr,osr,osrclr,st4,dst4,tgz,tg0,radlwg,lwgclr, |
520 |
. turbu,turbv,turbt,turbq,moistu,moistv,moistt,moistq, |
. turbu,turbv,turbt,turbq,moistu,moistv,moistt,moistq, |