Dear Vladimir, dear Burkhadt,
the lstomata paramater is related to variable RSMIN INT2LM: “RSMIN” > prs_min_lm, see "src_gribtabs.f90" CCLM: "RSMIN" > rsmin2d, see “src_setup_vartab.f90”
In both codes the sea-land-dependency-flag “l” is set, see the aformementioned modules.
In INT2LM, “prs_min_lm” is initialized with the value “zero”, and obviously also above water bodies where it is not used.
Question: might it be possible that the sea-land-dependency of
RSMIN
is not taken into account coorectly in
CCLM
?
I only found one place in the code (module src_soil_multlay.f90, abount line 2382) where
RSMIN
(=rsmin2d) is used if lstomate=.TRUE.,
and there we have it in the denominator of an expression; division by zero might lead to NaN.
However, the expression is excecuted only
- IF (lstomata) THEN
and
- IF (llandmask(i,j))
THEN
! land points only
and two further conditions.
Vladimir, you should check the distribution and values of
RSMIN
in your laf-file of your 0.12 simulation.
It must not be equal to zero at any gridpoint over land
Hans-Juergen
Dear Vladimir, dear Burkhadt,
the lstomata paramater is related to variable RSMIN INT2LM: “RSMIN” > prs_min_lm, see "src_gribtabs.f90" CCLM: "RSMIN" > rsmin2d, see “src_setup_vartab.f90”
In both codes the sea-land-dependency-flag “l” is set, see the aformementioned modules.
In INT2LM, “prs_min_lm” is initialized with the value “zero”, and obviously also above water bodies where it is not used.
Question: might it be possible that the sea-land-dependency of RSMIN is not taken into account coorectly in CCLM ?
I only found one place in the code (module src_soil_multlay.f90, abount line 2382) where RSMIN (=rsmin2d) is used if lstomate=.TRUE.,
and there we have it in the denominator of an expression; division by zero might lead to NaN.
However, the expression is excecuted only
- IF (lstomata) THEN
and
- IF (llandmask(i,j)) THEN ! land points only
and two further conditions.
Vladimir, you should check the distribution and values of RSMIN in your laf-file of your 0.12 simulation.
It must not be equal to zero at any gridpoint over land
Hans-Juergen