Changeset 4795 for LMDZ6/branches/LMDZ_DECOUPLE/libf/phylmd/lsc_scav.F90

Timestamp:

Feb 1, 2024, 4:53:24 PM (4 months ago)

Author:

nfevrier

Message:

First save of N. Février's modifications

File:

• LMDZ6/branches/LMDZ_DECOUPLE/libf/phylmd/lsc_scav.F90 (modified) (10 diffs)

LMDZ6/branches/LMDZ_DECOUPLE/libf/phylmd/lsc_scav.F90

@@ -7 +7 @@
                     oliq,flxr,flxs,rneb,beta_fisrt,  &
                     beta_v1,pplay,paprs,t,tr_seri,d_tr_insc,          &
-                    d_tr_bcscav,d_tr_evap,qPrls) 
+                    d_tr_bcscav,d_tr_evap,qPrls,beta_st,j) 
   USE ioipsl
   USE dimphy
@@ -13 +13 @@
   USE mod_phys_lmdz_para
   USE traclmdz_mod
-  USE infotrac_phy,ONLY : nbtr
+  USE infotrac_phy,ONLY : nbtr,solsym
   USE iophy
   IMPLICIT NONE 
@@ -47 +47 @@
   REAL,DIMENSION(klon,klev,nbtr),INTENT(OUT)  :: d_tr_bcscav  ! tendance de traceur
   REAL,DIMENSION(klon,klev,nbtr),INTENT(OUT)  :: d_tr_evap
+  ! Nicolas Février
+  REAL,DIMENSION(klon,klev,12),INTENT(OUT)       :: beta_st
+  INTEGER,INTENT(IN)                             :: j           ! indice parmi les 12 traceurs de WETDEP
+
   REAL,DIMENSION(klon,nbtr),INTENT(OUT)       :: qPrls      !jyg: concentration tra dans pluie LS a la surf.
   REAL :: dxin,dxev                              ! tendance temporaire de traceur incloud
@@ -111 +115 @@
       R_r=0.001            !  mean raindrop radius (m) 
       R_s=0.001            !  mean snow crystal radius (m)
-      frac_fine_scav=0.7
-      frac_coar_scav=0.7
+      !frac_fine_scav=0.7
+      !frac_coar_scav=0.7
 !     frac_aer=0.5 ~ droplet size shrinks by evap
       frac_aer=0.5
@@ -133 +137 @@
    print*,'R_r',R_r
    print*,'R_s',R_s
-   print*,'frac_fine_scav',frac_fine_scav
-   print*,'frac_coar_scav',frac_coar_scav
+   !print*,'frac_fine_scav',frac_fine_scav
+   !print*,'frac_coar_scav',frac_coar_scav
    print*,'frac_aer ev',frac_aer
 
@@ -142 +146 @@
 !
 ! initialization
+
+  ! Nicolas Février
+  ! Différentiation BC soluble/insoluble
+
+  IF ((TRIM(solsym(it))=='AIBCM') .OR. (TRIM(solsym(it))=='AIPOMM') .OR. (TRIM(solsym(it))=='CIDUSTM')) THEN
+     frac_fine_scav=0.2
+     frac_coar_scav=0.2
+     !PRINT *,'insoluble okidoki'
+  ELSE
+     frac_fine_scav=0.7
+     frac_coar_scav=0.7
+  ENDIF
+
   dxin=0.
   dxev=0.
@@ -223 +240 @@
       his_dh(i)=his_dh(i)-dxin*zmass(i,k)/pdtime !  kg/m2/s
       d_tr_insc(i,k,it)=dxin
+      ! Nicolas Février
+      beta_st(i,k,j)=beta_st(i,k,j)-(exp(-scav(i,k)*beta*pdtime)-1.)/pdtime
 
 !  below-cloud impaction
@@ -239 +258 @@
      his_dh(i)=his_dh(i)-dxbc(i,k)*zmass(i,k)/pdtime !  kg/m2/s
      d_tr_bcscav(i,k,it)=dxbc(i,k)
+     ! Nicolas Février
+     beta_st(i,k,j)=beta_st(i,k,j)+3./4.*(water+ice)
     ENDIF
 
@@ -279 +300 @@
    ENDDO  !! do k
 
+! Debug beta_st
+    print *,"BETA LSC_SCAV",beta_st(5000,35,7)
+
+
 !jyg (20130114)
    DO i = 1,klon
@@ -298 +323 @@
 !     ENDDO
 
+!print*, "TAGADATSOINTSOIN"
+
 END SUBROUTINE lsc_scav