Changeset 3118

Timestamp:

Nov 8, 2023, 4:34:28 PM (13 months ago)

Author:

llange

Message:

Mars PCM
Following -r 3113, fixing newstart & startarchive by implementing the changes in phyredem1 and phyetat0, i.e., adding qsoil and nqsoil in the arguments of the subroutine
I've tested newstart with existing files and it is okay. start2archive does not seems to work because of
an issue with perenial_ice. I'll look at it later.
LL

Location:

trunk/LMDZ.MARS

Files:

• changelog.txt (modified) (1 diff)
• libf/dynphy_lonlat/phymars/newstart.F (modified) (3 diffs)
• libf/dynphy_lonlat/phymars/start2archive.F (modified) (2 diffs)

trunk/LMDZ.MARS/changelog.txt

@@ -4317 +4317 @@
 == 08/11/2023 == JBC
 Correction of the case when startfiles1D = .true. and there is no "startfi.nc": the flags used in "init_testphys1d_mod.F90" were not correct so the model crashed and didn't create a "startfi.nc" from scratch as intended.
+
+== 08/11/2023 == LL
+Following -r 3113, fixing newstart & startarchive by implementing the changes
+in phyredem1 and phyetat0, i.e., adding qsoil and nqsoil in the arguments of
+the subroutine

trunk/LMDZ.MARS/libf/dynphy_lonlat/phymars/newstart.F

@@ -34 +34 @@
       use comsoil_h, only: inertiedat, inertiesoil,layer, mlayer,
      & nsoilmx,tsoil,ini_comsoil_h_slope_var, end_comsoil_h_slope_var,
-     & flux_geo
+     & flux_geo,qsoil,nqsoil
       use control_mod, only: day_step, iphysiq, anneeref, planet_type
       use geometry_mod, only: longitude,latitude,cell_area
@@ -465 +465 @@
         fichnom = 'startfi.nc'
         CALL phyetat0 (fichnom,tab0,Lmodif,nsoilmx,ngridmx,llm,nqtot,
-     &        day_ini,time,tsurf,tsoil,albedo,emis,
-     &        q2,qsurf,tauscaling,totcloudfrac,
+     &        nqsoil,day_ini,time,tsurf,tsoil,albedo,emis,
+     &        q2,qsurf,qsoil,tauscaling,totcloudfrac,
      &        wstar,watercap,perenial_co2ice,
      &        def_slope,def_slope_mean,subslope_dist)
@@ -2007 +2007 @@
      &              albfi,ithfi,def_slope,
      &              subslope_dist)
-      call physdem1("restartfi.nc",nsoilmx,ngridmx,llm,nqtot,
+      call physdem1("restartfi.nc",nsoilmx,ngridmx,llm,nqtot,nqsoil,
      &              dtphys,hour_ini,
-     &              tsurf,tsoil,inertiesoil,albedo,emis,q2,qsurf,
+     &              tsurf,tsoil,inertiesoil,albedo,
+     &              emis,q2,qsurf,qsoil,
      &              tauscaling,totcloudfrac,wstar,watercap,
      &              perenial_co2ice)

trunk/LMDZ.MARS/libf/dynphy_lonlat/phymars/start2archive.F

@@ -20 +20 @@
 
       use infotrac, only: infotrac_init, nqtot, tname
-      use comsoil_h, only: nsoilmx, inertiedat, inertiesoil
+      use comsoil_h, only: nsoilmx, inertiedat, inertiesoil,
+     &                     nqsoil, qsoil
       use surfdat_h, only: ini_surfdat_h, qsurf,watercaptag
       use comsoil_h, only: ini_comsoil_h
@@ -186 +187 @@
       allocate(subslope_distS(ip1jmp1,nslope))
 
-      CALL phyetat0 (fichnom,0,Lmodif,nsoilmx,ngridmx,llm,nqtot,
-     &      day_ini_fi,timefi,tsurf,tsoil,albedo,emis,q2,qsurf,
+      CALL phyetat0(fichnom,0,Lmodif,nsoilmx,ngridmx,llm,nqtot,nqsoil,
+     &      day_ini_fi,timefi,tsurf,tsoil,albedo,emis,q2,qsurf,qsoil,
      &      tauscaling,totcloudfrac,wstar,watercap,perenial_co2ice,
      &      def_slope, def_slope_mean,subslope_dist)