Changeset 2655 for trunk/LMDZ.GENERIC/libf/phystd/inifis_mod.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/inifis_mod.F90 (modified) (2 diffs)

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

@@ -259 +259 @@
 
      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
+       ": prohibit calculations outside CIA T grid?"
+     strictboundcia=.true. ! default value
+     call getin_p("strictboundcia",strictboundcia)
+     if (is_master) write(*,*) trim(rname)//&
+       ": strictboundcia = ",strictboundcia
 
      if (is_master) write(*,*) trim(rname)//&
@@ -304 +304 @@
      if (versH2H2cia.ne.2011 .and. versH2H2cia.ne.2018) then
         call abort_physic(rname,"Error: Choose a correct value (2011 or 2018) for versH2H2cia !",1)
-     endif 
+     endif
+     
+     if (is_master) write(*,*) trim(rname)//& 
+       ": CIA - normal or equilibrium ortho-para mixture for H2?"
+     H2orthopara_mixture = 'normal'
+     call getin_p("H2orthopara_mixture",H2orthopara_mixture)
+     if (is_master) write(*,*)trim(rname)//&
+       ": H2orthopara_mixture = ",trim(H2orthopara_mixture)
 
      if (is_master) write(*,*) trim(rname)//&