Changeset 2860 for trunk/LMDZ.GENERIC/libf

Timestamp:

Jan 4, 2023, 12:42:33 PM (2 years ago)

Author:

mturbet

Message:

add the CO2-H2 CIA opacity from Turbet2020 in the model

Location:

trunk/LMDZ.GENERIC/libf/phystd

Files:

• interpolateCO2H2.F90 (added)
• optci.F90 (modified) (1 diff)
• optcv.F90 (modified) (1 diff)
• sugas_corrk.F90 (modified) (1 diff)

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

@@ -200 +200 @@
                        call interpolateN2H2(wn_cont,T_cont,p_cross,p_cont,dtempc,.false.,interm)
                        indi(nw,igas,jgas) = interm
+                    elseif(jgas.eq.igas_CO2)then 
+                       interm = indi(nw,igas,jgas)
+                       call interpolateCO2H2(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)

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

@@ -208 +208 @@
                        indv(nw,igas,jgas) = interm
                        ! should be irrelevant in the visible
+                    elseif(jgas.eq.igas_CO2)then 
+                       interm = indv(nw,igas,jgas)
+                       call interpolateCO2H2(wn_cont,T_cont,p_cross,p_cont,dtempc,.false.,interm)
+                       indv(nw,igas,jgas) = interm
+                       ! might not be relevant in the visible
                     elseif(jgas.eq.igas_He)then 
                        interm = indv(nw,igas,jgas)

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

@@ -732 +732 @@
                   dummy = -9999
                   call interpolateN2H2(592.D+0,278.15D+0,200000.D+0,10000.D+0,testcont,.true.,dummy)
+               elseif (jgas .eq. igas_CO2) then
+                  dummy = -9999
+                  call interpolateCO2H2(592.D+0,278.15D+0,200000.D+0,10000.D+0,testcont,.true.,dummy)
                elseif (jgas .eq. igas_He) then
                   dummy = -9999