/[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.144 by jmc, Fri Jan 13 14:19:25 2012 UTC revision 1.145 by jmc, Fri Jan 13 21:09:35 2012 UTC
# Line 95  C     constants Line 95  C     constants
95        _RL TBC, ICE2SNOW        _RL TBC, ICE2SNOW
96        _RL QI, QS        _RL QI, QS
97    
98    C--   TmixLoc        :: ocean surface/mixed-layer temperature (in K)
99          _RL TmixLoc       (1:sNx,1:sNy)
100    
101  C     a_QbyATM_cover :: available heat (in W/m^2) due to the interaction of  C     a_QbyATM_cover :: available heat (in W/m^2) due to the interaction of
102  C             the atmosphere and the ocean surface - for ice covered water  C             the atmosphere and the ocean surface - for ice covered water
103  C     a_QbyATM_open  :: same but for open water  C     a_QbyATM_open  :: same but for open water
# Line 687  C determine available heat due to the at Line 690  C determine available heat due to the at
690  C ================================================================  C ================================================================
691    
692  C ocean surface/mixed layer temperature  C ocean surface/mixed layer temperature
693          DO J=1,sNy          DO j=1,sNy
694           DO I=1,sNx           DO i=1,sNx
695            TMIX(I,J,bi,bj)=theta(I,J,kSurface,bi,bj)+celsius2K            TmixLoc(i,j) = theta(i,j,kSurface,bi,bj)+celsius2K
696           ENDDO           ENDDO
697          ENDDO          ENDDO
698    
# Line 704  C wind speed from exf Line 707  C wind speed from exf
707  CADJ STORE qnet(:,:,bi,bj) = comlev1_bibj, key = iicekey,byte=isbyte  CADJ STORE qnet(:,:,bi,bj) = comlev1_bibj, key = iicekey,byte=isbyte
708  CADJ STORE qsw(:,:,bi,bj)  = comlev1_bibj, key = iicekey,byte=isbyte  CADJ STORE qsw(:,:,bi,bj)  = comlev1_bibj, key = iicekey,byte=isbyte
709  cCADJ STORE UG = comlev1_bibj, key = iicekey,byte=isbyte  cCADJ STORE UG = comlev1_bibj, key = iicekey,byte=isbyte
710  cCADJ STORE TMIX(:,:,bi,bj)  = comlev1_bibj, key = iicekey,byte=isbyte  cCADJ STORE TmixLoc = comlev1_bibj, key = iicekey,byte=isbyte
711  #endif /* ALLOW_AUTODIFF_TAMC */  #endif /* ALLOW_AUTODIFF_TAMC */
712    
713          CALL SEAICE_BUDGET_OCEAN(          CALL SEAICE_BUDGET_OCEAN(
714       I       UG,       I       UG,
715       U       TMIX,       I       TmixLoc,
716       O       a_QbyATM_open, a_QSWbyATM_open,       O       a_QbyATM_open, a_QSWbyATM_open,
717       I       bi, bj, myTime, myIter, myThid )       I       bi, bj, myTime, myIter, myThid )
718    
# Line 1128  c           heat flux overcoming potenti Line 1131  c           heat flux overcoming potenti
1131       &            (1.0 _d 0 - AREApreTH(i,J))       &            (1.0 _d 0 - AREApreTH(i,J))
1132  c           Penetrative shortwave flux beyond first layer  c           Penetrative shortwave flux beyond first layer
1133  c           that is therefore not available to ice growth/melt  c           that is therefore not available to ice growth/melt
1134              tmpscal2=SWFRACB * a_QSWbyATM_open(I,J)              tmpscal2=SWFracB * a_QSWbyATM_open(I,J)
1135  #ifdef SEAICE_DO_OPEN_WATER_MELT  #ifdef SEAICE_DO_OPEN_WATER_MELT
1136  C           allow not only growth but also melt by open ocean heat flux  C           allow not only growth but also melt by open ocean heat flux
1137              tmpscal3=MAX(tmpscal1-tmpscal2, -HEFF(I,J,bi,bj))              tmpscal3=MAX(tmpscal1-tmpscal2, -HEFF(I,J,bi,bj))
Legend:
Removed from v.1.144  
changed lines
  Added in v.1.145
  ViewVC Help
Powered by ViewVC 1.1.22