Changeset 4147 for trunk/LMDZ.COMMON/libf/evolution/pem.F90

Timestamp:

Mar 23, 2026, 2:55:10 PM (3 weeks ago)

Author:

jbclement

Message:

PEM:

• Simplification of subroutines to convert data between the physical and the dynamical/lon-lat grids + making them more robust.
• Correction for air mass to give back to the PCM. The variable is extensive so poles must be treated specifically.
• Making the PEM able to do 0 year.
• Explicit information about the frost values computed by the PEM + enforcing positivity of yearly minima.

JBC

File:

• trunk/LMDZ.COMMON/libf/evolution/pem.F90 (modified) (2 diffs)

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

@@ -19 +19 @@
 ! ------------
 ! Common modules
+use job_mod,        only: timelimit, antetime, timewall
 use parse_args_mod, only: parse_args
-use job_mod,        only: timelimit, antetime, timewall
 ! PEM modules
 use allocation,         only: ini_allocation, end_allocation
@@ -399 +399 @@
         end do
         totmass_ini = max(totmass_atmco2_ini + totmass_co2ice_ini + totmass_adsco2_ini,minieps_qp) ! To avoid division by 0
-        call print_msg(" > Relative total CO2 mass balance = "//real2str(100._qp*(totmass_atmco2 + totmass_co2ice + totmass_adsco2 - totmass_atmco2_ini - totmass_co2ice_ini - totmass_adsco2_ini)/totmass_ini)//' %',LVL_NFO)
+        call print_msg('> Relative total CO2 mass balance = '//real2str(100._qp*(totmass_atmco2 + totmass_co2ice + totmass_adsco2 - totmass_atmco2_ini - totmass_co2ice_ini - totmass_adsco2_ini)/totmass_ini)//' %',LVL_NFO)
         if (abs((totmass_atmco2 + totmass_co2ice + totmass_adsco2 - totmass_atmco2_ini - totmass_co2ice_ini - totmass_adsco2_ini)/totmass_ini) > 0.01_qp) then
             call print_msg('Mass balance is not conserved!',LVL_WRN)