/[MITgcm]/MITgcm/pkg/seaice/seaice_growth.F
ViewVC logotype

Diff of /MITgcm/pkg/seaice/seaice_growth.F

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

revision 1.54 by mlosch, Wed Mar 18 14:02:25 2009 UTC revision 1.55 by mlosch, Wed Mar 18 14:35:37 2009 UTC
# Line 787  C     use (abuse) GHEFF to diagnose the Line 787  C     use (abuse) GHEFF to diagnose the
787  C     convert snow to ice if submerged  C     convert snow to ice if submerged
788           DO J=1,sNy           DO J=1,sNy
789            DO I=1,sNx            DO I=1,sNx
790             hDraft     = (HSNOW(I,J,bi,bj)*330. _d 0             hDraft     = (HSNOW(I,J,bi,bj)*SEAICE_rhoSnow
791       &              +HEFF(I,J,1,bi,bj)*SEAICE_rhoIce)/1000. _d 0       &              +HEFF(I,J,1,bi,bj)*SEAICE_rhoIce)/1000. _d 0
792  C     here GEFF is the gain of ice due to flooding  C     here GEFF is the gain of ice due to flooding
793             GHEFF(I,J) = hDraft - MIN(hDraft,HEFF(I,J,1,bi,bj))             GHEFF(I,J) = hDraft - MIN(hDraft,HEFF(I,J,1,bi,bj))
# Line 816  C     turn GHEFF into a rate Line 816  C     turn GHEFF into a rate
816           DO J=1,sNy           DO J=1,sNy
817            DO I=1,sNx            DO I=1,sNx
818             sIceLoad(i,j,bi,bj) = HEFF(I,J,1,bi,bj)*SEAICE_rhoIce             sIceLoad(i,j,bi,bj) = HEFF(I,J,1,bi,bj)*SEAICE_rhoIce
819       &                         + HSNOW(I,J,bi,bj)* 330. _d 0       &                         + HSNOW(I,J,bi,bj)*SEAICE_rhoSnow
820            ENDDO            ENDDO
821           ENDDO           ENDDO
822          ENDIF          ENDIF
Legend:
Removed from v.1.54  
changed lines
  Added in v.1.55
  ViewVC Help
Powered by ViewVC 1.1.22