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

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

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

revision 1.2 by jmc, Wed Apr 7 23:40:34 2004 UTC revision 1.3 by jmc, Thu Jul 8 15:56:36 2004 UTC
# Line 24  C     === Global variables === Line 24  C     === Global variables ===
24  #include "EEPARAMS.h"  #include "EEPARAMS.h"
25  #include "PARAMS.h"  #include "PARAMS.h"
26  #include "GRID.h"  #include "GRID.h"
27    #include "SURFACE.h"
28  #include "DYNVARS.h"  #include "DYNVARS.h"
29  #include "FFIELDS.h"  #include "FFIELDS.h"
30  #include "THSICE_PARAMS.h"  #include "THSICE_PARAMS.h"
# Line 72  C---+----1----+----2----+----3----+----4 Line 73  C---+----1----+----2----+----3----+----4
73        DO bj=myByLo(myThid),myByHi(myThid)        DO bj=myByLo(myThid),myByHi(myThid)
74         DO bi=myBxLo(myThid),myBxHi(myThid)         DO bi=myBxLo(myThid),myBxHi(myThid)
75    
76    C--     Mixed layer thickness: take the 1rst layer
77    #ifdef NONLIN_FRSURF
78            IF ( staggerTimeStep .AND. nonlinFreeSurf.GT.0 ) THEN
79             IF ( select_rStar.GT.0 ) THEN
80              DO j = jMin, jMax
81               DO i = iMin, iMax
82                 hOceMxL(i,j,bi,bj) = drF(1)*h0FacC(i,j,1,bi,bj)
83         &                                  *rStarFacC(i,j,bi,bj)
84               ENDDO
85              ENDDO
86             ELSE
87              DO j = jMin, jMax
88               DO i = iMin, iMax
89                IF ( ksurfC(i,j,bi,bj).EQ.1 ) THEN
90                 hOceMxL(i,j,bi,bj) = drF(1)*hFac_surfC(i,j,bi,bj)
91                ELSE
92                 hOceMxL(i,j,bi,bj) = drF(1)*hfacC(i,j,1,bi,bj)
93                ENDIF
94               ENDDO
95              ENDDO
96             ENDIF
97            ELSE
98    #else /* ndef NONLIN_FRSURF */
99            IF (.TRUE.) THEN
100    #endif /* NONLIN_FRSURF */
101              DO j = jMin, jMax
102               DO i = iMin, iMax
103                 hOceMxL(i,j,bi,bj) = drF(1)*hfacC(i,j,1,bi,bj)
104               ENDDO
105              ENDDO
106            ENDIF
107    
108           DO j = jMin, jMax           DO j = jMin, jMax
109            DO i = iMin, iMax            DO i = iMin, iMax
            hOceMxL(i,j,bi,bj) = hfacC(i,j,1,bi,bj)*drF(1)  
110             tOceMxL(i,j,bi,bj) = theta(i,j,1,bi,bj)             tOceMxL(i,j,bi,bj) = theta(i,j,1,bi,bj)
111             sOceMxL(i,j,bi,bj) = salt (i,j,1,bi,bj)             sOceMxL(i,j,bi,bj) = salt (i,j,1,bi,bj)
112             v2ocMxL(i,j,bi,bj) =             v2ocMxL(i,j,bi,bj) =
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22