#199 closed Bug (fixed)
sbcblk_core.F90 : division by 0 in blk_ice_core
Reported by: | molines@… | Owned by: | nemo |
---|---|---|---|
Priority: | highest | Milestone: | |
Component: | OCE | Version: | v2 |
Severity: | Keywords: | LIM* OPA SBC v2 | |
Cc: |
Description
In the core bulk formula for ice, latent heat fluxuses an expression like
p_qla(ji,jj,jl) = MAX( 0.e0, rhoa * Ls * Cice * z_wnds_t(ji,jj) & & * ( 11637800. * EXP( -5897.8 / pst(ji,jj,jl) ) / rhoa - sf(jp_humi)%fnow(ji,jj) )
and its derivative with respect to temperature
p_dqla(ji,jj,jl) = zcoef_dqla * z_wnds_t(ji,jj) / ( zst2 ) * EXP( -5897.8 / pst(ji,jj,jl) )
Here pst is the sist (Sea Ice Surface Temperature) in K and zst2=pst*pst. Unfortunatly ..., pst can be 0 which produces a FPE (division by 0)
The fix : in sbcice_lim_2 (and probably sbcice_lim_3) the call to blk_ice_core passes sist as a local 'zsist'. I suggest to modified zsist, from
zsist (:,:,1) = sist (:,:)
to
zsist (:,:,1) = sist (:,:) + rt0 * ( 1-tmask(:,:,1))
Commit History (1)
Changeset | Author | Time | ChangeLog |
---|---|---|---|
1109 | ctlod | 2008-06-13T11:59:54+02:00 | trunk: ensure that surface sea-ice temperature is set to 273 degK over land area, see ticket: #199 |
Change History (7)
comment:1 Changed 16 years ago by ctlod
- Resolution set to fixed
- Status changed from new to closed
comment:2 Changed 8 years ago by nicolasmartin
- Milestone New Surface Module deleted
Milestone New Surface Module deleted
comment:3 Changed 8 years ago by nicolasmartin
- Keywords SBC LIM2 added; sbcice_lim_2 removed
comment:4 Changed 8 years ago by nicolasmartin
- Keywords LIM3 added; sbcice_lim_3 removed
comment:5 Changed 7 years ago by nemo
- Keywords LIM* added; LIM2 LIM3 removed
comment:6 Changed 2 years ago by nemo
- Keywords OPA r2 added
comment:7 Changed 2 years ago by nemo
- Keywords v2 added; r2 removed
Note: See
TracTickets for help on using
tickets.
it has also been done in sbcice_lim.F90 for the field t_su(:,:,:)