Changeset 709 for trunk/LMDZ.MARS/libf

Timestamp:

Jun 22, 2012, 3:45:44 PM (13 years ago)

Author:

aslmd

Message:

MESOSCALE: set default r_aspect_thermals to 1 also in mesoscale. updated DEF files for reference THARSIS case (latest settings by TN).

File:

• trunk/LMDZ.MARS/libf/phymars/calltherm_interface.F90 (modified) (2 diffs)

trunk/LMDZ.MARS/libf/phymars/calltherm_interface.F90

@@ -209 +209 @@
                                 ! very well with a sub-timestep of 10.
          IF (firstcall) THEN
+            r_aspect_thermals=1.  ! same value is OK for GCM and mesoscale
 #ifdef MESOSCALE
             !! valid for timesteps < 200s
             nsplit_thermals=4
-            r_aspect_thermals=0.7
 #else
             IF ((ptimestep .le. 3699.*24./96.) .and. (iradia .eq. 1)) THEN
@@ -219 +219 @@
                nsplit_thermals=35
             ENDIF
-            r_aspect_thermals=1.
 #endif
             call getin("nsplit_thermals",nsplit_thermals)