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

Timestamp:

Sep 19, 2018, 1:04:56 PM (6 years ago)

Author:

emillour

Message:

Mars GCM:

• Bug fix in newstart. lon/lat and cell_area were not correctly sent to the restartfi.nc file.

EM

File:

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

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

@@ -30 +30 @@
       use comsoil_h, only: inertiedat, layer, mlayer, nsoilmx, tsoil
       use control_mod, only: day_step, iphysiq, anneeref, planet_type
+      use geometry_mod, only: longitude,latitude,cell_area
       use phyetat0_mod, only: phyetat0
       use phyredem, only: physdem0, physdem1
@@ -108 +109 @@
       real ith(iip1,jjp1,nsoilmx),ithfi(ngridmx,nsoilmx) ! thermal inertia (3D)
       real surfith(iip1,jjp1),surfithfi(ngridmx) ! surface thermal inertia (2D)
-      REAL latfi(ngridmx),lonfi(ngridmx),airefi(ngridmx)
+!      REAL latfi(ngridmx),lonfi(ngridmx),airefi(ngridmx)
 
       INTEGER i,j,l,isoil,ig,idum
@@ -1646 +1647 @@
 C Ecriture etat initial physique
 C
-
-      call physdem0("restartfi.nc",lonfi,latfi,nsoilmx,ngridmx,llm,
-     .              nqtot,dtphys,real(day_ini),0.0,
-     .              airefi,albfi,ithfi,zmea,zstd,zsig,zgam,zthe,
-     .              hmons)
+      call physdem0("restartfi.nc",longitude,latitude,
+     &              nsoilmx,ngridmx,llm,
+     &              nqtot,dtphys,real(day_ini),0.0,cell_area,
+     &              albfi,ithfi,zmea,zstd,zsig,zgam,zthe,
+     &              hmons)
       call physdem1("restartfi.nc",nsoilmx,ngridmx,llm,nqtot,
      &              dtphys,hour_ini,