Skip to content

Commit

Permalink
Merge pull request #377 from FESOM/refactoring_bugfix_ice_cutoff
Browse files Browse the repository at this point in the history
 make sure when we cutoff a_ice=0.0 that we also cutoff m_ice and m_snow
  • Loading branch information
koldunovn authored Jul 26, 2024
2 parents f59e8c1 + 2eab17b commit d6cbab8
Show file tree
Hide file tree
Showing 2 changed files with 10 additions and 10 deletions.
12 changes: 6 additions & 6 deletions setups/test_pi/setup.yml
Original file line number Diff line number Diff line change
Expand Up @@ -57,12 +57,12 @@ namelist.io:
prec: 8

fcheck:
a_ice: 0.2692498167586123
salt: 23.944089812053637
sst: 8.526792796322084
temp: 1.7018189804265962
u: -0.001431070160064853
v: 0.00014314235109323336
a_ice: 0.2692496175281426
salt: 23.94409404180963
sst: 8.527030848425198
temp: 1.7018242410566116
u: -0.001431025773970824
v: 0.00014310903096321633



Expand Down
8 changes: 4 additions & 4 deletions src/ice_thermo_oce.F90
Original file line number Diff line number Diff line change
Expand Up @@ -100,20 +100,20 @@ subroutine cut_off(ice, partit, mesh)
if (a_ice(n) > 1.0_WP) a_ice(n)=1.0_WP
! upper cutoff: a_ice
if (a_ice(n) < .1e-8_WP) then
a_ice(n)=0.0_WP
#if defined (__oifs) || defined (__ifsinterface)
a_ice(n) =0.0_WP
m_ice(n) =0.0_WP
m_snow(n) =0.0_WP
#if defined (__oifs) || defined (__ifsinterface)
ice_temp(n)=273.15_WP
#endif /* (__oifs) */
end if
!___________________________________________________________________________
! lower cutoff: m_ice
if (m_ice(n) < .1e-8_WP) then
m_ice(n)=0.0_WP
#if defined (__oifs) || defined (__ifsinterface)
m_ice(n) =0.0_WP
m_snow(n) =0.0_WP
a_ice(n) =0.0_WP
#if defined (__oifs) || defined (__ifsinterface)
ice_temp(n)=273.15_WP
#endif /* (__oifs) */
end if
Expand Down

0 comments on commit d6cbab8

Please sign in to comment.