Changeset 3571 for trunk/LMDZ.COMMON/libf/evolution/recomp_tend_co2_mod.F90

Timestamp:

Jan 10, 2025, 5:45:03 PM (10 days ago)

Author:

jbclement

Message:

PEM:

• New way to manage the pressure: now the PEM manages only the average pressure and keeps the pressure deviation with the instantaneous pressure from the start to reconstruct the pressure at the end ('ps_avg = ps_start + ps_dev'). As a consequence, everything related to pressure in the PEM is modified accordingly.
• Surface temperatures management is now simpler. It follows the strategy for the pressure (and soil temperature) described above.
• Soil temperatures are now adapted to match the surface temperature changes occured during the PEM by modifying the soil temperature deviation at the end.
• Few simplifications/optimizations: notably, the two PCM years are now read in one go in 'read_data_PCM_mod.F90' and only the needed variables are extracted.
• Deletion of unused variables and unnecessary intermediate variables (memory saving and loop deletion in some cases).
• Renaming of variables and subroutines to make everything clearer. In particular, the suffixes: '_avg' = average, '_start' = PCM start file, '_dev' = deviation, '_ini' or '0' = initial, '_dyn' = dynamical grid, '_timeseries' = daily average of last PCM year.
• Cosmetic cleanings for readability.

JBC

File:

• trunk/LMDZ.COMMON/libf/evolution/recomp_tend_co2_mod.F90 (moved) (moved from trunk/LMDZ.COMMON/libf/evolution/recomp_tend_co2_slope_mod.F90) (5 diffs)

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

@@ -1 @@
-MODULE recomp_tend_co2_slope_mod
+MODULE recomp_tend_co2_mod
 
 implicit none
@@ -8 +8 @@
 
 SUBROUTINE recomp_tend_co2(ngrid,nslope,timelen,d_co2ice_phys,d_co2ice_ini,co2ice,emissivity, &
-                           vmr_co2_PCM,vmr_co2_PEM,ps_PCM_2,global_avg_press_PCM,global_avg_press_new)
+                           vmr_co2_PCM,vmr_co2_PEM,ps_PCM,ps_avg_global_ini,ps_avg_global)
 
 use constants_marspem_mod, only : alpha_clap_co2, beta_clap_co2, sigmaB, Lco2, sols_per_my, sec_per_sol
@@ -20 +20 @@
 !=======================================================================
 
-!   arguments:
-!   ----------
-!   INPUT
+! Inputs
+! ------
 integer,                        intent(in) :: timelen, ngrid, nslope
-real, dimension(ngrid,timelen), intent(in) :: vmr_co2_PCM          ! physical point field: Volume mixing ratio of co2 in the first layer
-real, dimension(ngrid,timelen), intent(in) :: vmr_co2_PEM          ! physical point field: Volume mixing ratio of co2 in the first layer
-real, dimension(ngrid,timelen), intent(in) :: ps_PCM_2             ! physical point field: Surface pressure in the PCM
-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) :: d_co2ice_ini         ! physical point field: Evolution of perennial ice over one year
-real, dimension(ngrid,nslope),  intent(in) :: co2ice               ! CO2 ice per mesh and sub-grid slope (kg/m^2)
-real, dimension(ngrid,nslope),  intent(in) :: emissivity           ! Emissivity per mesh and sub-grid slope(1)
-!   OUTPUT
+real, dimension(ngrid,timelen), intent(in) :: vmr_co2_PCM       ! physical point field: Volume mixing ratio of co2 in the first layer
+real, dimension(ngrid,timelen), intent(in) :: vmr_co2_PEM       ! physical point field: Volume mixing ratio of co2 in the first layer
+real, dimension(ngrid,timelen), intent(in) :: ps_PCM            ! physical point field: Surface pressure in the PCM
+real,                           intent(in) :: ps_avg_global_ini ! global averaged pressure at previous timestep
+real,                           intent(in) :: ps_avg_global     ! global averaged pressure at current timestep
+real, dimension(ngrid,nslope),  intent(in) :: d_co2ice_ini      ! physical point field: Evolution of perennial ice over one year
+real, dimension(ngrid,nslope),  intent(in) :: co2ice            ! CO2 ice per mesh and sub-grid slope (kg/m^2)
+real, dimension(ngrid,nslope),  intent(in) :: emissivity        ! Emissivity per mesh and sub-grid slope(1)
+! Outputs
+! -------
 real, dimension(ngrid,nslope), intent(inout) :: d_co2ice_phys ! physical point field: Evolution of perennial ice over one year
 
