Changeset 2633

Timestamp:

Mar 1, 2022, 11:56:54 PM (3 years ago)

Author:

dbardet

Message:

Adding of strictboundH2H2cia flag to allow the run to continue even for temperature bigger than 400K (that is out of the current H2H2 cia tables)

Location:

trunk/LMDZ.GENERIC/libf/phystd

Files:

• callkeys_mod.F90 (modified) (1 diff)
• inifis_mod.F90 (modified) (1 diff)
• interpolateH2H2.F90 (modified) (2 diffs)

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

@@ -14 +14 @@
       logical,save :: strictboundcorrk                                     
 !$OMP THREADPRIVATE(strictboundcorrk)
+      logical,save :: strictboundH2H2cia                                     
+!$OMP THREADPRIVATE(strictboundH2H2cia)
 
       logical,save :: enertest

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

@@ -266 +266 @@
 
      if (is_master) write(*,*) trim(rname)//&
+       ": prohibit calculations outside H2H2 CIA T grid?"
+     strictboundH2H2cia=.true. ! default value
+     call getin_p("strictboundH2H2cia",strictboundH2H2cia)
+     if (is_master) write(*,*) trim(rname)//&
+       ": strictboundH2H2cia = ",strictboundH2H2cia
+
+     if (is_master) write(*,*) trim(rname)//&
        ": Minimum atmospheric temperature for Planck function integration ?"
      tplanckmin=30.0 ! default value

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

@@ -17 +17 @@
 !==================================================================
 
-      use callkeys_mod, only: versH2H2cia
+      use callkeys_mod, only: versH2H2cia, strictboundH2H2cia
       use datafile_mod, only: datadir
 
@@ -62 +62 @@
 
       if(temp.gt.400)then
-         print*,'Your temperatures are too high for this H2-H2 CIA dataset. If you '
-         print*,'really want to run simulations with hydrogen at T > 400 K, contact'
-         print*,'Robin Wordsworth [rwordsworth@uchicago.edu].'
-         stop
+         if (strictboundH2H2cia) then
+            print*,'Your temperatures are too high for this H2-H2 CIA dataset. If you '
+            print*,'really want to run simulations with hydrogen at T > 400 K, contact'
+            print*,'Robin Wordsworth [rwordsworth@uchicago.edu].'
+            stop
+         else
+            print*,'Your temperatures are too high for this H2-H2 CIA dataset'
+            print*,'you have chosen strictboundH2H2cia = ', strictboundH2H2cia
+            print*,'*********************************************************'
+            print*,' we allow model to continue but with temp = 400          '
+            print*,'  ...       for H2-H2 CIA dataset          ...           '
+            print*,'  ... we assume we know what you are doing ...           '
+            print*,'*********************************************************'
+            temp = 400
+         endif
       endif