/[MITgcm]/MITgcm/pkg/fizhi/do_fizhi.F
ViewVC logotype

Diff of /MITgcm/pkg/fizhi/do_fizhi.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.29 by molod, Thu Jul 29 15:21:50 2004 UTC revision 1.41 by molod, Thu Oct 7 01:20:47 2004 UTC
# Line 10  C $Name$ Line 10  C $Name$
10       . albvisdr,albvisdf,ityp,chfr,alai,agrn,igrd,chlat,chlon,       . albvisdr,albvisdf,ityp,chfr,alai,agrn,igrd,chlat,chlon,
11       . tcanopy,tdeep,ecanopy,swetshal,swetroot,swetdeep,snodep,capac,       . tcanopy,tdeep,ecanopy,swetshal,swetroot,swetdeep,snodep,capac,
12       . o3,qstr,co2,cfc11,cfc12,cfc22,n2o,methane,       . o3,qstr,co2,cfc11,cfc12,cfc22,n2o,methane,
13         . iras,nlwcld,cldtot_lw,cldras_lw,cldlsp_lw,nlwlz,lwlz,
14         . nswcld,cldtot_sw,cldras_sw,cldlsp_sw,nswlz,swlz,
15         . imstturbsw,imstturblw,qliqavesw,qliqavelw,fccavesw,fccavelw,
16         . raincon,rainlsp,snowfall,
17       . duphy,dvphy,dthphy,dsphy)       . duphy,dvphy,dthphy,dsphy)
18  c-----------------------------------------------------------------------  c-----------------------------------------------------------------------
19  c Interface routine to calculate physics increments - calls fizhi_driver.  c Interface routine to calculate physics increments - calls fizhi_driver.
# Line 28  c--------------------------------------- Line 32  c---------------------------------------
32    
33  C Argument list declarations  C Argument list declarations
34        integer myid,im1,im2,jm1,jm2,idim1,idim2,jdim1,jdim2        integer myid,im1,im2,jm1,jm2,idim1,idim2,jdim1,jdim2
35        integer Nrphin,Nsxin,Nsyin,bi,bj,nchp,nchptot,nchpland        integer Nrphin,Nsxin,Nsyin,bi,bj,nchp
36          integer nchptot(Nsxin,Nsyin),nchpland(Nsxin,Nsyin)
37        _RL uphy(idim1:idim2,jdim1:jdim2,Nrphin,Nsxin,Nsyin)        _RL uphy(idim1:idim2,jdim1:jdim2,Nrphin,Nsxin,Nsyin)
38        _RL vphy(idim1:idim2,jdim1:jdim2,Nrphin,Nsxin,Nsyin)        _RL vphy(idim1:idim2,jdim1:jdim2,Nrphin,Nsxin,Nsyin)
39        _RL thphy(idim1:idim2,jdim1:jdim2,Nrphin,Nsxin,Nsyin)        _RL thphy(idim1:idim2,jdim1:jdim2,Nrphin,Nsxin,Nsyin)
# Line 60  C Argument list declarations Line 65  C Argument list declarations
65        _RL o3(im2,jm2,Nrphin,Nsxin,Nsyin)        _RL o3(im2,jm2,Nrphin,Nsxin,Nsyin)
66        _RL qstr(im2,jm2,Nrphin,Nsxin,Nsyin)        _RL qstr(im2,jm2,Nrphin,Nsxin,Nsyin)
67        _RL co2,cfc11,cfc12,cfc22,n2o(Nrphin),methane(Nrphin)        _RL co2,cfc11,cfc12,cfc22,n2o(Nrphin),methane(Nrphin)
68    
69          integer iras,nlwcld,nlwlz,nswcld,nswlz
70          integer imstturbsw,imstturblw
71          _RL cldtot_lw(idim1:idim2,jdim1:jdim2,Nrphin,Nsxin,Nsyin)
72          _RL cldras_lw(idim1:idim2,jdim1:jdim2,Nrphin,Nsxin,Nsyin)
73          _RL cldlsp_lw(idim1:idim2,jdim1:jdim2,Nrphin,Nsxin,Nsyin)
74          _RL lwlz(idim1:idim2,jdim1:jdim2,Nrphin,Nsxin,Nsyin)
75          _RL cldtot_sw(idim1:idim2,jdim1:jdim2,Nrphin,Nsxin,Nsyin)
76          _RL cldras_sw(idim1:idim2,jdim1:jdim2,Nrphin,Nsxin,Nsyin)
77          _RL cldlsp_sw(idim1:idim2,jdim1:jdim2,Nrphin,Nsxin,Nsyin)
78          _RL swlz(idim1:idim2,jdim1:jdim2,Nrphin,Nsxin,Nsyin)
79          _RL qliqavesw(idim1:idim2,jdim1:jdim2,Nrphin,Nsxin,Nsyin)
80          _RL qliqavelw(idim1:idim2,jdim1:jdim2,Nrphin,Nsxin,Nsyin)
81          _RL fccavesw(idim1:idim2,jdim1:jdim2,Nrphin,Nsxin,Nsyin)
82          _RL fccavelw(idim1:idim2,jdim1:jdim2,Nrphin,Nsxin,Nsyin)
83          _RL rainlsp(idim1:idim2,jdim1:jdim2,Nsxin,Nsyin)
84          _RL raincon(idim1:idim2,jdim1:jdim2,Nsxin,Nsyin)
85          _RL snowfall(idim1:idim2,jdim1:jdim2,Nsxin,Nsyin)
86    
87    
88        _RL duphy(idim1:idim2,jdim1:jdim2,Nrphin,Nsxin,Nsyin)        _RL duphy(idim1:idim2,jdim1:jdim2,Nrphin,Nsxin,Nsyin)
89        _RL dvphy(idim1:idim2,jdim1:jdim2,Nrphin,Nsxin,Nsyin)        _RL dvphy(idim1:idim2,jdim1:jdim2,Nrphin,Nsxin,Nsyin)
90        _RL dthphy(idim1:idim2,jdim1:jdim2,Nrphin,Nsxin,Nsyin)        _RL dthphy(idim1:idim2,jdim1:jdim2,Nrphin,Nsxin,Nsyin)
# Line 69  c Local Variables Line 94  c Local Variables
94        integer ptracer,ntracer        integer ptracer,ntracer
95        parameter (ptracer = 1)        parameter (ptracer = 1)
96        parameter (ntracer = 1)        parameter (ntracer = 1)
       integer iras,nlwcld,nlwlz,nswcld,nswlz  
       integer imstturbsw,imstturblw  
