Changeset 3086 for LMDZ6/trunk/libf/phylmd/regr_pr_comb_coefoz_m.F90

Timestamp:

Nov 22, 2017, 11:01:27 PM (7 years ago)

Author:

dcugnet

Message:

Modification of the ozone interpolation routine to ensure more realistic
profiles, in particular when the tropopause adjustment is activated, as
well as an improved accuracy on the chemical tropopause of streched ozone
field:

• stretching function has no longer a trapezoidal shape, but is a hat

function, which top is at the target (ldmz) tropopause, with lower /
upper bounds respectively lower / higher than the lowest / highest
tropopause (lmdz and file). This helps to avoid stretching the profile
too far from the strictly required zone.

• revert to a classical linear vertical interpolation on log-pressure

because the stairs-shape effect of the conservative interpolation can not
be avoided, even at second order method. Note that the second order
conservative method is still used for Cariolle linearized ozone
coefficients interpolation: the regr_pr_time routine now switches from
one method to the other depending on flag "Ploc" (pressures must be given
at layers interfaces - Ploc=='I' - for conservative method and layers
centers - Ploc=='C' - for the regular one).

• force the stretching function to be discretized on at least three

levels to avoid the effect of the stretching to be reduced to a single
point shifting of the ozone profile when discretized.

File:

• LMDZ6/trunk/libf/phylmd/regr_pr_comb_coefoz_m.F90 (modified) (1 diff)

LMDZ6/trunk/libf/phylmd/regr_pr_comb_coefoz_m.F90

@@ -129 +129 @@
 
     call regr_pr_time_av(ncid, (/"a2           ", "a4           ", "a6           ", &
-         "P_net_Mob    ", "r_Mob        ", "temp_Mob     ", "R_Het        "/), REAL(julien-1), &
-         press_in_edg, paprs, coefoz)
+         "P_net_Mob    ", "r_Mob        ", "temp_Mob     ", "R_Het        "/),      &
+         REAL(julien-1), 'I', press_in_edg, paprs, coefoz)
     a2 = coefoz(:, :, 1)
     a4_mass = coefoz(:, :, 2) * 48. / 29.