/[MITgcm]/MITgcm/pkg/gmredi/gmredi_k3d.F
ViewVC logotype

Diff of /MITgcm/pkg/gmredi/gmredi_k3d.F

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

revision 1.7 by jmc, Thu Aug 22 23:01:47 2013 UTC revision 1.8 by m_bates, Sat Aug 31 17:33:55 2013 UTC
# Line 64  C     kLow_V :: Local version of the ind Line 64  C     kLow_V :: Local version of the ind
64  C     N2loc  :: local N**2  C     N2loc  :: local N**2
65  C     slope  :: local slope  C     slope  :: local slope
66  C     Req    :: local equatorial deformation radius (m)  C     Req    :: local equatorial deformation radius (m)
 C     Rurms  :: a local mixing length used in calculation of urms (m)  
67  C     deltaH :: local thickness of Eady integration (m)  C     deltaH :: local thickness of Eady integration (m)
68  C     g_reciprho_sq :: (gravity*recip_rhoConst)**2  C     g_reciprho_sq :: (gravity*recip_rhoConst)**2
69  C     M4loc  :: local M**4  C     M4loc  :: local M**4
# Line 77  C     small  :: a small number (to avoid Line 76  C     small  :: a small number (to avoid
76        _RL N2loc        _RL N2loc
77        _RL slope        _RL slope
78        _RL Req        _RL Req
       _RL Rurms  
79        _RL deltaH        _RL deltaH
80        _RL g_reciprho_sq        _RL g_reciprho_sq
81        _RL M4loc        _RL M4loc
# Line 89  C     small  :: a small number (to avoid Line 87  C     small  :: a small number (to avoid
87  C     dfdy    :: gradient of the Coriolis paramter, df/dy, 1/(m*s)  C     dfdy    :: gradient of the Coriolis paramter, df/dy, 1/(m*s)
88  C     dfdx    :: gradient of the Coriolis paramter, df/dx, 1/(m*s)  C     dfdx    :: gradient of the Coriolis paramter, df/dx, 1/(m*s)
89  C     gradf   :: gradient of the Coriolis paramter at a cell centre, 1/(m*s)  C     gradf   :: gradient of the Coriolis paramter at a cell centre, 1/(m*s)
90    C     Rurms  :: a local mixing length used in calculation of urms (m)
91  C     RRhines :: The Rhines scale (m)  C     RRhines :: The Rhines scale (m)
92  C     Rmix    :: Mixing length  C     Rmix    :: Mixing length
93  C     N2      :: Square of the buoyancy frequency (1/s**2)  C     N2      :: Square of the buoyancy frequency (1/s**2)
# Line 103  C     Ubaro   :: Barotropic velocity (m/ Line 102  C     Ubaro   :: Barotropic velocity (m/
102        _RL dfdx(   1-Olx:sNx+Olx,1-Oly:sNy+Oly)        _RL dfdx(   1-Olx:sNx+Olx,1-Oly:sNy+Oly)
103        _RL gradf(  1-Olx:sNx+Olx,1-Oly:sNy+Oly)        _RL gradf(  1-Olx:sNx+Olx,1-Oly:sNy+Oly)
104        _RL dummy(  1-Olx:sNx+Olx,1-Oly:sNy+Oly)        _RL dummy(  1-Olx:sNx+Olx,1-Oly:sNy+Oly)
105          _RL Rurms(  1-Olx:sNx+Olx,1-Oly:sNy+Oly)
106        _RL RRhines(1-Olx:sNx+Olx,1-Oly:sNy+Oly)        _RL RRhines(1-Olx:sNx+Olx,1-Oly:sNy+Oly)
107        _RL Rmix(   1-Olx:sNx+Olx,1-Oly:sNy+Oly)        _RL Rmix(   1-Olx:sNx+Olx,1-Oly:sNy+Oly)
108        _RL N2(     1-Olx:sNx+Olx,1-Oly:sNy+Oly,Nr)        _RL N2(     1-Olx:sNx+Olx,1-Oly:sNy+Oly,Nr)
# Line 434  C     ================================== Line 434  C     ==================================
434        DO j=1-Oly+1,sNy+Oly        DO j=1-Oly+1,sNy+Oly
435         DO i=1-Olx+1,sNx+Olx-1         DO i=1-Olx+1,sNx+Olx-1
436  C       Calculate the Visbeck velocity  C       Calculate the Visbeck velocity
437          Rurms = MIN(Rdef(i,j,bi,bj),GM_K3D_maxLurms)          Rurms(i,j) = MIN(Rdef(i,j,bi,bj),GM_K3D_maxLurms)
438          urms(i,j,1) = GM_K3D_Lambda*eady(i,j)*Rurms          urms(i,j,1) = GM_K3D_Lambda*eady(i,j)*Rurms(i,j)
439  C       Set the bottom urms to zero  C       Set the bottom urms to zero
440          k=kLow_C(i,j)          k=kLow_C(i,j)
441          IF (k.GT.0) urms(i,j,k) = 0.0          IF (k.GT.0) urms(i,j,k) = 0.0
# Line 899  C     Diagnostics Line 899  C     Diagnostics
899          CALL DIAGNOSTICS_FILL(K3D,    'GM_K3D  ',0,Nr,0,1,1,myThid)          CALL DIAGNOSTICS_FILL(K3D,    'GM_K3D  ',0,Nr,0,1,1,myThid)
900          CALL DIAGNOSTICS_FILL(urms,   'GM_URMS ',0,Nr,0,1,1,myThid)          CALL DIAGNOSTICS_FILL(urms,   'GM_URMS ',0,Nr,0,1,1,myThid)
901          CALL DIAGNOSTICS_FILL(Rdef,   'GM_RDEF ',0, 1,0,1,1,myThid)          CALL DIAGNOSTICS_FILL(Rdef,   'GM_RDEF ',0, 1,0,1,1,myThid)
902          CALL DIAGNOSTICS_FILL(RRhines,'GM_LRHNS',0, 1,0,1,1,myThid)          CALL DIAGNOSTICS_FILL(Rurms,  'GM_RURMS',0, 1,0,1,1,myThid)
903            CALL DIAGNOSTICS_FILL(RRhines,'GM_RRHNS',0, 1,0,1,1,myThid)
904          CALL DIAGNOSTICS_FILL(Rmix,   'GM_RMIX ',0, 1,0,1,1,myThid)          CALL DIAGNOSTICS_FILL(Rmix,   'GM_RMIX ',0, 1,0,1,1,myThid)
905          CALL DIAGNOSTICS_FILL(supp,   'GM_SUPP ',0,Nr,0,1,1,myThid)          CALL DIAGNOSTICS_FILL(supp,   'GM_SUPP ',0,Nr,0,1,1,myThid)
906          CALL DIAGNOSTICS_FILL(Xix,    'GM_Xix  ',0,Nr,0,1,1,myThid)          CALL DIAGNOSTICS_FILL(Xix,    'GM_Xix  ',0,Nr,0,1,1,myThid)
Legend:
Removed from v.1.7  
changed lines
  Added in v.1.8
  ViewVC Help
Powered by ViewVC 1.1.22