Changeset 2520 for trunk/LMDZ.GENERIC/libf/phystd/optci.F90

Timestamp:

May 12, 2021, 2:45:26 PM (4 years ago)

Author:

Guillaume Chaverot

Message:

update of the H2O continuum

File:

• trunk/LMDZ.GENERIC/libf/phystd/optci.F90 (modified) (2 diffs)

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

@@ -14 +14 @@
   use gases_h, only: gfrac, ngasmx, igas_N2, igas_He, igas_H2O, igas_H2
   use comcstfi_mod, only: g, r, mugaz
-  use callkeys_mod, only: kastprof,continuum,graybody,H2Ocont_simple
+  use callkeys_mod, only: kastprof,continuum,graybody
   implicit none
 
@@ -205 +205 @@
                  enddo
 
-              elseif(igas.eq.igas_H2O.and.T_cont.gt.200.0)then
-
+              elseif(igas.eq.igas_H2O.and.T_cont.gt.100.0)then
+                 ! Compute self and foreign (with air) continuum of H2O 
                  p_air = dble(PMID(k)*scalep) - p_cont ! note assumes background is air!
-                 if(H2Ocont_simple)then
-                    call interpolateH2Ocont_PPC(wn_cont,T_cont,p_cont,p_air,dtemp,.false.)
-                 else
-                    interm = indi(nw,igas,igas)
-                    call interpolateH2Ocont_CKD(wn_cont,T_cont,p_cont,p_air,dtemp,.false.,interm)
-                    indi(nw,igas,igas) = interm
-                 endif
+                 interm = indi(nw,igas,igas)
+                 call interpolateH2O_self_foreign(wn_cont,T_cont,p_cont,p_air,dtemp,.false.,interm) ! MTCKD v3.3
+                 indi(nw,igas,igas) = interm
 
               endif