/[MITgcm]/MITgcm/pkg/thsice/thsice_get_ocean.F
ViewVC logotype

Diff of /MITgcm/pkg/thsice/thsice_get_ocean.F

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

revision 1.1 by jmc, Sat Mar 30 18:38:50 2013 UTC revision 1.2 by jmc, Mon Apr 1 18:18:15 2013 UTC
# Line 2  C $Header$ Line 2  C $Header$
2  C $Name$  C $Name$
3    
4  #include "THSICE_OPTIONS.h"  #include "THSICE_OPTIONS.h"
5    #ifdef ALLOW_SEAICE
6    # include "SEAICE_OPTIONS.h"
7    #endif /* ALLOW_SEAICE */
8    
9  CBOP  CBOP
10  C     !ROUTINE: THSICE_GET_OCEAN  C     !ROUTINE: THSICE_GET_OCEAN
# Line 29  C     === Global variables === Line 32  C     === Global variables ===
32  c#include "THSICE_PARAMS.h"  c#include "THSICE_PARAMS.h"
33  #include "THSICE_SIZE.h"  #include "THSICE_SIZE.h"
34  #include "THSICE_VARS.h"  #include "THSICE_VARS.h"
35    #ifdef ALLOW_SEAICE
36    # include "SEAICE_SIZE.h"
37    # include "SEAICE.h"
38    #endif /* ALLOW_SEAICE */
39    
40  C     !INPUT/OUTPUT PARAMETERS:  C     !INPUT/OUTPUT PARAMETERS:
41  C     === Routine arguments ===  C     === Routine arguments ===
# Line 47  C     !LOCAL VARIABLES: Line 54  C     !LOCAL VARIABLES:
54  C     === Local variables ===  C     === Local variables ===
55        INTEGER i,j        INTEGER i,j
56        INTEGER ks        INTEGER ks
57    #ifdef ALLOW_SEAICE
58          _RL uRel1, uRel2, vRel1, vRel2
59    #endif /* ALLOW_SEAICE */
60    
61  C--     Mixed layer thickness: take the 1rst layer  C--     Mixed layer thickness: take the 1rst layer
62          ks = 1          ks = 1
# Line 93  C--     Mixed layer thickness: take the Line 103  C--     Mixed layer thickness: take the
103            siceAlb (i,j,bi,bj) = 0. _d 0            siceAlb (i,j,bi,bj) = 0. _d 0
104           ENDDO           ENDDO
105          ENDDO          ENDDO
106          DO j=1-OLy,sNy+OLy-1          IF ( .NOT.useSEAICE ) THEN
107           DO i=1-OLx,sNx+OLx-1           DO j=1-OLy,sNy+OLy-1
108            v2ocMxL(i,j,bi,bj) =            DO i=1-OLx,sNx+OLx-1
109               v2ocMxL(i,j,bi,bj) =
110       &              ( uVel(i,j,ks,bi,bj) * uVel(i,j,ks,bi,bj)       &              ( uVel(i,j,ks,bi,bj) * uVel(i,j,ks,bi,bj)
111       &              + uVel(i+1,j,ks,bi,bj)*uVel(i+1,j,ks,bi,bj)       &              + uVel(i+1,j,ks,bi,bj)*uVel(i+1,j,ks,bi,bj)
112       &              + vVel(i,j+1,ks,bi,bj)*vVel(i,j+1,ks,bi,bj)       &              + vVel(i,j+1,ks,bi,bj)*vVel(i,j+1,ks,bi,bj)
# Line 103  C--     Mixed layer thickness: take the Line 114  C--     Mixed layer thickness: take the
114       &              )*0.5 _d 0       &              )*0.5 _d 0
115           ENDDO           ENDDO
116          ENDDO          ENDDO
117    #ifdef ALLOW_SEAICE
118            ELSE
119             DO j=1-OLy,sNy+OLy-1
120              DO i=1-OLx,sNx+OLx-1
121               uRel1 = uVel( i, j,ks,bi,bj)-uIce( i, j,bi,bj)
122               uRel2 = uVel(i+1,j,ks,bi,bj)-uIce(i+1,j,bi,bj)
123               vRel1 = vVel(i, j, ks,bi,bj)-vIce(i, j, bi,bj)
124               vRel2 = vVel(i,j+1,ks,bi,bj)-vIce(i,j+1,bi,bj)
125               v2ocMxL(i,j,bi,bj) =
126         &              ( ( uRel1*uRel1 + uRel2*uRel2 )
127         &              + ( vRel1*vRel1 + vRel2*vRel2 )
128         &              )*0.5 _d 0
129              ENDDO
130             ENDDO
131    #endif /* ALLOW_SEAICE */
132            ENDIF
133    
134  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
135  #endif  /*ALLOW_THSICE*/  #endif  /*ALLOW_THSICE*/
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22