/[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.15 by molod, Tue Jul 27 04:10:56 2004 UTC revision 1.19 by molod, Thu Jul 29 15:21:50 2004 UTC
# Line 60  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 108  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          _RL tempij(im,jm)
119    
120        istrip = im*jm        istrip = im*jm
121        npcs   = 1        npcs   = 1
122    
# Line 145  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 )
164    
165        ENDIF        ENDIF
166                                                                                                                                                                    
   
167  C LONG WAVE RADIATION  C LONG WAVE RADIATION
168  C ===================  C ===================
169        IF ( alarm('radlw') ) THEN        IF ( alarm('radlw') ) THEN
# Line 175  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,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    
       print *,' Done with radiation '  
       stop  
                                                                                   
190  C TURBULENCE  C TURBULENCE
191  C ==========  C ==========
192        IF ( alarm('turb') ) THEN        IF ( alarm('turb') ) THEN
Legend:
Removed from v.1.15  
changed lines
  Added in v.1.19
  ViewVC Help
Powered by ViewVC 1.1.22