Changeset 3114

Timestamp:

Nov 3, 2023, 5:47:39 PM (13 months ago)

Author:

jbclement

Message:

PEM:
Following r3113, addition of 'nqsoil' and 'qsoil' in the arguments of the subroutines 'phyetat0' and 'physdem1' to be able to compile.
JBC

Location:

trunk/LMDZ.COMMON/libf/evolution

Files:

• changelog.txt (modified) (1 diff)
• pem.F90 (modified) (6 diffs)

trunk/LMDZ.COMMON/libf/evolution/changelog.txt

@@ -127 +127 @@
 == 27/10/2023 == JBC
 Addition of the possibility to ask for a number of Mars years to be simulated in total (previously it was only in Earth years).
+
+== 03/11/2023 == JBC
+Following r3113, addition of 'nqsoil' and 'qsoil' in the arguments of the subroutines 'phyetat0' and 'physdem1' to be able to compile.

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

@@ -73 +73 @@
 
 #ifndef CPP_STD
-    use comsoil_h,          only: tsoil, nsoilmx, ini_comsoil_h, inertiedat, mlayer, volcapa, inertiesoil
+    use comsoil_h,          only: tsoil, nsoilmx, ini_comsoil_h, inertiedat, mlayer, volcapa, inertiesoil, nqsoil, qsoil
     use surfdat_h,          only: tsurf, emis, qsurf, watercap, ini_surfdat_h, &
                                   albedodat, zmea, zstd, zsig, zgam, zthe,     &
@@ -367 +367 @@
 ! First we read the initial state (starfi.nc)
 #ifndef CPP_STD
-    call phyetat0(FILE_NAME,0,0,nsoilmx,ngrid,nlayer,nq,day_ini,time_phys,tsurf, &
-                  tsoil,albedo,emis,q2,qsurf,tauscaling,totcloudfrac,wstar,      &
+    call phyetat0(FILE_NAME,0,0,nsoilmx,ngrid,nlayer,nq,nqsoil,day_ini,time_phys,tsurf, &
+                  tsoil,albedo,emis,q2,qsurf,qsoil,tauscaling,totcloudfrac,wstar,       &
                   watercap,perenial_co2ice,def_slope,def_slope_mean,subslope_dist)
 
@@ -383 +383 @@
     allocate(qsurf_read_generic(ngrid,nq))
     allocate(tsoil_read_generic(ngrid,nsoilmx))
+    allocate(qsoil_read_generic(ngrid,nsoilmx,nqsoil,nslope))
     allocate(emis_read_generic(ngrid))
     allocate(tsurf(ngrid,1))
@@ -393 +394 @@
     allocate(albedo(ngrid,2,1))
     allocate(inertiesoil(ngrid,nsoilmx,1))
-    call phyetat0(.true.,ngrid,nlayer,FILE_NAME,0,0,nsoilmx,nq,day_ini,time_phys, &
-                  tsurf_read_generic,tsoil_read_generic,emis_read_generic,q2,     &
-                  qsurf_read_generic,cloudfrac,totcloudfrac,hice,rnat,pctsrf_sic, &
-                  tslab, tsea_ice,sea_ice)
+    call phyetat0(.true.,ngrid,nlayer,FILE_NAME,0,0,nsoilmx,nq,nqsoil,day_ini,time_phys, &
+                  tsurf_read_generic,tsoil_read_generic,emis_read_generic,q2,            &
+                  qsurf_read_generic,qsoil_read_generic,cloudfrac,totcloudfrac,hice,     &
+                  rnat,pctsrf_sic,tslab,tsea_ice,sea_ice)
     call surfini(ngrid,nq,qsurf_read_generic,albedo_read_generic,albedo_bareground,albedo_snow_SPECTV,albedo_co2_ice_SPECTV)
 
@@ -1112 +1113 @@
                   inertiedat,def_slope,subslope_dist)
 
-    call physdem1("restartfi_evol.nc",nsoilmx,ngrid,nlayer,nq,  &
-                  ptimestep,ztime_fin,tsurf,tsoil,inertiesoil,  &
-                  albedo,emis,q2,qsurf,tauscaling,totcloudfrac, &
+    call physdem1("restartfi_evol.nc",nsoilmx,ngrid,nlayer,nq,nqsoil, &
+                  ptimestep,ztime_fin,tsurf,tsoil,inertiesoil,        &
+                  albedo,emis,q2,qsurf,qsoil,tauscaling,totcloudfrac, &
                   wstar,watercap,perenial_co2ice)
 #else
@@ -1121 +1122 @@
                   albedo_bareground,inertiedat,zmea,zstd,zsig,zgam,zthe)
 
-    call physdem1("restartfi_evol.nc",nsoilmx,ngrid,nlayer,nq,   &
-                  ptimestep,ztime_fin,tsurf,tsoil,emis,q2,qsurf, &
-                  cloudfrac,totcloudfrac,hice,rnat,pctsrf_sic,   &
-                  tslab,tsea_ice,sea_ice)
+    call physdem1("restartfi_evol.nc",nsoilmx,ngrid,nlayer,nq,nqsoil,  &
+                  ptimestep,ztime_fin,tsurf,tsoil,emis,q2,qsurf,qsoil, &
+                  cloudfrac,totcloudfrac,hice,rnat,pctsrf_sic,tslab,   &
+                  tsea_ice,sea_ice)
 #endif
 write(*,*) "restartfi_evol.nc has been written"