Changeset 2194

Timestamp:

Dec 13, 2019, 11:37:38 AM (5 years ago)

Author:

slebonnois

Message:

SL: ajustement commit precedent pour la microphysique de Sabrina

Location:

trunk/LMDZ.VENUS/libf/phyvenus

Files:

• physiq_mod.F (modified) (2 diffs)
• phytrac_chimie.F (modified) (2 diffs)

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

@@ -1113 +1113 @@
 !        tendency due to sedimentation 
 
-               do iq = nqmax - nmicro + 1,nqmax
-                  tr_seri(:,:,iq)  = tr_seri(:,:,iq) + d_tr_sed(:,:,iq)
+               do iq = nqmax-nmicro+1,nqmax
                   d_tr_sed(:,:,iq) = d_tr_sed(:,:,iq)/zctime
                end do
@@ -1130 +1129 @@
 !        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.
+!====================================================================
 
          end if     ! tr_scheme

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

@@ -269 +269 @@
       end do
 
-!     liquid phase
+!     liquid phase or moments
 
       if (ok_cloud .and. cl_scheme == 1) then
@@ -282 +282 @@
       end if
 
+
+      if (ok_cloud .and. cl_scheme == 2) then
+         do iq = nqmax-nmicro+1,nqmax
+            d_tr_chem(:,:,iq) = (ztrac(:,:,iq) - trac(:,:,iq))/pdtphys
+         end do
+      end if
+
       end subroutine phytrac_chimie