97    
98        _RL xlats(sNx,sNy),xlons(sNx,sNy),sea_ice(sNx,sNy)        _RL xlats(sNx,sNy),xlons(sNx,sNy),sea_ice(sNx,sNy)
99        _RL p(sNx,sNy,Nsx,Nsy)        _RL p(sNx,sNy,Nsx,Nsy)
# Line 99  c Local Variables Line 122  c Local Variables
122        _RL tg0(sNx,sNy,Nsx,Nsy),radlwg(sNx,sNy,Nsx,Nsy)        _RL tg0(sNx,sNy,Nsx,Nsy),radlwg(sNx,sNy,Nsx,Nsy)
123        _RL lwgclr(sNx,sNy,Nsx,Nsy),st4(sNx,sNy,Nsx,Nsy)        _RL lwgclr(sNx,sNy,Nsx,Nsy),st4(sNx,sNy,Nsx,Nsy)
124        _RL dst4(sNx,sNy,Nsx,Nsy),dlwdtg(sNx,sNy,Nrphys,Nsx,Nsy)        _RL dst4(sNx,sNy,Nsx,Nsy),dlwdtg(sNx,sNy,Nrphys,Nsx,Nsy)
       _RL rainlsp(sNx,sNy,Nsx,Nsy),raincon(sNx,sNy,Nsx,Nsy)  
       _RL snowfall(sNx,sNy,Nsx,Nsy)  
       _RL cldtot_lw(sNx,sNy,Nrphys,Nsx,Nsy)  
       _RL cldras_lw(sNx,sNy,Nrphys,Nsx,Nsy)  
       _RL cldlsp_lw(sNx,sNy,Nrphys,Nsx,Nsy)  
       _RL lwlz(sNx,sNy,Nrphys,Nsx,Nsy)  
       _RL cldtot_sw(sNx,sNy,Nrphys,Nsx,Nsy)  
       _RL cldras_sw(sNx,sNy,Nrphys,Nsx,Nsy)  
       _RL cldlsp_sw(sNx,sNy,Nrphys,Nsx,Nsy)  
       _RL swlz(sNx,sNy,Nrphys,Nsx,Nsy)  
       _RL qliqavesw(sNx,sNy,Nrphys,Nsx,Nsy)  
       _RL qliqavelw(sNx,sNy,Nrphys,Nsx,Nsy)  
       _RL fccavesw(sNx,sNy,Nrphys,Nsx,Nsy)  
       _RL fccavelw(sNx,sNy,Nrphys,Nsx,Nsy)  
