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

Timestamp:

Nov 7, 2024, 2:48:08 PM (3 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/recomp_tend_co2_slope_mod.F90 (modified) (4 diffs)

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

@@ -7 +7 @@
 !=======================================================================
 
-SUBROUTINE recomp_tend_co2_slope(ngrid,nslope,timelen,tendencies_co2_ice_phys,tendencies_co2_ice_phys_ini,co2ice_slope,emissivity_slope, &
+SUBROUTINE recomp_tend_co2_slope(ngrid,nslope,timelen,d_co2ice_phys,d_co2ice_phys_ini,co2ice_slope,emissivity_slope, &
                                  vmr_co2_PCM,vmr_co2_pem,ps_PCM_2,global_avg_press_PCM,global_avg_press_new)
 
@@ -29 +29 @@
 real,                           intent(in) :: global_avg_press_PCM        ! global averaged pressure at previous timestep
 real,                           intent(in) :: global_avg_press_new        ! global averaged pressure at current timestep
-real, dimension(ngrid,nslope),  intent(in) :: tendencies_co2_ice_phys_ini ! physical point field: Evolution of perennial ice over one year
+real, dimension(ngrid,nslope),  intent(in) :: d_co2ice_phys_ini ! physical point field: Evolution of perennial ice over one year
 real, dimension(ngrid,nslope),  intent(in) :: co2ice_slope                ! CO2 ice per mesh and sub-grid slope (kg/m^2)
 real, dimension(ngrid,nslope),  intent(in) :: emissivity_slope            ! Emissivity per mesh and sub-grid slope(1)
 !   OUTPUT
-real, dimension(ngrid,nslope), intent(inout) ::  tendencies_co2_ice_phys ! physical point field: Evolution of perennial ice over one year
+real, dimension(ngrid,nslope), intent(inout) ::  d_co2ice_phys ! physical point field: Evolution of perennial ice over one year
 
 !   local:
@@ -47 +47 @@
         coef = sols_per_my*sec_per_sol*emissivity_slope(i,islope)*sigmaB/Lco2
         ave = 0.
-        if (co2ice_slope(i,islope) > 1.e-4 .and. abs(tendencies_co2_ice_phys(i,islope)) > 1.e-5) then
+        if (co2ice_slope(i,islope) > 1.e-4 .and. abs(d_co2ice_phys(i,islope)) > 1.e-5) then
             do t=1,timelen
                 ave = ave + (beta_clap_co2/(alpha_clap_co2-log(vmr_co2_PCM(i,t)*ps_PCM_2(i,t)/100.)))**4 &
@@ -53 +53 @@
             enddo
             if (ave < 1e-4) ave = 0.
-            tendencies_co2_ice_phys(i,islope) = tendencies_co2_ice_phys_ini(i,islope) - coef*ave/timelen
+            d_co2ice_phys(i,islope) = d_co2ice_phys_ini(i,islope) - coef*ave/timelen
         endif
     enddo