Changeset 3803

Timestamp:

Jun 10, 2025, 2:28:53 PM (4 days ago)

Author:

slebonnois

Message:

SL, VENUS PCM : correction following Martin's correction, optcv

File:

• trunk/LMDZ.VENUS/libf/phyvenus/optcv.F90 (modified) (2 diffs)

trunk/LMDZ.VENUS/libf/phyvenus/optcv.F90

@@ -73 +73 @@
   real*8 DRAYAER
   double precision wn_cont, p_cont, p_air, T_cont, dtemp, dtempc
+  double precision p_cross
 
   ! variable species mixing ratio variables
@@ -174 +175 @@
                  p_cont  = dble(PMID(k)*scalep*gfrac(igas)*(1.-QVAR(k)))
               endif
-
-              dtemp=0.0
+ 
+              do jgas=1,ngasmx
+                if(gfrac(jgas).eq.-1)then ! variable
+                  p_cross  = dble(PMID(k)*scalep*QVAR(k)) ! qvar = mol/mol
+                else
+                  p_cross  = dble(PMID(k)*scalep*gfrac(jgas)*(1.-QVAR(k)))
+                endif
+               
+                dtemp=0.0
 
 ! For Venus: only H2O and CO2
-              do jgas=1,ngasmx
                 if ( ((igas .eq. igas_H2O) .and. (jgas .eq. igas_H2O)) .or.  &
-                     ((igas .eq. igas_H2O) .and. (jgas .eq. igas_CO2)) .or.   &
+                     ((igas .eq. igas_H2O) .and. (jgas .eq. igas_CO2)) .or.  &
                      ((igas .eq. igas_CO2) .and. (jgas .eq. igas_CO2)) ) then
-                  call interpolate_continuum('',igas,jgas,'VI',nw,T_cont,p_cont,p_cont,dtemp,.false.)
+                     
+                  call interpolate_continuum('',igas,jgas,'VI',nw,T_cont,p_cont,p_cross,dtemp,.false.)
+
                 endif
-              enddo
-
-              DCONT = DCONT + dtemp
-
-           enddo
+                
+                DCONT = DCONT + dtemp
+ 
+              enddo ! jgas=1,ngasmx
+           enddo ! igas=1,ngasmx
 
            DCONT = DCONT*dz(k)