Changeset 2837

Timestamp:

Dec 1, 2022, 10:55:44 AM (2 years ago)

Author:

aslmd

Message:

Changed the H2H2 CIA file in the Hot Jupiter case

Modified the default H2H2 CIA file. The file used now
(H2-H2_2011_extended.cia) is a combinaison of multiple files existing
for hot temperature. It is very similar to the H2H2 CIA used in Baudino
et al, 2015, but has been extended to 7000 K instead of 5000 K.

LT

File:

• trunk/LMDZ.GENERIC/libf/phystd/interpolateH2H2.F90 (modified) (3 diffs)

trunk/LMDZ.GENERIC/libf/phystd/interpolateH2H2.F90

@@ -63 +63 @@
       integer ind
       if ((H2orthopara_mixture .eq. "hot")) then
-        ! print*,"We're in the Hot Jupiter case"
-        if (temp .gt. 3000.) then
+        if (temp .gt. 7000.) then
           if (strictboundcia) then
             if (is_master) then
-              print*,'Your temperatures are too high for this H2-H2 CIA dataset (Hot Jupiter case). If you '
-              print*,'really want to run simulations with hydrogen at T > 400 K, contact'
-              print*,'Robin Wordsworth [rwordsworth@uchicago.edu].'
+              print*,'Your temperatures are too high for this H2-H2 CIA dataset (>7000K, Hot Jupiter case)'
             endif !is_master
             stop
@@ -77 +74 @@
               print*,'you have chosen strictboundcia = ', strictboundcia
               print*,'*********************************************************'
-              print*,' we allow model to continue but with temp = 3000          '
+              print*,' we allow model to continue but with temp = 7000          '
               print*,'  ...       for H2-H2 CIA dataset          ...           '
               print*,'  ... we assume we know what you are doing ...           '
               print*,'*********************************************************'
             endif !is_master
-            temp = 3000.
+            temp = 7000.
           endif !strictboundcia
-        endif !(temp .gt. 3000.)
-      else ! if not Hot Jupiter
+        endif !(temp .gt. 7000.)
+      else ! if not "hot"
         if(temp.gt.400)then
           if (strictboundcia) then
@@ -143 +140 @@
              dt_file=TRIM(datadir)//'/continuum_data/H2-H2_eq_2011.cia'
            else if (H2orthopara_mixture.eq."hot") then 
-            dt_file=TRIM(datadir)//'/continuum_data/H2-H2_2011.cia'
-            ns = 9981
+            dt_file=TRIM(datadir)//'/continuum_data/H2-H2_2011_extended.cia'
+            ns = 800
            endif
          else if (versH2H2cia.eq.2018) then