Changeset 2198

Timestamp:

Dec 13, 2019, 1:14:02 PM (5 years ago)

Author:

flefevre

Message:

reglages Diogo Quirino & Gabriella Gilli pour la haute atmosphere

Location:

trunk/LMDZ.VENUS/libf/phyvenus

Files:

• diffusion.h (modified) (1 diff)
• flott_gwd_ran.F90 (modified) (1 diff)
• nirco2abs.F (modified) (4 diffs)
• nlteparams.h (modified) (1 diff)

trunk/LMDZ.VENUS/libf/phyvenus/diffusion.h

@@ -11 +11 @@
       parameter (Pdiff=1.e-1)      ! pressure below which diffusion is computed
       parameter (tdiffmin=5.d0)
-      parameter (dzres=2.d0)     ! grid resolution (km) for diffusion
+      parameter (dzres=0.1d0)    ! grid resolution (km) for diffusion
       

trunk/LMDZ.VENUS/libf/phyvenus/flott_gwd_ran.F90

@@ -146 +146 @@
     endif
     CMIN = 1.           ! Min phase velocity
-    CMAX = 111.          ! Max phase speed velocity
+    CMAX = 61.          ! Max phase speed velocity
 !    XLAUNCH=0.6         ! Parameter that control launching altitude
     XLAUNCH=5e-3        ! Value for top of cloud convective region

trunk/LMDZ.VENUS/libf/phyvenus/nirco2abs.F

@@ -23 +23 @@
 c   Stephen Lewis 2000
 c
+c   apr 2019 d.quirino   Improving NLTE params, SOIR/SPICAV Temp comparison 
 c   oct 2014 g.gilli     Coupling with photochemical model
 C   jan 2014 g.gilli     Revision (following martian non-lte param)    
@@ -104 +105 @@
 
     
-cc "Nominal" values 
-       parameter (n_a = 18.13/86400.0)     !c     K/Eday  ---> K/sec   
-       parameter (p_ctop=13.2e2)
-c    -- NLTE Param v3  --
-       parameter (n_p0=0.008)  
+cc "Nominal" values used in Gilli+2'17
+c       parameter (n_a = 18.13/86400.0)     !c     K/Eday  ---> K/sec   
+c       parameter (p_ctop=13.2e2)
+c       parameter (n_p0=0.008) 
+
+cc "New" values used to improve SPICAV/SOIR Temperature comparision (D.Quirino)
+       parameter (n_a = 15.92/86400.0)     !c     K/Eday  ---> K/sec   
+       parameter (p_ctop=19.85e2)
+       parameter (n_p0=0.1)  
        parameter (n_b=1.362)
-
-cc   -- Varoxy5
-C       parameter (n_a = 20/86400.0)
-C       parameter (p_ctop=870)   ! [Pa] 
-C       parameter (n_b=1.98)
-C       parameter (n_p0=0.045) 
-
-c      parameter (n_p0=0.1) !!!!cccc test varoxy5mod
-c      parameter (n_b=0.9)
-
 
 c    -- NLTE Param v2  --
@@ -125 +120 @@
 c       parameter (n_b = 1.3)
    
-
 
 c     Variables added to implement NLTE correction factor (feb 2011)
@@ -141 +135 @@
 
 CCCCCC   TEST: reduce by X% nir Heating
-
 c      co2heat0  = co2heat0 * 0.8
-
 
 c----------------------------------------------------------------------

trunk/LMDZ.VENUS/libf/phyvenus/nlteparams.h

@@ -7 +7 @@
 
       real ptrans           ! central pressure for transition (Pa)
-      parameter (ptrans = 0.5)
+c     "Nominal" values in Gilli+2017
+c      parameter (ptrans = 0.5)
+c     "New" values by D.Quirino  (better agreement with SOIR/SPICAV temperature)
+      parameter (ptrans = 0.2)
       real zw               ! half-width for transition (scale heights)
       parameter (zw = 0.5)