Ignore:
Timestamp:
Aug 4, 2025, 3:03:07 PM (5 days ago)
Author:
aborella
Message:

Additional diags for contrails + simplified coupling between deep conv and cirrus clouds + small modifsin RRTM for RF of contrails alone

File:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ6/branches/contrails/libf/phylmd/phys_local_var_mod.F90

    r5794 r5796  
    722722      REAL, SAVE, ALLOCATABLE :: fsurv_contrails(:,:), section_contrails(:,:)
    723723      !$OMP THREADPRIVATE(fsurv_contrails, section_contrails)
     724      REAL, SAVE, ALLOCATABLE :: nice_ygcont(:,:), iwc_ygcont(:,:)
     725      !$OMP THREADPRIVATE(nice_ygcont, iwc_ygcont)
     726      REAL, SAVE, ALLOCATABLE :: rvol_ygcont(:,:), tau_ygcont(:,:)
     727      !$OMP THREADPRIVATE(rvol_ygcont, tau_ygcont)
     728      REAL, SAVE, ALLOCATABLE :: nice_cont(:,:), iwc_cont(:,:)
     729      !$OMP THREADPRIVATE(nice_cont, iwc_cont)
     730      REAL, SAVE, ALLOCATABLE :: rvol_cont(:,:), tau_cont(:,:)
     731      !$OMP THREADPRIVATE(rvol_cont, tau_cont)
    724732      REAL, SAVE, ALLOCATABLE :: qice_cont(:,:)
    725733      !$OMP THREADPRIVATE(qice_cont)
     
    738746      REAL, SAVE, ALLOCATABLE :: dnic_agg(:,:)
    739747      !$OMP THREADPRIVATE(dnic_agg)
    740       REAL, SAVE, ALLOCATABLE :: cldfra_nocont(:,:), cldtau_nocont(:,:), cldemi_nocont(:,:)
    741       !$OMP THREADPRIVATE(cldfra_nocont, cldtau_nocont, cldemi_nocont)
     748      REAL, SAVE, ALLOCATABLE :: cldfra_cont(:,:), cldtau_nocont(:,:), cldemi_nocont(:,:)
     749      !$OMP THREADPRIVATE(cldfra_cont, cldtau_nocont, cldemi_nocont)
    742750      REAL, SAVE, ALLOCATABLE :: cldh_nocont(:), contcov(:), conttau(:,:), contemi(:,:)
    743751      !$OMP THREADPRIVATE(cldh_nocont, contcov, conttau, contemi)
    744       REAL, SAVE, ALLOCATABLE :: fiwp_nocont(:), fiwc_nocont(:,:), ref_ice_nocont(:,:)
    745       !$OMP THREADPRIVATE(fiwp_nocont, fiwc_nocont, ref_ice_nocont)
     752      REAL, SAVE, ALLOCATABLE :: fiwp_cont(:), fiwc_cont(:,:), ref_ice_cont(:,:)
     753      !$OMP THREADPRIVATE(fiwp_cont, fiwc_cont, ref_ice_cont)
    746754      REAL, SAVE, ALLOCATABLE :: topsw_nocont(:), toplw_nocont(:)
    747755      !$OMP THREADPRIVATE(topsw_nocont, toplw_nocont)
     
    13331341      ALLOCATE(AEI_contrails(klon,klev), AEI_surv_contrails(klon,klev))
    13341342      ALLOCATE(fsurv_contrails(klon,klev), section_contrails(klon,klev))
     1343      ALLOCATE(nice_ygcont(klon,klev), iwc_ygcont(klon,klev))
     1344      ALLOCATE(rvol_ygcont(klon,klev), tau_ygcont(klon,klev))
     1345      ALLOCATE(nice_cont(klon,klev), iwc_cont(klon,klev))
     1346      ALLOCATE(rvol_cont(klon,klev), tau_cont(klon,klev))
    13351347      ALLOCATE(qice_cont(klon,klev))
    13361348      ALLOCATE(contfra(klon,klev), qradice_cont(klon,klev))
     
    13411353      ALLOCATE(dcfc_sed(klon,klev), dqic_sed(klon,klev), dqtc_sed(klon,klev), dnic_sed(klon,klev))
    13421354      ALLOCATE(dcfc_auto(klon,klev), dqic_auto(klon,klev), dqtc_auto(klon,klev), dnic_auto(klon,klev))
    1343       ALLOCATE(cldfra_nocont(klon,klev), cldtau_nocont(klon,klev), cldemi_nocont(klon,klev))
     1355      ALLOCATE(cldfra_cont(klon,klev), cldtau_nocont(klon,klev), cldemi_nocont(klon,klev))
    13441356      ALLOCATE(cldh_nocont(klon), contcov(klon), conttau(klon,klev), contemi(klon,klev))
    1345       ALLOCATE(fiwp_nocont(klon), fiwc_nocont(klon,klev), ref_ice_nocont(klon,klev))
     1357      ALLOCATE(fiwp_cont(klon), fiwc_cont(klon,klev), ref_ice_cont(klon,klev))
    13461358      ALLOCATE(topsw_nocont(klon), toplw_nocont(klon))
    13471359      ALLOCATE(solsw_nocont(klon), sollw_nocont(klon))
     
    17791791      DEALLOCATE(Tcritcont, qcritcont, potcontfraP, potcontfraNP)
    17801792      DEALLOCATE(AEI_contrails, AEI_surv_contrails, fsurv_contrails, section_contrails)
     1793      DEALLOCATE(nice_ygcont, iwc_ygcont, rvol_ygcont, tau_ygcont)
     1794      DEALLOCATE(nice_cont, iwc_cont, rvol_cont, tau_cont)
    17811795      DEALLOCATE(qice_cont, contfra, qradice_cont)
    17821796      DEALLOCATE(dcfc_ini, dqic_ini, dqtc_ini, dnic_ini)
     
    17861800      DEALLOCATE(dcfc_sed, dqic_sed, dqtc_sed, dnic_sed)
    17871801      DEALLOCATE(dcfc_auto, dqic_auto, dqtc_auto, dnic_auto)
    1788       DEALLOCATE(cldfra_nocont, cldtau_nocont, cldemi_nocont, conttau, contemi)
    1789       DEALLOCATE(cldh_nocont, contcov, fiwp_nocont, fiwc_nocont, ref_ice_nocont)
     1802      DEALLOCATE(cldfra_cont, cldtau_nocont, cldemi_nocont, conttau, contemi)
     1803      DEALLOCATE(cldh_nocont, contcov, fiwp_cont, fiwc_cont, ref_ice_cont)
    17901804      DEALLOCATE(topsw_nocont, toplw_nocont)
    17911805      DEALLOCATE(solsw_nocont, sollw_nocont)
Note: See TracChangeset for help on using the changeset viewer.