/[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.96 by gforget, Sun Oct 31 20:07:01 2010 UTC revision 1.97 by mlosch, Fri Nov 5 08:10:57 2010 UTC
# Line 70  c or a_QbyATM_cover, the naming conventi Line 70  c or a_QbyATM_cover, the naming conventi
70  c    The prefix 'a_' means available, the prefix 'd_' means delta  c    The prefix 'a_' means available, the prefix 'd_' means delta
71  c       (i.e. increment), and the prefix 'r_' means residual.  c       (i.e. increment), and the prefix 'r_' means residual.
72  c    The suffix '_cover' denotes a value for the ice covered fraction  c    The suffix '_cover' denotes a value for the ice covered fraction
73  c       of the grid cell, whereas '_cover' is for the open water fraction.  c       of the grid cell, whereas '_open' is for the open water fraction.
74  c    The main part of the name states what ice/snow stock is concerned  c    The main part of the name states what ice/snow stock is concerned
75  c       (e.g. QbyATM or HEFF), and how it is affected (e.g. d_HEFFbyATMonOCN  c       (e.g. QbyATM or HEFF), and how it is affected (e.g. d_HEFFbyATMonOCN
76  c       is the increment of HEFF due to the ATMosphere extracting heat from the  c       is the increment of HEFF due to the ATMosphere extracting heat from the
# Line 724  cgf warming conditions, the lab_sea resu Line 724  cgf warming conditions, the lab_sea resu
724            tmpscal2 = MAX(-HEFF(I,J,bi,bj),r_QbyATM_cover(I,J))            tmpscal2 = MAX(-HEFF(I,J,bi,bj),r_QbyATM_cover(I,J))
725            d_HEFFbyATMonOCN(I,J)=d_HEFFbyATMonOCN(I,J)+tmpscal2            d_HEFFbyATMonOCN(I,J)=d_HEFFbyATMonOCN(I,J)+tmpscal2
726            r_QbyATM_cover(I,J)=r_QbyATM_cover(I,J)-tmpscal2            r_QbyATM_cover(I,J)=r_QbyATM_cover(I,J)-tmpscal2
 c The following line is what I referred to as the 'vintage bug'  
 c If I mistook a desirable feature for a bug, then we will un-comment  
 c this line. Otherwise we may want to delete the comment at some point.  
 c     &                                            *AREApreTH(I,J)  
727            HEFF(I,J,bi,bj) = HEFF(I,J,bi,bj) + tmpscal2            HEFF(I,J,bi,bj) = HEFF(I,J,bi,bj) + tmpscal2
728           ENDDO           ENDDO
729          ENDDO          ENDDO
# Line 1098  CADJ &                         key = iic Line 1094  CADJ &                         key = iic
1094            AREA(I,J,bi,bj) = AREA(I,J,bi,bj)*HEFFM(I,J,bi,bj)            AREA(I,J,bi,bj) = AREA(I,J,bi,bj)*HEFFM(I,J,bi,bj)
1095            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)
1096  #ifdef SEAICE_CAP_HEFF  #ifdef SEAICE_CAP_HEFF
1097            HEFF(I,J,bi,bj)=MIN(MAX_HEFF,HEFF(I,J,bi,bj))  C     This is not energy conserving, but at least it conserves fresh water
1098              tmpscal0         = MAX(HEFF(I,J,bi,bj)-MAX_HEFF,0. _d 0)
1099              HEFF(I,J,bi,bj)  = HEFF(I,J,bi,bj) - tmpscal0
1100              EmPmR(I,J,bi,bj) = EmPmR(I,J,bi,bj)
1101         &         - tmpscal0*SEAICE_rhoIce/SEAICE_deltaTtherm
1102  #endif /* SEAICE_CAP_HEFF */  #endif /* SEAICE_CAP_HEFF */
1103            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)
1104           ENDDO           ENDDO
Legend:
Removed from v.1.96  
changed lines
  Added in v.1.97
  ViewVC Help
Powered by ViewVC 1.1.22