Changeset 2200 for trunk/LMDZ.VENUS/libf/phyvenus/cloudvenus

Timestamp:

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

Author:

flefevre

Message:

correction de bugs dans le supercycling de la chimie

File:

• trunk/LMDZ.VENUS/libf/phyvenus/cloudvenus/new_cloud_sedim.F (modified) (7 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