Changeset 2835 for trunk/LMDZ.COMMON/libf/evolution/compute_tendencies_mod_slope.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/compute_tendencies_mod_slope.F90 (modified) (3 diffs)

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

@@ -6 +6 @@
 
       IMPLICIT NONE
-
 
 !=======================================================================
@@ -34 +33 @@
 !=======================================================================
 
-
 !  We compute the difference
-!  tendencies_h2o_ice(:,:,:)=min_h2o_ice_Y2(:,:,:)-min_h2o_ice_Y1(:,:,:)
 
   DO j=1,jjm_input+1
@@ -46 +43 @@
   ENDDO
 
-     print *, "jjm_input+1", jjm_input+1
-     print *, "iim_input+1", iim_input+1
-     print *, "nslope+1", nslope+1
-
 !  If the difference is too small; there is no evolution
   DO j=1,jjm_input+1
     DO i = 1, iim_input
        DO islope = 1, nslope
-!         print *, "tendencies_h2o_ice(i,j,islope)LAAA", tendencies_h2o_ice(i,j,islope)
          if(abs(tendencies_h2o_ice(i,j,islope)).LT.1.0E-10) then
             tendencies_h2o_ice(i,j,islope)=0.
          endif
-!         print *, "tendencies_h2o_ice(i,j,islope)HERE", tendencies_h2o_ice(i,j,islope)
        enddo
     ENDDO
   ENDDO
 
-
-!  We reorganise the difference on the physical grid
-  tendencies_h2o_ice_phys(1,:)=tendencies_h2o_ice(1,1,:)
-
-  ig0 = 2
-  DO j=2,jjm_input
-    DO i = 1, iim_input
-       tendencies_h2o_ice_phys(ig0,:)  =tendencies_h2o_ice(i,j,:)
-       ig0= ig0 + 1
-    ENDDO
+  DO islope = 1,nslope
+    CALL gr_dyn_fi(1,iim_input+1,jjm_input+1,ngrid,tendencies_h2o_ice(:,:,islope),tendencies_h2o_ice_phys(:,islope))
   ENDDO
-
-  tendencies_h2o_ice_phys(ig0,:) = tendencies_h2o_ice(1,jjm_input+1,:)
-
-!  print *, "tendencies_h2o_ice_physze", tendencies_h2o_ice_phys(:,:)
-
 
 END SUBROUTINE compute_tendencies_slope
 
-
-
-
-