/[MITgcm]/MITgcm/pkg/aim_v23/phy_suflux_sice.F
ViewVC logotype

Diff of /MITgcm/pkg/aim_v23/phy_suflux_sice.F

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

revision 1.5 by jmc, Thu Jul 22 22:58:38 2004 UTC revision 1.6 by jmc, Mon Mar 13 03:58:32 2006 UTC
# Line 32  C     Resolution parameters Line 32  C     Resolution parameters
32  C-- size for MITgcm & Physics package :  C-- size for MITgcm & Physics package :
33  #include "AIM_SIZE.h"  #include "AIM_SIZE.h"
34  #include "EEPARAMS.h"  #include "EEPARAMS.h"
35    #include "PARAMS.h"
36    
37  C-- Physics package  C-- Physics package
38  #include "AIM_PARAMS.h"  #include "AIM_PARAMS.h"
# Line 91  CEOP Line 92  CEOP
92    
93  C-- Local variables:  C-- Local variables:
94  C    CDENVV :: surf. heat flux (sens.,lat.) coeff including stability effect  C    CDENVV :: surf. heat flux (sens.,lat.) coeff including stability effect
95  C    ALHevp :: Latent Heat of evaporation  C    ALHevp :: Latent Heat of evaporation
96        _RL CDENVV(NGP), RDTH, FSSICE        _RL CDENVV(NGP), RDTH, FSSICE
97        _RL ALHevp, Fstb0, dTstb, dFstb        _RL ALHevp, Fstb0, dTstb, dFstb
98        _RL QSAT0(NGP,2)        _RL QSAT0(NGP,2)
# Line 102  C---+----1----+----2----+----3----+----4 Line 103  C---+----1----+----2----+----3----+----4
103    
104        ALHevp = ALHC        ALHevp = ALHC
105  C     Evap of snow/ice: account for Latent Heat of freezing :  C     Evap of snow/ice: account for Latent Heat of freezing :
106        IF ( aim_energPrecip ) ALHevp = ALHC + ALHF        IF ( aim_energPrecip .OR. useThSIce ) ALHevp = ALHC + ALHF
107    
108  C     1.5 Define effective skin temperature to compensate for  C     1.5 Define effective skin temperature to compensate for
109  C         non-linearity of heat/moisture fluxes during the daily cycle  C         non-linearity of heat/moisture fluxes during the daily cycle
# Line 138  C   the derivative using the discrete fo Line 139  C   the derivative using the discrete fo
139         ENDDO         ENDDO
140  C-    deBug part:  C-    deBug part:
141  c      J = 6 + (17-1)*sNx  c      J = 6 + (17-1)*sNx
142  c      IF ( bi.EQ.3 .AND. J.LE.NGP )  c      IF ( bi.EQ.3 .AND. J.LE.NGP )
143  c    &  WRITE(6,1020)'SUFLUX_SICE: Stab=',Shf0(J),CDENVV(J),dShf(J)  c    &  WRITE(6,1020)'SUFLUX_SICE: Stab=',Shf0(J),CDENVV(J),dShf(J)
144        ENDIF        ENDIF
145    
# Line 176  C-    account for stability function der Line 177  C-    account for stability function der
177       &            + dShf(J)*CP*(TSKIN(J)-T0(J))       &            + dShf(J)*CP*(TSKIN(J)-T0(J))
178          dShf(J) = MAX( dShf(J), 0. _d 0 )          dShf(J) = MAX( dShf(J), 0. _d 0 )
179  C--   do not allow negative derivative vs Ts of Sensible+Latent H.flux:  C--   do not allow negative derivative vs Ts of Sensible+Latent H.flux:
180  C     a) quiet unrealistic ;  C     a) quiet unrealistic ;
181  C     b) garantee positive deriv. of total H.flux (needed for implicit solver)  C     b) garantee positive deriv. of total H.flux (needed for implicit solver)
182          dEvp(J) = MAX( dEvp(J), -dShf(J)/ALHevp )          dEvp(J) = MAX( dEvp(J), -dShf(J)/ALHevp )
183         ENDDO         ENDDO
Legend:
Removed from v.1.5  
changed lines
  Added in v.1.6
  ViewVC Help
Powered by ViewVC 1.1.22