Changeset 3793 for trunk/LMDZ.GENERIC/libf/phystd/optcv.F90

Timestamp:

Jun 5, 2025, 1:10:40 PM (3 days ago)

Author:

mturbet

Message:

fix a bug and update new continuum interpolation routine in optci and optcv

File:

• trunk/LMDZ.GENERIC/libf/phystd/optcv.F90 (modified) (3 diffs)

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

@@ -195 +195 @@
   
   !     we ignore K=1...
+  
   do K=2,L_LEVELS
 
      do NW=1,L_NSPECTV
-
+     
         DRAYAER = TRAY(K,NW)
         !     DRAYAER is Tau RAYleigh scattering, plus AERosol opacity
@@ -222 +223 @@
               endif
              
-              dtemp=0.0
-              
-              if (igas .eq. igas_N2) then
-               call interpolate_continuum('',igas_N2,igas_N2,'VI',nw,T_cont,p_cont,p_cont,dtemp,.false.)
-               do jgas=1,ngasmx
-                if (jgas .eq. igas_H2) then
-                 call interpolate_continuum('',igas_N2,igas_H2,'VI',nw,T_cont,p_cont,p_cont,dtemp,.false.)
-                elseif (jgas .eq. igas_O2) then
-                 call interpolate_continuum('',igas_N2,igas_O2,'VI',nw,T_cont,p_cont,p_cont,dtemp,.false.)
-                elseif (jgas .eq. igas_CH4) then
-                 call interpolate_continuum('',igas_N2,igas_CH4,'VI',nw,T_cont,p_cont,p_cont,dtemp,.false.)
+              do jgas=1,ngasmx
+                if(gfrac(jgas).eq.-1)then ! variable
+                  p_cross  = dble(PMID(k)*scalep*QVAR(k)) ! qvar = mol/mol
+                elseif(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
-               enddo
-              elseif (igas .eq. igas_O2) then
-               call interpolate_continuum('',igas_O2,igas_O2,'VI',nw,T_cont,p_cont,p_cont,dtemp,.false.)
-               do jgas=1,ngasmx
-                if (jgas .eq. igas_CO2) then
-                 call interpolate_continuum('',igas_CO2,igas_O2,'VI',nw,T_cont,p_cont,p_cont,dtemp,.false.)
+               
+                dtemp=0.0
+
+                if ( ((igas .eq. igas_N2) .and. (jgas .eq. igas_N2)) .or.    &
+                     ((igas .eq. igas_N2) .and. (jgas .eq. igas_H2)) .or.    &
+                     ((igas .eq. igas_N2) .and. (jgas .eq. igas_O2)) .or.    &
+                     ((igas .eq. igas_N2) .and. (jgas .eq. igas_CH4)) .or.   &
+                     ((igas .eq. igas_O2) .and. (jgas .eq. igas_O2)) .or.    &
+                     ((igas .eq. igas_CO2) .and. (jgas .eq. igas_O2)) .or.   &
+                     ((igas .eq. igas_H2) .and. (jgas .eq. igas_H2)) .or.    &
+                     ((igas .eq. igas_H2) .and. (jgas .eq. igas_CH4)) .or.   &
+                     ((igas .eq. igas_H2) .and. (jgas .eq. igas_He)) .or.    &
+                     ((igas .eq. igas_CH4) .and. (jgas .eq. igas_CH4)) .or.  &
+                     ((igas .eq. igas_H2O) .and. (jgas .eq. igas_H2O)) .or.  &
+                     ((igas .eq. igas_H2O) .and. (jgas .eq. igas_N2)) .or.   &
+                     ((igas .eq. igas_H2O) .and. (jgas .eq. igas_O2)) .or.   &
+                     ((igas .eq. igas_H2O) .and. (jgas .eq. igas_CO2)) .or.  &
+                     ((igas .eq. igas_CO2) .and. (jgas .eq. igas_CO2)) .or.  &
+                     ((igas .eq. igas_CO2) .and. (jgas .eq. igas_H2)) .or.   &
+                     ((igas .eq. igas_CO2) .and. (jgas .eq. igas_CH4)) ) then
+
+                  call interpolate_continuum('',igas,jgas,'VI',nw,T_cont,p_cont,p_cross,dtemp,.false.)
+
                 endif
-               enddo
-              elseif (igas .eq. igas_H2) then
-               call interpolate_continuum('',igas_H2,igas_H2,'VI',nw,T_cont,p_cont,p_cont,dtemp,.false.)
-               do jgas=1,ngasmx
-                if (jgas .eq. igas_CH4) then
-                 call interpolate_continuum('',igas_H2,igas_CH4,'VI',nw,T_cont,p_cont,p_cont,dtemp,.false.)
-                elseif (jgas .eq. igas_He) then
-                 call interpolate_continuum('',igas_H2,igas_He,'VI',nw,T_cont,p_cont,p_cont,dtemp,.false.)
-                endif
-               enddo     
-              elseif (igas .eq. igas_CH4) then
-               call interpolate_continuum('',igas_CH4,igas_CH4,'VI',nw,T_cont,p_cont,p_cont,dtemp,.false.)
-              elseif (igas .eq. igas_H2O) then
-               call interpolate_continuum('',igas_H2O,igas_H2O,'VI',nw,T_cont,p_cont,p_cont,dtemp,.false.)
-               do jgas=1,ngasmx
-                if (jgas .eq. igas_N2) then
-                 call interpolate_continuum('',igas_H2O,igas_N2,'VI',nw,T_cont,p_cont,p_cont,dtemp,.false.)
-                elseif (jgas .eq. igas_O2) then
-                 call interpolate_continuum('',igas_H2O,igas_O2,'VI',nw,T_cont,p_cont,p_cont,dtemp,.false.)
-                elseif (jgas .eq. igas_CO2) then
-                 call interpolate_continuum('',igas_H2O,igas_CO2,'VI',nw,T_cont,p_cont,p_cont,dtemp,.false.)
-                endif
-               enddo     
-              elseif (igas .eq. igas_CO2) then
-               call interpolate_continuum('',igas_CO2,igas_CO2,'VI',nw,T_cont,p_cont,p_cont,dtemp,.false.)
-               do jgas=1,ngasmx
-                if (jgas .eq. igas_H2) then
-                 call interpolate_continuum('',igas_CO2,igas_H2,'VI',nw,T_cont,p_cont,p_cont,dtemp,.false.)
-                elseif (jgas .eq. igas_CH4) then
-                 call interpolate_continuum('',igas_CO2,igas_CH4,'VI',nw,T_cont,p_cont,p_cont,dtemp,.false.)
-                endif
-               enddo
-              endif
-               
-              DCONT = DCONT + dtemp
+                
+                DCONT = DCONT + dtemp
+                
+              enddo ! jgas=1,ngasmx
                
             enddo ! igas=1,ngasmx
@@ -456 +439 @@
   end do
 
-
   !=======================================================================
   !     Now the full treatment for the layers, where besides the opacity