/[MITgcm]/MITgcm/pkg/fizhi/fizhi_driver.F
ViewVC logotype

Diff of /MITgcm/pkg/fizhi/fizhi_driver.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.3 by molod, Thu Jun 17 16:07:06 2004 UTC revision 1.6 by molod, Mon Jun 28 16:12:17 2004 UTC
# Line 2  C $Header$ Line 2  C $Header$
2  C $Name$  C $Name$
3                                                                                                                                                                    
4  #include "CPP_OPTIONS.h"  #include "CPP_OPTIONS.h"
5        subroutine fizhi_driver (myid,im,jm,lm,ptracer,ntracer,        subroutine fizhi_driver (myid,im,jm,lm,bi,bj,ptracer,ntracer,
6       . xlats,xlons,       . xlats,xlons,
7       . p,u,v,t,q,pl,ple,dpres,pkht,pkl,fracland,landtype,radswt,       . p,u,v,t,q,pl,ple,dpres,pkht,pkl,fracland,landtype,radswt,
8       . phis_var,tgz,sea_ice,       . phis_var,tgz,sea_ice,
# Line 13  C $Name$ Line 13  C $Name$
13       . chemq,chemo3,co2,cfc11,cfc12,cfc22,methane,n2o,       . chemq,chemo3,co2,cfc11,cfc12,cfc22,methane,n2o,
14       . lwdt,lwdtclr,swdt,swdtclr,turbu,turbv,turbt,turbq,       . lwdt,lwdtclr,swdt,swdtclr,turbu,turbv,turbt,turbq,
15       . moistu,moistv,moistt,moistq,       . moistu,moistv,moistt,moistq,
16       . radswg,swgclr,albedo,fdirpar,fdifpar,osr,osrclr,tg0,radlwg,       . radswg,swgclr,fdirpar,fdifpar,osr,osrclr,tg0,radlwg,
17       . st4,dst4,dlwdtg,rainlsp,raincon,snowfall,iras,nlwcld,       . st4,dst4,dlwdtg,rainlsp,raincon,snowfall,iras,nlwcld,
18       . cldtot_lw,clras_lw,cldlsp_lw,nlwlz,lwlz,nswcld,cldtot_sw,       . cldtot_lw,clras_lw,cldlsp_lw,nlwlz,lwlz,nswcld,cldtot_sw,
19       . clras_sw,cldlsp_sw,nswlz,swlz,imstturbsw,imstturblw,qliqavesw,       . clras_sw,cldlsp_sw,nswlz,swlz,imstturbsw,imstturblw,qliqavesw,
20       . qliqavelw,fccavesw,fccavelw,qq,u2m,v2m,t2m,q2m,u10m,v10m,t10m,       . qliqavelw,fccavesw,fccavelw,qq)
      . q10m,gwet,snow)  
