/[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.19 by jmc, Tue Apr 28 18:21:18 2009 UTC revision 1.20 by jmc, Thu Jun 11 19:23:54 2009 UTC
# Line 111  C--     Mixed layer thickness: take the Line 111  C--     Mixed layer thickness: take the
111  #ifdef NONLIN_FRSURF  #ifdef NONLIN_FRSURF
112          IF ( staggerTimeStep .AND. nonlinFreeSurf.GT.0 ) THEN          IF ( staggerTimeStep .AND. nonlinFreeSurf.GT.0 ) THEN
113           IF ( select_rStar.GT.0 ) THEN           IF ( select_rStar.GT.0 ) THEN
114            DO j = jMin, jMax            DO j=1-OLy,sNy+OLy
115             DO i = iMin, iMax             DO i=1-OLx,sNx+OLx
116               hOceMxL(i,j,bi,bj) = drF(1)*h0FacC(i,j,1,bi,bj)               hOceMxL(i,j,bi,bj) = drF(1)*h0FacC(i,j,1,bi,bj)
117       &                                  *rStarFacC(i,j,bi,bj)       &                                  *rStarFacC(i,j,bi,bj)
118             ENDDO             ENDDO
119            ENDDO            ENDDO
120           ELSE           ELSE
121            DO j = jMin, jMax            DO j=1-OLy,sNy+OLy
122             DO i = iMin, iMax             DO i=1-OLx,sNx+OLx
123              IF ( ksurfC(i,j,bi,bj).EQ.1 ) THEN              IF ( ksurfC(i,j,bi,bj).EQ.1 ) THEN
124               hOceMxL(i,j,bi,bj) = drF(1)*hFac_surfC(i,j,bi,bj)               hOceMxL(i,j,bi,bj) = drF(1)*hFac_surfC(i,j,bi,bj)
125              ELSE              ELSE
# Line 132  C--     Mixed layer thickness: take the Line 132  C--     Mixed layer thickness: take the
132  #else /* ndef NONLIN_FRSURF */  #else /* ndef NONLIN_FRSURF */
133          IF (.TRUE.) THEN          IF (.TRUE.) THEN
134  #endif /* NONLIN_FRSURF */  #endif /* NONLIN_FRSURF */
135            DO j = jMin, jMax            DO j=1-OLy,sNy+OLy
136             DO i = iMin, iMax             DO i=1-OLx,sNx+OLx
137               hOceMxL(i,j,bi,bj) = drF(1)*hFacC(i,j,1,bi,bj)               hOceMxL(i,j,bi,bj) = drF(1)*hFacC(i,j,1,bi,bj)
138             ENDDO             ENDDO
139            ENDDO            ENDDO
Legend:
Removed from v.1.19  
changed lines
  Added in v.1.20
  ViewVC Help
Powered by ViewVC 1.1.22