/[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.24 by molod, Thu Aug 26 14:30:20 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 71  c ---------------- Line 72  c ----------------
72        _RL lwdt(im,jm,lm),lwdtclr(im,jm,lm)        _RL lwdt(im,jm,lm),lwdtclr(im,jm,lm)
73        _RL swdt(im,jm,lm),swdtclr(im,jm,lm)        _RL swdt(im,jm,lm),swdtclr(im,jm,lm)
74        _RL turbu(im,jm,lm),turbv(im,jm,lm),turbt(im,jm,lm)        _RL turbu(im,jm,lm),turbv(im,jm,lm),turbt(im,jm,lm)
75        _RL turbq(im,jm,lm)        _RL turbq(im,jm,lm,ntracer)
76        _RL moistu(im,jm,lm),moistv(im,jm,lm),moistt(im,jm,lm)        _RL moistu(im,jm,lm),moistv(im,jm,lm),moistt(im,jm,lm)
77        _RL moistq(im,jm,lm)        _RL moistq(im,jm,lm,ntracer)
78        _RL chemo3(im,jm,lm),chemq(im,jm,lm)        _RL chemo3(im,jm,lm),chemq(im,jm,lm)
79        _RL albvisdr(im,jm),albvisdf(im,jm),albnirdr(im,jm)        _RL albvisdr(im,jm),albvisdf(im,jm),albnirdr(im,jm)
80        _RL albnirdf(im,jm),emiss(im,jm)        _RL albnirdf(im,jm),emiss(im,jm,10)
81        _RL alai(nchp),agrn(nchp)        _RL alai(nchp),agrn(nchp)
82        _RL radswg(im,jm),swgclr(im,jm)        _RL radswg(im,jm),swgclr(im,jm)
83        _RL fdirpar(im,jm),fdifpar(im,jm),osr(im,jm),osrclr(im,jm)        _RL fdirpar(im,jm),fdifpar(im,jm),osr(im,jm),osrclr(im,jm)
# 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,
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
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.15  
changed lines
  Added in v.1.24
  ViewVC Help
Powered by ViewVC 1.1.22