Changeset 2224 for trunk/LMDZ.MARS/libf/phymars

Timestamp:

Jan 24, 2020, 3:15:03 PM (5 years ago)

Author:

abierjon

Message:

Mars GCM:
Bug fix in physiq_mod for 1D runs, where the routine getslopes reinitialized the slope inclination and orientation to 0, while they were assigned in the input files.
BR+JN+AB

File:

• trunk/LMDZ.MARS/libf/phymars/physiq_mod.F (modified) (1 diff)

trunk/LMDZ.MARS/libf/phymars/physiq_mod.F

@@ -610 +610 @@
 
 #ifndef MESOSCALE
-         if (callslope) call getslopes(ngrid,phisfi)
-
-         if (ngrid.ne.1) then ! no need to create a restart file in 1d
-         call physdem0("restartfi.nc",longitude,latitude,
-     &                 nsoilmx,ngrid,nlayer,nq,
-     &                 ptimestep,pday,time_phys,cell_area,
-     &                 albedodat,inertiedat,zmea,zstd,zsig,zgam,zthe,
-     &                 hmons,summit,base)
+
+         if (ngrid.ne.1) then
+           ! no need to compute slopes when in 1D; it is an input
+           if (callslope) call getslopes(ngrid,phisfi)
+           ! no need to create a restart file in 1d
+           call physdem0("restartfi.nc",longitude,latitude,
+     &                   nsoilmx,ngrid,nlayer,nq,
+     &                   ptimestep,pday,time_phys,cell_area,
+     &                   albedodat,inertiedat,zmea,zstd,zsig,zgam,zthe,
+     &                   hmons,summit,base)
          endif