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

Timestamp:

Nov 8, 2023, 4:34:28 PM (14 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/libf/dynphy_lonlat/phymars

Files:

• newstart.F (modified) (3 diffs)
• start2archive.F (modified) (2 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)

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)