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

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

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

revision 1.6 by molod, Mon Jul 26 18:45:17 2004 UTC revision 1.11 by molod, Sun Aug 29 19:43:43 2004 UTC
# Line 18  C*************************************** Line 18  C***************************************
18    
19        integer myThid,im1,im2,jm1,jm2,Nrphys,Nbi,Nbj,bi,bj,ntracer        integer myThid,im1,im2,jm1,jm2,Nrphys,Nbi,Nbj,bi,bj,ntracer
20        _RL p(im2,jm2,Nbi,Nbj)        _RL p(im2,jm2,Nbi,Nbj)
21        _RL uphy(im2,jm2,Nrphys,Nbi,Nbj),vphy(im2,jm2,Nrphys,Nbi,Nbj)        _RL uphy(im2,jm2,Nrphys,Nbi,Nbj)
22        _RL thphy(im2,jm2,Nrphys,Nbi,Nbj),sphy(im2,jm2,Nrphys,Nbi,Nbj)        _RL vphy(im2,jm2,Nrphys,Nbi,Nbj)
23          _RL thphy(im2,jm2,Nrphys,Nbi,Nbj)
24          _RL sphy(im2,jm2,Nrphys,Nbi,Nbj)
25        _RL qq(im2,jm2,Nrphys),pk(im2,jm2,Nrphys,Nbi,Nbj)        _RL qq(im2,jm2,Nrphys),pk(im2,jm2,Nrphys,Nbi,Nbj)
26        _RL dp(im2,jm2,Nrphys,Nbi,Nbj)        _RL dp(im2,jm2,Nrphys,Nbi,Nbj)
27        _RL radswt(im2,jm2,Nbi,Nbj),radswg(im2,jm2,Nbi,Nbj)        _RL radswt(im2,jm2,Nbi,Nbj),radswg(im2,jm2,Nbi,Nbj)
# Line 28  C*************************************** Line 30  C***************************************
30        _RL dst4(im2,jm2,Nbi,Nbj),tgz(im2,jm2,Nbi,Nbj)        _RL dst4(im2,jm2,Nbi,Nbj),tgz(im2,jm2,Nbi,Nbj)
31        _RL tg0(im2,jm2,Nbi,Nbj),radlwg(im2,jm2,Nbi,Nbj)        _RL tg0(im2,jm2,Nbi,Nbj),radlwg(im2,jm2,Nbi,Nbj)
32        _RL lwgclr(im2,jm2,Nbi,Nbj)        _RL lwgclr(im2,jm2,Nbi,Nbj)
33        _RL turbu(im2,jm2,Nrphys,Nbi,Nbj),turbv(im2,jm2,Nrphys,Nbi,Nbj)        _RL turbu(im2,jm2,Nrphys,Nbi,Nbj)
34          _RL turbv(im2,jm2,Nrphys,Nbi,Nbj)
35        _RL turbt(im2,jm2,Nrphys,Nbi,Nbj)        _RL turbt(im2,jm2,Nrphys,Nbi,Nbj)
36        _RL turbq(im2,jm2,Nrphys,ntracer,Nbi,Nbj)        _RL turbq(im2,jm2,Nrphys,ntracer,Nbi,Nbj)
37        _RL moistu(im2,jm2,Nrphys,Nbi,Nbj),moistv(im2,jm2,Nrphys,Nbi,Nbj)        _RL moistu(im2,jm2,Nrphys,Nbi,Nbj)
38          _RL moistv(im2,jm2,Nrphys,Nbi,Nbj)
39        _RL moistt(im2,jm2,Nrphys,Nbi,Nbj)        _RL moistt(im2,jm2,Nrphys,Nbi,Nbj)
40        _RL moistq(im2,jm2,Nrphys,ntracer,Nbi,Nbj)        _RL moistq(im2,jm2,Nrphys,ntracer,Nbi,Nbj)
41        _RL lwdt(im2,jm2,Nrphys,Nbi,Nbj),swdt(im2,jm2,Nrphys,Nbi,Nbj)        _RL lwdt(im2,jm2,Nrphys,Nbi,Nbj)
42          _RL swdt(im2,jm2,Nrphys,Nbi,Nbj)
43        _RL lwdtclr(im2,jm2,Nrphys,Nbi,Nbj)        _RL lwdtclr(im2,jm2,Nrphys,Nbi,Nbj)
44        _RL swdtclr(im2,jm2,Nrphys,Nbi,Nbj)        _RL swdtclr(im2,jm2,Nrphys,Nbi,Nbj)
45        _RL dlwdtg(im2,jm2,Nrphys,Nbi,Nbj)        _RL dlwdtg(im2,jm2,Nrphys,Nbi,Nbj)
# Line 44  C*************************************** Line 49  C***************************************
49    
50  C **********************************************************************          C **********************************************************************        
51    
52    #ifdef ALLOW_DIAGNOSTICS
53        do j=jm1,jm2        do j=jm1,jm2
54        do i=im1,im2        do i=im1,im2
55        pinv(i,j) = 1.0 / p(i,j,bi,bj)        pinv(i,j) = 1.0 / p(i,j,bi,bj)
# Line 140  c -------------------------------------- Line 146  c --------------------------------------
146        enddo        enddo
147        endif        endif
148                                                                                                                                                                    
149          if( (bi.eq.1) .and. (bj.eq.1) ) then
150        nradswt = nradswt + 1        nradswt = nradswt + 1
151        nradswg = nradswg + 1        nradswg = nradswg + 1
152        nswgclr = nswgclr + 1        nswgclr = nswgclr + 1
# Line 148  c -------------------------------------- Line 155  c --------------------------------------
155        nradlwg = nradlwg + 1        nradlwg = nradlwg + 1
156        nlwgclr = nlwgclr + 1        nlwgclr = nlwgclr + 1
157        nlwgup  = nlwgup  + 1        nlwgup  = nlwgup  + 1
158          endif
159                                                                                                                                                                    
160  C **********************************************************************          C **********************************************************************        
161        do L=1,Nrphys        do L=1,Nrphys
# Line 213  c -------------------------- Line 221  c --------------------------
221        enddo        enddo
222        enddo        enddo
223        endif        endif
224                                                                                    
225  c Longwave Heating Clear-Sky (deg/day)  c Longwave Heating Clear-Sky (deg/day)
226  c ------------------------------------  c ------------------------------------
227        if (ilwclr.ne.0) then                                                            if (ilwclr.ne.0) then                                                    
# Line 307  c ----------------------- Line 315  c -----------------------
315                                                                                                                                                                    
316        enddo        enddo
317    
318          if( (bi.eq.1) .and. (bj.eq.1) ) then
319    
320        ndiabu = ndiabu + 1        ndiabu = ndiabu + 1
321        ndiabv = ndiabv + 1        ndiabv = ndiabv + 1
322        ndiabt = ndiabt + 1        ndiabt = ndiabt + 1
# Line 321  c ----------------------- Line 331  c -----------------------
331        ntke   = ntke   + 1        ntke   = ntke   + 1
332        nsphu  = nsphu  + 1        nsphu  = nsphu  + 1
333    
334          endif
335    
336  C **********************************************************************          C **********************************************************************        
337    
338  c Vertically Averaged Moist-T Increment (K/day)  c Vertically Averaged Moist-T Increment (K/day)
# Line 420  c -------------------------------------- Line 432  c --------------------------------------
432        enddo        enddo
433        endif        endif
434    
435          if( (bi.eq.1) .and. (bj.eq.1) ) then
436        nvdtmoist = nvdtmoist + 1        nvdtmoist = nvdtmoist + 1
437        nvdtturb  = nvdtturb  + 1        nvdtturb  = nvdtturb  + 1
438        nvdtradlw = nvdtradlw + 1        nvdtradlw = nvdtradlw + 1
439        nvdtradsw = nvdtradsw + 1        nvdtradsw = nvdtradsw + 1
440          endif
441    
442    #endif
443        return        return
444        end        end
Legend:
Removed from v.1.6  
changed lines
  Added in v.1.11
  ViewVC Help
Powered by ViewVC 1.1.22