Changeset 3608

Timestamp:

Feb 3, 2025, 3:58:30 PM (6 hours ago)

Author:

jbclement

Message:

PEM:
Rescaling the pressure deviation to avoid disproportionate oscillations in case of huge average pressure drop.
JBC

Location:

trunk/LMDZ.COMMON/libf/evolution

Files:

• changelog.txt (modified) (1 diff)
• pem.F90 (modified) (1 diff)
• soil_thermalproperties_mod.F90 (modified) (1 diff)

trunk/LMDZ.COMMON/libf/evolution/changelog.txt

@@ -566 +566 @@
 - Bug correction about the pressure/teta reconstruction for the PCM (mismatch between the physical and dynamical grid).
 - Improvement of messages giving info about the PEM workflow in the terminal.
+
+== 03/02/2025 == JBC
+Rescaling the pressure deviation to avoid disproportionate oscillations in case of huge average pressure drop.

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

@@ -1088 +1088 @@
 write(*,*) '> Reconstructing the pressure for the PCM'
 allocate(ps_start(ngrid))
-ps_start = ps_avg + ps_dev
+! The pressure deviation is rescaled to avoid disproportionate oscillations in case of huge average pressure drop
+ps_start = ps_avg + ps_dev*ps_avg_global_new/ps_avg_global_ini
 deallocate(ps_avg,ps_dev)
 

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

@@ -94 +94 @@
 real    :: ice_inertia                            ! Inertia of water ice [SI]
 
-write(*,*) "Update soil properties"
+write(*,*) "> Updating soil properties"
 
 ! 1. Modification of the regolith thermal inertia.