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

Timestamp:

Jul 16, 2025, 3:25:48 PM (3 days ago)

Author:

jbclement

Message:

PEM:

• Making the computation of CO2 mass balance more robust, especially regarding 'CO2cond_ps'.
• Small correction about the dust tendency management for the layering algorithm.
• Small improvement of the visualization in "visu_evol_layering.py".
• File "log_launchPEM.txt" renamed into "launchPEM.log".

JBC

File:

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

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

@@ -291 +291 @@
 ! Loop variables
 integer :: i, l, ig, nnq, t, islope, ig_loop, islope_loop, isoil, icap
+real    :: totmass_ini
 logical :: num_str
 
@@ -995 +996 @@
 
     ! Checking mass balance for CO2
-    totmass_co2ice = 0.
-    totmass_atmco2 = 0.
-    do ig = 1,ngrid
-        totmass_atmco2 = totmass_atmco2 + cell_area(ig)*ps_avg(ig)/g
-        do islope = 1,nslope
-            totmass_co2ice = totmass_co2ice + co2_ice(ig,islope)*cell_area(ig)*subslope_dist(ig,islope)/cos(pi*def_slope_mean(islope)/180.)
-        enddo
-    enddo
-    write(*,'(a,f8.3,a)') " > Relative total CO2 mass balance = ", 100.*(totmass_atmco2 + totmass_co2ice + totmass_adsco2 - totmass_atmco2_ini - totmass_co2ice_ini - totmass_adsco2_ini)/(totmass_atmco2_ini + totmass_co2ice_ini + totmass_adsco2_ini), ' %'
-    if ((totmass_atmco2 + totmass_co2ice + totmass_adsco2 - totmass_atmco2_ini - totmass_co2ice_ini - totmass_adsco2_ini)/(totmass_atmco2_ini + totmass_co2ice_ini + totmass_adsco2_ini) > 0.01) then
-        write(*,*) '  /!\ Warning: mass balance is not conseved!'
-        write(*,'(a,f8.3,a)') '       Atmospheric CO2 mass balance = ', 100.*(totmass_atmco2 - totmass_atmco2_ini)/totmass_atmco2_ini, ' %'
-        write(*,'(a,f8.3,a)') '       CO2 ice mass balance         = ', 100.*(totmass_co2ice - totmass_co2ice_ini)/totmass_co2ice_ini, ' %'
-        write(*,'(a,f8.3,a)') '       Adsorbed CO2 mass balance    = ', 100.*(totmass_adsco2 - totmass_adsco2_ini)/totmass_adsco2_ini, ' %'
+    if (abs(CO2cond_ps - 1.) < 1.e-10) then
+        totmass_co2ice = 0.
+        totmass_atmco2 = 0.
+        do ig = 1,ngrid
+            totmass_atmco2 = totmass_atmco2 + cell_area(ig)*ps_avg(ig)/g
+            do islope = 1,nslope
+                totmass_co2ice = totmass_co2ice + co2_ice(ig,islope)*cell_area(ig)*subslope_dist(ig,islope)/cos(pi*def_slope_mean(islope)/180.)
+            enddo
+        enddo
+        totmass_ini = max(totmass_atmco2_ini + totmass_co2ice_ini + totmass_adsco2_ini,1.e-10)
+        write(*,'(a,f8.3,a)') " > Relative total CO2 mass balance = ", 100.*(totmass_atmco2 + totmass_co2ice + totmass_adsco2 - totmass_atmco2_ini - totmass_co2ice_ini - totmass_adsco2_ini)/totmass_ini, ' %'
+        if ((totmass_atmco2 + totmass_co2ice + totmass_adsco2 - totmass_atmco2_ini - totmass_co2ice_ini - totmass_adsco2_ini)/totmass_ini > 0.01) then
+            write(*,*) '  /!\ Warning: mass balance is not conseved!'
+            totmass_ini = max(totmass_atmco2_ini,1.e-10)
+            write(*,'(a,f8.3,a)') '       Atmospheric CO2 mass balance = ', 100.*(totmass_atmco2 - totmass_atmco2_ini)/totmass_ini, ' %'
+            totmass_ini = max(totmass_co2ice_ini,1.e-10)
+            write(*,'(a,f8.3,a)') '       CO2 ice mass balance         = ', 100.*(totmass_co2ice - totmass_co2ice_ini)/totmass_ini, ' %'
+            totmass_ini = max(totmass_adsco2_ini,1.e-10)
+            write(*,'(a,f8.3,a)') '       Adsorbed CO2 mass balance    = ', 100.*(totmass_adsco2 - totmass_adsco2_ini)/totmass_ini, ' %'
+        endif
     endif