Changeset 3118 for trunk/LMDZ.MARS/libf/dynphy_lonlat/phymars/newstart.F

Timestamp:

Nov 8, 2023, 4:34:28 PM (16 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

File:

• trunk/LMDZ.MARS/libf/dynphy_lonlat/phymars/newstart.F (modified) (3 diffs)

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)