125        _RL qq(sNx,sNy,Nrphys,Nsx,Nsy)        _RL qq(sNx,sNy,Nrphys,Nsx,Nsy)
   
126        integer i,j,L        integer i,j,L
127        _RL getcon, kappa, p0kappa, s0, ra        _RL getcon, kappa, p0kappa, s0, ra
128        _RL cosz(sNx,sNy)        _RL cosz(sNx,sNy)
# Line 124  c Local Variables Line 132  c Local Variables
132        logical alarm        logical alarm
133        external alarm        external alarm
134    
135        save lwdt,lwdtclr,swdt,swdtclr,turbu,turbv,turbt,turbq        common /saver/ lwdt,lwdtclr,swdt,swdtclr,turbu,turbv,turbt,turbq
136        save moistu,moistv,moistt,moistq        common /saver/ moistu,moistv,moistt,moistq
137        save radswg,swgclr,fdirpar,fdifpar,osr,osrclr,tg0,radlwg        common /saver/ radswg,swgclr,fdirpar,fdifpar,osr,osrclr,tg0,radlwg
138        save st4,dst4,dlwdtg,rainlsp,raincon,snowfall,iras        common /saver/ st4,dst4,dlwdtg
139        save nlwcld,cldtot_lw,cldras_lw,cldlsp_lw,nlwlz,lwlz        common /saver/ qq
140        save nswcld,cldtot_sw,cldras_sw,cldlsp_sw,nswlz,swlz        common /saver/ pl,ple,dpres,pkle,pkl
141        save imstturbsw,imstturblw,qliqavesw,qliqavelw,fccavesw,fccavelw  
       save qq  
       save pl,ple,dpres,pkle,pkl  
   
142  C***********************************************************************  C***********************************************************************
143  C Unshadow input arrays (and make 'fizhi theta' from true theta)  C Unshadow input arrays (and make 'fizhi theta' from true theta)
144  C***********************************************************************  C***********************************************************************
145    
146        if( (nhms.eq.nhms0) .and. (nymd.eq.nymd0) ) then        if( (nhms.eq.nhms0) .and. (nymd.eq.nymd0) ) then
147         if(myid.eq.1) print *,' nhms, nhms0 ',nhms,' ',nhms0         _BEGIN_MASTER(myid)
        if(myid.eq.1) print *,' nymd, nymd0 ',nymd,' ',nymd0  
