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

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

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

revision 1.4 by molod, Fri Oct 22 17:56:03 2004 UTC revision 1.7 by molod, Tue Dec 14 19:56:45 2004 UTC
# Line 37  C ************************************** Line 37  C **************************************
37    
38  #ifdef ALLOW_DIAGNOSTICS  #ifdef ALLOW_DIAGNOSTICS
39  #include "SIZE.h"  #include "SIZE.h"
40  #include "diagnostics_SIZE.h"  #include "DIAGNOSTICS_SIZE.h"
41  #include "diagnostics.h"  #include "DIAGNOSTICS.h"
42  #endif  #endif
43    
44        integer myid,im,jm,lm,bi,bj        integer myid,im,jm,lm,bi,bj
# Line 71  C ************************************** Line 71  C **************************************
71         enddo         enddo
72        enddo        enddo
73    
       do L=1,lm  
        do j=1,jm  
        do i=1,im  
         dz(i,j,L) = cpog * (pekap(i,j,L+1)-pekap(i,j,L)) * t(i,j,L) *  
      .                (1.+virtcon*s(i,j,L))  
        enddo  
        enddo  
       enddo  
   
74        do j=1,jm        do j=1,jm
75        do i=1,im        do i=1,im
76         z(i,j,lm) = zsurf(i,j) +  0.5 * dz(i,j,lm)         z(i,j,lm) = zsurf(i,j) +  0.5 * dz(i,j,lm)
# Line 97  C ************************************** Line 88  C **************************************
88        do L=1,lm        do L=1,lm
89         do j=1,jm         do j=1,jm
90         do i=1,im         do i=1,im
91          rf(i,j,L) = dampcoef*(1+tanh((z(i,j,L)-80000.)/5000.))/86400.          rf(i,j,L) = dampcoef*(1+tanh((z(i,j,L)-50000.)/5000.))/86400.
92          rfu(i,j,L) = - rf(i,j,L) * u(i,j,L)          rfu(i,j,L) = - rf(i,j,L) * u(i,j,L)
93          rfv(i,j,L) = - rf(i,j,L) * v(i,j,L)          rfv(i,j,L) = - rf(i,j,L) * v(i,j,L)
94          rft(i,j,L) = -(u(i,j,L)*rfu(i,j,L) + v(i,j,L)*rfv(i,j,L) )*cpinv          rft(i,j,L) = -(u(i,j,L)*rfu(i,j,L) + v(i,j,L)*rfv(i,j,L) )*cpinv
Legend:
Removed from v.1.4  
changed lines
  Added in v.1.7
  ViewVC Help
Powered by ViewVC 1.1.22