Changeset 4055 for trunk/LMDZ.GENERIC/libf/phygeneric/optcv.F90

Timestamp:

Feb 6, 2026, 10:20:21 PM (3 weeks ago)

Author:

mturbet

Message:

Cleaning of continuum routines following 04-02-26 hackathon. RIP our beloved bilinear and bilinearbig routines. You will be deeply missed...

File:

• trunk/LMDZ.GENERIC/libf/phygeneric/optcv.F90 (modified) (5 diffs)

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

@@ -10 +10 @@
 
   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 radcommon_h, only: gasv, tlimit, wrefVAR, Cmk, tgasref, pfgasref,wnov,scalep,glat_ig
   use gases_h, only: gfrac, ngasmx, igas_H2, igas_H2O, igas_He, igas_N2, &
                      igas_CH4, igas_CO2, igas_O2
   use comcstfi_mod, only: g, r, mugaz
   use callkeys_mod, only: kastprof,continuum,graybody,callgasvis,varspec, &
-                          generic_continuum_database,rayleigh
+                          rayleigh
   use recombin_corrk_mod, only: corrk_recombin, gasv_recomb
   use tpindex_mod, only: tpindex
@@ -107 +107 @@
   integer igas, jgas
 
-  integer interm
-
   logical :: firstcall=.true.
 !$OMP THREADPRIVATE(firstcall)
@@ -120 +118 @@
     firstcall=.false.
   endif ! of if (firstcall)
-
-  !! AS: to save time in computing continuum (see bilinearbig)
-  IF (.not.ALLOCATED(indv)) THEN
-      ALLOCATE(indv(L_NSPECTV,ngasmx,ngasmx))
-      indv = -9999 ! this initial value means "to be calculated"
-  ENDIF
 
   !=======================================================================
@@ -211 +203 @@
            ! include continua if necessary
            
-          if(generic_continuum_database)then
             T_cont  = dble(TMID(k))
             do igas=1,ngasmx
@@ -261 +252 @@
                
             enddo ! igas=1,ngasmx
-            
-          else ! generic_continuum_database
-           
-           
-           wn_cont = dble(wnov(nw))
-           T_cont  = dble(TMID(k))
-           do igas=1,ngasmx
-
-              if(gfrac(igas).eq.-1)then ! variable
-                 p_cont  = dble(PMID(k)*scalep*QVAR(k)) ! qvar = mol/mol
-              elseif(varspec) then
-                 p_cont  = dble(PMID(k)*scalep*FRACVAR(igas,k)*(1.-QVAR(k)))
-              else
-                 p_cont  = dble(PMID(k)*scalep*gfrac(igas)*(1.-QVAR(k)))
-              endif
-
-              dtemp=0.0
-              if(igas.eq.igas_N2)then
-
-                 interm = indv(nw,igas,igas)
-!                 call interpolateN2N2(wn_cont,T_cont,p_cont,dtemp,.false.,interm)
-                 indv(nw,igas,igas) = interm
-                 ! only goes to 500 cm^-1, so unless we're around a cold brown dwarf, this is irrelevant in the visible
-
-              elseif(igas.eq.igas_H2)then
-
-                 ! first do self-induced absorption
-                 interm = indv(nw,igas,igas)
-                 call interpolateH2H2(wn_cont,T_cont,p_cont,dtemp,.false.,interm)
-                 indv(nw,igas,igas) = interm
-
-                 ! then cross-interactions with other gases
-                 do jgas=1,ngasmx
-                    if(varspec) then
-                      p_cross = dble(PMID(k)*scalep*FRACVAR(jgas,k)*(1.-QVAR(k)))
-                    else
-                      p_cross = dble(PMID(k)*scalep*gfrac(jgas)*(1.-QVAR(k)))
-                    endif
-                    dtempc  = 0.0
-                    if(jgas.eq.igas_N2)then 
-                       interm = indv(nw,igas,jgas)
-                       call interpolateN2H2(wn_cont,T_cont,p_cross,p_cont,dtempc,.false.,interm)
-                       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)
-                       call interpolateH2He(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_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
-                    if(varspec) then
-                      p_cross = dble(PMID(k)*scalep*FRACVAR(jgas,k)*(1.-QVAR(k)))
-                    else
-                      p_cross = dble(PMID(k)*scalep*gfrac(jgas)*(1.-QVAR(k)))
-                    endif
-                    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_CO2)then 
-                       interm = indv(nw,igas,jgas)
-                       call interpolateCO2CH4(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)
-                       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
-                 ! Compute self and foreign (with air) continuum of H2O 
-                 p_air = dble(PMID(k)*scalep) - p_cont ! note assumes background is air!
-                 interm = indv(nw,igas,igas)
-                 call interpolateH2O_self_foreign(wn_cont,T_cont,p_cont,p_air,dtemp,.false.,interm) ! MTCKD v3.3
-                 indv(nw,igas,igas) = interm
-
-              endif
-
-              DCONT = DCONT + dtemp
-
-           enddo
-
-          endif ! generic_continuum_database
           
           DCONT = DCONT*dz(k)