105 |
real pcheckpress(numpcheck) |
real pcheckpress(numpcheck) |
106 |
data pcheckpress/950.,750.,700.,400.,10./ |
data pcheckpress/950.,750.,700.,400.,10./ |
107 |
|
|
108 |
|
integer low_level,mid_level,nltop,nsubmin,nsubmax,Lup |
109 |
integer ndmoist,ndturb,ndlw,ndsw |
integer ndmoist,ndturb,ndlw,ndsw |
110 |
integer istrip,npcs |
integer istrip,npcs |
111 |
integer i,j,L,n |
integer i,j,L,n |
134 |
C Fill array of model levels closest to a given pressure value |
C Fill array of model levels closest to a given pressure value |
135 |
|
|
136 |
call getpwhere(myid,numpcheck,pcheckpress,pchecklevs) |
call getpwhere(myid,numpcheck,pcheckpress,pchecklevs) |
137 |
|
low_level = pchecklevs(3) |
138 |
|
mid_level = pchecklevs(4) |
139 |
|
nltop = pchecklevs(5) |
140 |
|
nsubmin = pchecklevs(1) |
141 |
|
nsubmax = pchecklevs(2) |
142 |
|
Lup = pchecklevs(3) |
143 |
|
|
144 |
C ********************************************************************** |
C ********************************************************************** |
145 |
C **** Call Physics Mini-Drivers **** |
C **** Call Physics Mini-Drivers **** |
150 |
IF ( alarm('radsw') ) THEN |
IF ( alarm('radsw') ) THEN |
151 |
|
|
152 |
call swrio ( nymd,nhms,bi,bj,ndsw,myid,istrip,npcs, |
call swrio ( nymd,nhms,bi,bj,ndsw,myid,istrip,npcs, |
153 |
|
. low_level,mid_level, |
154 |
. p,pl,ple,dpres,pkht,pkl,t,chemq,chemo3,co2, |
. p,pl,ple,dpres,pkht,pkl,t,chemq,chemo3,co2, |
155 |
. albvisdr,albvisdf,albnirdr,albnirdf,swdt,swdtclr, |
. albvisdr,albvisdf,albnirdr,albnirdf,swdt,swdtclr, |
156 |
. radswg,swgclr,fdifpar,fdirpar,osr,osrclr,im,jm,lm, |
. radswg,swgclr,fdifpar,fdirpar,osr,osrclr,im,jm,lm, |
173 |
enddo |
enddo |
174 |
|
|
175 |
call lwrio ( nymd,nhms,bi,bj,istrip,npcs, |
call lwrio ( nymd,nhms,bi,bj,istrip,npcs, |
176 |
|
. low_level,mid_level, |
177 |
. p,pl,ple,dpres,pkht,pkl,t,chemq,pkht,chemo3,co2, |
. p,pl,ple,dpres,pkht,pkl,t,chemq,pkht,chemo3,co2, |
178 |
. cfc11,cfc12,cfc22,methane,n2o,emiss,tgz,radlwg,st4,dst4, |
. cfc11,cfc12,cfc22,methane,n2o,emiss,tgz,radlwg,st4,dst4, |
179 |
. lwdt,dlwdtg,lwdtclr,lwgclr,im,jm,lm,ptop, |
. lwdt,dlwdtg,lwdtclr,lwgclr,im,jm,lm,ptop, |
213 |
C =============== |
C =============== |
214 |
if ( alarm('moist') ) then |
if ( alarm('moist') ) then |
215 |
call moistio (ndmoist,istrip,npcs, |
call moistio (ndmoist,istrip,npcs, |
216 |
|
. low_level,mid_level,nltop,nsubmin,nsubmax,Lup, |
217 |
. p,pl,ple,dpres,pkht,pkl,t,q,bi,bj,ntracer,ptracer,qq, |
. p,pl,ple,dpres,pkht,pkl,t,q,bi,bj,ntracer,ptracer,qq, |
218 |
. moistu,moistv,moistt,moistq,im,jm,lm,ptop,iras, |
. moistu,moistv,moistt,moistq,im,jm,lm,ptop,iras, |
219 |
. rainlsp,rainconv,snowfall, |
. rainlsp,rainconv,snowfall, |