Changeset 2860


Ignore:
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
Files:
1 added
4 edited

Legend:

Unmodified
Added
Removed
  • trunk/LMDZ.GENERIC/README

    r2831 r2860  
    17641764aerovenus3, aerovenusUV flags which may be specified in callphys.def).
    17651765
     1766== 04/01/2023 == MT
     1767- Add the possibility to include CO2-H2 CIA opacity (based on Turbet et al. 2020, Icarus, Volume 346, article id. 113762)
     1768- Add interpolateCO2H2 routine + updated associated RT routines
     1769(Be careful, if you use the CO2-H2 CIA, to use the associated correlated-k table, with CO2 and H2, and not just CO2)
     1770
  • trunk/LMDZ.GENERIC/libf/phystd/optci.F90

    r2655 r2860  
    200200                       call interpolateN2H2(wn_cont,T_cont,p_cross,p_cont,dtempc,.false.,interm)
    201201                       indi(nw,igas,jgas) = interm
     202                    elseif(jgas.eq.igas_CO2)then
     203                       interm = indi(nw,igas,jgas)
     204                       call interpolateCO2H2(wn_cont,T_cont,p_cross,p_cont,dtempc,.false.,interm)
     205                       indi(nw,igas,jgas) = interm
    202206                    elseif(jgas.eq.igas_He)then
    203207                       interm = indi(nw,igas,jgas)
  • trunk/LMDZ.GENERIC/libf/phystd/optcv.F90

    r2655 r2860  
    208208                       indv(nw,igas,jgas) = interm
    209209                       ! should be irrelevant in the visible
     210                    elseif(jgas.eq.igas_CO2)then
     211                       interm = indv(nw,igas,jgas)
     212                       call interpolateCO2H2(wn_cont,T_cont,p_cross,p_cont,dtempc,.false.,interm)
     213                       indv(nw,igas,jgas) = interm
     214                       ! might not be relevant in the visible
    210215                    elseif(jgas.eq.igas_He)then
    211216                       interm = indv(nw,igas,jgas)
  • trunk/LMDZ.GENERIC/libf/phystd/sugas_corrk.F90

    r2655 r2860  
    732732                  dummy = -9999
    733733                  call interpolateN2H2(592.D+0,278.15D+0,200000.D+0,10000.D+0,testcont,.true.,dummy)
     734               elseif (jgas .eq. igas_CO2) then
     735                  dummy = -9999
     736                  call interpolateCO2H2(592.D+0,278.15D+0,200000.D+0,10000.D+0,testcont,.true.,dummy)
    734737               elseif (jgas .eq. igas_He) then
    735738                  dummy = -9999
Note: See TracChangeset for help on using the changeset viewer.