Changeset 2835 for trunk/LMDZ.COMMON/libf/evolution/update_soil.F90

Timestamp:

Nov 30, 2022, 11:29:29 AM (2 years ago)

Author:

romain.vande

Message:

Mars PEM:
Introduction of the possibility to follow an orbital forcing.
Introduction of new control parameters.
Cleaning of the PEM (removing unused files, add comments and new files)

A file named run_PEM.def can be added to the run.def. It contains the following variables:

_ evol_orbit_pem: Boolean. Do you want to follow an orbital forcing predefined (read in ob_ex_lsp.asc for example)? (default=false)
_ year_bp_ini: Integer. Number of year before present to start the pem run if evol_orbit_pem=.true. , default=0
_ Max_iter_pem: Integer. Maximal number of iteration if none of the stopping criterion is reached and if evol_orbit_pem=.false., default=99999999
_ dt_pem: Integer. Time step of the PEM in year, default=1
_ alpha_criterion: Real. Acceptance rate of sublimating ice surface change, default=0.2
_ soil_pem: Boolean. Do you want to run with subsurface physical processes in the PEM? default=.true.

RV

File:

• trunk/LMDZ.COMMON/libf/evolution/update_soil.F90 (modified) (8 diffs)

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

@@ -1 @@
-   SUBROUTINE update_soil(ngrid,nslope,nsoil_PEM,tend_h2oglaciers,tend_co2glaciers,co2ice,waterice,ps_new,&
-                          cellarea,ice_depth,TI_PEM)
-
+   SUBROUTINE update_soil(ngrid,nslope,nsoil_PEM,tend_h2oglaciers,tend_co2glaciers,co2ice,waterice,p_avg_new,&
+                          ice_depth,TI_PEM)
 
  USE comsoil_h, only:  inertiedat, volcapa
@@ -10 +9 @@
  INTEGER,INTENT(IN) ::  ngrid, nslope, nsoil_PEM
  REAL,INTENT(IN) :: tend_h2oglaciers(ngrid,nslope),tend_co2glaciers(ngrid,nslope)
- REAL,INTENT(IN) :: ps_new(ngrid)
- REAL,INTENT(IN) :: cellarea(ngrid)
+ REAL,INTENT(IN) :: p_avg_new
  REAL,INTENT(IN) :: co2ice(ngrid,nslope)
  REAL,INTENT(IN) :: waterice(ngrid,nslope)
  REAL,INTENT(in) :: ice_depth(ngrid,nslope)
- 
 
 ! Outputs:
@@ -44 +41 @@
  REAL :: regolith_inertia(ngrid,nslope) ! TI of the regolith 
  REAL :: d(ngrid,nsoil_PEM,nslope)
- REAL :: p_avg_new
  REAL :: Total_surface
  INTEGER :: ispermanent_co2glaciers(ngrid,nslope)
  INTEGER :: ispermanent_h2oglaciers(ngrid,nslope)
 
-
 ! 0. Initialisation
-
-
-
-
-  Total_surface = sum(cellarea)
-  p_avg_new = 0.
-  do ig = 1,ngrid
-    p_avg_new = p_avg_new + ps_new(ig)*cellarea(ig)/Total_surface
-  enddo
-
-
 
   do ig = 1,ngrid
@@ -80 +64 @@
   enddo
 
-
  ispermanent_co2glaciers(:,:) = 0
  ispermanent_h2oglaciers(:,:) = 0
 
-
 !  1.Ice TI feedback
 
-!  do ig=1,ngrid
-!    do islope=1,nslope
-!      if ((ispermanent_co2glaciers(ig,islope).eq.1).or.(ispermanent_h2oglaciers(ig,islope).eq.1)) then
-!       do iloop = 1,n_1km
-!         TI_PEM(ig,iloop,islope) = surfaceice_inertia
-!       enddo
-!      endif
-!    enddo
-!   enddo
     do islope = 1,nslope
      call soil_TIfeedback_PEM(ngrid,nsoil_PEM,waterice(:,islope),   TI_PEM(:,:,islope))
@@ -101 +74 @@
 
 ! 2. Modification of the regolith thermal inertia.
-
-
 
   do ig=1,ngrid
@@ -115 +86 @@
    enddo
 enddo
-
-
-
-
-
 
 !  3. Build new TI for the PEM
@@ -142 +108 @@
          exit
         endif
-       
       enddo
-
-
 
       ! 4.2 Build the new ti
@@ -173 +136 @@
     enddo !ig
 
-
-
-
-
-
-
 !=======================================================================
       RETURN