21  C***********************************************************************          C***********************************************************************        
22  C  Purpose                                                                        C  Purpose                                                                      
23  C  -------                                                                        C  -------                                                                      
# Line 47  c ------------- Line 46  c -------------
46    
47  c Input Parameters  c Input Parameters
48  c ----------------  c ----------------
49        integer myid,im,jm,lm,ptracer,ntracer        integer myid,im,jm,lm,bi,bj,ptracer,ntracer
50        integer nchp,chlat,chlon,igrd,nchpland,chfr,ityp        integer nchp,chlat,chlon,igrd,nchpland,chfr,ityp
51        integer nwatlevs,nwatlats        integer nwatlevs,nwatlats
52        integer nozolevs,nozolats        integer nozolevs,nozolats
# Line 81  c ---------------- Line 80  c ----------------
80        real albvisdr(im,jm),albvisdf(im,jm),albnirdr(im,jm)        real albvisdr(im,jm),albvisdf(im,jm),albnirdr(im,jm)
81        real albnirdf(im,jm),emiss(im,jm)        real albnirdf(im,jm),emiss(im,jm)
82        real alai(nchp),agrn(nchp)        real alai(nchp),agrn(nchp)
83        real radswg(im,jm),swgclr(im,jm),albedo(im,jm)        real radswg(im,jm),swgclr(im,jm)
84        real fdirpar(im,jm),fdifpar(im,jm),osr(im,jm),osrclr(im,jm)        real fdirpar(im,jm),fdifpar(im,jm),osr(im,jm),osrclr(im,jm)
85        real tg0(im,jm),radlwg(im,jm),st4(im,jm),dst4(im,jm)        real tg0(im,jm),radlwg(im,jm),st4(im,jm),dst4(im,jm)
86        real dlwdtg(im,jm,lm)        real dlwdtg(im,jm,lm)
# Line 93  c ---------------- Line 92  c ----------------
92        real qliqavesw(im,jm,lm),qliqavelw(im,jm,lm)        real qliqavesw(im,jm,lm),qliqavelw(im,jm,lm)
93        real fccavesw(im,jm,lm),fccavelw(im,jm,lm)        real fccavesw(im,jm,lm),fccavelw(im,jm,lm)
94        real qq(im,jm,lm)        real qq(im,jm,lm)
95        real u2m(im,jm),v2m(im,jm),t2m(im,jm),q2m(im,jm),u10m(im,jm)        real tgz(im,jm)
       real v10m(im,jm),t10m(im,jm),q10m(im,jm),tgz(im,jm)  
       real gwet(im,jm),snow(im,jm)  
96    
97  c Local Variables  c Local Variables
98  c ---------------  c ---------------
99        logical     alarm        logical     alarm
100        external    alarm        external    alarm
101          real massinlayer(im,jm,lm)
102    
103        integer ndmoist,ndturb,ndlw,ndsw        integer ndmoist,ndturb,ndlw,ndsw
104        integer istrip,npcs        integer istrip,npcs
# Line 127  C ************************************** Line 125  C **************************************
125        call get_alarm ( 'pnt',ndum,ndum,ndpnt,ndum )        call get_alarm ( 'pnt',ndum,ndum,ndpnt,ndum )
126        lpnt = ndpnt.ne.0        lpnt = ndpnt.ne.0
127    
128  C Compute pl,ple,pkht,pkl,dpres   and pcheck for use by mini-drivers  C pchecks for use by mini-drivers
   
 C Make 'model theta' out of theta, ie, divide by p0**kappa  
