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

Timestamp:

Dec 24, 2019, 1:57:50 PM (5 years ago)

Author:

flefevre

Message:

• correction du supercycling de la chimie
• changements cosmetiques dans la sedimentation

File:

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

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

@@ -300 +300 @@
       REAL :: d_drop_sed(klev),d_ccn_sed(klev,2),d_liq_sed(klev,2)
       REAL :: aer_flux(klev)
-      REAL :: d_tr_ssed(klon)
 c
 c Variables du changement
@@ -1015 +1014 @@
      $                              d_tr_chem,
      $                              d_tr_sed(:,:,1:2),
-     $                              d_tr_ssed,
      $                              nqmax,
      $                              Fsedim)
@@ -1035 +1033 @@
                end do
 
-!        tendency due to sedimentation 
-
-               d_tr_sed(:,:,:) = d_tr_sed(:,:,:)/zctime
+!        tendency due to condensation and sedimentation 
+
+               d_tr_sed(:,:,1:2) = d_tr_sed(:,:,1:2)/zctime
                Fsedim(:,1:klev) = Fsedim(:,1:klev)/zctime
                Fsedim(:,klev+1) = 0.
@@ -1119 +1117 @@
 #endif
             end if  ! cl_scheme
+
+!        update gaseous tracers (chemistry)
+
+            do iq = 1, nqmax - nmicro
+               tr_seri(:,:,iq) = tr_seri(:,:,iq) 
+     $                         + d_tr_chem(:,:,iq)*zctime
+            end do
+
+!        update condensed tracers (condensation + sedimentation)
+
+            if (cl_scheme == 1) then 
+               tr_seri(:,:,i_h2so4liq) = max(tr_seri(:,:,i_h2so4liq) 
+     $                                 + d_tr_sed(:,:,1)*zctime, 1.e-30)
+               tr_seri(:,:,i_h2oliq)   = max(tr_seri(:,:,i_h2oliq) 
+     $                                 + d_tr_sed(:,:,2)*zctime, 1.e-30)
+            else if (cl_scheme == 2) then
+               do iq = nqmax-nmicro+1,nqmax
+                  tr_seri(:,:,iq) = tr_seri(:,:,iq) 
+     $                            + d_tr_sed(:,:,iq)*zctime
+               end do
+            end if  ! cl_scheme
+
          end if     ! ok_sedim
          end if     ! mod(itap,chempas)  <------- end of chemistry supercycling
 
-!        update tracers (chemistry)
-
-         do iq = 1, nqmax - nmicro
-            tr_seri(:,:,iq) = tr_seri(:,:,iq) + d_tr_chem(:,:,iq)*dtime
-         end do
-
-!        update tracers (sedimentation)
-
-         if (ok_sedim) then 
-            if (cl_scheme == 1) then 
-         tr_seri(:,:,i_h2so4liq) = tr_seri(:,:,i_h2so4liq) 
-     $                           + d_tr_sed(:,:,1)*dtime
-         tr_seri(:,:,i_h2oliq)   = tr_seri(:,:,i_h2oliq) 
-     $                           + d_tr_sed(:,:,2)*dtime
-            else if (cl_scheme == 2) then
-         do iq = nqmax-nmicro+1,nqmax
-            tr_seri(:,:,iq) = tr_seri(:,:,iq) + d_tr_sed(:,:,iq)*dtime
-         end do
-            end if  ! cl_scheme
-         end if     ! ok_sedim
 !====================================================================
 ! End Case 3: Full chemistry and/or clouds.