-!   local:
-!   ------
+! Local:
+! ------
 integer :: i, t, islope
-real    :: coef, ave
+real    :: coef, avg
 
 write(*,*) "Update of the CO2 tendency from the current pressure"
@@ -46 +46 @@
     do islope = 1,nslope
         coef = sols_per_my*sec_per_sol*emissivity(i,islope)*sigmaB/Lco2
-        ave = 0.
+        avg = 0.
         if (co2ice(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 &
-                      - (beta_clap_co2/(alpha_clap_co2-log(vmr_co2_PEM(i,t)*ps_PCM_2(i,t)*(global_avg_press_new/global_avg_press_PCM)/100.)))**4
+                avg = avg + (beta_clap_co2/(alpha_clap_co2-log(vmr_co2_PCM(i,t)*ps_PCM(i,t)/100.)))**4 &
+                      - (beta_clap_co2/(alpha_clap_co2-log(vmr_co2_PEM(i,t)*ps_PCM(i,t)*(ps_avg_global/ps_avg_global_ini)/100.)))**4
             enddo
-            if (ave < 1.e-4) ave = 0.
-            d_co2ice_phys(i,islope) = d_co2ice_ini(i,islope) - coef*ave/timelen
+            if (avg < 1.e-4) avg = 0.
+            d_co2ice_phys(i,islope) = d_co2ice_ini(i,islope) - coef*avg/timelen
         endif
     enddo
@@ -59 +59 @@
 
 END SUBROUTINE recomp_tend_co2
+!=======================================================================
 
-END MODULE recomp_tend_co2_slope_mod
+SUBROUTINE recomp_tend_h2o(ngrid,nslope,timelen,d_h2oice,PCM_temp,PEM_temp)
+
+implicit none
+
+!=======================================================================
+!
+!  Routine that compute the evolution of the tendencie for h2o ice
+!
+!=======================================================================
+
+! Inputs
+! ------
+integer,            intent(in) :: timelen, ngrid, nslope
+real, dimension(:), intent(in) :: PCM_temp, PEM_temp
+! Outputs
+! -------
+real, dimension(ngrid,nslope), intent(inout) :: d_h2oice ! physical point field: Evolution of perennial ice over one year
+
+! Local:
+! ------
+
+write(*,*) "Update of the H2O tendency due to lag layer"
+
+! Flux correction due to lag layer
+!~ Rz_old = h2oice_depth_old*0.0325/4.e-4              ! resistance from PCM
+!~ Rz_new = h2oice_depth_new*0.0325/4.e-4              ! new resistance based on new depth
+!~ R_dec = (1./Rz_old)/(1./Rz_new)                     ! decrease because of resistance
+!~ soil_psv_old = psv(max(PCM_temp(h2oice_depth_old))) ! the maxmimum annual mean saturation vapor pressure at the temperature of the GCM run temperature at the old ice location
+!~ soil_psv_new = psv(max(PEM_temp(h2oice_depth_new))) ! the maxmimum annual mean saturation vapor pressure at the temperature of the PEM run temperature at the new ice location
+!~ hum_dec = soil_psv_old/soil_psv_new                 ! decrease because of lower water vapor pressure at the new depth
+!~ d_h2oice = d_h2oice*R_dec*hum_dec                   ! decrease of flux
+
+END SUBROUTINE recomp_tend_h2o
+
+END MODULE recomp_tend_co2_mod