/[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.34 by molod, Tue Dec 14 19:56:44 2004 UTC revision 1.41 by molod, Thu Jun 29 18:17:26 2006 UTC
# Line 107  c --------------- Line 107  c ---------------
107        integer pchecklevs(numpcheck)        integer pchecklevs(numpcheck)
108        _RL pcheckpress(numpcheck)        _RL pcheckpress(numpcheck)
109  C     data pcheckpress/950.,750.,700.,400.,10./  C     data pcheckpress/950.,750.,700.,400.,10./
110        data pcheckpress/950.,900.,700.,400.,150./        data pcheckpress/950.,850.,700.,400.,10./
111    
112        integer low_level,mid_level,nltop,nsubmin,nsubmax,Lup        integer low_level,mid_level,nltop,nsubmin,nsubmax,Lup
113        integer ndmoist,ndturb,ndlw,ndsw        integer ndmoist,ndturb,ndlw,ndsw
# Line 117  C     data pcheckpress/950.,750.,700.,40 Line 117  C     data pcheckpress/950.,750.,700.,40
117        _RL akap,getcon        _RL akap,getcon
118        _RL ptop        _RL ptop
119        logical lpnt,cumfric        logical lpnt,cumfric
120          integer imglobal
121    
122        istrip = im*jm        istrip = im*jm
123        npcs   = 1        npcs   = 1
# Line 146  C Fill array of model levels closest to Line 147  C Fill array of model levels closest to
147        nsubmax = pchecklevs(2)        nsubmax = pchecklevs(2)
148        Lup = pchecklevs(3)        Lup = pchecklevs(3)
149    
150    CCC   cumfric = .true.
151        cumfric = .false.        cumfric = .false.
152    
153  C **********************************************************************          C **********************************************************************        
# Line 179  c -------------------------------- Line 181  c --------------------------------
181        enddo        enddo
182        enddo        enddo
183    
184        call lwrio ( nymd,nhms,bi,bj,istrip,npcs,        call lwrio ( nymd,nhms,bi,bj,myid,istrip,npcs,
185       .   low_level,mid_level,im,jm,lm,       .   low_level,mid_level,im,jm,lm,
186       .   p,pl,ple,dpres,pkht,pkl,t,chemq,chemo3,co2,       .   p,pl,ple,dpres,pkht,pkl,t,chemq,chemo3,co2,
187       .   cfc11,cfc12,cfc22,methane,n2o,emiss,tgz,radlwg,st4,dst4,       .   cfc11,cfc12,cfc22,methane,n2o,emiss,tgz,radlwg,st4,dst4,
# Line 194  C TURBULENCE Line 196  C TURBULENCE
196  C ==========  C ==========
197        IF ( alarm('turb') ) THEN        IF ( alarm('turb') ) THEN
198    
199        call turbio (im,jm,lm,istrip,nymd,nhms,bi,bj,ndturb,        call turbio (im,jm,lm,istrip,nymd,nhms,bi,bj,ndturb,nltop,
200       .  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,
201       .  ctmt,xxmt,yymt,zetamt,xlmt,khmt,tke,       .  ctmt,xxmt,yymt,zetamt,xlmt,khmt,tke,
202       .  tgz,fracland,landtype,       .  tgz,fracland,landtype,
# Line 208  C ========== Line 210  C ==========
210    
211  c Add Gravity-Wave Drag Tendency  c Add Gravity-Wave Drag Tendency
212  c ------------------------------  c ------------------------------
213  C Comment this out for now - run tropospheric resolutions only.  C Comment this out for now
 C Need to see how to get zonal wavelengths on the cube.  
214    
215  c     call gwdrag (myid,p,pl,ple,dpres,pkht,pkl,u,v,t,q,phis_var,        imglobal = 128
216  c    .   turbu,turbv,turbt,im,jm,lm,ptop,istrip,npcs,        call gwdrag (myid,p,pl,ple,dpres,pkl,u,v,t,q,phis_var,
217  c    .   imglobal)       .   turbu,turbv,turbt,im,jm,lm,bi,bj,istrip,npcs,imglobal)
218    
219  c Add Rayleigh Friction Damping Above 70 Km  c Add Rayleigh Friction Damping Above 70 Km
220  c -----------------------------------------  c -----------------------------------------
Legend:
Removed from v.1.34  
changed lines
  Added in v.1.41
  ViewVC Help
Powered by ViewVC 1.1.22