148         if(myid.eq.1) print *,' Initializing fizhi arrays '         if(myid.eq.1) print *,' Initializing fizhi arrays '
149           _END_MASTER(myid)
150         imstturblw = 0         imstturblw = 0
151         imstturbsw = 0         imstturbsw = 0
152         iras = 0         iras = 0
# Line 184  C*************************************** Line 189  C***************************************
189          snowfall(i,j,bi,bj) = 0.          snowfall(i,j,bi,bj) = 0.
190         enddo         enddo
191         enddo         enddo
192    
193        endif        endif
194    
195        kappa = getcon('KAPPA')        kappa = getcon('KAPPA')
# Line 226  C compute pressures - all pressure are c Line 232  C compute pressures - all pressure are c
232         pl(i,j,L,bi,bj) = (pephy(i,j,L,bi,bj)+pephy(i,j,L+1,bi,bj))/200.         pl(i,j,L,bi,bj) = (pephy(i,j,L,bi,bj)+pephy(i,j,L+1,bi,bj))/200.
233         dpres(i,j,L,bi,bj)=(pephy(i,j,L+1,bi,bj)-pephy(i,j,L,bi,bj))/100.         dpres(i,j,L,bi,bj)=(pephy(i,j,L+1,bi,bj)-pephy(i,j,L,bi,bj))/100.
234         ple(i,j,L,bi,bj) = pephy(i,j,L,bi,bj)/100.         ple(i,j,L,bi,bj) = pephy(i,j,L,bi,bj)/100.
235         pkle(i,j,L,bi,bj) = ple(i,j,L,bi,bj) **kappa         if (ple(i,j,L,bi,bj).gt.0.) then
236            pkle(i,j,L,bi,bj) = ple(i,j,L,bi,bj) **kappa
237           else
238            pkle(i,j,L,bi,bj) = 0.
239           endif
240        enddo        enddo
241        enddo        enddo
242        enddo        enddo
# Line 234  C compute pressures - all pressure are c Line 244  C compute pressures - all pressure are c
244        call pkappa (im2,jm2,Nrphys,ple(1,1,1,bi,bj),pkle(1,1,1,bi,bj),        call pkappa (im2,jm2,Nrphys,ple(1,1,1,bi,bj),pkle(1,1,1,bi,bj),
245       .                                                 pkl(1,1,1,bi,bj))       .                                                 pkl(1,1,1,bi,bj))
246    
247          if(2.eq.1 )then
248          print *,' In do fizhi, before fizhi driver - bi = ',bi
249          do L = 1,Nrphys+1
250           do j = jm1,jm2
251           do i = im1,im2
252            tempij(i,j) = pephy(i,j,L,bi,bj)
253           enddo
254           enddo
255    c      print *,' pephy at level ',l,' ',tempij
256          enddo
257          do L = 1,Nrphys
258           do j = jm1,jm2
259           do i = im1,im2
260            tempij(i,j) = u(i,j,L)
261           enddo
262           enddo
263    c      print *,' u at level ',l,' ',tempij
264          enddo
265          do L = 1,Nrphys
266           do j = jm1,jm2
267           do i = im1,im2
268            tempij(i,j) = v(i,j,L)
269           enddo
270           enddo
271    c      print *,' v at level ',l,' ',tempij
272          enddo
273          do L = 1,Nrphys
274           do j = jm1,jm2
275           do i = im1,im2
276            tempij(i,j) = t(i,j,L)
277           enddo
278           enddo
279           print *,' t at level ',l,' ',tempij
280          enddo
281          do L = 1,Nrphys
282           do j = jm1,jm2
283           do i = im1,im2
284            tempij(i,j) = q(i,j,L,1)
285           enddo
286           enddo
287           print *,' q at level ',l,' ',tempij
288          enddo
289          do L = 1,Nrphys
290           do j = jm1,jm2
291           do i = im1,im2
292            tempij(i,j) = qstr(i,j,L,bi,bj)
293           enddo
294           enddo
295           print *,' radiation q at level ',l,' ',tempij
296          enddo
297          do L = 1,Nrphys
298           do j = jm1,jm2
299           do i = im1,im2
300            tempij(i,j) = pl(i,j,L,bi,bj)
301           enddo
302           enddo
303    c      print *,' pl at level ',l,' ',tempij
304          enddo
305          do L = 1,Nrphys+1
306           do j = jm1,jm2
307           do i = im1,im2
308            tempij(i,j) = ple(i,j,L,bi,bj)
309           enddo
310           enddo
311    c      print *,' ple at level ',l,' ',tempij
312          enddo
313          do L = 1,Nrphys
314           do j = jm1,jm2
315           do i = im1,im2
316            tempij(i,j) = pkl(i,j,L,bi,bj)
317           enddo
318           enddo
319    c      print *,' pkl at level ',l,' ',tempij
320          enddo
321          do L = 1,Nrphys+1
322           do j = jm1,jm2
323           do i = im1,im2
324            tempij(i,j) = pkle(i,j,L,bi,bj)
325           enddo
326           enddo
327    c      print *,' pkle at level ',l,' ',tempij
328          enddo
329          do L = 1,Nrphys
330           do j = jm1,jm2
331           do i = im1,im2
332            tempij(i,j) = dpres(i,j,L,bi,bj)
333           enddo
334           enddo
335    c      print *,' dpres at level ',l,' ',tempij
336          enddo
337          do j = jm1,jm2
338          do i = im1,im2
339           tempij(i,j) = tgz(i,j,bi,bj)
340          enddo
341          enddo
342    c     print *,' tgz ',tempij
343    
344          print *,' Just before fizhi driver call '
345    
346          endif
347    
348          
349    
350        call fizhi_driver(myid,im2,jm2,Nrphys,bi,bj,ptracer,ntracer,xlats,        call fizhi_driver(myid,im2,jm2,Nrphys,bi,bj,ptracer,ntracer,xlats,
351       . xlons,p(1,1,bi,bj),u,v,t,q,pl(1,1,1,bi,bj),ple(1,1,1,bi,bj),       . xlons,p(1,1,bi,bj),u,v,t,q,pl(1,1,1,bi,bj),ple(1,1,1,bi,bj),
352       . dpres(1,1,1,bi,bj),pkle(1,1,1,bi,bj),pkl(1,1,1,bi,bj),       . dpres(1,1,1,bi,bj),pkle(1,1,1,bi,bj),pkl(1,1,1,bi,bj),
353       . fracland(1,1,bi,bj),landtype(1,1,bi,bj),radswt(1,1,bi,bj),       . fracland(1,1,bi,bj),landtype(1,1,bi,bj),radswt(1,1,bi,bj),
354       . phis_var(1,1,bi,bj),tgz(1,1,bi,bj),sea_ice,nchp,chlat(1,bi,bj),       . phis_var(1,1,bi,bj),tgz(1,1,bi,bj),sea_ice,nchp,chlat(1,bi,bj),
355       . chlon(1,bi,bj),igrd(1,bi,bj),nchptot,nchpland,chfr(1,bi,bj),       . chlon(1,bi,bj),igrd(1,bi,bj),nchptot(bi,bj),nchpland(bi,bj),
356       . ityp(1,bi,bj),tcanopy(1,bi,bj),tdeep(1,bi,bj),ecanopy(1,bi,bj),       . chfr(1,bi,bj),ityp(1,bi,bj),tcanopy(1,bi,bj),tdeep(1,bi,bj),
357       . swetshal(1,bi,bj),swetroot(1,bi,bj),swetdeep(1,bi,bj),       . ecanopy(1,bi,bj),swetshal(1,bi,bj),swetroot(1,bi,bj),
358       . capac(1,bi,bj),snodep(1,bi,bj),       . swetdeep(1,bi,bj),capac(1,bi,bj),snodep(1,bi,bj),
359       . ctmt(1,bi,bj),xxmt(1,bi,bj),yymt(1,bi,bj),zetamt(1,bi,bj),       . ctmt(1,bi,bj),xxmt(1,bi,bj),yymt(1,bi,bj),zetamt(1,bi,bj),
360       . xlmt(1,1,bi,bj),khmt(1,1,bi,bj),tke(1,1,bi,bj),       . xlmt(1,1,bi,bj),khmt(1,1,bi,bj),tke(1,1,bi,bj),
361       . albvisdr(1,1,bi,bj),albvisdf(1,1,bi,bj),albnirdr(1,1,bi,bj),       . albvisdr(1,1,bi,bj),albvisdf(1,1,bi,bj),albnirdr(1,1,bi,bj),
# Line 267  C compute pressures - all pressure are c Line 380  C compute pressures - all pressure are c
380       . qliqavelw(1,1,1,bi,bj),fccavesw(1,1,1,bi,bj),       . qliqavelw(1,1,1,bi,bj),fccavesw(1,1,1,bi,bj),
381       . fccavelw(1,1,1,bi,bj),qq(1,1,1,bi,bj))       . fccavelw(1,1,1,bi,bj),qq(1,1,1,bi,bj))
382    
383        endif        if(2.eq.1)then
384          print *,' In do fizhi, after fizhi driver - bi = ',bi
       print *,' In do fizhi, bi = ',bi  
