Changeset 2870 for trunk/LMDZ.GENERIC/libf/phystd/interpolateN2H2.F90

Timestamp:

Jan 16, 2023, 4:47:10 PM (23 months ago)

Author:

jleconte

Message:

Removed too many prints in CIA interpolation + correction for hot H2H2

File:

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

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

@@ -14 +14 @@
   !==================================================================
 
+  use callkeys_mod, only: strictboundcia
   use datafile_mod, only: datadir
   use mod_phys_lmdz_para, only : is_master
@@ -54 +55 @@
 
   character*20 bleh
-  double precision blah, Ttemp
+  double precision blah, Ttemp, ztemp
   integer nres
   integer ind
 
-  if(temp.gt.400)then
-     if (is_master) then
-       print*,'Your temperatures are too high for this N2-H2 CIA dataset.'
-       print*,'Please run mixed N2-H2 atmospheres below T = 400 K.'      
-     endif
-     stop
+  ztemp=temp
+
+  if(ztemp.gt.400)then
+    if (strictboundcia) then
+      if (is_master) then
+        print*,'Your temperatures are too high for this N2-H2 CIA dataset. If you '
+        print*,'really want to run simulations with hydrogen at T > 400 K,'
+        print*,'find relevant data.'
+      endif !is_master
+      stop
+    else
+      !if (is_master) then
+      !  print*,'Your temperatures are too high for this N2-H2 CIA dataset'
+      !  print*,'you have chosen strictboundcia = ', strictboundcia
+      !  print*,'*********************************************************'
+      !  print*,' we allow model to continue but with temp = 4000          '
+      !  print*,'  ...       for N2-H2 CIA dataset          ...           '
+      !  print*,'  ... we assume we know what you are doing ...           '
+      !  print*,'*********************************************************'
+      !endif !is_master
+      ztemp = 400.
+    endif !strictboundcia
   endif
 
@@ -121 +138 @@
   endif
 
-  call bilinearbig(nS,nT,wn_arr,temp_arr,abs_arr,wn,temp,abcoef,ind)
+  call bilinearbig(nS,nT,wn_arr,temp_arr,abs_arr,wn,ztemp,abcoef,ind)
 
 !     print*,'the absorption is ',abcoef,' cm^5 molecule^-2'