Changeset 2522

Timestamp:

May 15, 2021, 9:36:14 PM (4 years ago)

Author:

jnaar

Message:

Added flag for temperature-dependant water contact parameter. Default is false for now !
JN

Location:

trunk/LMDZ.MARS

Files:

• README (modified) (1 diff)
• libf/phymars/callkeys.h (modified) (2 diffs)
• libf/phymars/conf_phys.F (modified) (1 diff)
• libf/phymars/nuclea.F (modified) (3 diffs)

trunk/LMDZ.MARS/README

@@ -3379 +3379 @@
 Correction of the threshold value used to prevent from negative values in watercloud_mod.F: the threshold is set to 1e-16 (=qperemin, an already existing parameter for the water cycle), instead of 1e-8. This former value was at the origin of a strong and unrealistic drying out of the upper atmosphere, when the photochemistry was active, during the first part of the martian year (months 3 and 4). (bug reported by Francisco Gonzalez-Galindo).
 
+== 15/05/2021 == JN
+Added a logical flag allowing for a temperature-dependant parametrization of the water contact parameter : "temp_dependant_m". Default value is false for now.
+
+

trunk/LMDZ.MARS/libf/phymars/callkeys.h

@@ -17 +17 @@
      &   ,CLFvarying,satindexco2,rdstorm,slpwind,calllott_nonoro        &
      &   ,latentheat_surfwater,gwd_convective_source,startphy_file      &
-     &   ,hdo,hdofrac,cap_albedo
+     &   ,hdo,hdofrac,cap_albedo,temp_dependant_m
      
       COMMON/callkeys_i/iradia,iaervar,ilwd,ilwb,ilwn,ncouche           &
@@ -70 +70 @@
       logical latentheat_surfwater ! latent heat release from ground water ice sublimation/condensation
       logical cap_albedo ! polar cap albedo remains unchanged by water frost deposition
+      logical temp_dependant_m ! temperature-dependant water contact parameter
       logical sedimentation
       logical activice,tifeedback,supersat,caps

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

@@ -706 +706 @@
 ! microphysical parameter contact       
         write(*,*) "water contact parameter ?"
-        mteta  = 0.95
-        call getin_p("mteta",mteta)
-        write(*,*) "mteta = ", mteta
+        mteta  = 0.95 ! default value
+        temp_dependant_m  = .false. ! default value
+        call getin_p("temp_dependant_m",temp_dependant_m)
+        if (temp_dependant_m) then
+           print*,'You have chosen a temperature-dependant water'
+           print*,'contact parameter ! From Maattanen et al. 2014'
+        else if (.not.temp_dependant_m) then
+           print*,'Water contact parameter is constant'
+           call getin_p("mteta",mteta)
+           write(*,*) "mteta = ", mteta
+        endif
         
 ! scavenging

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

@@ -16 +16 @@
 
 #include "microphys.h"
+      include "callkeys.h"
 
 c     Inputs
@@ -55 +56 @@
       mtetalocal = mteta  !! use mtetalocal for better performance
 
+      IF (temp_dependant_m) THEN
+         mtetalocal = min(0.0044*temp + 0.1831,0.97)
+      ENDIF ! (temp_dependant_m) THEN
 cccccccccccccccccccccccccccccccccccccccccccccccccc
 ccccccccccc ESSAIS TN MTETA = F (T) cccccccccccccc
@@ -70 +74 @@
 cccccccccccccccccccccccccccccccccccccccccccccccccc
 cccccccccccccccccccccccccccccccccccccccccccccccccc 
-      IF (firstcall) THEN
+      IF (firstcall.and.temp_dependant_m) THEN
+          print*, ' '  
+          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*, ' '  
+         firstcall=.false.
+      ELSE IF (firstcall.and.(.not.(temp_dependant_m))) THEN
           print*, ' '  
           print*, 'dear user, please keep in mind that'
           print*, 'contact parameter IS constant'
-          !print*, 'contact parameter IS NOT constant:'
-          !print*, 'max(mteta - 6005*exp(-0.065*temp),0.1)'
-          !print*, 'max(mteta - 6005*exp(-0.068*temp),0.1)'
           print*, ' '  
          firstcall=.false.