385        do L = 1,Nrphys        do L = 1,Nrphys
386         do j = jm1,jm2         do j = jm1,jm2
387         do i = im1,im2         do i = im1,im2
388          tempij(i,j) = turbu(i,j,L,bi,bj)          tempij(i,j) = turbu(i,j,L,bi,bj)
389         enddo         enddo
390         enddo         enddo
391         print *,' turbu at level ',l,' ',tempij  c      print *,' turbu at level ',l,' ',tempij
392        enddo        enddo
393        do L = 1,Nrphys        do L = 1,Nrphys
394         do j = jm1,jm2         do j = jm1,jm2
# Line 284  C compute pressures - all pressure are c Line 396  C compute pressures - all pressure are c
396          tempij(i,j) = turbv(i,j,L,bi,bj)          tempij(i,j) = turbv(i,j,L,bi,bj)
397         enddo         enddo
398         enddo         enddo
399         print *,' turbv at level ',l,' ',tempij  c      print *,' turbv at level ',l,' ',tempij
400        enddo        enddo
401        do L = 1,Nrphys        do L = 1,Nrphys
402         do j = jm1,jm2         do j = jm1,jm2
# Line 300  C compute pressures - all pressure are c Line 412  C compute pressures - all pressure are c
412          tempij(i,j) = turbq(i,j,L,1,bi,bj)/p(i,j,bi,bj)          tempij(i,j) = turbq(i,j,L,1,bi,bj)/p(i,j,bi,bj)
413         enddo         enddo
414         enddo         enddo
415         print *,' turbq at level ',l,' ',tempij  c      print *,' turbq at level ',l,' ',tempij
416        enddo        enddo
417        do L = 1,Nrphys        do L = 1,Nrphys
418         do j = jm1,jm2         do j = jm1,jm2
# Line 308  C compute pressures - all pressure are c Line 420  C compute pressures - all pressure are c
420          tempij(i,j) = moistu(i,j,L,bi,bj)          tempij(i,j) = moistu(i,j,L,bi,bj)
421         enddo         enddo
422         enddo         enddo
423         print *,' moistu at level ',l,' ',tempij  c      print *,' moistu at level ',l,' ',tempij
424        enddo        enddo
425        do L = 1,Nrphys        do L = 1,Nrphys
426         do j = jm1,jm2         do j = jm1,jm2
# Line 316  C compute pressures - all pressure are c Line 428  C compute pressures - all pressure are c
428          tempij(i,j) = moistv(i,j,L,bi,bj)          tempij(i,j) = moistv(i,j,L,bi,bj)
429         enddo         enddo
430         enddo         enddo
431         print *,' moistv at level ',l,' ',tempij  c      print *,' moistv at level ',l,' ',tempij
432        enddo        enddo
433        do L = 1,Nrphys        do L = 1,Nrphys
434         do j = jm1,jm2         do j = jm1,jm2
# Line 332  C compute pressures - all pressure are c Line 444  C compute pressures - all pressure are c
444          tempij(i,j) = moistq(i,j,L,1,bi,bj)/p(i,j,bi,bj)          tempij(i,j) = moistq(i,j,L,1,bi,bj)/p(i,j,bi,bj)
445         enddo         enddo
446         enddo         enddo
447         print *,' moistq at level ',l,' ',tempij  c      print *,' moistq at level ',l,' ',tempij
448        enddo        enddo
449        do L = 1,Nrphys        do L = 1,Nrphys
450         do j = jm1,jm2         do j = jm1,jm2
# Line 346  C compute pressures - all pressure are c Line 458  C compute pressures - all pressure are c
458        do L = 1,Nrphys        do L = 1,Nrphys
459         do j = jm1,jm2         do j = jm1,jm2
460         do i = im1,im2         do i = im1,im2
461            tempij(i,j) = lwdt(i,j,L,bi,bj)
462           enddo
463           enddo
464           print *,' lwdt alone at level ',l,' ',tempij
465          enddo
466          do L = 1,Nrphys
467           do j = jm1,jm2
468           do i = im1,im2
469          tempij(i,j) = (lwdt(i,j,L,bi,bj)+          tempij(i,j) = (lwdt(i,j,L,bi,bj)+
470       .         dlwdtg(i,j,L,bi,bj)*(tgz(i,j,bi,bj)-tg0(i,j,bi,bj)))       .         dlwdtg(i,j,L,bi,bj)*(tgz(i,j,bi,bj)-tg0(i,j,bi,bj)))
471       .                 *p0kappa/p(i,j,bi,bj)       .                 *p0kappa/p(i,j,bi,bj)
472         enddo         enddo
473         enddo         enddo
474         print *,' lwdt at level ',l,' ',tempij         print *,' net lwdt at level ',l,' ',tempij
475        enddo        enddo
476          endif
477    
478        do L = 1,Nrphys        do L = 1,Nrphys
479        do j = jm1,jm2        do j = jm1,jm2
# Line 369  C compute pressures - all pressure are c Line 490  C compute pressures - all pressure are c
490        enddo        enddo
491        enddo        enddo
492    
493          if(2.eq.1 )then
494          print *,' In do fizhi, computed fizhi tendencies ',bi
495          do L = 1,Nrphys
496           do j = jm1,jm2
497           do i = im1,im2
498            tempij(i,j) = duphy(i,j,L,bi,bj)
499           enddo
500           enddo
501    c      print *,' duphy at level ',l,' ',tempij
502          enddo
503          do L = 1,Nrphys
504           do j = jm1,jm2
505           do i = im1,im2
506            tempij(i,j) = dvphy(i,j,L,bi,bj)
507           enddo
508           enddo
509    c      print *,' dvphy at level ',l,' ',tempij
510          enddo
511          do L = 1,Nrphys
512           do j = jm1,jm2
513           do i = im1,im2
514            tempij(i,j) = dthphy(i,j,L,bi,bj)
515           enddo
516           enddo
517           print *,' dthphy at level ',l,' ',tempij
518          enddo
519          do L = 1,Nrphys
520           do j = jm1,jm2
521           do i = im1,im2
522            tempij(i,j) = dsphy(i,j,L,bi,bj)
523           enddo
524           enddo
525           print *,' dsphy at level ',l,' ',tempij
526          enddo
527          endif
528    
529          endif
530    
531        call fizhi_step_diag(myid,p,uphy,vphy,thphy,sphy,qq,pkl,dpres,        call fizhi_step_diag(myid,p,uphy,vphy,thphy,sphy,qq,pkl,dpres,
532       .  radswt,radswg,swgclr,osr,osrclr,st4,dst4,tgz,tg0,radlwg,lwgclr,       .  radswt,radswg,swgclr,osr,osrclr,st4,dst4,tgz,tg0,radlwg,lwgclr,
533       .  turbu,turbv,turbt,turbq,moistu,moistv,moistt,moistq,       .  turbu,turbv,turbt,turbq,moistu,moistv,moistt,moistq,
Legend:
Removed from v.1.29  
changed lines
  Added in v.1.41
  ViewVC Help
Powered by ViewVC 1.1.22