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

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

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

revision 1.39 by molod, Wed Oct 6 23:23:34 2004 UTC revision 1.40 by molod, Thu Oct 7 00:12:28 2004 UTC
# Line 10  C $Name$ Line 10  C $Name$
10       . albvisdr,albvisdf,ityp,chfr,alai,agrn,igrd,chlat,chlon,       . albvisdr,albvisdf,ityp,chfr,alai,agrn,igrd,chlat,chlon,
11       . tcanopy,tdeep,ecanopy,swetshal,swetroot,swetdeep,snodep,capac,       . tcanopy,tdeep,ecanopy,swetshal,swetroot,swetdeep,snodep,capac,
12       . o3,qstr,co2,cfc11,cfc12,cfc22,n2o,methane,       . o3,qstr,co2,cfc11,cfc12,cfc22,n2o,methane,
13         . iras,nlwcld,cldtot_lw,cldras_lw,cldlsp_lw,nlwlz,lwlz,
14         . nswcld,cldtot_sw,cldras_sw,cldlsp_sw,nswlz,swlz,
15         . imstturbsw,imstturblw,qliqavesw,qliqavelw,fccavesw,fccavelw,
16       . duphy,dvphy,dthphy,dsphy)       . duphy,dvphy,dthphy,dsphy)
17  c-----------------------------------------------------------------------  c-----------------------------------------------------------------------
18  c Interface routine to calculate physics increments - calls fizhi_driver.  c Interface routine to calculate physics increments - calls fizhi_driver.
# Line 61  C Argument list declarations Line 64  C Argument list declarations
64        _RL o3(im2,jm2,Nrphin,Nsxin,Nsyin)        _RL o3(im2,jm2,Nrphin,Nsxin,Nsyin)
65        _RL qstr(im2,jm2,Nrphin,Nsxin,Nsyin)        _RL qstr(im2,jm2,Nrphin,Nsxin,Nsyin)
66        _RL co2,cfc11,cfc12,cfc22,n2o(Nrphin),methane(Nrphin)        _RL co2,cfc11,cfc12,cfc22,n2o(Nrphin),methane(Nrphin)
67    
68          integer iras,nlwcld,nlwlz,nswcld,nswlz
69          integer imstturbsw,imstturblw
70          _RL cldtot_lw(idim1:idim2,jdim1:jdim2,Nrphin,Nsxin,Nsyin)
71          _RL cldras_lw(idim1:idim2,jdim1:jdim2,Nrphin,Nsxin,Nsyin)
72          _RL cldlsp_lw(idim1:idim2,jdim1:jdim2,Nrphin,Nsxin,Nsyin)
73          _RL lwlz(idim1:idim2,jdim1:jdim2,Nrphin,Nsxin,Nsyin)
74          _RL cldtot_sw(idim1:idim2,jdim1:jdim2,Nrphin,Nsxin,Nsyin)
75          _RL cldras_sw(idim1:idim2,jdim1:jdim2,Nrphin,Nsxin,Nsyin)
76          _RL cldlsp_sw(idim1:idim2,jdim1:jdim2,Nrphin,Nsxin,Nsyin)
77          _RL swlz(idim1:idim2,jdim1:jdim2,Nrphin,Nsxin,Nsyin)
78          _RL qliqavesw(idim1:idim2,jdim1:jdim2,Nrphin,Nsxin,Nsyin)
79          _RL qliqavelw(idim1:idim2,jdim1:jdim2,Nrphin,Nsxin,Nsyin)
80          _RL fccavesw(idim1:idim2,jdim1:jdim2,Nrphin,Nsxin,Nsyin)
81          _RL fccavelw(idim1:idim2,jdim1:jdim2,Nrphin,Nsxin,Nsyin)
82    
83        _RL duphy(idim1:idim2,jdim1:jdim2,Nrphin,Nsxin,Nsyin)        _RL duphy(idim1:idim2,jdim1:jdim2,Nrphin,Nsxin,Nsyin)
84        _RL dvphy(idim1:idim2,jdim1:jdim2,Nrphin,Nsxin,Nsyin)        _RL dvphy(idim1:idim2,jdim1:jdim2,Nrphin,Nsxin,Nsyin)
85        _RL dthphy(idim1:idim2,jdim1:jdim2,Nrphin,Nsxin,Nsyin)        _RL dthphy(idim1:idim2,jdim1:jdim2,Nrphin,Nsxin,Nsyin)
# Line 70  c Local Variables Line 89  c Local Variables
89        integer ptracer,ntracer        integer ptracer,ntracer
90        parameter (ptracer = 1)        parameter (ptracer = 1)
91        parameter (ntracer = 1)        parameter (ntracer = 1)
       integer iras,nlwcld,nlwlz,nswcld,nswlz  
       integer imstturbsw,imstturblw  
