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

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

@@ -12 +12 @@
                       L_NLEVRAD, L_REFVAR, naerkind
   use radcommon_h, only: gasi,tlimit,wrefVAR,Cmk,tgasref,pfgasref,wnoi,scalep,indi,glat_ig
-  use gases_h, only: gfrac, ngasmx, igas_N2, igas_He, igas_H2O, igas_H2
+  use gases_h, only: gfrac, ngasmx, igas_N2, igas_He, igas_H2O, igas_H2, igas_CH4
   use comcstfi_mod, only: g, r, mugaz
   use callkeys_mod, only: kastprof,continuum,graybody
@@ -208 +208 @@
                  enddo
 
+              elseif(igas.eq.igas_CH4)then
+
+                 ! first do self-induced absorption
+                 interm = indi(nw,igas,igas)
+                 call interpolateCH4CH4(wn_cont,T_cont,p_cont,dtemp,.false.,interm)
+                 indi(nw,igas,igas) = interm
+
+                 ! then cross-interactions with other gases
+                 do jgas=1,ngasmx
+                    p_cross = dble(PMID(k)*scalep*gfrac(jgas)*(1.-QVAR(k)))
+                    dtempc  = 0.0d0
+                    if(jgas.eq.igas_H2)then 
+                       interm = indi(nw,igas,jgas)
+                       call interpolateH2CH4(wn_cont,T_cont,p_cross,p_cont,dtempc,.false.,interm)
+                       indi(nw,igas,jgas) = interm
+                    elseif(jgas.eq.igas_He)then 
+                       interm = indi(nw,igas,jgas)
+                       call interpolateHeCH4(wn_cont,T_cont,p_cross,p_cont,dtempc,.false.,interm)
+                       indi(nw,igas,jgas) = interm
+                    endif
+                    dtemp = dtemp + dtempc
+                 enddo
+
               elseif(igas.eq.igas_H2O.and.T_cont.gt.100.0)then
                  ! Compute self and foreign (with air) continuum of H2O