Changeset 4055 for trunk/LMDZ.GENERIC/libf/phygeneric/sugas_corrk.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/sugas_corrk.F90 (modified) (4 diffs)

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

@@ -34 +34 @@
       use ioipsl_getin_p_mod, only: getin_p
       use callkeys_mod, only: varactive,varfixed,graybody,callgasvis,&
-                continuum,generic_continuum_database
+                continuum
       use tracer_h, only : nqtot, moderntracdef, is_recomb, noms
       use recombin_corrk_mod, only: su_recombin,        &
@@ -62 +62 @@
 
       double precision testcont ! for continuum absorption initialisation
-
-      integer :: dummy
 
       if (.not. moderntracdef) use_premix=.true. ! Added by JVO for compatibility with 'old' traceur.def
@@ -717 +715 @@
 !     Initialise the continuum absorption data
       if(continuum)then
-      if(generic_continuum_database)then
+      
        do igas=1,ngasmx ! we loop on all pairs of molecules that have data available
        ! data can be downloaded from https://web.lmd.jussieu.fr/~lmdz/planets/generic/datagcm/continuum_data/
@@ -806 +804 @@
        enddo ! igas=1,ngasmx
        
-      else ! generic_continuum_database tag
-        do igas=1,ngasmx
-         if (igas .eq. igas_N2) then
-
-            dummy = -9999
-            call interpolateN2N2(100.D+0,250.D+0,17500.D+0,testcont,.true.,dummy)
-
-         elseif (igas .eq. igas_H2) then
-
-            ! first do self-induced absorption
-            dummy = -9999
-            call interpolateH2H2(500.D+0,250.D+0,17500.D+0,testcont,.true.,dummy)
-            ! then cross-interactions with other gases
-            do jgas=1,ngasmx
-               if (jgas .eq. igas_N2) then
-                  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
-                  call interpolateH2He(500.D+0,250.D+0,200000.D+0,10000.D+0,testcont,.true.,dummy)
-               endif
-            enddo
-               
-         
-         elseif (igas .eq. igas_CH4) then
-         
-            ! first do self-induced absorption
-            dummy = -9999
-            call interpolateCH4CH4(600.0,66.7,400.0,testcont,.true.,dummy)
-            ! then cross-interactions with other gases
-            do jgas=1,ngasmx
-               if (jgas .eq. igas_H2) then
-                  dummy = -9999
-                  call interpolateH2CH4(500.D+0,250.D+0,200000.D+0,10000.D+0,testcont,.true.,dummy)
-               elseif (jgas .eq. igas_CO2) then
-                  dummy = -9999
-                  call interpolateCO2CH4(592.D+0,278.15D+0,200000.D+0,10000.D+0,testcont,.true.,dummy)
-               elseif (jgas .eq. igas_He) then
-                  dummy = -9999
-                  call interpolateHeCH4(500.D+0,250.D+0,200000.D+0,10000.D+0,testcont,.true.,dummy)
-               endif
-            enddo
-            
-
-         elseif (igas .eq. igas_H2O) then
-
-            ! Compute self and foreign (with air) continuum of H2O 
-            dummy = -9999
-            call interpolateH2O_self_foreign(990.D+0,296.D+0,683.2D+0*2,0.D+0,testcont,.true.,dummy)  
-
-         endif
-
-      enddo ! igas=1,ngasmx
-      endif ! generic_continuum_database tag
       endif ! continuum flag