179 |
enddo |
enddo |
180 |
enddo |
enddo |
181 |
|
|
182 |
call lwrio ( nymd,nhms,bi,bj,istrip,npcs, |
call lwrio ( nymd,nhms,bi,bj,myid,istrip,npcs, |
183 |
. low_level,mid_level,im,jm,lm, |
. low_level,mid_level,im,jm,lm, |
184 |
. p,pl,ple,dpres,pkht,pkl,t,chemq,chemo3,co2, |
. p,pl,ple,dpres,pkht,pkl,t,chemq,chemo3,co2, |
185 |
. cfc11,cfc12,cfc22,methane,n2o,emiss,tgz,radlwg,st4,dst4, |
. cfc11,cfc12,cfc22,methane,n2o,emiss,tgz,radlwg,st4,dst4, |
194 |
C ========== |
C ========== |
195 |
IF ( alarm('turb') ) THEN |
IF ( alarm('turb') ) THEN |
196 |
|
|
197 |
call turbio (im,jm,lm,istrip,nymd,nhms,bi,bj,ndturb, |
call turbio (im,jm,lm,istrip,nymd,nhms,bi,bj,ndturb,nltop, |
198 |
. 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, |
199 |
. ctmt,xxmt,yymt,zetamt,xlmt,khmt,tke, |
. ctmt,xxmt,yymt,zetamt,xlmt,khmt,tke, |
200 |
. tgz,fracland,landtype, |
. tgz,fracland,landtype, |
208 |
|
|
209 |
c Add Gravity-Wave Drag Tendency |
c Add Gravity-Wave Drag Tendency |
210 |
c ------------------------------ |
c ------------------------------ |
211 |
C Comment this out for now - run tropospheric resolutions only. |
C Comment this out for now |
|
C Need to see how to get zonal wavelengths on the cube. |
|
212 |
|
|
213 |
c call gwdrag (myid,p,pl,ple,dpres,pkht,pkl,u,v,t,q,phis_var, |
c call gwdrag (myid,p,pl,ple,dpres,pkl,u,v,t,q,phis_var, |
214 |
c . turbu,turbv,turbt,im,jm,lm,ptop,istrip,npcs, |
c . turbu,turbv,turbt,im,jm,lm,bi,bj,istrip,npcs,imglobal) |
|
c . imglobal) |
|
215 |
|
|
216 |
c Add Rayleigh Friction Damping Above 70 Km |
c Add Rayleigh Friction Damping Above 70 Km |
217 |
c ----------------------------------------- |
c ----------------------------------------- |