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

Timestamp:

Feb 7, 2024, 12:14:38 PM (11 months ago)

Author:

jbclement

Message:

Mars PCM:

• Addition of the possibility to use subslopes parametization in 1D. To do so, put 'nslope=x' in "run.def" where x (1, 3, 5 or 7) is the number of slopes you want to, or modify it in an appropriate "startfi.nc". Then, default subslopes definition and distribution are used by the model. The already existing case of using one slope whose inclination and orientation are defined in "run.def" is still possible.
• Some cleanings throughout the code, in particular related to unused variables.

JBC

File:

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

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

@@ -1761 +1761 @@
       select case(nslope)
           case(1)
+              default_def_slope(1) = -90.
+              default_def_slope(2) = 90.
+          case(3)
               default_def_slope(1) = -50.
-              default_def_slope(2) = 50.
+              default_def_slope(2) = -3.
+              default_def_slope(3) = 3.
+              default_def_slope(4) = 50.
           case(5)
               default_def_slope(1) = -43.
@@ -1781 +1786 @@
           case default
               write(*,*) 'Number of slopes not possible: nslope should 
-     & be 1, 5 or 7!'
+     & be 1, 3, 5 or 7!'
               call abort
       end select
@@ -1801 +1806 @@
        ENDDO
 
-       call subslope_mola(ngridmx,nslope_new,def_slope,subslope_dist)
+       if (ngridmx /= 1) then
+         call subslope_mola(ngridmx,nslope_new,def_slope,subslope_dist)
+       endif
 
 ! Surfdat related stuff