Changeset 3850 for trunk/LMDZ.COMMON/libf/evolution/pemetat0.F90

Timestamp:

Jul 15, 2025, 2:43:41 PM (11 days ago)

Author:

jbclement

Message:

PEM:

• Handling the case where the new date is positive (present-day climate) in "recomp_orb_param_mod.F90", especially because of rounding errors.
• Correction in "visu_evol_layering.py" to get the proper extent of the stratication + replacing 'imshow' by 'pcolormesh' to better process the data.
• Addition of a graphic to plot the changes in obliquity directly onto the evolution of the stratification over time for better analysis. Obliquity curve is colored according to stratification gain/loss.
• Correction of numbering in "lib_launchPEM.sh" when requesting a relaunch.
• Update of "PEMrun.job" due to new parsing of arguments introduced in r3836/r3837.

JBC

File:

• trunk/LMDZ.COMMON/libf/evolution/pemetat0.F90 (modified) (1 diff)

trunk/LMDZ.COMMON/libf/evolution/pemetat0.F90

@@ -393 +393 @@
                 do islope = 1,nslope
                     call ini_layering(layerings_map(ig,islope))
-                    print*, 'coucou', 1.05*ini_huge_h2oice/rho_h2oice, layerings_map(ig,islope)%top%top_elevation
                     call add_stratum(layerings_map(ig,islope),1.05*ini_huge_h2oice/rho_h2oice,0.,ini_huge_h2oice/rho_h2oice,0.05*ini_huge_h2oice/rho_h2oice,0.,0.)
                 enddo