Changeset 4160 for trunk/LMDZ.MARS/libf/dynphy_lonlat/phymars/newstart.F

Timestamp:

Mar 31, 2026, 3:38:07 PM (12 days ago)

Author:

jbclement

Message:

Mars PCM:
Giving standard explicit names to variables related to "paleoclimate" module + adding yearly average of flux exchanged with subsurface ice in XIOS output file intended for the PEM.
JBC

File:

• trunk/LMDZ.MARS/libf/dynphy_lonlat/phymars/newstart.F (modified) (6 diffs)

trunk/LMDZ.MARS/libf/dynphy_lonlat/phymars/newstart.F

@@ -60 +60 @@
       use comslope_mod, ONLY: nslope,def_slope,def_slope_mean,
      &             subslope_dist,end_comslope_h,ini_comslope_h
-      use paleoclimate_mod, only: h2o_ice_depth, lag_co2_ice, d_coef,
-     &             ini_paleoclimate_h, end_paleoclimate_h
+      use paleoclimate_mod, only: h2oice_depth, co2ice_depth, coef_ssdif,
+     &             ini_paleoclimate, end_paleoclimate
       use subslope_mola_mod, ONLY: subslope_mola
       use parse_args_mod, only: parse_args
@@ -208 +208 @@
       real, allocatable, dimension(:,:,:) :: albedo_old_slope        ! Surface albedo in each solar band
       real, allocatable, dimension(:,:)   :: flux_geo_old_slope
-      real, allocatable, dimension(:,:)   :: h2o_ice_depth_old_slope
-      real, allocatable, dimension(:,:)   :: lag_co2_ice_old_slope
-      real, allocatable, dimension(:,:)   :: d_coef_old_slope
+      real, allocatable, dimension(:,:)   :: h2oice_depth_old_slope
+      real, allocatable, dimension(:,:)   :: co2ice_depth_old_slope
+      real, allocatable, dimension(:,:)   :: coef_ssdif_old_slope
       integer :: iflat, nslope_old, nslope_new
 
@@ -1850 +1850 @@
 
 ! Paleoclimate related stuff
-        allocate(h2o_ice_depth_old_slope(ngridmx,nslope_old))
-        allocate(lag_co2_ice_old_slope(ngridmx,nslope_old))
-        allocate(d_coef_old_slope(ngridmx,nslope_old))
-        h2o_ice_depth_old_slope = h2o_ice_depth
-        lag_co2_ice_old_slope = lag_co2_ice
-        d_coef_old_slope = d_coef
-        call end_paleoclimate_h
-        call ini_paleoclimate_h(ngridmx,nslope_new)
+        allocate(h2oice_depth_old_slope(ngridmx,nslope_old))
+        allocate(co2ice_depth_old_slope(ngridmx,nslope_old))
+        allocate(coef_ssdif_old_slope(ngridmx,nslope_old))
+        h2oice_depth_old_slope = h2oice_depth
+        co2ice_depth_old_slope = co2ice_depth
+        coef_ssdif_old_slope = coef_ssdif
+        call end_paleoclimate
+        call ini_paleoclimate(ngridmx,nslope_new)
 
 ! Update of the variables with the new "slopes" situation according to the old one
@@ -1872 +1872 @@
              inertiesoil(:,:,islope) = inertiesoil_old_slope(:,:,1)
              flux_geo(:,islope) = flux_geo_old_slope(:,1)
-             h2o_ice_depth(:,islope) = h2o_ice_depth_old_slope(:,1)
-             lag_co2_ice(:,islope) = lag_co2_ice_old_slope(:,1)
-             d_coef(:,islope) = d_coef_old_slope(:,1)
+             h2oice_depth(:,islope) = h2oice_depth_old_slope(:,1)
+             co2ice_depth(:,islope) = co2ice_depth_old_slope(:,1)
+             coef_ssdif(:,islope) = coef_ssdif_old_slope(:,1)
           enddo
         else
@@ -1888 +1888 @@
              inertiesoil(:,:,islope) = inertiesoil_old_slope(:,:,iflat)
              flux_geo(:,islope) = flux_geo_old_slope(:,iflat)
-             h2o_ice_depth(:,islope) = h2o_ice_depth_old_slope(:,iflat)
-             lag_co2_ice(:,islope) = lag_co2_ice_old_slope(:,iflat)
-             d_coef(:,islope) = d_coef_old_slope(:,iflat)
+             h2oice_depth(:,islope) = h2oice_depth_old_slope(:,iflat)
+             co2ice_depth(:,islope) = co2ice_depth_old_slope(:,iflat)
+             coef_ssdif(:,islope) = coef_ssdif_old_slope(:,iflat)
           enddo
         endif
@@ -1898 +1898 @@
         deallocate(perennial_co2_old_slope,tsoil_old_slope)
         deallocate(inertiesoil_old_slope,flux_geo_old_slope)
-        deallocate(albedo_old_slope,h2o_ice_depth_old_slope)
-        deallocate(lag_co2_ice_old_slope,d_coef_old_slope)
+        deallocate(albedo_old_slope,h2oice_depth_old_slope)
+        deallocate(co2ice_depth_old_slope,coef_ssdif_old_slope)
 
         endif !nslope=nslope_new