Custom Query (2350 matches)

Filters
 
Or
 
  
 
Columns

Show under each result:


Results (4 - 6 of 2350)

1 2 3 4 5 6 7 8 9 10 11 12
Ticket Resolution Summary Owner Reporter
#2549 fixed wrong diagnotic of error handling in stp_ctl systeam cbricaud
Description

Context

Some mpi splitting can lead to have land over all the inner local domain of a processor and ocean only on the halo bands.

stp_ctl compute MAXVAL of different fields to chec model values:

IF sn_cfctl%l_runstat= .True. the MAXVAL values are computed first on each local domain and then over all the processor , in stp_ctl.

IF sn_cfctl%l_runstat= .False. the MAXVAL values are only computed on each local domain.

When sn_cfctl%l_runstat= .False., the model stops on processors where the inner local domain have only land with values such as: -1.797693134862316E+308 e+308

Analysis

In stp_ctl, MAXVAL are computed with a mask: zmax(3) = MAXVAL( -ts(:,:,:,jp_sal,Kmm), mask = llmsk )

llmsk is false over the halo band

When inner local domain have only land, all llmsk is set to false

MAXVAL documentation (https://gcc.gnu.org/onlinedocs/gfortran/MAXVAL.html) says that, if all of the elements of MASK are .FALSE., then the result is -HUGE(ARRAY)

Fix

check tmask values in the inner domain and set zmax to zero if the inner local domain have only land?

#2547 fixed in tools/DOMAINcfg/src/dommsk.F90: wrong arguments when calling usr_def_fmask systeam cbricaud
Description

Context

in tools/DOMAINcfg/src/dommsk.F90, usr_def_fmask with cn_cfg and nn_cfg as arguments

Analysis

These variables are not filled.

Fix

They shloud be replaced by cp_cfg and jp_cfg: replace CALL usr_def_fmask( cn_cfg, nn_cfg, fmask ) by CALL usr_def_fmask( cp_cfg, jp_cfg, fmask )

Declarations of cn_cfg and nn_cfg in par_oce.f90 should be removed.

#2545 fixed bug in heat budget for SI3 systeam clem
Description

Context

they are several flaws in the way SI3 treats its heat exchanges with the ocean.

1) heat conservation is not ensured while landfast is activated because we switched off the production of ice in the leads (in icethd_do.F90) without setting qlead to 0 when ice is thick enough to have friction on the sea floor (tau_icebfr/=0, icethd.F90). Also, switching off ice formation as soon as there is friction whatever the ice velocity or the ice concentration is not very smart. It was done to avoid a continuous formation of ice when ice is grounded.

=⇒ solution 1: set qlead to 0 and no ice production in the leads where there is friction (tau_icebfr/=0) and ice concentration is close to its maximum value (rn_amax)

=⇒ solution 2 (better?): set qlead to 0 and no ice production in the leads where ice velocity drops below 5.e-04 m/s (commonly accepted value for landfast) and ice concentration is close to its max value (rn_amax)

2) There is a mixed up between solar and non-solar fluxes reaching the ocean after ice thermodynamics, though the total heat is conserved (see iceupdate.F90). Both solar and non solar fluxes reaching the leads are used to melt sea ice, so this energy is somewhat lost for the ocean to increase its temperature. However, at the end, the ocean will still see the total solar flux that reached the leads while the non-solar is adapted to conserve heat.

=⇒ solution: consider that the solar flux which has reached the 1st level of the ocean is entirely consumed by the ice when the ocean is warming up and there is still some ice left.

3) The parameterization of lateral melting (icethd_da.F90) does not depend on any heat balance (contrary to basal and surface melting). It only depdends on the difference of temperature between the ocean and ice. Therefore, we can temporarily end up with a local supercooling.

=⇒ solution: none.

At least the first 2 defects must imperatively be corrected in 4.0 and the trunk

1 2 3 4 5 6 7 8 9 10 11 12
Note: See TracQuery for help on using queries.