/[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.13 by molod, Mon Jul 26 18:45:17 2004 UTC revision 1.19 by molod, Thu Jul 29 15:21:50 2004 UTC
# Line 37  C*************************************** Line 37  C***************************************
37  c Diagnostic Common  c Diagnostic Common
38  c -----------------  c -----------------
39  #ifdef ALLOW_DIAGNOSTICS  #ifdef ALLOW_DIAGNOSTICS
40  #include 'diagnostics.h'  #include "SIZE.h"
41    #include "diagnostics_SIZE.h"
42    #include "diagnostics.h"
43  #endif  #endif
44    
45  c Timers Common  c Timers Common
# Line 58  c ---------------- Line 60  c ----------------
60        _RL pl(im,jm,lm),ple(im,jm,lm+1),dpres(im,jm,lm)        _RL pl(im,jm,lm),ple(im,jm,lm+1),dpres(im,jm,lm)
61        _RL pkht(im,jm,lm)        _RL pkht(im,jm,lm)
62        _RL pkl(im,jm,lm)        _RL pkl(im,jm,lm)
63        _RL radswt(im,jm), fracland(im,jm), landtype(im,jm)        _RL radswt(im,jm), fracland(im,jm)
64          integer landtype(im,jm)
65        _RL phis_var(im,jm), sea_ice(im,jm)        _RL phis_var(im,jm), sea_ice(im,jm)
66        _RL chlat(nchp),chlon(nchp),chfr(nchp)        _RL chlat(nchp),chlon(nchp),chfr(nchp)
67        _RL tcanopy(nchp),tdeep(nchp),ecanopy(nchp),swetshal(nchp)        _RL tcanopy(nchp),tdeep(nchp),ecanopy(nchp),swetshal(nchp)
# Line 106  c --------------- Line 109  c ---------------
109        integer low_level,mid_level,nltop,nsubmin,nsubmax,Lup        integer low_level,mid_level,nltop,nsubmin,nsubmax,Lup
110        integer ndmoist,ndturb,ndlw,ndsw        integer ndmoist,ndturb,ndlw,ndsw
111        integer istrip,npcs        integer istrip,npcs
112        integer i,j        integer i,j,L
113        integer ndum,ndpnt        integer ndum,ndpnt
114        _RL akap,getcon        _RL akap,getcon
115        _RL ptop        _RL ptop
116        logical lpnt        logical lpnt
117    
118        istrip = min( im*jm,71 )        _RL tempij(im,jm)
119        npcs   = ((im*jm-1)/istrip) + 1  
120          istrip = im*jm
121          npcs   = 1
122    
123        ptop   =  0.        ptop   =  0.
124        akap = getcon('KAPPA')        akap = getcon('KAPPA')
# Line 132  C ************************************** Line 137  C **************************************
137  C Fill array of model levels closest to a given pressure value  C Fill array of model levels closest to a given pressure value
138    
139        call getpwhere(myid,numpcheck,pcheckpress,pchecklevs)        call getpwhere(myid,numpcheck,pcheckpress,pchecklevs)
140    
141        low_level = pchecklevs(3)        low_level = pchecklevs(3)
142        mid_level = pchecklevs(4)        mid_level = pchecklevs(4)
143        nltop = pchecklevs(5)        nltop = pchecklevs(5)
# Line 142  C Fill array of model levels closest to Line 148  C Fill array of model levels closest to
148  C **********************************************************************          C **********************************************************************        
149  C ****                  Call Physics Mini-Drivers                   ****          C ****                  Call Physics Mini-Drivers                   ****        
150  C **********************************************************************          C **********************************************************************        
151                                                                                    
152  C SHORT WAVE RADIATION  C SHORT WAVE RADIATION
153  C ====================  C ====================
154        IF ( alarm('radsw') ) THEN        IF ( alarm('radsw') ) THEN
155    
156        call swrio ( nymd,nhms,bi,bj,ndsw,myid,istrip,npcs,        call swrio ( nymd,nhms,bi,bj,ndsw,myid,istrip,npcs,
157       .   low_level,mid_level,       .   low_level,mid_level,im,jm,lm,
158       .   p,pl,ple,dpres,pkht,pkl,t,chemq,chemo3,co2,       .   p,pl,ple,dpres,pkht,pkl,t,chemq,chemo3,co2,
159       .   albvisdr,albvisdf,albnirdr,albnirdf,swdt,swdtclr,       .   albvisdr,albvisdf,albnirdr,albnirdf,swdt,swdtclr,
160       .   radswg,swgclr,fdifpar,fdirpar,osr,osrclr,im,jm,lm,       .   radswg,swgclr,fdifpar,fdirpar,osr,osrclr,
161       .   ptop,nswcld,cldtot_sw,cldras_sw,nswlz,swlz,       .   ptop,nswcld,cldtot_sw,cldras_sw,nswlz,swlz,
162       .   .false.,imstturbsw,qliqavesw,       .   .false.,imstturbsw,qliqavesw,
163       .   fccavesw,landtype,xlats,xlons )       .   fccavesw,landtype,xlats,xlons )
# Line 171  c -------------------------------- Line 177  c --------------------------------
177        enddo        enddo
178    
179        call lwrio ( nymd,nhms,bi,bj,istrip,npcs,        call lwrio ( nymd,nhms,bi,bj,istrip,npcs,
180       .   low_level,mid_level,       .   low_level,mid_level,im,jm,lm,
181       .   p,pl,ple,dpres,pkht,pkl,t,chemq,pkht,chemo3,co2,       .   p,pl,ple,dpres,pkht,pkl,t,chemq,chemo3,co2,
182       .   cfc11,cfc12,cfc22,methane,n2o,emiss,tgz,radlwg,st4,dst4,       .   cfc11,cfc12,cfc22,methane,n2o,emiss,tgz,radlwg,st4,dst4,
183       .   lwdt,dlwdtg,lwdtclr,lwgclr,im,jm,lm,ptop,       .   lwdt,dlwdtg,lwdtclr,lwgclr,ptop,
184       .   nlwcld,cldtot_lw,cldras_lw,nlwlz,lwlz,       .   nlwcld,cldtot_lw,cldras_lw,nlwlz,lwlz,
185       .   .false.,imstturblw,qliqavelw,       .   .false.,imstturblw,qliqavelw,
186       .   fccavelw,landtype )       .   fccavelw,landtype )
187    
188        ENDIF        ENDIF
189                                                                                    
190  C TURBULENCE  C TURBULENCE
191  C ==========  C ==========
192        IF ( alarm('turb') ) THEN        IF ( alarm('turb') ) THEN
Legend:
Removed from v.1.13  
changed lines
  Added in v.1.19
  ViewVC Help
Powered by ViewVC 1.1.22