/[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.76 by gforget, Fri Oct 8 21:17:41 2010 UTC revision 1.77 by gforget, Fri Oct 8 22:08:14 2010 UTC
# Line 460  CADJ &                          key = ii Line 460  CADJ &                          key = ii
460             IF ( theta(I,J,kSurface,bi,bj) .GE. TBC ) THEN             IF ( theta(I,J,kSurface,bi,bj) .GE. TBC ) THEN
461                a_QbyICE(i,j) = SEAICE_availHeatFrac                a_QbyICE(i,j) = SEAICE_availHeatFrac
462       &             * (theta(I,J,kSurface,bi,bj)-TBC) * dRf(kSurface)       &             * (theta(I,J,kSurface,bi,bj)-TBC) * dRf(kSurface)
463       &             * hFacC(i,j,kSurface,bi,bj) / 72.0764 _d 0       &             * hFacC(i,j,kSurface,bi,bj) *
464       &             * ( - QI / SEAICE_deltaTtherm )       &             (- HeatCapacity_Cp*rhoConst/SEAICE_deltaTtherm)
465             ELSE             ELSE
466                a_QbyICE(i,j) = SEAICE_availHeatFracFrz                a_QbyICE(i,j) = SEAICE_availHeatFracFrz
467       &             * (theta(I,J,kSurface,bi,bj)-TBC) * dRf(kSurface)       &             * (theta(I,J,kSurface,bi,bj)-TBC) * dRf(kSurface)
468       &             * hFacC(i,j,kSurface,bi,bj) / 72.0764 _d 0       &             * hFacC(i,j,kSurface,bi,bj) *
469       &             * ( - QI / SEAICE_deltaTtherm )       &             (- HeatCapacity_Cp*rhoConst/SEAICE_deltaTtherm)
470             ENDIF             ENDIF
471            ELSE            ELSE
472             a_QbyICE(i,j) = 0.             a_QbyICE(i,j) = 0.
473            ENDIF            ENDIF
474  cgf heat and water conservation: NOT OK -- the "QI" is not  cgf heat and water conservation: ok -- since rid of 72.0764 factor
 cgf consistent with the value involved in the 72.0764 factor  
 cgf (see comment in seaice_growth.F v1.70) -- but the error  
 cgf is at least made consistently when 72.0764 is used  
475           ENDDO           ENDDO
476          ENDDO          ENDDO
477    
# Line 904  C ====================================== Line 901  C ======================================
901    
902          DO J=1,sNy          DO J=1,sNy
903           DO I=1,sNx           DO I=1,sNx
           tmpscal1=  
      &         -maskC(I,J,kSurface,bi,bj)  
      &         *(HeatCapacity_Cp*rUnit2mass/QI)  
      &         *72.0764 _d 0  
904            QNET(I,J,bi,bj) = QNET(I,J,bi,bj)            QNET(I,J,bi,bj) = QNET(I,J,bi,bj)
905       &         - ( d_QbyICE(I,J) + d_QbySNW(I,J) )       &         + ( d_QbyICE(I,J) + d_QbySNW(I,J) )
906       &         * tmpscal1       &         * maskC(I,J,kSurface,bi,bj)
907  cgf heat and water conservation: NOT OK -- the "QI" is not  cgf heat and water conservation: ok -- since rid of 72.0764 factor
 cgf consistent with the value involved in the 72.0764 factor  
 cgf (see comment in seaice_growth.F v1.70) -- but the error  
 cgf is at least made consistently when 72.0764 is used  
908           ENDDO           ENDDO
909          ENDDO          ENDDO
910    
Legend:
Removed from v.1.76  
changed lines
  Added in v.1.77
  ViewVC Help
Powered by ViewVC 1.1.22