/[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.99 by gforget, Fri Nov 12 23:10:07 2010 UTC revision 1.100 by mlosch, Fri Nov 19 14:32:11 2010 UTC
# Line 117  c Line 117  c
117        _RL d_HEFFbyOCNonICE    (1:sNx,1:sNy)        _RL d_HEFFbyOCNonICE    (1:sNx,1:sNy)
118        _RL d_HEFFbyATMonOCN    (1:sNx,1:sNy)        _RL d_HEFFbyATMonOCN    (1:sNx,1:sNy)
119        _RL d_HEFFbyFLOODING    (1:sNx,1:sNy)        _RL d_HEFFbyFLOODING    (1:sNx,1:sNy)
120    #ifdef SEAICE_CAP_HEFF
121          _RL d_HEFFbyCAPPING     (1:sNx,1:sNy)
122    #endif /* SEAICE_CAP_HEFF */
123  c  c
124        _RL d_HEFFbyATMonOCN_open(1:sNx,1:sNy)        _RL d_HEFFbyATMonOCN_open(1:sNx,1:sNy)
125  c  c
# Line 262  c Line 265  c
265            d_HEFFbyOCNonICE(I,J)      = 0.0 _d 0            d_HEFFbyOCNonICE(I,J)      = 0.0 _d 0
266            d_HEFFbyATMonOCN(I,J)      = 0.0 _d 0            d_HEFFbyATMonOCN(I,J)      = 0.0 _d 0
267            d_HEFFbyFLOODING(I,J)      = 0.0 _d 0            d_HEFFbyFLOODING(I,J)      = 0.0 _d 0
268    #ifdef SEAICE_CAP_HEFF
269              d_HEFFbyCAPPING(I,J)       = 0.0 _d 0
270    #endif /* SEAICE_CAP_HEFF */
271  c  c
272            d_HEFFbyATMonOCN_open(I,J) = 0.0 _d 0            d_HEFFbyATMonOCN_open(I,J) = 0.0 _d 0
273  c  c
# Line 1088  CADJ &                         key = iic Line 1094  CADJ &                         key = iic
1094            HEFF(I,J,bi,bj) = HEFF(I,J,bi,bj)*HEFFM(I,J,bi,bj)            HEFF(I,J,bi,bj) = HEFF(I,J,bi,bj)*HEFFM(I,J,bi,bj)
1095  #ifdef SEAICE_CAP_HEFF  #ifdef SEAICE_CAP_HEFF
1096  C     This is not energy conserving, but at least it conserves fresh water  C     This is not energy conserving, but at least it conserves fresh water
1097            tmpscal0         = MAX(HEFF(I,J,bi,bj)-MAX_HEFF,0. _d 0)            d_HEFFbyCAPPING(I,J) = -MAX(HEFF(I,J,bi,bj)-MAX_HEFF,0. _d 0)
1098            HEFF(I,J,bi,bj)  = HEFF(I,J,bi,bj) - tmpscal0            HEFF(I,J,bi,bj)  = HEFF(I,J,bi,bj) + d_HEFFbyCAPPING(I,J)
           EmPmR(I,J,bi,bj) = EmPmR(I,J,bi,bj)  
      &         - tmpscal0*SEAICE_rhoIce/SEAICE_deltaTtherm  
1099  #endif /* SEAICE_CAP_HEFF */  #endif /* SEAICE_CAP_HEFF */
1100            HSNOW(I,J,bi,bj)  = HSNOW(I,J,bi,bj)*HEFFM(I,J,bi,bj)            HSNOW(I,J,bi,bj)  = HSNOW(I,J,bi,bj)*HEFFM(I,J,bi,bj)
1101           ENDDO           ENDDO
# Line 1271  C ====================================== Line 1275  C ======================================
1275       &            +d_HEFFbyATMonOCN(I,J)       &            +d_HEFFbyATMonOCN(I,J)
1276       &            +d_HEFFbyNEG(I,J)       &            +d_HEFFbyNEG(I,J)
1277       &            +d_HSNWbyNEG(I,J)/ICE2SNOW       &            +d_HSNWbyNEG(I,J)/ICE2SNOW
1278    #ifdef SEAICE_CAP_HEFF
1279         &            +d_HEFFbyCAPPING(I,J)
1280    #endif /* SEAICE_CAP_HEFF */
1281            EmPmR(I,J,bi,bj)  = maskC(I,J,kSurface,bi,bj)*(            EmPmR(I,J,bi,bj)  = maskC(I,J,kSurface,bi,bj)*(
1282       &         ( EVAP(I,J,bi,bj)-PRECIP(I,J,bi,bj) )       &         ( EVAP(I,J,bi,bj)-PRECIP(I,J,bi,bj) )
1283       &         * ( ONE - AREApreTH(I,J) )       &         * ( ONE - AREApreTH(I,J) )
Legend:
Removed from v.1.99  
changed lines
  Added in v.1.100
  ViewVC Help
Powered by ViewVC 1.1.22