92    
93        _RL xlats(sNx,sNy),xlons(sNx,sNy),sea_ice(sNx,sNy)        _RL xlats(sNx,sNy),xlons(sNx,sNy),sea_ice(sNx,sNy)
94        _RL p(sNx,sNy,Nsx,Nsy)        _RL p(sNx,sNy,Nsx,Nsy)
# Line 100  c Local Variables Line 117  c Local Variables
117        _RL tg0(sNx,sNy,Nsx,Nsy),radlwg(sNx,sNy,Nsx,Nsy)        _RL tg0(sNx,sNy,Nsx,Nsy),radlwg(sNx,sNy,Nsx,Nsy)
118        _RL lwgclr(sNx,sNy,Nsx,Nsy),st4(sNx,sNy,Nsx,Nsy)        _RL lwgclr(sNx,sNy,Nsx,Nsy),st4(sNx,sNy,Nsx,Nsy)
119        _RL dst4(sNx,sNy,Nsx,Nsy),dlwdtg(sNx,sNy,Nrphys,Nsx,Nsy)        _RL dst4(sNx,sNy,Nsx,Nsy),dlwdtg(sNx,sNy,Nrphys,Nsx,Nsy)
120          _RL qq(sNx,sNy,Nrphys,Nsx,Nsy)
121        _RL rainlsp(sNx,sNy,Nsx,Nsy),raincon(sNx,sNy,Nsx,Nsy)        _RL rainlsp(sNx,sNy,Nsx,Nsy),raincon(sNx,sNy,Nsx,Nsy)
122        _RL snowfall(sNx,sNy,Nsx,Nsy)        _RL snowfall(sNx,sNy,Nsx,Nsy)
       _RL cldtot_lw(sNx,sNy,Nrphys,Nsx,Nsy)  
       _RL cldras_lw(sNx,sNy,Nrphys,Nsx,Nsy)  
       _RL cldlsp_lw(sNx,sNy,Nrphys,Nsx,Nsy)  
       _RL lwlz(sNx,sNy,Nrphys,Nsx,Nsy)  
       _RL cldtot_sw(sNx,sNy,Nrphys,Nsx,Nsy)  
       _RL cldras_sw(sNx,sNy,Nrphys,Nsx,Nsy)  
       _RL cldlsp_sw(sNx,sNy,Nrphys,Nsx,Nsy)  
       _RL swlz(sNx,sNy,Nrphys,Nsx,Nsy)  
       _RL qliqavesw(sNx,sNy,Nrphys,Nsx,Nsy)  
       _RL qliqavelw(sNx,sNy,Nrphys,Nsx,Nsy)  
       _RL fccavesw(sNx,sNy,Nrphys,Nsx,Nsy)  
       _RL fccavelw(sNx,sNy,Nrphys,Nsx,Nsy)  
       _RL qq(sNx,sNy,Nrphys,Nsx,Nsy)  
123    
124        integer i,j,L        integer i,j,L
125        _RL getcon, kappa, p0kappa, s0, ra        _RL getcon, kappa, p0kappa, s0, ra
# Line 129  c Local Variables Line 134  c Local Variables
134        common /saver/ moistu,moistv,moistt,moistq        common /saver/ moistu,moistv,moistt,moistq
135        common /saver/ radswg,swgclr,fdirpar,fdifpar,osr,osrclr,tg0,radlwg        common /saver/ radswg,swgclr,fdirpar,fdifpar,osr,osrclr,tg0,radlwg
136        common /saver/ st4,dst4,dlwdtg,rainlsp,raincon,snowfall        common /saver/ st4,dst4,dlwdtg,rainlsp,raincon,snowfall
       common /saver/ cldtot_lw,cldras_lw,cldlsp_lw,lwlz  
       common /saver/ cldtot_sw,cldras_sw,cldlsp_sw,swlz  
       common /saver/ imstturbsw,imstturblw,qliqavesw,qliqavelw,fccavesw  
       common /saver/ fccavelw  
137        common /saver/ qq        common /saver/ qq
138        common /saver/ pl,ple,dpres,pkle,pkl        common /saver/ pl,ple,dpres,pkle,pkl
       common /saver/ nlwcld,nlwlz  
       common /saver/ nswcld,nswlz  
       common /saver/ iras  
139    
140  C***********************************************************************  C***********************************************************************
141  C Unshadow input arrays (and make 'fizhi theta' from true theta)  C Unshadow input arrays (and make 'fizhi theta' from true theta)
# Line 145  C*************************************** Line 143  C***************************************
143    
144        if( (nhms.eq.nhms0) .and. (nymd.eq.nymd0) ) then        if( (nhms.eq.nhms0) .and. (nymd.eq.nymd0) ) then
145         _BEGIN_MASTER(myid)         _BEGIN_MASTER(myid)
146         if(bi.eq.1) print *,' Initializing fizhi arrays '         if(myid.eq.1) print *,' Initializing fizhi arrays '
147         _END_MASTER(myid)         _END_MASTER(myid)
148         imstturblw = 0         imstturblw = 0
149         imstturbsw = 0         imstturbsw = 0
# Line 189  C*************************************** Line 187  C***************************************
187          snowfall(i,j,bi,bj) = 0.          snowfall(i,j,bi,bj) = 0.
188         enddo         enddo
189         enddo         enddo
190    
191        endif        endif
192    
193        kappa = getcon('KAPPA')        kappa = getcon('KAPPA')
Legend:
Removed from v.1.39  
changed lines
  Added in v.1.40
  ViewVC Help
Powered by ViewVC 1.1.22