Changeset 4087

Timestamp:

Feb 27, 2026, 8:22:38 AM (2 weeks ago)

Author:

streel

Message:

First follow-up of the microphysics scheme commit :

• Removal of wsa_new_bis.f90 (was not used)
• Correction of change_wsa call in mad_muphy + cleaner definition of the r_merge
• Correction of the nucleation variable status
• Correction of a bug in physiq_mod concerning the non-activated CCN + cleaning

NS

Location:

trunk/LMDZ.VENUS/libf/phyvenus

Files:

• cloudvenus/mad_muphy.f90 (modified) (2 diffs)
• cloudvenus/nucleations.f90 (modified) (3 diffs)
• cloudvenus/wsa_new_bis.f90 (deleted)
• physiq_mod.F (modified) (3 diffs)

trunk/LMDZ.VENUS/libf/phyvenus/cloudvenus/mad_muphy.f90

@@ -204 +204 @@
   write(*,*)'RHOpi et RHOeqi', RHOSAp, RHOSA   
   ENDIF
-  call change_wsa(TAIR,PAIR,MRSAv,MRSAss,WSAv)
+  call change_wsa(TAIR,PAIR,MRSAv,MRSAss)
   
   IF (dbg_glob) then
@@ -854 +854 @@
      !Computation of the 9th decile
 
-     r_merge = exp(r2 + sigma2 * -1.2816)
+     r_merge = exp(r2 - sigma2 * 1.2816)
 
      !Mode merging 2->1  

trunk/LMDZ.VENUS/libf/phyvenus/cloudvenus/nucleations.f90

@@ -134 +134 @@
 
 !*****************************************************************************
-SUBROUTINE newnuklefit(t,p,rhoa,jnuc,rc,ntotal)
+SUBROUTINE newnuklefit(temp,p,rhoacid,jnuc,rc,ntotal)
 
   !    Fortran 90 subroutine binapara
@@ -165 +165 @@
 
   ! Inputs
-  real,intent(inout) :: t        ! temperature in K (165-400 K)
-  real,intent(inout) :: p        ! pressure in Pa
-  real,intent(inout) :: rhoa     ! sulfuric acid concentration in #/m3 (10^10 - 10^19) 
+  real, intent(in) :: temp     ! temperature in K (165-400 K)
+  real, intent(in) :: p        ! pressure in Pa
+  real, intent(in) :: rhoacid  ! sulfuric acid concentration in #/m3 (10^10 - 10^19) 
 
   ! Outputs
@@ -177 +177 @@
   real :: x          ! mole fraction of H2SO4 in the critical cluster     
   real :: nac        ! local variable for number of acid molecules
+  real :: t          ! local variable for temperature
+  real :: rhoa       ! local variable for acid concentration
 
   ! Initialization
+  rhoa = rhoacid
+  t = temp
   jnuc = 0.e0
   ntotal = 0.e0

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

@@ -1194 +1194 @@
 !        mode 2
 
-               d_tr_sed(:,:,:) = 0.
                call sedim_mad_m2(klon, klev, zctime, pplay, paprs,
      $                           rho,zzlev,t_seri,d_tr_sed(:,:,1:5),
@@ -1279 +1278 @@
                do i = 1, klon
                   if (tr_seri(i,k,i_m3_mode1sa) < 1e-18) then
-                     tr_seri(:,:,i_m0_aer) = tr_seri(:,:,i_m0_aer) 
+                     tr_seri(i,k,i_m0_aer) = tr_seri(i,k,i_m0_aer) 
      $                                     + tr_seri(i,k,i_m0_mode1ccn)
 
-                     tr_seri(:,:,i_m3_aer) = tr_seri(:,:,i_m3_aer) 
+                     tr_seri(i,k,i_m3_aer) = tr_seri(i,k,i_m3_aer) 
      $                                     + tr_seri(i,k,i_m3_mode1ccn)
 
@@ -1293 +1292 @@
        
                   if (tr_seri(i,k,i_m3_mode2sa) < 1e-18) then
-                     tr_seri(:,:,i_m0_aer) = tr_seri(:,:,i_m0_aer)
+                     tr_seri(i,k,i_m0_aer) = tr_seri(i,k,i_m0_aer)
      $                                     + tr_seri(i,k,i_m0_mode2ccn)
 
-                     tr_seri(:,:,i_m3_aer) = tr_seri(:,:,i_m3_aer)
+                     tr_seri(i,k,i_m3_aer) = tr_seri(i,k,i_m3_aer)
      $                                     + tr_seri(i,k,i_m3_mode2ccn)