Changeset 2242 for trunk/LMDZ.TITAN/README

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/README (modified) (1 diff)

trunk/LMDZ.TITAN/README

@@ -1563 +1563 @@
  but would require to be able to have writediag in 5D)
  EDIT (22/05/19) : To have correct calculations, output is now exp(-tau), you need to postproc it to have tau.
+
+== 17/02/2020 == JVO (r2242)
+- 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.
+- TBD ? Add a sanity check within YAMMS filtering high unphysical radii ?
+- 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.