Changeset 5473

Timestamp:

Jan 13, 2025, 11:22:23 AM (7 hours ago)

Author:

jyg

Message:

output the total wet deposit of tracers (convective + large scale)

Location:

LMDZ6/trunk/libf/phylmd

Files:

• phys_output_ctrlout_mod.F90 (modified) (1 diff)
• phys_output_mod.F90 (modified) (2 diffs)
• phys_output_write_mod.F90 (modified) (3 diffs)
• phytrac_mod.f90 (modified) (6 diffs)

LMDZ6/trunk/libf/phylmd/phys_output_ctrlout_mod.F90

@@ -2001 +2001 @@
   TYPE(ctrl_out), SAVE, ALLOCATABLE :: o_dtr_sat(:)
   TYPE(ctrl_out), SAVE, ALLOCATABLE :: o_dtr_uscav(:)
-  TYPE(ctrl_out), SAVE, ALLOCATABLE :: o_dtr_wet_con(:)
+  TYPE(ctrl_out), SAVE, ALLOCATABLE :: o_dtr_wet_cv(:)
+  TYPE(ctrl_out), SAVE, ALLOCATABLE :: o_dtr_wet(:)
   TYPE(ctrl_out), SAVE, ALLOCATABLE :: o_dtr_dry(:)
 

LMDZ6/trunk/libf/phylmd/phys_output_mod.F90

@@ -172 +172 @@
     ALLOCATE(o_dtr_evapls(nqtot),o_dtr_ls(nqtot),o_dtr_trsp(nqtot))
     ALLOCATE(o_dtr_sscav(nqtot),o_dtr_sat(nqtot),o_dtr_uscav(nqtot))
-    ALLOCATE(o_dtr_wet_con(nqtot))
+    ALLOCATE(o_dtr_wet_cv(nqtot), o_dtr_wet(nqtot))
     ALLOCATE(o_dtr_dry(nqtot),o_dtr_vdf(nqtot))
 IF (CPPKEY_STRATAER) THEN
@@ -542 +542 @@
 
             lnam = 'tracer convective wet deposition'//TRIM(tracers(iq)%longName)
-            tnam = TRIM(dn)//'wet_con';       o_dtr_wet_con       (itr) = ctrl_out(flag, tnam, lnam, "-", [('',i=1,nfiles)])
+            tnam = TRIM(dn)//'wet_cv';       o_dtr_wet_cv       (itr) = ctrl_out(flag, tnam, lnam, "-", [('',i=1,nfiles)])
+            lnam = 'tracer total wet deposition'//TRIM(tracers(iq)%longName)
+            tnam = TRIM(dn)//'wet';       o_dtr_wet       (itr) = ctrl_out(flag, tnam, lnam, "-", [('',i=1,nfiles)])
             lnam = 'tracer tendency dry deposition'//TRIM(tracers(iq)%longName)
             tnam = 'cum'//TRIM(dn)//'dry';  o_dtr_dry       (itr) = ctrl_out(flag, tnam, lnam, "-", [('',i=1,nfiles)])

LMDZ6/trunk/libf/phylmd/phys_output_write_mod.F90

@@ -6 +6 @@
   USE phytrac_mod, ONLY : d_tr_cl, d_tr_th, d_tr_cv, d_tr_lessi_impa, &
        d_tr_lessi_nucl, d_tr_insc, d_tr_bcscav, d_tr_evapls, d_tr_ls,  &
-       d_tr_trsp, d_tr_sscav, d_tr_sat, d_tr_uscav, flux_tr_wet, flux_tr_dry
+       d_tr_trsp, d_tr_sscav, d_tr_sat, d_tr_uscav,  &
+       flux_tr_wet_cv, flux_tr_wet, flux_tr_dry
 
   ! Author: Abderrahmane IDELKADI (original include file)
@@ -189 +190 @@
          o_dtr_insc, o_dtr_bcscav, o_dtr_evapls, &
          o_dtr_ls, o_dtr_trsp, o_dtr_sscav, o_dtr_dry, &
-         o_dtr_sat, o_dtr_uscav, o_dtr_wet_con, & 
+         o_dtr_sat, o_dtr_uscav, o_dtr_wet_cv, o_dtr_wet, & 
          o_trac_cum, o_du_gwd_rando, o_dv_gwd_rando, &
          o_ustr_gwd_hines,o_vstr_gwd_hines,o_ustr_gwd_rando,o_vstr_gwd_rando, &
@@ -2857 +2858 @@
              CALL histwrite_phy(o_dtr_uscav(itr),d_tr_uscav(:,:,itr))
             !--2D fields
-             CALL histwrite_phy(o_dtr_wet_con(itr), flux_tr_wet(:,itr))
+             CALL histwrite_phy(o_dtr_wet_cv(itr), flux_tr_wet_cv(:,itr))
+             CALL histwrite_phy(o_dtr_wet(itr), flux_tr_wet(:,itr))
              CALL histwrite_phy(o_dtr_dry(itr), flux_tr_dry(:,itr))
              zx_tmp_fi2d=0.

LMDZ6/trunk/libf/phylmd/phytrac_mod.f90

