/[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.17 by molod, Fri Jun 17 01:04:24 2005 UTC revision 1.18 by molod, Tue Jan 17 15:59:44 2006 UTC
# Line 39  C*************************************** Line 39  C***************************************
39        _RL dlwdtg(im2,jm2,Nrphys,Nbi,Nbj)        _RL dlwdtg(im2,jm2,Nrphys,Nbi,Nbj)
40    
41        integer  i,j,L        integer  i,j,L
42          _RL getcon, gravity
43        _RL pinv(im2,jm2), qbar(im2,jm2),tmpdiag(im2,jm2)        _RL pinv(im2,jm2), qbar(im2,jm2),tmpdiag(im2,jm2)
44  #ifdef ALLOW_DIAGNOSTICS  #ifdef ALLOW_DIAGNOSTICS
45        logical  diagnostics_is_on        logical  diagnostics_is_on
# Line 426  c -------------------------------------- Line 427  c --------------------------------------
427         call diagnostics_fill(tmpdiag,'VDTRADSW',0,1,3,bi,bj,myid)         call diagnostics_fill(tmpdiag,'VDTRADSW',0,1,3,bi,bj,myid)
428        endif        endif
429    
430    c Total Precipitable Water (g/cm^2)
431    c ---------------------------------------------
432          if(diagnostics_is_on('TPW     ',myid) ) then
433           gravity = getcon('GRAVITY')
434           do j=jm1,jm2
435           do i=im1,im2
436           qbar(i,j) = 0.0
437           enddo
438           enddo
439           do L=1,Nrphys
440           do j=jm1,jm2
441           do i=im1,im2
442           qbar(i,j) = qbar(i,j) +
443         .             sphy(i,j,L,bi,bj)*dp(i,j,L,bi,bj)
444           enddo
445           enddo
446           enddo
447           do j=jm1,jm2
448           do i=im1,im2
449           tmpdiag(i,j) = qbar(i,j)*10. _d 0 /gravity
450           enddo
451           enddo
452           call diagnostics_fill(tmpdiag,'TPW     ',0,1,3,bi,bj,myid)
453          endif
454  #endif  #endif
455        return        return
456        end        end
Legend:
Removed from v.1.17  
changed lines
  Added in v.1.18
  ViewVC Help
Powered by ViewVC 1.1.22