Changeset 3384 for trunk/LMDZ.COMMON/libf

Timestamp:

Jun 18, 2024, 5:42:24 PM (6 months ago)

Author:

jbclement

Message:

PEM:
Correction related to r3330 for the initialization of the global average pressure taken from the PCM: it has to be done before the call of "pemetat0".
JBC

Location:

trunk/LMDZ.COMMON/libf/evolution

Files:

• changelog.txt (modified) (1 diff)
• pem.F90 (modified) (4 diffs)

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

@@ -371 +371 @@
 
 == 17/06/2024 == JBC
-Small correction for the initialization value of dtphys + Corrections for the launching script.
+Small correction for the initialization value of 'dtphys' + Corrections for the launching script.
+
+== 18/06/2024 == JBC
+Correction related to r3330 for the initialization of the global average pressure taken from the PCM: it has to be done before the call of "pemetat0".

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

@@ -7 +7 @@
 !    I_e Initialization of the PEM variable and soil
 !    I_f Compute tendencies
-!    I_g Read the "startpem.nc"
-!    I_h Save the initial situation
+!    I_g Save the initial situation
+!    I_h Read the "startpem.nc"
 !    I_i Compute orbit criterion
 
@@ -598 +598 @@
 !------------------------
 ! I   Initialization
-!    I_g Read the "startpem.nc"
+!    I_g Save the initial situation
+!------------------------
+allocate(zplev_PCM(ngrid,nlayer + 1))
+Total_surface = 0.
+do ig = 1,ngrid
+    Total_surface = Total_surface + cell_area(ig)
+    zplev_PCM(ig,:) = ap + bp*ps_start_PCM(ig)
+enddo
+global_avg_press_old = sum(cell_area*ps_start_PCM)/Total_surface
+global_avg_press_PCM = global_avg_press_old
+global_avg_press_new = global_avg_press_old
+write(*,*) "Total surface of the planet =", Total_surface
+write(*,*) "Initial global average pressure =", global_avg_press_PCM
+
+!------------------------
+! I   Initialization
+!    I_h Read the "startpem.nc"
 !------------------------
 allocate(co2_ice(ngrid,nslope),h2o_ice(ngrid,nslope))
@@ -616 +632 @@
               watersurf_density_avg,watersoil_density_PEM_avg,co2_adsorbded_phys,delta_co2_adsorbded,                &
               h2o_adsorbded_phys,delta_h2o_adsorbded,stratif)
+
+! We save the places where h2o ice is sublimating
+! We compute the surface of h2o ice sublimating
+allocate(ini_co2ice_sublim(ngrid,nslope),ini_h2oice_sublim(ngrid,nslope),is_co2ice_ini(ngrid,nslope))
+co2ice_ini_surf = 0.
+h2oice_ini_surf = 0.
+ini_co2ice_sublim = .false.
+ini_h2oice_sublim = .false.
+is_co2ice_ini = .false.
+co2ice_disappeared = .false.
+do i = 1,ngrid
+    do islope = 1,nslope
+        if (co2_ice(i,islope) > 0.) is_co2ice_ini(i,islope) = .true.
+        if (tend_co2_ice(i,islope) < 0. .and. co2_ice(i,islope) > 0.) then
+            ini_co2ice_sublim(i,islope) = .true.
+            co2ice_ini_surf = co2ice_ini_surf + cell_area(i)*subslope_dist(i,islope)
+        endif
+        if (tend_h2o_ice(i,islope) < 0. .and. h2o_ice(i,islope) > 0.) then
+            ini_h2oice_sublim(i,islope) = .true.
+            h2oice_ini_surf = h2oice_ini_surf + cell_area(i)*subslope_dist(i,islope)
+        endif
+    enddo
+enddo
+write(*,*) "Total initial surface of co2 ice sublimating (slope) =", co2ice_ini_surf
+write(*,*) "Total initial surface of h2o ice sublimating (slope) =", h2oice_ini_surf
 
 delta_h2o_icetablesublim = 0.
@@ -643 +684 @@
 endif ! adsorption
 deallocate(tsurf_avg_yr1)
-
-!------------------------
-! I   Initialization
-!    I_h Save the initial situation
-!------------------------
-! We save the places where h2o ice is sublimating
-! We compute the surface of h2o ice sublimating
-allocate(ini_co2ice_sublim(ngrid,nslope),ini_h2oice_sublim(ngrid,nslope),is_co2ice_ini(ngrid,nslope))
-co2ice_ini_surf = 0.
-h2oice_ini_surf = 0.
-Total_surface = 0.
-ini_co2ice_sublim = .false.
-ini_h2oice_sublim = .false.
-is_co2ice_ini = .false.
-co2ice_disappeared = .false.
-do i = 1,ngrid
-    Total_surface = Total_surface + cell_area(i)
-    do islope = 1,nslope
-        if (co2_ice(i,islope) > 0.) is_co2ice_ini(i,islope) = .true.
-        if (tend_co2_ice(i,islope) < 0. .and. co2_ice(i,islope) > 0.) then
-            ini_co2ice_sublim(i,islope) = .true.
-            co2ice_ini_surf = co2ice_ini_surf + cell_area(i)*subslope_dist(i,islope)
-        endif
-        if (tend_h2o_ice(i,islope) < 0. .and. h2o_ice(i,islope) > 0.) then
-            ini_h2oice_sublim(i,islope) = .true.
-            h2oice_ini_surf = h2oice_ini_surf + cell_area(i)*subslope_dist(i,islope)
-        endif
-    enddo
-enddo
-
-write(*,*) "Total initial surface of co2 ice sublimating (slope) =", co2ice_ini_surf
-write(*,*) "Total initial surface of h2o ice sublimating (slope) =", h2oice_ini_surf
-write(*,*) "Total surface of the planet =", Total_surface
-allocate(zplev_PCM(ngrid,nlayer + 1))
-
-do ig = 1,ngrid
-    zplev_PCM(ig,:) = ap + bp*ps_start_PCM(ig)
-enddo
-
-global_avg_press_old = sum(cell_area*ps_start_PCM)/Total_surface
-global_avg_press_PCM = global_avg_press_old
-global_avg_press_new = global_avg_press_old
-write(*,*) "Initial global average pressure =", global_avg_press_PCM
 
 !------------------------