Changeset 3174

Timestamp:

Jan 11, 2024, 3:15:26 PM (11 months ago)

Author:

jbclement

Message:

Mars PCM:
Correction of a bug: 'perennial_co2ice' was not initialized if there was no "startfi.nc" file.
JBC

Location:

trunk/LMDZ.MARS

Files:

• changelog.txt (modified) (1 diff)
• libf/phymars/dyn1d/init_testphys1d_mod.F90 (modified) (3 diffs)
• libf/phymars/phyetat0_mod.F90 (modified) (2 diffs)

trunk/LMDZ.MARS/changelog.txt

@@ -4430 +4430 @@
 Minor fix for XIOS output of a scalar in write_output_mod.F90.
 Also added more helpful error messages in xios_output_mod.F90.
+
+== 11/01/2024 == JBC
+Correction of a bug: 'perennial_co2ice' was not initialized if there was no "startfi.nc" file.

trunk/LMDZ.MARS/libf/phymars/dyn1d/init_testphys1d_mod.F90

@@ -15 +15 @@
 use surfdat_h,                only: albedodat, z0_default, z0, emissiv, emisice, albedice, iceradius, dtemisice,      &
                                     zmea, zstd, zsig, zgam, zthe, hmons, summit, base, phisfi, watercaptag, watercap, &
-                                    tsurf, emis, qsurf
+                                    tsurf, emis, qsurf, perennial_co2ice
 use infotrac,                 only: nqtot, tname, nqperes, nqfils
 use read_profile_mod,         only: read_profile
@@ -40 +40 @@
 use nonoro_gwd_ran_mod,       only: du_nonoro_gwd, dv_nonoro_gwd
 use conf_phys_mod,            only: conf_phys
+use paleoclimate_mod,         only: paleoclimate
 ! Mostly for XIOS outputs:
 use mod_const_mpi,            only: COMM_LMDZ
@@ -652 +653 @@
 write(*,*) " watercaptag = ",watercaptag
 
+! Initialize perennial_co2ice
+! ---------------------------
+if (.not. therestartfi .and. paleoclimate) perennial_co2ice = 0.
+
 ! Check if the atmospheric water profile is specified
 ! ---------------------------------------------------

trunk/LMDZ.MARS/libf/phymars/phyetat0_mod.F90

@@ -853 +853 @@
 
 ! Perennial CO2 ice
-    perennial_co2ice(:,:) = 0.
+    perennial_co2ice = 0.
     call get_field("perennial_co2ice",perennial_co2ice,found,indextime)
     if (.not. found) then
@@ -859 +859 @@
                    "<perennial_co2ice> is set as 10m at the South Pole"
         if (abs(latitude(ngrid) - (-pi/2.)) < 1.e-5) then
-            do islope = 1,nslope
-                perennial_co2ice(ngrid,islope) = 10*1.6e3 ! 10m which is convert to kg/m^2
-               qsurf(ngrid,igcm_co2_tmp,islope) = qsurf(ngrid - 1,igcm_co2_tmp,islope) + perennial_co2ice(ngrid,islope) ! perennial ice + frost
-            enddo
+            perennial_co2ice(ngrid,:) = 10*1.6e3 ! 10m which is convert to kg/m^2
+            qsurf(ngrid,igcm_co2_tmp,:) = qsurf(ngrid - 1,igcm_co2_tmp,:) + perennial_co2ice(ngrid,:) ! perennial ice + frost
         endif
     endif ! not found
   else ! no startfiphyle
-    h2o_ice_depth(:,:) = -1.
-    lag_co2_ice(:,:) = -1.
-    d_coef(:,:)= 4.e-4
-    perennial_co2ice(:,:) = 0.
+    h2o_ice_depth = -1.
+    lag_co2_ice = -1.
+    d_coef = 4.e-4
+    perennial_co2ice = 0.
     if (abs(latitude(ngrid) - (-pi/2.)) < 1.e-5) then
-        do islope = 1,nslope
-            perennial_co2ice(ngrid,islope) = 10*1.6e3 ! 10m which is convert to kg/m^2
-            qsurf(ngrid,igcm_co2_tmp,islope) = qsurf(ngrid - 1,igcm_co2_tmp,islope) + perennial_co2ice(ngrid,islope) ! perennial ice + frost
-        enddo
+        perennial_co2ice(ngrid,:) = 10*1.6e3 ! 10m which is convert to kg/m^2
+        qsurf(ngrid,igcm_co2_tmp,:) = qsurf(ngrid - 1,igcm_co2_tmp,:) + perennial_co2ice(ngrid,:) ! perennial ice + frost
     endif
   endif !startphy_file
 else
-   h2o_ice_depth(:,:) = -1.
-   lag_co2_ice(:,:) = -1.
-   d_coef(:,:)= 4.e-4
-   perennial_co2ice(:,:) = 0.
+   h2o_ice_depth = -1.
+   lag_co2_ice = -1.
+   d_coef = 4.e-4
+   perennial_co2ice = 0.
 endif !paleoclimate