/[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.9 by molod, Tue Aug 10 15:13:47 2004 UTC revision 1.13 by molod, Tue Dec 14 19:56:45 2004 UTC
# Line 12  C*************************************** Line 12  C***************************************
12    
13  #ifdef ALLOW_DIAGNOSTICS  #ifdef ALLOW_DIAGNOSTICS
14  #include "SIZE.h"  #include "SIZE.h"
15  #include "diagnostics_SIZE.h"  #include "DIAGNOSTICS_SIZE.h"
16  #include "diagnostics.h"  #include "DIAGNOSTICS.h"
17  #endif  #endif
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
# Line 221  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 240  c --------------------------------- Line 240  c ---------------------------------
240        if (iradsw.ne.0) then        if (iradsw.ne.0) then
241        do j=jm1,jm2        do j=jm1,jm2
242        do i=im1,im2        do i=im1,im2
243        qdiag(i,j,iradsw+l-1,bi,bj) = qdiag(i,j,iradsw+l-1,bi,bj) +        qdiag(i,j,iradsw+l-1,bi,bj) = qdiag(i,j,iradsw+l-1,bi,bj)
244       .  + swdt(i,j,l,bi,bj)*radswt(i,j,bi,bj)*       .  + swdt(i,j,l,bi,bj)*radswt(i,j,bi,bj)*
245       .                   pk(i,j,l,bi,bj)*pinv(i,j)*86400       .                   pk(i,j,l,bi,bj)*pinv(i,j)*86400
246        enddo        enddo
# Line 316  c ----------------------- Line 316  c -----------------------
316        enddo        enddo
317    
318        if( (bi.eq.1) .and. (bj.eq.1) ) then        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 329  c ----------------------- Line 330  c -----------------------
330        ntmpu  = ntmpu  + 1                                  ntmpu  = ntmpu  + 1                          
331        ntke   = ntke   + 1        ntke   = ntke   + 1
332        nsphu  = nsphu  + 1        nsphu  = nsphu  + 1
333    
334        endif        endif
335    
336  C **********************************************************************          C **********************************************************************        
Legend:
Removed from v.1.9  
changed lines
  Added in v.1.13
  ViewVC Help
Powered by ViewVC 1.1.22