Changeset 2861 for trunk/LMDZ.GENERIC


Ignore:
Timestamp:
Jan 4, 2023, 1:01:09 PM (2 years ago)
Author:
mturbet
Message:

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

Location:
trunk/LMDZ.GENERIC
Files:
1 added
5 edited

Legend:

Unmodified
Added
Removed

  • trunk/LMDZ.GENERIC/README

    r2860 r2861  
    17691769(Be careful, if you use the CO2-H2 CIA, to use the associated correlated-k table, with CO2 and H2, and not just CO2)
    17701770
     1771== 04/01/2023 bis == MT
     1772- Add the possibility to include CO2-CH4 CIA opacity (based on Turbet et al. 2020, Icarus, Volume 346, article id. 113762)
     1773- Add interpolateCO2CH4 routine + updated associated RT routines

  • trunk/LMDZ.GENERIC/libf/phystd/interpolateCO2H2.F90

    r2860 r2861  
    1010  !     Authors
    1111  !     -------
    12   !     M. Turbet (2011)
     12  !     M. Turbet (2023)
    1313  !     
    1414  !==================================================================
    … …  
    7272  if(firstcall)then ! called by sugas_corrk only
    7373     print*,'----------------------------------------------------'
    74      print*,'Initialising CO2-H2 continuum from HITRAN database...'
     74     print*,'Initialising CO2-H2 continuum from Turbet et al. 2020'
    7575
    7676     !     1.1 Open the ASCII files

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

    r2860 r2861  
    227227                       call interpolateH2CH4(wn_cont,T_cont,p_cross,p_cont,dtempc,.false.,interm)
    228228                       indi(nw,igas,jgas) = interm
     229                    elseif(jgas.eq.igas_CO2)then
     230                       interm = indi(nw,igas,jgas)
     231                       call interpolateCO2CH4(wn_cont,T_cont,p_cross,p_cont,dtempc,.false.,interm)
     232                       indi(nw,igas,jgas) = interm
    229233                    elseif(jgas.eq.igas_He)then
    230234                       interm = indi(nw,igas,jgas)

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

    r2860 r2861  
    236236                       call interpolateH2CH4(wn_cont,T_cont,p_cross,p_cont,dtempc,.false.,interm)
    237237                       indv(nw,igas,jgas) = interm
     238                    elseif(jgas.eq.igas_CO2)then
     239                       interm = indv(nw,igas,jgas)
     240                       call interpolateCO2CH4(wn_cont,T_cont,p_cross,p_cont,dtempc,.false.,interm)
     241                       indv(nw,igas,jgas) = interm
     242                       ! might not be relevant in the visible
    238243                    elseif(jgas.eq.igas_He)then
    239244                       interm = indv(nw,igas,jgas)

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

    r2860 r2861  
    752752                  dummy = -9999
    753753                  call interpolateH2CH4(500.D+0,250.D+0,200000.D+0,10000.D+0,testcont,.true.,dummy)
     754               elseif (jgas .eq. igas_CO2) then
     755                  dummy = -9999
     756                  call interpolateCO2CH4(592.D+0,278.15D+0,200000.D+0,10000.D+0,testcont,.true.,dummy)
    754757               elseif (jgas .eq. igas_He) then
    755758                  dummy = -9999
Note: See TracChangeset for help on using the changeset viewer.