Changeset 2655 for trunk/LMDZ.GENERIC/libf/phystd/interpolateH2H2.F90

Timestamp:

Mar 30, 2022, 12:19:20 PM (3 years ago)

Author:

gmilcareck

Message:

Major changes to CIA interpolation:
1) Add contribution from CH4 (H2-CH4,He-CH4,CH4-CH4) ;
2) Add some tests before interpolation for H2, He and CH4 ;
3) Add the possibility to choose between a normal or equilibrium ortho:para
fraction for CIA H2;
4) Change "strictboundH2H2cia" to the generic "strictboundcia" for H2,He,CH4.
It can be added for others CIA (N2,H,CO2...) if you want.

File:

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

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

@@ +1 @@
+
      subroutine interpolateH2H2(wn,temp,pres,abcoef,firstcall,ind)
 
@@ -17 +18 @@
 !==================================================================
 
-      use callkeys_mod, only: versH2H2cia, strictboundH2H2cia
+      use callkeys_mod, only: versH2H2cia, strictboundcia, H2orthopara_mixture
       use datafile_mod, only: datadir
 
@@ -60 +61 @@
 
       integer ind
-
+      
       if(temp.gt.400)then
-         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
+        if (strictboundcia) 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 strictboundcia = ', strictboundcia
+          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
+      elseif(temp.lt.40)then     
+        if (strictboundcia) then
+          print*,'Your temperatures are too low for this H2-H2 CIA dataset. If you '
+          print*,'really want to run simulations with hydrogen at T < 40 K, contact'
+          print*,'Robin Wordsworth [rwordsworth@uchicago.edu].'
+          stop
+        else
+          print*,'Your temperatures are too low for this H2-H2 CIA dataset'
+          print*,'you have chosen strictboundcia = ', strictboundcia
+          print*,'*********************************************************'
+          print*,' we allow model to continue but with temp = 40           '
+          print*,'  ...       for H2-H2 CIA dataset          ...           '
+          print*,'  ... we assume we know what you are doing ...           '
+          print*,'*********************************************************'
+          temp = 40
+        endif            
       endif
 
@@ -88 +105 @@
          ! Only two possible versions for now : 2011 or 2018 (sanity check in inifis_mod)
          if (versH2H2cia.eq.2011) then
-           dt_file=TRIM(datadir)//'/continuum_data/H2-H2_norm_2011.cia'
            nS = 2428
+           if (H2orthopara_mixture.eq."normal") then
+             dt_file=TRIM(datadir)//'/continuum_data/H2-H2_norm_2011.cia'
+           else if (H2orthopara_mixture.eq."equilibrium") then
+             dt_file=TRIM(datadir)//'/continuum_data/H2-H2_eq_2011.cia'
+           endif
          else if (versH2H2cia.eq.2018) then
-           dt_file=TRIM(datadir)//'/continuum_data/H2-H2_norm_2018.cia'
            nS = 9600
+           if (H2orthopara_mixture.eq."normal") then
+             dt_file=TRIM(datadir)//'/continuum_data/H2-H2_norm_2018.cia'
+           else if (H2orthopara_mixture.eq."equilibrium") then
+             dt_file=TRIM(datadir)//'/continuum_data/H2-H2_eq_2018.cia'
+           endif
          endif
 
@@ -106 +131 @@
            write(*,*) 'is correct. You can change it in callphys.def with:'
            write(*,*) 'datadir = /absolute/path/to/datagcm'
-           write(*,*) 'Also check that the continuum data continuum_data/H2-H2_norm_2011.cia or H2-H2_norm_2018.cia is there.'
+           write(*,*) 'Also check that the continuum data is there.'
            call abort
          else
@@ -113 +138 @@
            write(*,*) '... You are using H2-H2 CIA from 2011 but you should use more recent data available on HITRAN !'
            write(*,*) '... (Especially if you are running a giant planet atmosphere)'
-           write(*,*) '... Just find out the H2-H2_norm_2018.cia, put it in your datadir and have a look at interpolateH2H2.F90 ! .'
+           write(*,*) '... Just find out the H2-H2 CIA from 2018, put it in your datadir and have a look at interpolateH2H2.F90 ! .'
          endif
 
@@ -149 +174 @@
 
       endif
-
+      
+            
+            
          call bilinearbig(nS,nT,wn_arr,temp_arr,abs_arr,wn,temp,abcoef,ind)