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

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

@@ -11 +11 @@
   use radinc_h, only: L_NLAYRAD, L_NLEVRAD, L_LEVELS, L_NSPECTV, L_NGAUSS, L_REFVAR, NAERKIND
   use radcommon_h, only: gasv, tlimit, wrefVAR, Cmk, tgasref, pfgasref,wnov,scalep,indv,glat_ig
-  use gases_h, only: gfrac, ngasmx, igas_H2, igas_H2O, igas_He, igas_N2
+  use gases_h, only: gfrac, ngasmx, igas_H2, igas_H2O, igas_He, igas_N2, igas_CH4
   use comcstfi_mod, only: g, r, mugaz
   use callkeys_mod, only: kastprof,continuum,graybody,callgasvis
@@ -215 +215 @@
                     dtemp = dtemp + dtempc
                  enddo
+                 
+              elseif(igas.eq.igas_CH4)then
+
+                 ! first do self-induced absorption
+                 interm = indv(nw,igas,igas)
+                 call interpolateCH4CH4(wn_cont,T_cont,p_cont,dtemp,.false.,interm)
+                 indv(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 = indv(nw,igas,jgas)
+                       call interpolateH2CH4(wn_cont,T_cont,p_cross,p_cont,dtempc,.false.,interm)
+                       indv(nw,igas,jgas) = interm
+                    elseif(jgas.eq.igas_He)then 
+                       interm = indv(nw,igas,jgas)
+                       call interpolateHeCH4(wn_cont,T_cont,p_cross,p_cont,dtempc,.false.,interm)
+                       indv(nw,igas,jgas) = interm
+                    endif
+                    dtemp = dtemp + dtempc
+                 enddo
 
               elseif(igas.eq.igas_H2O.and.T_cont.gt.100.0)then