@@ -35 +35 @@
   REAL,DIMENSION(:,:,:),ALLOCATABLE,SAVE :: d_tr_sat
   REAL,DIMENSION(:,:,:),ALLOCATABLE,SAVE :: d_tr_uscav
-  REAL,DIMENSION(:,:),ALLOCATABLE,SAVE   :: flux_tr_wet ! tracer wet deposit (surface)                    jyg
+  REAL,DIMENSION(:,:),ALLOCATABLE,SAVE   :: flux_tr_wet    ! tracer wet deposit (surface)                    jyg
+  REAL,DIMENSION(:,:),ALLOCATABLE,SAVE   :: flux_tr_wet_cv ! tracer convective wet deposit (surface)         jyg
   REAL,DIMENSION(:,:,:),ALLOCATABLE,SAVE :: qPr,qDi ! concentration tra dans pluie,air descente insaturee
   REAL,DIMENSION(:,:,:),ALLOCATABLE,SAVE :: qPa,qMel
@@ -48 +49 @@
 
 !$OMP THREADPRIVATE(qPa,qMel,qTrdi,dtrcvMA,d_tr_th,d_tr_lessi_impa,d_tr_lessi_nucl)
-!$OMP THREADPRIVATE(d_tr_trsp,d_tr_sscav,d_tr_sat,d_tr_uscav,flux_tr_wet,qPr,qDi)
+!$OMP THREADPRIVATE(d_tr_trsp,d_tr_sscav,d_tr_sat,d_tr_uscav)
+!$OMP THREADPRIVATE(flux_tr_wet,flux_tr_wet_cv,qPr,qDi)
 !$OMP THREADPRIVATE(d_tr_insc,d_tr_bcscav,d_tr_evapls,d_tr_ls,qPrls)
 !$OMP THREADPRIVATE(d_tr_cl,d_tr_dry,flux_tr_dry,d_tr_dec,d_tr_cv)
@@ -69 +71 @@
     ALLOCATE(d_tr_sscav(klon,klev,nbtr),d_tr_sat(klon,klev,nbtr))
     ALLOCATE(d_tr_uscav(klon,klev,nbtr),qPr(klon,klev,nbtr),qDi(klon,klev,nbtr))
-    ALLOCATE(flux_tr_wet(klon,nbtr))
+    ALLOCATE(flux_tr_wet(klon,nbtr),flux_tr_wet_cv(klon,nbtr))
     ALLOCATE(qPa(klon,klev,nbtr),qMel(klon,klev,nbtr))
     ALLOCATE(qTrdi(klon,klev,nbtr),dtrcvMA(klon,klev,nbtr))
@@ -411 +413 @@
           flux_tr_dry(i,it)=0.
           flux_tr_wet(i,it)=0.
+          flux_tr_wet_cv(i,it)=0.
        ENDDO
     ENDDO
@@ -700 +703 @@
                 !--with the full array tr_seri even if only item it is processed
 
-                CALL cvltr_scav(pdtphys, da, phi,phi2,d1a,dam, mp,ep,                &
-                     sigd,sij,wght_cvfd,clw,elij,epmlmMm,eplaMm,                     &     
-                     pmflxr,pmflxs,evap,t_seri,wdtrainA,wdtrainM,                    &   
-                     paprs,it,tr_seri,upwd,dnwd,itop_con,ibas_con,                   &
-                     ccntrAA_3d,ccntrENV_3d,coefcoli_3d,                             &
-                     d_tr_cv,d_tr_trsp,d_tr_sscav,d_tr_sat,d_tr_uscav,flux_tr_wet,   &
-                     qDi,qPr,                                                        &
-                     qPa,qMel,qTrdi,dtrcvMA,Mint,                                    &
+                CALL cvltr_scav(pdtphys, da, phi,phi2,d1a,dam, mp,ep,                 &
+                     sigd,sij,wght_cvfd,clw,elij,epmlmMm,eplaMm,                      &    
+                     pmflxr,pmflxs,evap,t_seri,wdtrainA,wdtrainM,                     &  
+                     paprs,it,tr_seri,upwd,dnwd,itop_con,ibas_con,                    &
+                     ccntrAA_3d,ccntrENV_3d,coefcoli_3d,                              &
+                     d_tr_cv,d_tr_trsp,d_tr_sscav,d_tr_sat,d_tr_uscav,flux_tr_wet_cv, &
+                     qDi,qPr,                                                         &
+                     qPa,qMel,qTrdi,dtrcvMA,Mint,                                     &
                      zmfd1a,zmfphi2,zmfdam)
 
@@ -923 +926 @@
                            beta_v1,pplay,paprs,t_seri,tr_seri,d_tr_insc,d_tr_bcscav,d_tr_evapls,qPrls)
 
+             !total wet deposit = large scale wet deposit + convective wet deposit
+             DO i = 1, klon
+               flux_tr_wet(i, it) = flux_tr_wet_cv(i, it) + &
+                                    qPrls(i, it)*(prfl(i, 1)+psfl(i, 1))*pdtphys
+             ENDDO  ! i = 1, klon
+
              !large scale scavenging tendency
              DO k = 1, klev