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

Timestamp:

Sep 28, 2023, 3:23:58 PM (22 months ago)

Author:

jnaar

Message:

Adaptative timestep and temperature-dependent contact parameter now required together to run. Actualized relation and coefficients corresponding to the latest tuning of the global water cycle. JN

Location:

trunk/LMDZ.MARS/libf/phymars

Files:

• conf_phys.F (modified) (2 diffs)
• improvedclouds_mod.F (modified) (1 diff)
• nuclea.F (modified) (2 diffs)

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

@@ -748 +748 @@
         temp_dependant_m  = .false. ! default value
         call getin_p("temp_dependant_m",temp_dependant_m)
-        if (temp_dependant_m) then
+        if (temp_dependant_m) then !(JN 2023)
            print*,'You have chosen a temperature-dependant water'
            print*,'contact parameter ! From Maattanen et al. 2014'
@@ -763 +763 @@
          call getin_p("cloud_adapt_ts",cloud_adapt_ts)
          write(*,*)"cloud_adapt_ts= ",cloud_adapt_ts
+
+! Test of incompatibility:
+! If you want the adaptative timestep, you should use the
+! temperature dependent contact parameter (otherwise the 
+! global water cycle is nuts !)
         
+        if ((cloud_adapt_ts.and..not.temp_dependant_m).or.
+     &       (temp_dependant_m.and..not.cloud_adapt_ts)) then
+           print*,'Water cycle v6 : cloud_adapt_ts and '
+           print*,'temp_dependant_m must be used together !'
+           print*,'Otherwise water cycle is off. At least one '
+           print*,'of these flags is turned off. '
+           call abort_physic(modname,
+     &          "water cycle v6 needs cloud_adapt_ts and temp_dep_m",1)
+        endif
 ! scavenging
         write(*,*)"Dust scavenging by H2O/CO2 snowfall ?"

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

@@ -740 +740 @@
       integer,intent(out) :: zimicro ! number of ptimestep division
 
-c       zimicro = 30
-c      Coefficients good enough for present-day Mars :
-c       alpha = 1.87485684e+09
-c       beta = 1.45655856e+00
-c      Coefficients covering high obliquity scenarios :
-       alpha=3.36711332e+15
-       beta=1.98636414e+00
-c       nimicro=min(max(alpha*abs(potcond)**beta,5.),7000.)
-c       zimicro=ceiling((ptimestep/timeleft)*nimicro)
-c       zimicro=ceiling((timeleft/ptimestep)*nimicro)
+c      Conservative coefficients :
+       alpha=1.81846504e+11
+       beta=1.54550140e+00
        zimicro=ceiling(min(max(alpha*abs(potcond)**beta,5.),7000.))
 

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

@@ -60 +60 @@
 
       IF (temp_dependant_m) THEN
+c     J.Naar - sep 2023 :
 c     Simple linear parametrisation from Maattaanen 2014
-c     Smectite sample
 c     Maxed out at 0.97 for physical realism
-         mtetalocal = min(0.0044*temp + 0.1831,0.97)
+c     Used to tune the WC without adaptative ts of microphy
+c         mtetalocal = min(0.0044*temp + 0.1831,0.97)
+c     With adapt_ts need to use this relation (same paper, tanh):
+         mtetalocal = 0.469+((0.972-0.469)*tanh((temp/158.282)**4.244))
       ENDIF ! (temp_dependant_m) THEN
 cccccccccccccccccccccccccccccccccccccccccccccccccc
@@ -84 +87 @@
           print*, 'dear user, please keep in mind that'
           print*, 'contact parameter IS NOT constant ;'
-          print*, 'Using the following linear fit from'
-          print*, 'Maattanen et al. 2014 (SM linear fit) :'
-          print*, 'min(0.0044*temp + 0.1831,0.97)'
+          print*, 'Using the tanh fit n2 (jsc sample) from'
+          print*, 'Maattanen et al. 2014 (10.1016/j.grj.2014.09.002) :'
+          print*, 'm=0.469+((0.972-0.469)*tanh((temp/158.282)**4.244))'
           print*, ' '  
          firstcall=.false.