Adapt "time"-variables from ECHAM accumulation to ICON accumulation
Posted by David Neubauer 11/02/2023 https://redmine.hammoz.ethz.ch/issues/890
In the Lohmann 2-moment cloud microphysics scheme are several "time"-variables used to track how often clouds, cloud liquid, cloud ice etc. are present. I.e. cloud_time, burden_time, cliwc_time, burdic_time, reffl_time, reffi_time, conv_time
. These are used to compute e.g. in-cloud variables.
In ECHAM-HAM accumulated variables were divided by the total number of seconds of the output interval, that's why these "time"-variables are diagnosed as the sum of the length of the timestep in seconds - zdtime
, over all timesteps where a condition is fullfilled e.g. presence of liquid cloud water. However, in ICON the time accumulation is done differently and no division by the total number of seconds of the output interval is done automatically. Therefore instead of summing up zdtime
, in ICON-HAM, 1._dp
should be summed up over all timesteps the specific condition is fullfilled.
For example replace:
ztmp1(jl) = zdtime
pcloud_time(jl) = MERGE(ztmp1(jl), 0._dp, ld_liqcl(jl))
with:
ztmp1(jl) = 1._dp
pcloud_time(jl) = MERGE(ztmp1(jl), 0._dp, ld_liqcl(jl))
and similar for other "time"-variables.