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

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

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

revision 1.1 by jmc, Wed Apr 7 23:40:34 2004 UTC revision 1.2 by jmc, Thu Jun 3 16:34:36 2004 UTC
# Line 78  C--    add restoring (backward) toward c Line 78  C--    add restoring (backward) toward c
78         DO j=1,sNy         DO j=1,sNy
79          DO i=1,sNx          DO i=1,sNx
80           IF ( hOceMxL(i,j,bi,bj).NE.0. _d 0 ) THEN           IF ( hOceMxL(i,j,bi,bj).NE.0. _d 0 ) THEN
81              sOceMxL(i,j,bi,bj) =
82         &         (sOceMxL(i,j,bi,bj) + dtFac*sMxL_default)*dtFacR
83              oceTfreez = -mu_Tf*sOceMxL(i,j,bi,bj)
84            locTemp = ( aim_sWght0*aim_sst0(i,j,bi,bj)            locTemp = ( aim_sWght0*aim_sst0(i,j,bi,bj)
85       &              + aim_sWght1*aim_sst1(i,j,bi,bj)       &              + aim_sWght1*aim_sst1(i,j,bi,bj)
86       &              ) - celsius2K       &              ) - celsius2K
87            locTemp = MAX( locTemp , oceTfreez )            locTemp = MAX( locTemp , oceTfreez )
88            tOceMxL(i,j,bi,bj) =            tOceMxL(i,j,bi,bj) =
89       &         (tOceMxL(i,j,bi,bj) + dtFac*locTemp)*dtFacR       &         (tOceMxL(i,j,bi,bj) + dtFac*locTemp)*dtFacR
           sOceMxL(i,j,bi,bj) =  
      &         (sOceMxL(i,j,bi,bj) + dtFac*sMxL_default)*dtFacR  
90           ENDIF           ENDIF
91          ENDDO          ENDDO
92         ENDDO         ENDDO
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22