Changeset 3611 for trunk/LMDZ.COMMON

Timestamp:

Feb 7, 2025, 12:14:55 PM (7 days ago)

Author:

jbclement

Message:

PEM:
Few bug corrections:

• global surface pressure was badlu updated (r3610);
• XIOS slope variables were not consistent;
• missing initialization for slope number (r3609);
• taking into account that there may be no "startfi.nc" at the beginning of the PEM simulation.

JBC

Location:

trunk/LMDZ.COMMON/libf/evolution

Files:

• changelog.txt (modified) (1 diff)
• deftank/lib_launchPEM.sh (modified) (1 diff)
• pem.F90 (modified) (3 diffs)
• read_data_PCM_mod.F90 (modified) (2 diffs)

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

@@ -577 +577 @@
 - The sublimation flux can now be modified by the new function 'recomp_tend_h2o' to account for the growth of a dust lag layer (see Eran Vos's note for the formula).
 - Addition of a function 'itp_tsoil' to compute the soil temperature at any point in the soil profile.
+
+== 07/02/2025 == JBC
+Few bug corrections:
+    - global surface pressure was badlu updated (r3610);
+    - XIOS slope variables were not consistent;
+    - missing initialization for slope number (r3609);
+    - taking into account that there may be no "startfi.nc" at the beginning of the PEM simulation.

trunk/LMDZ.COMMON/libf/evolution/deftank/lib_launchPEM.sh

@@ -144 +144 @@
     iPEM=1
     iPCM=1
-    cp startfi.nc starts/
+    if [ -f "startfi.nc" ]; then
+        cp startfi.nc starts/
+    fi
     if [ -f "start.nc" ]; then
         cp start.nc starts/

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

@@ -727 +727 @@
     do i = 1,ngrid
         do islope = 1,nslope
-            ps_avg_global_new = ps_avg_global_old - CO2cond_ps*g*cell_area(i)*d_co2ice(i,islope)*subslope_dist(i,islope)/cos(pi*def_slope_mean(islope)/180.)/total_surface
+            ps_avg_global_new = ps_avg_global_new - CO2cond_ps*g*cell_area(i)*d_co2ice(i,islope)*subslope_dist(i,islope)/cos(pi*def_slope_mean(islope)/180.)/total_surface
         enddo
     enddo
@@ -733 +733 @@
     if (adsorption_pem) then
         do i = 1,ngrid
-            ps_avg_global_new = ps_avg_global_old - g*cell_area(i)*delta_co2_adsorbed(i)/total_surface
+            ps_avg_global_new = ps_avg_global_new - g*cell_area(i)*delta_co2_adsorbed(i)/total_surface
         enddo
     endif
@@ -1251 +1251 @@
 call info_PEM(i_myear_leg,stopPEM,i_myear,n_myear)
 
-write(*,*) ">> The PEM has run for", i_myear_leg, "Martian years."
-write(*,*) ">> The chained simulation has run for", i_myear, "Martian years =", i_myear*convert_years, "Earth years."
-write(*,*) ">> The reached date is now", (year_bp_ini + i_myear)*convert_years, "Earth years."
-write(*,*) ">> PEM: so far, so good!"
+write(*,*) '************************************************************************'
+write(*,*) "> The PEM has run for", i_myear_leg, "Martian years."
+write(*,*) "> The chained simulation has run for", i_myear, "Martian years =", i_myear*convert_years, "Earth years."
+write(*,*) "> The reached date is now", (year_bp_ini + i_myear)*convert_years, "Earth years."
+write(*,*) "> PEM: so far, so good!"
+write(*,*) '************************************************************************'
 
 if (layering_algo) then

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

@@ -76 +76 @@
 do islope = 1,nslope
     if (nslope /= 1) then
+        num='  '
         write(num,'(i2.2)') islope
         num = '_slope'//num
@@ -196 +197 @@
 do islope = 1,nslope
     if (nslope /= 1) then
+        num='  '
         write(num,'(i2.2)') islope
         num = '_slope'//num