129    
130  C **********************************************************************          C **********************************************************************        
131  C ****                  Call Physics Mini-Drivers                   ****          C ****                  Call Physics Mini-Drivers                   ****        
# Line 139  C SHORT WAVE RADIATION Line 135  C SHORT WAVE RADIATION
135  C ====================  C ====================
136        IF ( alarm('radsw') ) THEN        IF ( alarm('radsw') ) THEN
137    
138        call swrio ( nymd,nhms,ndsw,myid,istrip,npcs,        call swrio ( nymd,nhms,bi,bj,ndsw,myid,istrip,npcs,
139       .   p,pl,ple,dpres,pkht,pkl,t,chemq,chemo3,co2,       .   p,pl,ple,dpres,pkht,pkl,t,chemq,chemo3,co2,
140       .   albvisdr,albvisdf,albnirdr,albnirdf,swdt,swdtclr,       .   albvisdr,albvisdf,albnirdr,albnirdf,swdt,swdtclr,
141       .   radswg,swgclr,albedo,fdifpar,fdirpar,osr,osrclr,im,jm,lm,       .   radswg,swgclr,fdifpar,fdirpar,osr,osrclr,im,jm,lm,
142       .   ptop,nswcld,cldtot_sw,cldras_sw,nswlz,swlz,       .   ptop,nswcld,cldtot_sw,cldras_sw,nswlz,swlz,
143       .   .false.,imstturbsw,qliqavesw,       .   .false.,imstturbsw,qliqavesw,
144       .   fccavesw,landtype,xlats,xlons )       .   fccavesw,landtype,xlats,xlons )
# Line 161  c -------------------------------- Line 157  c --------------------------------
157        enddo        enddo
158        enddo        enddo
159    
160        call lwrio ( nymd,nhms,istrip,npcs,        call lwrio ( nymd,nhms,bi,bj,istrip,npcs,
161       .   p,pl,ple,dpres,pkht,pkl,t,chemq,pkht,chemo3,co2,       .   p,pl,ple,dpres,pkht,pkl,t,chemq,pkht,chemo3,co2,
162       .   cfc11,cfc12,cfc22,methane,n2o,emiss,tgz,radlwg,st4,dst4,       .   cfc11,cfc12,cfc22,methane,n2o,emiss,tgz,radlwg,st4,dst4,
163       .   lwdt,dlwdtg,lwdtclr,lwgclr,im,jm,lm,ptop,       .   lwdt,dlwdtg,lwdtclr,lwgclr,im,jm,lm,ptop,
# Line 175  C TURBULENCE Line 171  C TURBULENCE
171  C ==========  C ==========
172        IF ( alarm('turb') ) THEN        IF ( alarm('turb') ) THEN
173    
174        call turbio (im,jm,lm,1,nymd,nhms,ndturb,        call turbio (im,jm,lm,1,nymd,nhms,bi,bj,ndturb,
175       .   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,
176       .   ctmt,xxmt,yymt,zetamt,xlmt,khmt,tke,       .   ctmt,xxmt,yymt,zetamt,xlmt,khmt,tke,
177       .   tgz,gwet,snow,fracland,landtype,       .   tgz,fracland,landtype,
178       .   tcanopy,ecanopy,tdeep,swetshal,swetroot,swetdeep,snodep,capac,       .   tcanopy,ecanopy,tdeep,swetshal,swetroot,swetdeep,snodep,capac,
179       .   nchp,nchpland,chfr,chlt,chlon,igrd,ityp,alai,agrn,sea_ice,lpnt,       .   nchp,nchpland,chfr,chlt,chlon,igrd,ityp,alai,agrn,sea_ice,lpnt,
180       .   turbu,turbv,turbt,turbq,radlwg,st4,dst4,radswg,radswt,       .   turbu,turbv,turbt,turbq,radlwg,st4,dst4,radswg,radswt,
181       .   fdifpar,fdirpar,rainlsp,rainconv,snowfall,tg0,       .   fdifpar,fdirpar,rainlsp,rainconv,snowfall,tg0,
182       .   imstturblw,imstturbsw,qliqavelw,qliqavesw,fccavelw,fccavesw,qq,       .   imstturblw,imstturbsw,qliqavelw,qliqavesw,fccavelw,fccavesw,qq,
183       .   u2m ,v2m ,t2m ,q2m,u10m,v10m,t10m,q10m,myid,comm )       .   myid,comm )
184    
185  c Add Gravity-Wave Drag Tendency  c Add Gravity-Wave Drag Tendency
186  c ------------------------------  c ------------------------------
# Line 201  C MOIST PROCESSES Line 197  C MOIST PROCESSES
197  C ===============  C ===============
198        if ( alarm('moist') ) then        if ( alarm('moist') ) then
199        call moistio (ndmoist,istrip,npcs,        call moistio (ndmoist,istrip,npcs,
200       .    p,pl,ple,dpres,pkht,pkl,t,q,ntracer,ptracer,qq,       .    p,pl,ple,dpres,pkht,pkl,t,q,bi,bj,ntracer,ptracer,qq,
201       .    moistu,moistv,moistt,moistq,im,jm,lm,ptop,iras,       .    moistu,moistv,moistt,moistq,im,jm,lm,ptop,iras,
202       .    rainlsp,rainconv,snowfall,       .    rainlsp,rainconv,snowfall,
203       .    nswcld,cldtot_sw,cldras_sw,cldlsp_sw,nswlz,swlz,       .    nswcld,cldtot_sw,cldras_sw,cldlsp_sw,nswlz,swlz,
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.6
  ViewVC Help
Powered by ViewVC 1.1.22