Changeset 4108 for trunk/LMDZ.VENUS/libf/phyvenus/physiq_mod.F

Timestamp:

Mar 6, 2026, 3:14:35 PM (5 weeks ago)

Author:

mlefevre

Message:

Correction to previous commits regardion of the polysulfur species in the Venus PCM. ML

File:

• trunk/LMDZ.VENUS/libf/phyvenus/physiq_mod.F (modified) (4 diffs)

trunk/LMDZ.VENUS/libf/phyvenus/physiq_mod.F

@@ -79 +79 @@
       use compo_hedin83_mod2, only: compo_hedin83_init2
       use compo_hedin83_mod2, only: compo_hedin83_mod
-      use moldiff_mpf_mod, only: moldiff_mpf
       use nirdata_mod, only: nir_leedat
       use radlwsw_newtoncool_mod, only: radlwsw_newtoncool
@@ -240 +239 @@
 c
       REAL Fsedim(klon,klev+1)  ! Flux de sedimentation (kg.m-2)
+      REAL Fsedims2(klon,klev+1) ! S2 Flux de sedimentation (kg.m-2)
+      REAL Fsedims3(klon,klev+1) ! S3 Flux de sedimentation (kg.m-2)
+      REAL Fsedims4(klon,klev+1) ! S4 Flux de sedimentation (kg.m-2)
+      REAL Fsedims8(klon,klev+1) ! S8 Flux de sedimentation (kg.m-2)      
 
 c
@@ -1168 +1171 @@
 #endif
             end if  ! cl_scheme
+            IF (1.eq.1) then
+               call polysulfur_sedim(klon,
+     $                               nlev,
+     $                               zctime,
+     $                               pplay,
+     $                               paprs,
+     $                               t_seri,
+     $                               tr_seri,
+     $                               d_tr_chem,
+     $                               d_tr_sedsn(:,:,1:4),
+     $                               nqmax,
+     $                               Fsedims2,
+     $                               Fsedims3,
+     $                               Fsedims4,
+     $                               Fsedims8)
+
+               d_tr_sedsn(:,:,1:4) = d_tr_sedsn(:,:,1:4)/zctime
+               Fsedims2(:,1:klev) = Fsedims2(:,1:klev)/zctime
+               Fsedims2(:,klev+1) = 0.
+               Fsedims3(:,1:klev) = Fsedims3(:,1:klev)/zctime
+               Fsedims3(:,klev+1) = 0.
+               Fsedims4(:,1:klev) = Fsedims4(:,1:klev)/zctime
+               Fsedims4(:,klev+1) = 0.
+               Fsedims8(:,1:klev) = Fsedims8(:,1:klev)/zctime
+               Fsedims8(:,klev+1) = 0.
+            end if            
          end if     ! ok_sedim
 
@@ -1275 +1304 @@
 
          end if     ! cl_scheme          <------- end of microphysics option
+         if (1.eq.1) then
+               tr_seri(:,:,i_s2con) = max(tr_seri(:,:,i_s2con)
+     $                              + d_tr_sedsn(:,:,1)*zctime, 1.e-30)
+               tr_seri(:,:,i_s3con) = max(tr_seri(:,:,i_s3con)
+     $                              + d_tr_sedsn(:,:,2)*zctime, 1.e-30)
+               tr_seri(:,:,i_s4con) = max(tr_seri(:,:,i_s4con)
+     $                              + d_tr_sedsn(:,:,3)*zctime, 1.e-30)
+               tr_seri(:,:,i_s8con) = max(tr_seri(:,:,i_s8con)
+     $                              + d_tr_sedsn(:,:,4)*zctime, 1.e-30)
+         end if          
          end if     ! mod(itap,chempas)  <------- end of chemistry supercycling