/[MITgcm]/MITgcm/pkg/thsice/thsice_get_exf.F
ViewVC logotype

Diff of /MITgcm/pkg/thsice/thsice_get_exf.F

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

revision 1.3 by jmc, Tue Jun 6 22:27:18 2006 UTC revision 1.4 by jmc, Sun Jun 25 22:22:38 2006 UTC
# Line 59  CEOP Line 59  CEOP
59    
60  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
61  C     === Local variables ===  C     === Local variables ===
62    C     hsLocal, hlLocal :: sensible & latent heat flux over sea-ice
63    
64        _RL     aln        _RL     aln
65    
66          _RL     hsLocal, hlLocal
67        integer iter        integer iter
68        _RL     delq        _RL     delq
69        _RL     deltap        _RL     deltap
# Line 213  C     coefficients. Line 215  C     coefficients.
215             tau    = tau*us(i,j,bi,bj)/sh(i,j,bi,bj)             tau    = tau*us(i,j,bi,bj)/sh(i,j,bi,bj)
216    
217             evapLoc       = -tau*qstar/ustar             evapLoc       = -tau*qstar/ustar
218             hl(i,j,bi,bj) = -lath*evapLoc             hlLocal = -lath*evapLoc
219             hs(i,j,bi,bj) = atmcp*tau*tstar/ustar             hsLocal = atmcp*tau*tstar/ustar
220  #ifndef EXF_READ_EVAP  #ifndef EXF_READ_EVAP
221  cdm             evap(i,j,bi,bj)    = tau*qstar/ustar  cdm             evap(i,j,bi,bj)    = tau*qstar/ustar
222  cdm !!! need to change sign and to convert from kg/m^2/s to m/s !!!  cdm !!! need to change sign and to convert from kg/m^2/s to m/s !!!
223             evap(i,j,bi,bj)    = -recip_rhonil*evapLoc  C jmc: do not reset evap which contains evaporation over ice-free ocean fraction
224    c          evap(i,j,bi,bj)    = -recip_rhonil*evapLoc
225  #endif  #endif
226    
227  C---  surf.Temp derivative of turbulent Fluxes  C---  surf.Temp derivative of turbulent Fluxes
# Line 241  C--   Total derivative with respect to s Line 244  C--   Total derivative with respect to s
244             df0dT = -dflwupdT+dfshdT+dflhdT             df0dT = -dflwupdT+dfshdT+dflhdT
245    
246             flwNet_dwn  = lwdown(i,j,bi,bj) - flwup             flwNet_dwn  = lwdown(i,j,bi,bj) - flwup
247             flxExceptSw = flwNet_dwn + hs(i,j,bi,bj) + hl(i,j,bi,bj)             flxExceptSw = flwNet_dwn + hsLocal + hlLocal
248    
249            endif            endif
250    
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.22