Changeset 2200 for trunk/LMDZ.VENUS

Timestamp:

Dec 16, 2019, 8:55:01 AM (5 years ago)

Author:

flefevre

Message:

correction de bugs dans le supercycling de la chimie

Location:

trunk/LMDZ.VENUS/libf/phyvenus

Files:

• cloudvenus/new_cloud_sedim.F (modified) (7 diffs)
• physiq_mod.F (modified) (1 diff)
• phytrac_chimie.F (modified) (2 diffs)

trunk/LMDZ.VENUS/libf/phyvenus/cloudvenus/new_cloud_sedim.F

@@ -1 @@
-      SUBROUTINE new_cloud_sedim(n_lon,n_lev,ptimestep,
-     &                pmidlay,pbndlay, 
-     &                pt,
-     &                pq, pdqsed,pdqs_sed,nq,F_sed)
+      subroutine new_cloud_sedim(n_lon, n_lev, ptimestep,
+     $                           pmidlay, pbndlay, pt, pq,
+     $                           d_tr_chem, pdqsed, pdqs_sed,
+     $                           nq, F_sed)
 
       USE ioipsl
@@ -8 +8 @@
       USE chemparam_mod
       IMPLICIT NONE
-
-c=======================================================================
-c      
-c=======================================================================
 
 c-----------------------------------------------------------------------
@@ -21 +17 @@
 c#include "tracer.h"
 c#include "callkeys.h"
-        
 c
 c   arguments:
@@ -36 +31 @@
       integer nq                    ! number of tracers
       real pq(n_lon,n_lev,nq)       ! tracers (kg/kg)
-c      real pdqfi(n_lon,n_lev,nq)   ! tendency before sedimentation (kg/kg.s-1)
       real pdqsed(n_lon,n_lev,2)    ! tendency due to sedimentation (kg/kg)
+      real d_tr_chem(n_lon,n_lev,nq)! tendency due to chemistry and clouds (kg/kg)
       real pdqs_sed(n_lon)          ! surface density (Flux if /ptimestep) at surface due to sedimentation (kg.m-2)
       
@@ -73 +68 @@
       REAL, PARAMETER :: rho_core = 2500.0
 
-      REAL, DIMENSION(n_lon,n_lev+1) ::
-     + wgt_SA                         ! Fraction of H2SO4 in droplet local
+      REAL, DIMENSION(n_lon,n_lev+1) :: wgt_SA   ! Fraction of H2SO4 in droplet local
 
 c     Stokes speed and sedimentation flux variable
@@ -87 +81 @@
      + F_sed                          ! Flux de sedimentation (kg.m-2.s-1 puis en output kg.m-2)
       
-      
       REAL :: R_mode0                 ! Rayon mode 0 (m), rayon le plus frequent
-
-
 
 !      PRINT*,'RHO_DROPLET new_cloud_sedim.F'
@@ -100 +91 @@
 c    -----------------
     
-c     Updating the droplet mass mixing ratio with the partition H2O/H2SO4
-c     ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
-
-      do l=1,n_lev
-         do ig=1,n_lon
-         zqi_wv(ig,l) = pq(ig,l,i_h2oliq)
-         zqi_sa(ig,l) = pq(ig,l,i_h2so4liq)
-         wgt_SA(ig,l) = WH2SO4(ig,l)
-         enddo
-      enddo 
+!     Updating the droplet mass mixing ratio with the partition H2O/H2SO4
+
+      zqi_wv(:,:) = pq(:,:,i_h2oliq) + d_tr_chem(:,:,i_h2oliq)*ptimestep
+      zqi_sa(:,:) = pq(:,:,i_h2so4liq) 
+     $            + d_tr_chem(:,:,i_h2so4liq)*ptimestep
+      wgt_SA(:,:) = wh2so4(:,:)
 
 c     Init F_sed

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

@@ -1013 +1013 @@
      $                              t_seri,
      $                              tr_seri,
+     $                              d_tr_chem,
      $                              d_tr_sed(:,:,1:2),
      $                              d_tr_ssed,

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

@@ -270 +270 @@
 
       do iq = 1,nqmax - nmicro
-         ztrac(:,:,iq) = ztrac(:,:,iq)*m_tr(iq)/mmean(:,:)
+         ztrac(:,:,iq) = max(ztrac(:,:,iq)*m_tr(iq)/mmean(:,:), 
+     $                       1.e-30)
          d_tr_chem(:,:,iq) = (ztrac(:,:,iq) - trac(:,:,iq))/pdtphys
       end do
@@ -277 +278 @@
 
       if (ok_cloud .and. cl_scheme == 1) then
-         ztrac(:,:,i_h2so4liq) = ztrac(:,:,i_h2so4liq)*m_tr(i_h2so4liq)
-     $                           /mmean(:,:)
-         ztrac(:,:,i_h2oliq) = ztrac(:,:,i_h2oliq)*m_tr(i_h2oliq)
-     $                         /mmean(:,:)
+         ztrac(:,:,i_h2so4liq) = max(ztrac(:,:,i_h2so4liq)
+     $                               *m_tr(i_h2so4liq)/mmean(:,:),
+     $                               1.e-30)
+         ztrac(:,:,i_h2oliq)   = max(ztrac(:,:,i_h2oliq)
+     $                               *m_tr(i_h2oliq)/mmean(:,:),
+     $                               1.e-30)
          d_tr_chem(:,:,i_h2so4liq) = (ztrac(:,:,i_h2so4liq) 
      $                              - trac(:,:,i_h2so4liq))/pdtphys