Changeset 3411 for LMDZ6/branches/DYNAMICO-conv/libf/phylmd/phytrac_mod.F90
- Timestamp:
- Nov 5, 2018, 3:24:59 PM (6 years ago)
- Location:
- LMDZ6/branches/DYNAMICO-conv
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/DYNAMICO-conv
- Property svn:mergeinfo changed
/LMDZ6/trunk removed
- Property svn:mergeinfo changed
-
LMDZ6/branches/DYNAMICO-conv/libf/phylmd/phytrac_mod.F90
r3356 r3411 49 49 !$OMP THREADPRIVATE(d_tr_trsp,d_tr_sscav,d_tr_sat,d_tr_uscav,qPr,qDi) 50 50 !$OMP THREADPRIVATE(d_tr_insc,d_tr_bcscav,d_tr_evapls,d_tr_ls,qPrls) 51 !$OMP THREADPRIVATE(d_tr _cl,d_tr_dry,flux_tr_dry,d_tr_dec,d_tr_cv)51 !$OMP THREADPRIVATE(d_tr,d_tr_cl,d_tr_dry,flux_tr_dry,d_tr_dec,d_tr_cv) 52 52 53 53 … … 104 104 USE traccoag_mod 105 105 USE phys_local_var_mod, ONLY: mdw 106 USE phys_local_var_mod, ONLY: budg_dep_dry_ocs, budg_dep_wet_ocs107 USE phys_local_var_mod, ONLY: budg_dep_dry_so2, budg_dep_wet_so2108 106 USE phys_local_var_mod, ONLY: budg_dep_dry_h2so4, budg_dep_wet_h2so4 109 107 USE phys_local_var_mod, ONLY: budg_dep_dry_part, budg_dep_wet_part 110 USE infotrac, ONLY: nbtr_sulgas, id_ OCS_strat, id_SO2_strat, id_H2SO4_strat108 USE infotrac, ONLY: nbtr_sulgas, id_SO2_strat, id_H2SO4_strat 111 109 USE aerophys 112 110 #endif … … 509 507 ENDIF 510 508 511 CASE('repr')512 flag_cvltr(it)=.false.513 514 509 CASE('inca') 515 510 ! IF ((it.EQ.id_Rn222) .OR. ((it.GE.id_SO2) .AND. (it.LE.id_NH3)) ) THEN … … 706 701 IF (type_trac=='coag') THEN 707 702 ! initialize wet deposition flux of sulfur 708 budg_dep_wet_ocs(:)=0.0709 budg_dep_wet_so2(:)=0.0710 703 budg_dep_wet_h2so4(:)=0.0 711 704 budg_dep_wet_part(:)=0.0 … … 716 709 DO it = 1, nbtr 717 710 !do not include SO2 because most of it comes trom the troposphere 718 IF (it==id_OCS_strat) THEN 719 budg_dep_wet_ocs(i)=budg_dep_wet_ocs(i)+d_tr_cv(i,k,it)*(mSatom/mOCSmol) & 720 & *(paprs(i,k)-paprs(i,k+1))/RG/pdtphys 721 ELSEIF (it==id_SO2_strat) THEN 722 budg_dep_wet_so2(i)=budg_dep_wet_so2(i)+d_tr_cv(i,k,it)*(mSatom/mSO2mol) & 723 & *(paprs(i,k)-paprs(i,k+1))/RG/pdtphys 724 ELSEIF (it==id_H2SO4_strat) THEN 711 IF (it==id_H2SO4_strat) THEN 725 712 budg_dep_wet_h2so4(i)=budg_dep_wet_h2so4(i)+d_tr_cv(i,k,it)*(mSatom/mH2SO4mol) & 726 713 & *(paprs(i,k)-paprs(i,k+1))/RG/pdtphys … … 747 734 d_tr_th(i,k,it)=0. 748 735 tr_seri(i,k,it)=MAX(tr_seri(i,k,it),0.) 749 ! the next safeguard causes some problem for stratospheric aerosol tracers (particle number) 750 ! and there is little justification for it so it is commented out (4 December 2017) by OB 751 ! if reinstated please keep the ifndef CPP_StratAer 752 !#ifndef CPP_StratAer 753 ! tr_seri(i,k,it)=MIN(tr_seri(i,k,it),1.e10) 754 !#endif 736 tr_seri(i,k,it)=MIN(tr_seri(i,k,it),1.e10) 755 737 END DO 756 738 END DO … … 788 770 789 771 ! initialize dry deposition flux of sulfur 790 budg_dep_dry_ocs(:)=0.0791 budg_dep_dry_so2(:)=0.0792 772 budg_dep_dry_h2so4(:)=0.0 793 773 budg_dep_dry_part(:)=0.0 … … 824 804 IF (type_trac=='coag') THEN 825 805 ! compute dry deposition flux of sulfur (sum over gases and particles) 826 IF (it==id_OCS_strat) THEN 827 budg_dep_dry_ocs(:)=budg_dep_dry_ocs(:)-source(:,it)*(mSatom/mOCSmol) 828 ELSEIF (it==id_SO2_strat) THEN 829 budg_dep_dry_so2(:)=budg_dep_dry_so2(:)-source(:,it)*(mSatom/mSO2mol) 830 ELSEIF (it==id_H2SO4_strat) THEN 806 IF (it==id_H2SO4_strat) THEN 831 807 budg_dep_dry_h2so4(:)=budg_dep_dry_h2so4(:)-source(:,it)*(mSatom/mH2SO4mol) 832 808 ELSEIF (it.GT.nbtr_sulgas) THEN … … 914 890 DO k = 1, klev 915 891 DO it = 1, nbtr 916 IF (it==id_OCS_strat) THEN 917 budg_dep_wet_ocs(i)=budg_dep_wet_ocs(i)+d_tr_ls(i,k,it)*(mSatom/mOCSmol) & 918 & *(paprs(i,k)-paprs(i,k+1))/RG/pdtphys 919 ELSEIF (it==id_SO2_strat) THEN 920 budg_dep_wet_so2(i)=budg_dep_wet_so2(i)+d_tr_ls(i,k,it)*(mSatom/mSO2mol) & 921 & *(paprs(i,k)-paprs(i,k+1))/RG/pdtphys 922 ELSEIF (it==id_H2SO4_strat) THEN 892 IF (it==id_H2SO4_strat) THEN 923 893 budg_dep_wet_h2so4(i)=budg_dep_wet_h2so4(i)+d_tr_ls(i,k,it)*(mSatom/mH2SO4mol) & 924 894 & *(paprs(i,k)-paprs(i,k+1))/RG/pdtphys
Note: See TracChangeset
for help on using the changeset viewer.