/[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.22 by molod, Tue Aug 10 15:13:47 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
193    
194        call turbio (im,jm,lm,1,nymd,nhms,bi,bj,ndturb,        if(2.eq.1)then
195          print *,' In fizhi driver, before turb - bi = ',bi
196          do j = 1,jm
197          do i = 1,im
198           tempij(i,j) = radswt(i,j)
199          enddo
200          enddo
201          print *,' radswt ',tempij
202          do j = 1,jm
203          do i = 1,im
204           tempij(i,j) = radswg(i,j)
205          enddo
206          enddo
207          print *,' radswg ',tempij
208          do j = 1,jm
209          do i = 1,im
210           tempij(i,j) = st4(i,j)
211          enddo
212          enddo
213          print *,' st4 ',tempij
214          do j = 1,jm
215          do i = 1,im
216           tempij(i,j) = dst4(i,j)
217          enddo
218          enddo
219          print *,' dst4 ',tempij
220          do j = 1,jm
221          do i = 1,im
222           tempij(i,j) = radlwg(i,j)
223          enddo
224          enddo
225          print *,' radlwg ',tempij
226          do j = 1,jm
227          do i = 1,im
228           tempij(i,j) = fdifpar(i,j)
229          enddo
230          enddo
231          print *,' fdifpar ',tempij
232          do j = 1,jm
233          do i = 1,im
234           tempij(i,j) = fdirpar(i,j)
235          enddo
236          enddo
237          print *,' fdirpar ',tempij
238          do j = 1,jm
239          do i = 1,im
240           tempij(i,j) = tgz(i,j)
241          enddo
242          enddo
243          print *,' tgz ',tempij
244          do j = 1,jm
245          do i = 1,im
246           tempij(i,j) = tg0(i,j)
247          enddo
248          enddo
249          print *,' tg0 ',tempij
250          endif
251    
252          call turbio (im,jm,lm,istrip,nymd,nhms,bi,bj,ndturb,
253       .  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,
254       .  ctmt,xxmt,yymt,zetamt,xlmt,khmt,tke,       .  ctmt,xxmt,yymt,zetamt,xlmt,khmt,tke,
255       .  tgz,fracland,landtype,       .  tgz,fracland,landtype,
Legend:
Removed from v.1.13  
changed lines
  Added in v.1.22
  ViewVC Help
Powered by ViewVC 1.1.22