/[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.18 by molod, Tue Jan 17 15:59:44 2006 UTC revision 1.19 by molod, Thu Mar 9 00:00:36 2006 UTC
# Line 12  C*************************************** Line 12  C***************************************
12    
13        integer myid,im1,im2,jm1,jm2,Nrphys,Nbi,Nbj,bi,bj,ntracer        integer myid,im1,im2,jm1,jm2,Nrphys,Nbi,Nbj,bi,bj,ntracer
14        _RL p(im2,jm2,Nbi,Nbj)        _RL p(im2,jm2,Nbi,Nbj)
15        _RL uphy(im2,jm2,Nrphys,Nbi,Nbj)        _RL uphy(im2,jm2,Nrphys)
16        _RL vphy(im2,jm2,Nrphys,Nbi,Nbj)        _RL vphy(im2,jm2,Nrphys)
17        _RL thphy(im2,jm2,Nrphys,Nbi,Nbj)        _RL thphy(im2,jm2,Nrphys)
18        _RL sphy(im2,jm2,Nrphys,Nbi,Nbj)        _RL sphy(im2,jm2,Nrphys)
19        _RL qq(im2,jm2,Nrphys,Nbi,Nbj),pk(im2,jm2,Nrphys,Nbi,Nbj)        _RL qq(im2,jm2,Nrphys,Nbi,Nbj),pk(im2,jm2,Nrphys,Nbi,Nbj)
20        _RL dp(im2,jm2,Nrphys,Nbi,Nbj)        _RL dp(im2,jm2,Nrphys,Nbi,Nbj)
21        _RL radswt(im2,jm2,Nbi,Nbj),radswg(im2,jm2,Nbi,Nbj)        _RL radswt(im2,jm2,Nbi,Nbj),radswg(im2,jm2,Nbi,Nbj)
# Line 275  c ------------------------ Line 275  c ------------------------
275        if(diagnostics_is_on('UWND    ',myid) ) then        if(diagnostics_is_on('UWND    ',myid) ) then
276         do j=jm1,jm2         do j=jm1,jm2
277         do i=im1,im2         do i=im1,im2
278          tmpdiag(i,j) = uphy(i,j,L,bi,bj)          tmpdiag(i,j) = uphy(i,j,L)
279         enddo         enddo
280         enddo         enddo
281         call diagnostics_fill(tmpdiag,'UWND    ',L,1,3,bi,bj,myid)         call diagnostics_fill(tmpdiag,'UWND    ',L,1,3,bi,bj,myid)
# Line 286  c ------------------------ Line 286  c ------------------------
286        if(diagnostics_is_on('VWND    ',myid) ) then        if(diagnostics_is_on('VWND    ',myid) ) then
287         do j=jm1,jm2         do j=jm1,jm2
288         do i=im1,im2         do i=im1,im2
289          tmpdiag(i,j) = vphy(i,j,L,bi,bj)          tmpdiag(i,j) = vphy(i,j,L)
290         enddo         enddo
291         enddo         enddo
292         call diagnostics_fill(tmpdiag,'VWND    ',L,1,3,bi,bj,myid)         call diagnostics_fill(tmpdiag,'VWND    ',L,1,3,bi,bj,myid)
# Line 297  c ---------------------- Line 297  c ----------------------
297        if(diagnostics_is_on('TMPU    ',myid) ) then        if(diagnostics_is_on('TMPU    ',myid) ) then
298         do j=jm1,jm2         do j=jm1,jm2
299         do i=im1,im2         do i=im1,im2
300          tmpdiag(i,j) = thphy(i,j,L,bi,bj)*pk(i,j,L,bi,bj)          tmpdiag(i,j) = thphy(i,j,L)*pk(i,j,L,bi,bj)
301         enddo         enddo
302         enddo         enddo
303         call diagnostics_fill(tmpdiag,'TMPU    ',L,1,3,bi,bj,myid)         call diagnostics_fill(tmpdiag,'TMPU    ',L,1,3,bi,bj,myid)
# Line 319  c ----------------------- Line 319  c -----------------------
319        if(diagnostics_is_on('SPHU    ',myid) ) then        if(diagnostics_is_on('SPHU    ',myid) ) then
320         do j=jm1,jm2         do j=jm1,jm2
321         do i=im1,im2         do i=im1,im2
322          tmpdiag(i,j) = sphy(i,j,L,bi,bj) * 1000.          tmpdiag(i,j) = sphy(i,j,L) * 1000.
323         enddo         enddo
324         enddo         enddo
325         call diagnostics_fill(tmpdiag,'SPHU    ',L,1,3,bi,bj,myid)         call diagnostics_fill(tmpdiag,'SPHU    ',L,1,3,bi,bj,myid)
# Line 440  c -------------------------------------- Line 440  c --------------------------------------
440         do j=jm1,jm2         do j=jm1,jm2
441         do i=im1,im2         do i=im1,im2
442         qbar(i,j) = qbar(i,j) +         qbar(i,j) = qbar(i,j) +
443       .             sphy(i,j,L,bi,bj)*dp(i,j,L,bi,bj)       .             sphy(i,j,L)*dp(i,j,L,bi,bj)
444         enddo         enddo
445         enddo         enddo
446         enddo         enddo
Legend:
Removed from v.1.18  
changed lines
  Added in v.1.19
  ViewVC Help
Powered by ViewVC 1.1.22