Index: LMDZ5/branches/testing/libf/phylmd/StratAer/aer_sedimnt.F90 =================================================================== --- LMDZ5/branches/testing/libf/phylmd/StratAer/aer_sedimnt.F90 (revision 2720) +++ LMDZ5/branches/testing/libf/phylmd/StratAer/aer_sedimnt.F90 (revision 2787) @@ -14,5 +14,5 @@ !----------------------------------------------------------------------- - USE phys_local_var_mod, ONLY: mdw, sfluxaer, DENSO4, f_r_wet, vsed_aer + USE phys_local_var_mod, ONLY: mdw, budg_sed_part, DENSO4, f_r_wet, vsed_aer USE dimphy, ONLY : klon,klev USE infotrac @@ -106,10 +106,10 @@ !---ZAERONWM1 now contains the surface concentration at the new timestep !---PFLUXAER in unit of xx m-2 s-1 -sfluxaer(:)=0.0 +budg_sed_part(:)=0.0 DO JL=1,klon ZRHO=pplay(JL,1)/(RD*t_seri(JL,1)) DO nb=1,nbtr_bin - !compute sfluxaer as sum over bins in kg(S)/m2/s - sfluxaer(JL)=sfluxaer(JL)+ZRHO*ZAERONWM1(JL,nb)*ZVAER(JL,1,nb)*(mSatom/mH2SO4mol) & + !compute budg_sed_part as sum over bins in kg(S)/m2/s + budg_sed_part(JL)=budg_sed_part(JL)+ZRHO*ZAERONWM1(JL,nb)*ZVAER(JL,1,nb)*(mSatom/mH2SO4mol) & & *dens_aer_dry*4./3.*RPI*(mdw(nb)/2.)**3 ENDDO Index: LMDZ5/branches/testing/libf/phylmd/StratAer/aerophys.F90 =================================================================== --- LMDZ5/branches/testing/libf/phylmd/StratAer/aerophys.F90 (revision 2720) +++ LMDZ5/branches/testing/libf/phylmd/StratAer/aerophys.F90 (revision 2787) @@ -6,5 +6,6 @@ REAL,PARAMETER :: dens_aer_dry=1848.682308 ! dry aerosol particle mass density at T_0=293K[kg/m3] REAL,PARAMETER :: dens_aer_ref=1483.905336 ! aerosol particle mass density with 75% H2SO4 at T_0=293K[kg/m3] - REAL,PARAMETER :: mdwmin=0.002e-6 ! dry diameter of smallest aerosol particles [m] +! REAL,PARAMETER :: mdwmin=0.002e-6 ! dry diameter of smallest aerosol particles [m] + REAL,PARAMETER :: mdwmin=0.2e-6 ! dry diameter of smallest aerosol particles [m] !--testing REAL,PARAMETER :: V_rat=2.0 ! volume ratio of neighboring size bins REAL,PARAMETER :: mfrac_H2SO4=0.75 ! default mass fraction of H2SO4 in the aerosol Index: LMDZ5/branches/testing/libf/phylmd/StratAer/interp_sulf_input.F90 =================================================================== --- LMDZ5/branches/testing/libf/phylmd/StratAer/interp_sulf_input.F90 (revision 2720) +++ LMDZ5/branches/testing/libf/phylmd/StratAer/interp_sulf_input.F90 (revision 2787) @@ -8,5 +8,5 @@ USE mod_phys_lmdz_mpi_data, ONLY : is_mpi_root USE mod_phys_lmdz_omp_data, ONLY : is_omp_root - USE phys_local_var_mod, ONLY : OCS_backgr_tend, SO2_backgr_tend + USE phys_local_var_mod, ONLY : budg_3D_backgr_ocs, budg_3D_backgr_so2 USE phys_local_var_mod, ONLY : OCS_lifetime, SO2_lifetime USE mod_phys_lmdz_para @@ -34,4 +34,5 @@ REAL OCS_tmp, SO2_tmp INTEGER, SAVE :: mth_pre +!$OMP THREADPRIVATE(mth_pre) ! Champs reconstitues @@ -244,6 +245,6 @@ tr_seri(i,k,id_SO2_strat)=SO2_clim(i,k) ENDIF - OCS_backgr_tend(i,k)=tr_seri(i,k,id_OCS_strat)-OCS_tmp - SO2_backgr_tend(i,k)=tr_seri(i,k,id_SO2_strat)-SO2_tmp + budg_3D_backgr_ocs(i,k)=tr_seri(i,k,id_OCS_strat)-OCS_tmp + budg_3D_backgr_so2(i,k)=tr_seri(i,k,id_SO2_strat)-SO2_tmp ENDDO ENDDO @@ -252,6 +253,6 @@ DO i=1, klon DO k=1, klev - SO2_backgr_tend(i,k)=SO2_backgr_tend(i,k)*mSatom/mSO2mol*(paprs(i,k)-paprs(i,k+1))/RG/pdtphys - OCS_backgr_tend(i,k)=OCS_backgr_tend(i,k)*mSatom/mOCSmol*(paprs(i,k)-paprs(i,k+1))/RG/pdtphys + budg_3D_backgr_ocs(i,k)=budg_3D_backgr_ocs(i,k)*mSatom/mOCSmol*(paprs(i,k)-paprs(i,k+1))/RG/pdtphys + budg_3D_backgr_so2(i,k)=budg_3D_backgr_so2(i,k)*mSatom/mSO2mol*(paprs(i,k)-paprs(i,k+1))/RG/pdtphys ENDDO ENDDO Index: LMDZ5/branches/testing/libf/phylmd/StratAer/micphy_tstep.F90 =================================================================== --- LMDZ5/branches/testing/libf/phylmd/StratAer/micphy_tstep.F90 (revision 2720) +++ LMDZ5/branches/testing/libf/phylmd/StratAer/micphy_tstep.F90 (revision 2787) @@ -4,5 +4,5 @@ USE aerophys USE infotrac - USE phys_local_var_mod, ONLY: mdw, sulf_nucl, sulf_cond_evap, R2SO4, DENSO4, f_r_wet + USE phys_local_var_mod, ONLY: mdw, budg_3D_nucl, budg_3D_cond_evap, R2SO4, DENSO4, f_r_wet USE nucleation_tstep_mod USE cond_evap_tstep_mod @@ -68,6 +68,6 @@ IF (is_strato(ilon,ilev)) THEN ! initialize sulfur fluxes - sulf_nucl(ilon,ilev)=0.0 - sulf_cond_evap(ilon,ilev)=0.0 + budg_3D_nucl(ilon,ilev)=0.0 + budg_3D_cond_evap(ilon,ilev)=0.0 H2SO4_init=tr_seri(ilon,ilev,id_H2SO4_strat) ! adaptive timestep for nucleation and condensation @@ -109,7 +109,7 @@ CALL nucleation_part(nucl_rate,ntot,x,dt,Vbin,tr_seri(ilon,ilev,:)) ! compute fluxes as diagnostic in [kg(S)/m2/layer/s] (now - for evap and + for cond) - sulf_cond_evap(ilon,ilev)=sulf_cond_evap(ilon,ilev)+mSatom/mH2SO4mol & + budg_3D_cond_evap(ilon,ilev)=budg_3D_cond_evap(ilon,ilev)+mSatom/mH2SO4mol & & *cond_evap_rate*(paprs(ilon,ilev)-paprs(ilon,ilev+1))/RG*dt/pdtphys - sulf_nucl(ilon,ilev)=sulf_nucl(ilon,ilev)+mSatom/mH2SO4mol & + budg_3D_nucl(ilon,ilev)=budg_3D_nucl(ilon,ilev)+mSatom/mH2SO4mol & & *nucl_rate*(paprs(ilon,ilev)-paprs(ilon,ilev+1))/RG*dt/pdtphys ! update time step @@ -139,5 +139,5 @@ CALL cond_evap_part(pdtphys,FL,ASO4,f_r_wet(ilon,ilev),RRSI,Vbin,tr_seri(ilon,ilev,:)) ! compute fluxes as diagnostic in [kg(S)/m2/layer/s] (now - for evap and + for cond) - sulf_cond_evap(ilon,ilev)=sulf_cond_evap(ilon,ilev)+mSatom/mH2SO4mol & + budg_3D_cond_evap(ilon,ilev)=budg_3D_cond_evap(ilon,ilev)+mSatom/mH2SO4mol & & *evap_rate*(paprs(ilon,ilev)-paprs(ilon,ilev+1))/RG ENDIF Index: LMDZ5/branches/testing/libf/phylmd/StratAer/ocs_to_so2.F90 =================================================================== --- LMDZ5/branches/testing/libf/phylmd/StratAer/ocs_to_so2.F90 (revision 2720) +++ LMDZ5/branches/testing/libf/phylmd/StratAer/ocs_to_so2.F90 (revision 2787) @@ -1,3 +1,3 @@ -SUBROUTINE ocs_to_so2(pdtphys,tr_seri,t_seri,pplay,paprs,sh,is_strato) +SUBROUTINE ocs_to_so2(pdtphys,tr_seri,t_seri,pplay,paprs,is_strato) USE dimphy, ONLY : klon,klev @@ -5,5 +5,5 @@ USE infotrac USE YOMCST, ONLY : RG - USE phys_local_var_mod, ONLY : OCS_lifetime, ocs_convert + USE phys_local_var_mod, ONLY : OCS_lifetime, budg_3D_ocs_to_so2, budg_ocs_to_so2 IMPLICIT NONE @@ -16,5 +16,4 @@ REAL,DIMENSION(klon,klev),INTENT(IN) :: pplay ! pression pour le mileu de chaque couche (en Pa) REAL,DIMENSION(klon,klev+1),INTENT(IN) :: paprs ! pression pour chaque inter-couche (en Pa) - REAL,DIMENSION(klon,klev),INTENT(IN) :: sh ! humidite specifique LOGICAL,DIMENSION(klon,klev),INTENT(IN) :: is_strato @@ -23,5 +22,7 @@ !--convert OCS to SO2 - ocs_convert(:,:)=0.0 + budg_3D_ocs_to_so2(:,:)=0.0 + budg_ocs_to_so2(:)=0.0 + DO ilon=1, klon DO ilev=1, klev @@ -29,10 +30,11 @@ IF (is_strato(ilon,ilev)) THEN IF (OCS_lifetime(ilon,ilev).GT.0.0) THEN - ocs_convert(ilon,ilev)=tr_seri(ilon,ilev,id_OCS_strat)*(1.0-exp(-pdtphys/OCS_lifetime(ilon,ilev))) + budg_3D_ocs_to_so2(ilon,ilev)=tr_seri(ilon,ilev,id_OCS_strat)*(1.0-exp(-pdtphys/OCS_lifetime(ilon,ilev))) ENDIF - tr_seri(ilon,ilev,id_OCS_strat)=tr_seri(ilon,ilev,id_OCS_strat) - ocs_convert(ilon,ilev) - tr_seri(ilon,ilev,id_SO2_strat)=tr_seri(ilon,ilev,id_SO2_strat) + mSO2mol/mOCSmol*ocs_convert(ilon,ilev) - !convert ocs_convert from kg(OCS)/kgA to kg(S)/m2/layer/s for saving as diagnostic - ocs_convert(ilon,ilev)=ocs_convert(ilon,ilev)*mSatom/mOCSmol*(paprs(ilon,ilev)-paprs(ilon,ilev+1))/RG/pdtphys + tr_seri(ilon,ilev,id_OCS_strat)=tr_seri(ilon,ilev,id_OCS_strat) - budg_3D_ocs_to_so2(ilon,ilev) + tr_seri(ilon,ilev,id_SO2_strat)=tr_seri(ilon,ilev,id_SO2_strat) + mSO2mol/mOCSmol*budg_3D_ocs_to_so2(ilon,ilev) + !convert budget from kg(OCS)/kgA to kg(S)/m2/layer/s for saving as diagnostic + budg_3D_ocs_to_so2(ilon,ilev)=budg_3D_ocs_to_so2(ilon,ilev)*mSatom/mOCSmol*(paprs(ilon,ilev)-paprs(ilon,ilev+1))/RG/pdtphys + budg_ocs_to_so2(ilon)=budg_ocs_to_so2(ilon)+budg_3D_ocs_to_so2(ilon,ilev) ENDIF ENDDO Index: LMDZ5/branches/testing/libf/phylmd/StratAer/so2_to_h2so4.F90 =================================================================== --- LMDZ5/branches/testing/libf/phylmd/StratAer/so2_to_h2so4.F90 (revision 2720) +++ LMDZ5/branches/testing/libf/phylmd/StratAer/so2_to_h2so4.F90 (revision 2787) @@ -1,3 +1,3 @@ -SUBROUTINE SO2_TO_H2SO4(pdtphys,tr_seri,t_seri,pplay,paprs,sh,is_strato) +SUBROUTINE SO2_TO_H2SO4(pdtphys,tr_seri,t_seri,pplay,paprs,is_strato) USE dimphy, ONLY : klon,klev @@ -5,5 +5,5 @@ USE infotrac USE YOMCST, ONLY : RG - USE phys_local_var_mod, ONLY : SO2_lifetime, sulf_convert + USE phys_local_var_mod, ONLY : SO2_lifetime, budg_3D_so2_to_h2so4, budg_so2_to_h2so4 IMPLICIT NONE @@ -17,5 +17,4 @@ REAL,DIMENSION(klon,klev),INTENT(IN) :: pplay ! pression pour le mileu de chaque couche (en Pa) REAL,DIMENSION(klon,klev+1),INTENT(IN) :: paprs ! pression pour chaque inter-couche (en Pa) - REAL,DIMENSION(klon,klev),INTENT(IN) :: sh ! humidite specifique LOGICAL,DIMENSION(klon,klev),INTENT(IN) :: is_strato ! stratospheric flag @@ -24,5 +23,7 @@ !--convert SO2 to H2SO4 - sulf_convert(:,:)=0.0 + budg_3D_so2_to_h2so4(:,:)=0.0 + budg_so2_to_h2so4(:)=0.0 + DO ilon=1, klon DO ilev=1, klev @@ -30,10 +31,11 @@ IF (is_strato(ilon,ilev)) THEN IF (SO2_lifetime(ilon,ilev).GT.0.0) THEN - sulf_convert(ilon,ilev)=tr_seri(ilon,ilev,id_SO2_strat)*(1.0-exp(-pdtphys/SO2_lifetime(ilon,ilev))) + budg_3D_so2_to_h2so4(ilon,ilev)=tr_seri(ilon,ilev,id_SO2_strat)*(1.0-exp(-pdtphys/SO2_lifetime(ilon,ilev))) ENDIF - tr_seri(ilon,ilev,id_SO2_strat)=tr_seri(ilon,ilev,id_SO2_strat) - sulf_convert(ilon,ilev) - tr_seri(ilon,ilev,id_H2SO4_strat)=tr_seri(ilon,ilev,id_H2SO4_strat) + mH2SO4mol/mSO2mol*sulf_convert(ilon,ilev) - !convert sulf_convert from kg(SO2)/kgA to kg(S)/m2/layer/s for saving as diagnostic - sulf_convert(ilon,ilev)=sulf_convert(ilon,ilev)*mSatom/mSO2mol*(paprs(ilon,ilev)-paprs(ilon,ilev+1))/RG/pdtphys + tr_seri(ilon,ilev,id_SO2_strat)=tr_seri(ilon,ilev,id_SO2_strat) - budg_3D_so2_to_h2so4(ilon,ilev) + tr_seri(ilon,ilev,id_H2SO4_strat)=tr_seri(ilon,ilev,id_H2SO4_strat) + mH2SO4mol/mSO2mol*budg_3D_so2_to_h2so4(ilon,ilev) + !convert budget from kg(SO2)/kgA to kg(S)/m2/layer/s for saving as diagnostic + budg_3D_so2_to_h2so4(ilon,ilev)=budg_3D_so2_to_h2so4(ilon,ilev)*mSatom/mSO2mol*(paprs(ilon,ilev)-paprs(ilon,ilev+1))/RG/pdtphys + budg_so2_to_h2so4(ilon)=budg_so2_to_h2so4(ilon)+budg_3D_so2_to_h2so4(ilon,ilev) ENDIF ENDDO Index: LMDZ5/branches/testing/libf/phylmd/StratAer/traccoag_mod.F90 =================================================================== --- LMDZ5/branches/testing/libf/phylmd/StratAer/traccoag_mod.F90 (revision 2720) +++ LMDZ5/branches/testing/libf/phylmd/StratAer/traccoag_mod.F90 (revision 2787) @@ -8,10 +8,8 @@ SUBROUTINE traccoag(pdtphys, gmtime, debutphy, julien, & presnivs, xlat, xlon, pphis, pphi, & - t_seri, pplay, paprs, sh, rh , & - tr_seri) - - USE phys_local_var_mod, ONLY: mdw, sulf_convert, sulf_nucl, sulf_cond_evap, & - & sfluxaer, ocs_convert, R2SO4, DENSO4, f_r_wet, SO2_backgr_tend, OCS_backgr_tend, & - & OCS_lifetime, SO2_lifetime, surf_PM25_sulf + t_seri, pplay, paprs, sh, rh, tr_seri) + + USE phys_local_var_mod, ONLY: mdw, R2SO4, DENSO4, f_r_wet, surf_PM25_sulf, & + & budg_emi_ocs, budg_emi_so2, budg_emi_h2so4, budg_emi_part USE dimphy @@ -71,5 +69,5 @@ REAL,PARAMETER :: altemiss_sai=17.e3 ! emission altitude in m REAL,PARAMETER :: sigma_alt_sai=1.e3 ! standard deviation of emission altitude in m - REAL,PARAMETER :: xlat_sai=0.0 ! latitude of SAI in degree + REAL,PARAMETER :: xlat_sai=0.01 ! latitude of SAI in degree REAL,PARAMETER :: xlon_sai=120.35 ! longitude of SAI in degree @@ -90,6 +88,8 @@ REAL :: zrho ! Density of air [kg/m3] REAL :: zdz ! thickness of atm. model layer in m + REAL,DIMENSION(klev) :: zdm ! mass of atm. model layer in kg REAL,DIMENSION(klon,klev) :: dens_aer ! density of aerosol particles [kg/m3 aerosol] with default H2SO4 mass fraction REAL :: dlat, dlon ! d latitude and d longitude of grid in degree + REAL :: emission ! emission IF (is_mpi_root) THEN @@ -141,17 +141,23 @@ DO ilon=1, klon DO ilev=1, klev - m_air_gridbox(ilon,ilev)=(paprs(ilon,ilev)-paprs(ilon,ilev+1)) / RG * cell_area(ilon) - ENDDO - ENDDO - - IF (debutphy) THEN - CALL gather(tr_seri, tr_seri_glo) - IF (MAXVAL(tr_seri_glo).LT.1.e-30) THEN + m_air_gridbox(ilon,ilev)=(paprs(ilon,ilev)-paprs(ilon,ilev+1))/RG*cell_area(ilon) + ENDDO + ENDDO + +! IF (debutphy) THEN +! CALL gather(tr_seri, tr_seri_glo) +! IF (MAXVAL(tr_seri_glo).LT.1.e-30) THEN !--initialising tracer concentrations to zero - DO it=1, nbtr - tr_seri(:,:,it)=0.0 - ENDDO - ENDIF - ENDIF +! DO it=1, nbtr +! tr_seri(:,:,it)=0.0 +! ENDDO +! ENDIF +! ENDIF + +!--initialise emission diagnostics + budg_emi_ocs(:)=0.0 + budg_emi_so2(:)=0.0 + budg_emi_h2so4(:)=0.0 + budg_emi_part(:)=0.0 !--sulfur emission, depending on chosen scenario (flag_sulf_emit) @@ -170,10 +176,13 @@ IF ( xlat(i).GE.xlat_vol-dlat .AND. xlat(i).LT.xlat_vol+dlat .AND. & xlon(i).GE.xlon_vol-dlon .AND. xlon(i).LT.xlon_vol+dlon ) THEN +! + PRINT *,'coordinates of volcanic injection point=',xlat(i), xlon(i), day_cur, mth_cur, year_cur ! compute altLMDz altLMDz(:)=0.0 DO k=1, klev - zrho=pplay(i,k)/t_seri(i,k)/RD !air density in kg/m3 - zdz=(paprs(i,k)-paprs(i,k+1))/zrho/RG !thickness of layer in m - altLMDz(k+1)=altLMDz(k)+zdz + zrho=pplay(i,k)/t_seri(i,k)/RD !air density in kg/m3 + zdm(k)=(paprs(i,k)-paprs(i,k+1))/RG !mass of layer in kg + zdz=zdm(k)/zrho !thickness of layer in m + altLMDz(k+1)=altLMDz(k)+zdz !altitude of interface ENDDO !compute distribution of emission to vertical model layers (based on Gaussian peak in altitude) @@ -194,7 +203,8 @@ !vertically distributed emission DO k=1, klev - tr_seri(i,k,id_SO2_strat)=tr_seri(i,k,id_SO2_strat)+ & - & m_aer_emiss_vol*(mSO2mol/mSatom)/m_air_gridbox(i,k)*f_lay_emiss(k) & - & /(1.*86400./pdtphys) ! stretch emission over one day of Pinatubo eruption + ! stretch emission over one day of Pinatubo eruption + emission=m_aer_emiss_vol*(mSO2mol/mSatom)/m_air_gridbox(i,k)*f_lay_emiss(k)/1./86400. + tr_seri(i,k,id_SO2_strat)=tr_seri(i,k,id_SO2_strat)+emission*pdtphys + budg_emi_so2(i)=budg_emi_so2(i)+emission*zdm(k)*mSatom/mSO2mol ENDDO ENDIF ! emission grid cell @@ -211,10 +221,13 @@ IF ( xlat(i).GE.xlat_sai-dlat .AND. xlat(i).LT.xlat_sai+dlat .AND. & & xlon(i).GE.xlon_sai-dlon .AND. xlon(i).LT.xlon_sai+dlon ) THEN +! + PRINT *,'coordinates of SAI point=',xlat(i), xlon(i), day_cur, mth_cur, year_cur ! compute altLMDz altLMDz(:)=0.0 DO k=1, klev - zrho=pplay(i,k)/t_seri(i,k)/RD !air density in kg/m3 - zdz=(paprs(i,k)-paprs(i,k+1))/zrho/RG !thickness of layer in m - altLMDz(k+1)=altLMDz(k)+zdz + zrho=pplay(i,k)/t_seri(i,k)/RD !air density in kg/m3 + zdm(k)=(paprs(i,k)-paprs(i,k+1))/RG !mass of layer in kg + zdz=zdm(k)/zrho !thickness of layer in m + altLMDz(k+1)=altLMDz(k)+zdz !altitude of interface ENDDO !compute distribution of emission to vertical model layers (based on Gaussian peak in altitude) @@ -235,17 +248,16 @@ !vertically distributed emission DO k=1, klev - tr_seri(i,k,id_SO2_strat)=tr_seri(i,k,id_SO2_strat)+ & - & m_aer_emiss_sai*(mSO2mol/mSatom)/m_air_gridbox(i,k)*f_lay_emiss(k) & - & /(360.*86400./pdtphys) ! stretch emission over whole year (360d) -! & /(60.*86400./pdtphys) ! stretch emission over 2 months (seasonal emission) -! & /7. ! distribute equally over 7 emission grid points + ! stretch emission over whole year (360d) + emission=m_aer_emiss_sai*(mSO2mol/mSatom)/m_air_gridbox(i,k)*f_lay_emiss(k)/360./86400. + tr_seri(i,k,id_SO2_strat)=tr_seri(i,k,id_SO2_strat)+emission*pdtphys + budg_emi_so2(i)=budg_emi_so2(i)+emission*zdm(k)*mSatom/mSO2mol ENDDO ! !emission as monodisperse particles with 0.1um dry radius (BIN21) ! !vertically distributed emission ! DO k=1, klev -! tr_seri(i,k,id_BIN01_strat+20)=tr_seri(i,k,id_BIN01_strat+20)+ & -! & m_aer_emiss*(mH2SO4mol/mSatom)/m_part_dry(21)/m_air_gridbox(i,k)*f_lay_emiss(k) & -! & /(360.*86400./pdtphys) & ! stretch emission over whole year (360d) -! & /7. ! distribute equally over 7 emission grid points +! ! stretch emission over whole year (360d) +! emission=m_aer_emiss*(mH2SO4mol/mSatom)/m_part_dry(21)/m_air_gridbox(i,k)*f_lay_emiss(k)/360./86400 +! tr_seri(i,k,id_BIN01_strat+20)=tr_seri(i,k,id_BIN01_strat+20)+emission*pdtphys +! budg_emi_part(i)=budg_emi_part(i)+emission*zdm(k)*mSatom/mH2SO4mol ! ENDDO ENDIF ! emission grid cell @@ -259,8 +271,8 @@ !--convert OCS to SO2 in the stratosphere - CALL ocs_to_so2(pdtphys,tr_seri,t_seri,pplay,paprs,sh,is_strato) + CALL ocs_to_so2(pdtphys,tr_seri,t_seri,pplay,paprs,is_strato) !--convert SO2 to H2SO4 - CALL so2_to_h2so4(pdtphys,tr_seri,t_seri,pplay,paprs,sh,is_strato) + CALL so2_to_h2so4(pdtphys,tr_seri,t_seri,pplay,paprs,is_strato) !--common routine for nucleation and condensation/evaporation with adaptive timestep @@ -282,5 +294,5 @@ surf_PM25_sulf(i)=surf_PM25_sulf(i)+tr_seri(i,1,it+nbtr_sulgas) & & *132./98.*dens_aer_dry*4./3.*RPI*(mdw(it)/2.)**3 & - & *pplay(i,1)/t_seri(i,1)/RD*1e9 + & *pplay(i,1)/t_seri(i,1)/RD*1.e9 ENDIF ENDDO