Changeset 2242 for trunk/LMDZ.TITAN/libf/muphytitan

Timestamp:

Feb 17, 2020, 5:44:12 PM (5 years ago)

Author:

jvatant

Message:

+ Update the interface with YAMMS so it now correctly handles the small values of the moments,
requiring dynamics to have a threshold quite low (set to 1D-200) and a sanity check in calmufi
corresponding to this value. Thus it removes 'most' of the unphysical radius obtained in
YAMMS. There are still some, but at least there is no more problem of model stability for the moments
and the code can run.
Still, take care the day you want to calculate opacities from the radii and not the moments.
Although, note that there are some negative values, in the output files, but theses negatives are

harmless, as they are only present in output files, dynamics reseting to epsilon after.

+ Given theses pbs with the radii, also update the optics so that it computes the opacities in
a simpler way, directly for M3, through look-up tables, M3 being a good proxy.
--JVO

File:

• trunk/LMDZ.TITAN/libf/muphytitan/mm_globals.f90 (modified) (1 diff)

trunk/LMDZ.TITAN/libf/muphytitan/mm_globals.f90

@@ -1066 +1066 @@
     ! il faudrait peut etre revoir la gestion des zeros ici et la 
     ! remplacer par une valeur seuil des moments.
+    !
+    !-> JVO 19 : Done. Zero threshold set at espilon from dynamics on the
+    ! input moments in calmufi (safer than here). Might still be some unphysical
+    ! values after the dynamics near the threshold. Could be a could idea to add
+    ! a sanity check filtering too high values of radii.
+    !
+    ! TBD : Add a sanity check for high radii ????
     WHERE(mm_m3aer_s > 0._mm_wp .AND. mm_m0aer_s > 0._mm_wp)
       mm_rcs = mm_get_rcs(mm_m0aer_s,mm_m3aer_s)