1 |
C $Header$ |
C $Header$ |
2 |
C $Name$ |
C $Name$ |
3 |
|
|
4 |
|
#include "CPP_OPTIONS.h" |
5 |
subroutine turbio (im,jm,nlay,maxtypes,istrip,nymd,nhms |
subroutine turbio (im,jm,nlay,maxtypes,istrip,nymd,nhms |
6 |
1 ,ndturb,ptop, pz, uz, vz, tz, qz, ntracers,ptracers, |
1 ,ndturb,ptop, pz, uz, vz, tz, qz, ntracers,ptracers, |
7 |
2 ,plz,plze,dpres,pkht,pkz,ctmt,xxmt,yymt,zetamt,xlmt,khmt,tke |
2 ,plz,plze,dpres,pkht,pkz,ctmt,xxmt,yymt,zetamt,xlmt,khmt,tke |
9 |
4 ,tcanopy,ecanopy,tdeep,swetshal,swetroot,swetdeep,snodep,capac |
4 ,tcanopy,ecanopy,tdeep,swetshal,swetroot,swetdeep,snodep,capac |
10 |
5 ,nchp,nchplnd,chfr,chlt,chlon,igrd,ityp,alai,agrn |
5 ,nchp,nchplnd,chfr,chlt,chlon,igrd,ityp,alai,agrn |
11 |
6 ,surftype,tilefrac,thkz,tprof |
6 ,surftype,tilefrac,thkz,tprof |
|
7 ,qdiag,ndiagsiz,qdiaglsm,ndlsm |
|
12 |
8 ,duturb, dvturb, dtturb,dqturb,radlwg,st4,dst4,radswg,radswt |
8 ,duturb, dvturb, dtturb,dqturb,radlwg,st4,dst4,radswg,radswt |
13 |
9 ,fdifpar,fdirpar,rainlsp,rainconv,snowfall,tempref |
9 ,fdifpar,fdirpar,rainlsp,rainconv,snowfall,tempref |
14 |
1 ,imstturblw,imstturbsw,qliqavelw,qliqavesw,fccavelw,fccavesw |
1 ,imstturblw,imstturbsw,qliqavelw,qliqavesw,fccavelw,fccavesw |
82 |
c swetdeep- deep soil level moisture field capacity fraction real[tiles] |
c swetdeep- deep soil level moisture field capacity fraction real[tiles] |
83 |
c snodep - depth of snow pack in cm liquid water equiv real[tiles] |
c snodep - depth of snow pack in cm liquid water equiv real[tiles] |
84 |
c capac - leaf canopy water reservoir in cm real[tiles] |
c capac - leaf canopy water reservoir in cm real[tiles] |
|
c qdiag - diagnostics array - real[im,jm,ndiagsiz] |
|
|
c qdiaglsm- tile diagnostics array - real[nchp,ndlsm] |
|
85 |
c output: |
c output: |
86 |
c duturb - change in zonal wind component due to turbulent processes |
c duturb - change in zonal wind component due to turbulent processes |
87 |
c per unit time in m/sec**2 - real[lon,lat,levels] |
c per unit time in m/sec**2 - real[lon,lat,levels] |
109 |
implicit none |
implicit none |
110 |
|
|
111 |
#include "mpif.h" |
#include "mpif.h" |
112 |
|
#ifdef ALLOW_DIAGNOSTICS |
113 |
#include "diagnostics.h" |
#include "diagnostics.h" |
114 |
|
#endif |
115 |
|
|
116 |
integer n,im,jm,nlay,maxtypes,istrip,nymd,nhms,ndturb,ndiagsiz |
integer n,im,jm,nlay,maxtypes,istrip,nymd,nhms,ndturb |
|
integer ndlsm |
|
117 |
real ptop |
real ptop |
118 |
integer nchp,nchplnd,landtype(im,jm),igrd(nchp) |
integer nchp,nchplnd,landtype(im,jm),igrd(nchp) |
119 |
integer ityp(nchp) |
integer ityp(nchp) |
148 |
real thkz(im,jm) |
real thkz(im,jm) |
149 |
|
|
150 |
logical tprof |
logical tprof |
|
real qdiag(im,jm,ndiagsiz) |
|
|
real qdiaglsm(nchp,ndlsm) |
|
151 |
|
|
152 |
real duturb(im,jm,nlay),dvturb(im,jm,nlay) |
real duturb(im,jm,nlay),dvturb(im,jm,nlay) |
153 |
real dtturb(im,jm,nlay),dqturb(im,jm,nlay,ntracers) |
real dtturb(im,jm,nlay),dqturb(im,jm,nlay,ntracers) |
266 |
real qliq(istrip,nlay),turbfcc(istrip,nlay) |
real qliq(istrip,nlay),turbfcc(istrip,nlay) |
267 |
real qliqmsc(nchp,nlay),fccmsc(nchp,nlay) |
real qliqmsc(nchp,nlay),fccmsc(nchp,nlay) |
268 |
|
|
269 |
|
integer ndlsm |
270 |
|
parameter ( ndlsm = 1 ) |
271 |
|
real qdiaglsm(nchp,ndlsm) |
272 |
|
|
273 |
integer nsecf,nmonf,ndayf |
integer nsecf,nmonf,ndayf |
274 |
nsecf(n) = n/10000*3600 + mod(n,10000)/100* 60 + mod(n,100) |
nsecf(n) = n/10000*3600 + mod(n,10000)/100* 60 + mod(n,100) |
275 |
nmonf(n) = mod(n,10000)/100 |
nmonf(n) = mod(n,10000)/100 |