Changeset 3002 for trunk/LMDZ.COMMON/libf/evolution/glaciers_mod.F90

Timestamp:

Jul 20, 2023, 11:01:17 AM (16 months ago)

Author:

llange

Message:

MARS PEM

• Fix a bug in conf_pem ("Flux_geo" was expected in the .def while it should be 'fluxgeo', corrected)
• Fix a bug in glaciers modules (the length for 'name_ice' was too large)
• Fix a bug when writing tsoil in the PEM (ngrid x nsoil_PEM x nslope was given, while it was expecting a ngrid x nsoil_GCM x nslope)

PEM runs correctly now
LL

File:

• trunk/LMDZ.COMMON/libf/evolution/glaciers_mod.F90 (modified) (6 diffs)

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

@@ -57 +57 @@
 
       call transfer_ice_duringflow(ngrid,nslope,iflat, subslope_dist,def_slope_mean,hmax,Tcond,"co2",co2ice,flag_co2flow,flag_co2flow_mesh)
-   RETURN
+   RETURN   
 end subroutine
 
@@ -96 +96 @@
 
       call compute_hmaxglaciers(ngrid,nslope,iflat,def_slope_mean,Tice,"h2o",hmax)
-
       call transfer_ice_duringflow(ngrid,nslope,iflat, subslope_dist,def_slope_mean,hmax,Tice,"h2o",h2oice,flag_h2oflow,flag_h2oflow_mesh)
 
@@ -127 +126 @@
       REAL,INTENT(IN) :: def_slope_mean(nslope) ! Slope field: Values of the subgrid slope angles [deg]
       REAL,INTENT(IN) :: Tice(ngrid,nslope)     ! Physical field:  ice temperature [K]
-      character(len=30), INTENT(IN) :: name_ice ! Nature of the ice
+      character(len=3), INTENT(IN) :: name_ice ! Nature of the ice
 ! Outputs
       REAL,INTENT(OUT) :: hmax(ngrid,nslope) ! Physical grid x Slope field: maximum  thickness before flaw [m]
@@ -135 +134 @@
       INTEGER :: ig,islope ! loop variables
       REAL :: slo_angle
-
 
 ! 1. Compute rho
@@ -198 +196 @@
       REAL, INTENT(IN) :: hmax(ngrid,nslope)            ! maximum height of the  glaciers before initiating flow [m]
       REAL, INTENT(IN) :: Tice(ngrid,nslope)            ! Ice temperature[K]
-      character(len=30), INTENT(IN) :: name_ice              ! Nature of the ice
+      character(len=3), INTENT(IN) :: name_ice              ! Nature of the ice
 
 ! Outputs
@@ -208 +206 @@
       REAL rho(ngrid,nslope) ! density of ice, temperature dependant [kg/m^3]
       INTEGER iaval ! ice will be transfered here
-
-
 
 ! 0. Compute rho