Changeset 3498 for trunk/LMDZ.COMMON/libf/evolution/evol_ice_mod.F90

Timestamp:

Nov 7, 2024, 2:48:08 PM (2 weeks ago)

Author:

jbclement

Message:

PEM:

• Correction of the variable name for the ice table depth in "pemetat0.F90". So it is now got as intended from the "startpem.nc" file;
• Renaming of the tendencies in the PEM with the prefix 'd_' instead of 'tend_';
• Modification of the PEM time step type from integer to real. As a consequence, all time variables are now of real type. This change adds the possibility to consider fractions of year as time step.

JBC

File:

• trunk/LMDZ.COMMON/libf/evolution/evol_ice_mod.F90 (modified) (12 diffs)

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

@@ -7 +7 @@
 !=======================================================================
 
-SUBROUTINE evol_co2_ice(ngrid,nslope,co2_ice,tend_co2_ice)
+SUBROUTINE evol_co2_ice(ngrid,nslope,co2_ice,d_co2ice)
 
-use time_evol_mod, only: dt_pem
+use time_evol_mod, only: dt
 
 implicit none
@@ -25 +25 @@
 !   OUTPUT
 real, dimension(ngrid,nslope), intent(inout) :: co2_ice      ! Previous and actual density of CO2 ice
-real, dimension(ngrid,nslope), intent(inout) :: tend_co2_ice ! Evolution of perennial ice over one year
+real, dimension(ngrid,nslope), intent(inout) :: d_co2ice ! Evolution of perennial ice over one year
 
 !   local:
@@ -35 +35 @@
 
 co2_ice_old = co2_ice
-co2_ice = co2_ice + tend_co2_ice*dt_pem
+co2_ice = co2_ice + d_co2ice*dt
 where (co2_ice < 0.)
     co2_ice = 0.
-    tend_co2_ice = -co2_ice_old/dt_pem
+    d_co2ice = -co2_ice_old/dt
 end where
 
@@ -45 +45 @@
 !=======================================================================
 
-SUBROUTINE evol_h2o_ice(ngrid,nslope,cell_area,delta_h2o_adsorbded,delta_h2o_icetablesublim,h2o_ice,tend_h2o_ice,stopPEM)
+SUBROUTINE evol_h2o_ice(ngrid,nslope,cell_area,delta_h2o_adsorbded,delta_h2o_icetablesublim,h2o_ice,d_h2oice,stopPEM)
 
-use time_evol_mod,          only: dt_pem
+use time_evol_mod,          only: dt
 use comslope_mod,           only: subslope_dist, def_slope_mean
 
@@ -73 +73 @@
 !   OUTPUT
 real, dimension(ngrid,nslope), intent(inout) :: h2o_ice      ! Previous and actual density of h2o ice (kg/m^2)
-real, dimension(ngrid,nslope), intent(inout) :: tend_h2o_ice ! Evolution of perennial ice over one year (kg/m^2/year)
+real, dimension(ngrid,nslope), intent(inout) :: d_h2oice ! Evolution of perennial ice over one year (kg/m^2/year)
 integer,                       intent(inout) :: stopPEM      ! Stopping criterion code
 
@@ -80 +80 @@
 integer                       :: i, islope                                           ! Loop variables
 real                          :: pos_tend, neg_tend, real_coefficient, negative_part ! Variable to conserve h2o
-real, dimension(ngrid,nslope) :: new_tendencies                                      ! Tendencies computed in order to conserve h2o ice on the surface, only exchange between surface are done
+real, dimension(ngrid,nslope) :: new_tend                                      ! Tendencies computed in order to conserve h2o ice on the surface, only exchange between surface are done
 !=======================================================================
 if (ngrid /= 1) then ! Not in 1D
@@ -99 +99 @@
         do islope = 1,nslope
             if (h2o_ice(i,islope) > 0.) then
-                if (tend_h2o_ice(i,islope) > 0.) then
-                    pos_tend = pos_tend + tend_h2o_ice(i,islope)*cell_area(i)*subslope_dist(i,islope)/cos(def_slope_mean(islope)*pi/180.)
+                if (d_h2oice(i,islope) > 0.) then
+                    pos_tend = pos_tend + d_h2oice(i,islope)*cell_area(i)*subslope_dist(i,islope)/cos(def_slope_mean(islope)*pi/180.)
                 else
-                    neg_tend = neg_tend - tend_h2o_ice(i,islope)*cell_area(i)*subslope_dist(i,islope)/cos(def_slope_mean(islope)*pi/180.)
+                    neg_tend = neg_tend - d_h2oice(i,islope)*cell_area(i)*subslope_dist(i,islope)/cos(def_slope_mean(islope)*pi/180.)
                 endif
             endif
@@ -108 +108 @@
     enddo
 
-    new_tendencies = 0.
+    new_tend = 0.
     if (abs(pos_tend) < 1.e-10 .or. abs(neg_tend) < 1.e-10) then
         write(*,*) "Reason of stopping: there is no sublimating or condensing h2o ice!"
@@ -118 +118 @@
         ! We adapt the tendencies to conserve h2o and do only exchange between grid points
         if (neg_tend > pos_tend .and. pos_tend > 0.) then ! More sublimation on the planet than condensation
-            where (tend_h2o_ice < 0.) ! We lower the sublimating rate by a coefficient
-                new_tendencies = tend_h2o_ice*pos_tend/neg_tend
+            where (d_h2oice < 0.) ! We lower the sublimating rate by a coefficient
+                new_tend = d_h2oice*pos_tend/neg_tend
             elsewhere ! We don't change the condensing rate
-                new_tendencies = tend_h2o_ice
+                new_tend = d_h2oice
             end where
         else if (neg_tend < pos_tend .and. neg_tend > 0.) then ! More condensation on the planet than sublimation
-            where (tend_h2o_ice < 0.) ! We don't change the sublimating rate
-                new_tendencies = tend_h2o_ice
+            where (d_h2oice < 0.) ! We don't change the sublimating rate
+                new_tend = d_h2oice
             elsewhere ! We lower the condensing rate by a coefficient
-                new_tendencies = tend_h2o_ice*neg_tend/pos_tend
+                new_tend = d_h2oice*neg_tend/pos_tend
             end where
         endif
@@ -133 +133 @@
 
     ! Evolution of the h2o ice for each physical point
-    h2o_ice = h2o_ice + new_tendencies*dt_pem
+    h2o_ice = h2o_ice + new_tend*dt
 
     ! We compute the amount of h2o that is sublimated in excess
@@ -142 +142 @@
                 negative_part = negative_part - h2o_ice(i,islope)*cell_area(i)*subslope_dist(i,islope)/cos(def_slope_mean(islope)*pi/180.)
                 h2o_ice(i,islope) = 0.
-                tend_h2o_ice(i,islope) = 0.
+                d_h2oice(i,islope) = 0.
             endif
         enddo
@@ -157 +157 @@
     do islope = 1,nslope
         do i = 1,ngrid
-             if (new_tendencies(i,islope) > 0.) h2o_ice(i,islope) = h2o_ice(i,islope) - new_tendencies(i,islope)*real_coefficient*dt_pem*cos(def_slope_mean(islope)*pi/180.)
+             if (new_tend(i,islope) > 0.) h2o_ice(i,islope) = h2o_ice(i,islope) - new_tend(i,islope)*real_coefficient*dt*cos(def_slope_mean(islope)*pi/180.)
         enddo
      enddo
 else ! ngrid == 1, i.e. in 1D
-    h2o_ice = h2o_ice + tend_h2o_ice*dt_pem
+    h2o_ice = h2o_ice + d_h2oice*dt
 endif