Changeset 5447 for LMDZ6/trunk/libf


Ignore:
Timestamp:
Dec 22, 2024, 5:48:23 PM (15 hours ago)
Author:
jyg
Message:

output the convective wet deposit of tracers

Location:
LMDZ6/trunk/libf/phylmd
Files:
5 edited

Legend:

Unmodified
Added
Removed
  • LMDZ6/trunk/libf/phylmd/cvltr_scav.f90

    r5292 r5447  
    33!
    44SUBROUTINE cvltr_scav(pdtime, da, phi,phi2,d1a,dam, mpIN,epIN, &
    5      sigd,sij,wght_cvfd,clw,elij,epmlmMm,eplaMm,    &
    6      pmflxrIN,pmflxsIN,ev,te,wdtrainA,wdtrainM,     &
    7      paprs,it,tr,upd,dnd,inb,icb,                   &
    8      ccntrAA_3d,ccntrENV_3d,coefcoli_3d,            &
    9      dtrcv,trsptd,dtrSscav,dtrsat,dtrUscav,qDi,qPr, &
    10      qPa,qMel,qTrdi,dtrcvMA,Mint,                   &
     5     sigd,sij,wght_cvfd,clw,elij,epmlmMm,eplaMm,        &
     6     pmflxrIN,pmflxsIN,ev,te,wdtrainA,wdtrainM,         &
     7     paprs,it,tr,upd,dnd,inb,icb,                       &
     8     ccntrAA_3d,ccntrENV_3d,coefcoli_3d,                &
     9     dtrcv,trsptd,dtrSscav,dtrsat,dtrUscav,flux_tr_wet, &
     10     qDi,qPr,                                           &
     11     qPa,qMel,qTrdi,dtrcvMA,Mint,                       &
    1112     zmfd1a,zmfphi2,zmfdam)
    1213  !
     
    7273  REAL,DIMENSION(klon,klev,nbtr),INTENT(OUT)     :: dtrsat    ! tendance trsp+sat scav
    7374  REAL,DIMENSION(klon,klev,nbtr),INTENT(OUT)     :: dtrUscav  ! tendance du lessivage courant unsat
     75  REAL,DIMENSION(klon,nbtr),     INTENT(OUT)     :: flux_tr_wet ! wet deposit
    7476  !
    7577  ! Variables locales
     
    685687     ENDDO
    686688  ENDDO
     689  DO i=1, klon
     690     flux_tr_wet(i,it) = (pmflxr(i,1)+pmflxs(i,1))*qPr(i,1,it)*pdtime    ! wet deposit
     691  ENDDO
    687692
    688693  ! test de conservation du traceur
  • LMDZ6/trunk/libf/phylmd/phys_output_ctrlout_mod.F90

    r5383 r5447  
    20012001  TYPE(ctrl_out), SAVE, ALLOCATABLE :: o_dtr_sat(:)
    20022002  TYPE(ctrl_out), SAVE, ALLOCATABLE :: o_dtr_uscav(:)
     2003  TYPE(ctrl_out), SAVE, ALLOCATABLE :: o_dtr_wet_con(:)
    20032004  TYPE(ctrl_out), SAVE, ALLOCATABLE :: o_dtr_dry(:)
    20042005
  • LMDZ6/trunk/libf/phylmd/phys_output_mod.F90

    r5394 r5447  
    172172    ALLOCATE(o_dtr_evapls(nqtot),o_dtr_ls(nqtot),o_dtr_trsp(nqtot))
    173173    ALLOCATE(o_dtr_sscav(nqtot),o_dtr_sat(nqtot),o_dtr_uscav(nqtot))
     174    ALLOCATE(o_dtr_wet_con(nqtot))
    174175    ALLOCATE(o_dtr_dry(nqtot),o_dtr_vdf(nqtot))
    175176IF (CPPKEY_STRATAER) THEN
     
    540541            tnam = TRIM(dn)//'uscav';       o_dtr_uscav     (itr) = ctrl_out(flag, tnam, lnam, "-", [('',i=1,nfiles)])
    541542
     543            lnam = 'tracer convective wet deposition'//TRIM(tracers(iq)%longName)
     544            tnam = TRIM(dn)//'wet_con';       o_dtr_wet_con       (itr) = ctrl_out(flag, tnam, lnam, "-", [('',i=1,nfiles)])
    542545            lnam = 'tracer tendency dry deposition'//TRIM(tracers(iq)%longName)
    543546            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

    r5384 r5447  
    66  USE phytrac_mod, ONLY : d_tr_cl, d_tr_th, d_tr_cv, d_tr_lessi_impa, &
    77       d_tr_lessi_nucl, d_tr_insc, d_tr_bcscav, d_tr_evapls, d_tr_ls,  &
    8        d_tr_trsp, d_tr_sscav, d_tr_sat, d_tr_uscav, flux_tr_dry
     8       d_tr_trsp, d_tr_sscav, d_tr_sat, d_tr_uscav, flux_tr_wet, flux_tr_dry
    99
    1010  ! Author: Abderrahmane IDELKADI (original include file)
     
    189189         o_dtr_insc, o_dtr_bcscav, o_dtr_evapls, &
    190190         o_dtr_ls, o_dtr_trsp, o_dtr_sscav, o_dtr_dry, &
    191          o_dtr_sat, o_dtr_uscav, o_trac_cum, o_du_gwd_rando, o_dv_gwd_rando, &
     191         o_dtr_sat, o_dtr_uscav, o_dtr_wet_con, &
     192         o_trac_cum, o_du_gwd_rando, o_dv_gwd_rando, &
    192193         o_ustr_gwd_hines,o_vstr_gwd_hines,o_ustr_gwd_rando,o_vstr_gwd_rando, &
    193194         o_ustr_gwd_front,o_vstr_gwd_front, &
     
    28562857             CALL histwrite_phy(o_dtr_uscav(itr),d_tr_uscav(:,:,itr))
    28572858            !--2D fields
     2859             CALL histwrite_phy(o_dtr_wet_con(itr), flux_tr_wet(:,itr))
    28582860             CALL histwrite_phy(o_dtr_dry(itr), flux_tr_dry(:,itr))
    28592861             zx_tmp_fi2d=0.
  • LMDZ6/trunk/libf/phylmd/phytrac_mod.f90

    r5330 r5447  
    3535  REAL,DIMENSION(:,:,:),ALLOCATABLE,SAVE :: d_tr_sat
    3636  REAL,DIMENSION(:,:,:),ALLOCATABLE,SAVE :: d_tr_uscav
     37  REAL,DIMENSION(:,:),ALLOCATABLE,SAVE   :: flux_tr_wet ! tracer wet deposit (surface)                    jyg
    3738  REAL,DIMENSION(:,:,:),ALLOCATABLE,SAVE :: qPr,qDi ! concentration tra dans pluie,air descente insaturee
    3839  REAL,DIMENSION(:,:,:),ALLOCATABLE,SAVE :: qPa,qMel
     
    4748
    4849!$OMP THREADPRIVATE(qPa,qMel,qTrdi,dtrcvMA,d_tr_th,d_tr_lessi_impa,d_tr_lessi_nucl)
    49 !$OMP THREADPRIVATE(d_tr_trsp,d_tr_sscav,d_tr_sat,d_tr_uscav,qPr,qDi)
     50!$OMP THREADPRIVATE(d_tr_trsp,d_tr_sscav,d_tr_sat,d_tr_uscav,flux_tr_wet,qPr,qDi)
    5051!$OMP THREADPRIVATE(d_tr_insc,d_tr_bcscav,d_tr_evapls,d_tr_ls,qPrls)
    5152!$OMP THREADPRIVATE(d_tr_cl,d_tr_dry,flux_tr_dry,d_tr_dec,d_tr_cv)
     
    6869    ALLOCATE(d_tr_sscav(klon,klev,nbtr),d_tr_sat(klon,klev,nbtr))
    6970    ALLOCATE(d_tr_uscav(klon,klev,nbtr),qPr(klon,klev,nbtr),qDi(klon,klev,nbtr))
     71    ALLOCATE(flux_tr_wet(klon,nbtr))
    7072    ALLOCATE(qPa(klon,klev,nbtr),qMel(klon,klev,nbtr))
    7173    ALLOCATE(qTrdi(klon,klev,nbtr),dtrcvMA(klon,klev,nbtr))
     
    408410          d_tr_dry(i,it)=0.
    409411          flux_tr_dry(i,it)=0.
     412          flux_tr_wet(i,it)=0.
    410413       ENDDO
    411414    ENDDO
     
    697700                !--with the full array tr_seri even if only item it is processed
    698701
    699                 CALL cvltr_scav(pdtphys, da, phi,phi2,d1a,dam, mp,ep,         &
    700                      sigd,sij,wght_cvfd,clw,elij,epmlmMm,eplaMm,              &     
    701                      pmflxr,pmflxs,evap,t_seri,wdtrainA,wdtrainM,             &   
    702                      paprs,it,tr_seri,upwd,dnwd,itop_con,ibas_con,            &
    703                      ccntrAA_3d,ccntrENV_3d,coefcoli_3d,                      &
    704                      d_tr_cv,d_tr_trsp,d_tr_sscav,d_tr_sat,d_tr_uscav,qDi,qPr,&
    705                      qPa,qMel,qTrdi,dtrcvMA,Mint,                             &
     702                CALL cvltr_scav(pdtphys, da, phi,phi2,d1a,dam, mp,ep,                &
     703                     sigd,sij,wght_cvfd,clw,elij,epmlmMm,eplaMm,                     &     
     704                     pmflxr,pmflxs,evap,t_seri,wdtrainA,wdtrainM,                    &   
     705                     paprs,it,tr_seri,upwd,dnwd,itop_con,ibas_con,                   &
     706                     ccntrAA_3d,ccntrENV_3d,coefcoli_3d,                             &
     707                     d_tr_cv,d_tr_trsp,d_tr_sscav,d_tr_sat,d_tr_uscav,flux_tr_wet,   &
     708                     qDi,qPr,                                                        &
     709                     qPa,qMel,qTrdi,dtrcvMA,Mint,                                    &
    706710                     zmfd1a,zmfphi2,zmfdam)
    707711
Note: See TracChangeset for help on using the changeset viewer.