Changeset 2695 for LMDZ5/trunk/libf/phylmd/StratAer/so2_to_h2so4.F90
- Timestamp:
- Nov 1, 2016, 11:19:45 AM (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ5/trunk/libf/phylmd/StratAer/so2_to_h2so4.F90
r2690 r2695 1 SUBROUTINE SO2_TO_H2SO4(pdtphys,tr_seri,t_seri,pplay,paprs,sh, SO2_lifetime,is_strato,sulf_convert_sub)1 SUBROUTINE SO2_TO_H2SO4(pdtphys,tr_seri,t_seri,pplay,paprs,sh,is_strato) 2 2 3 3 USE dimphy, ONLY : klon,klev … … 5 5 USE infotrac 6 6 USE YOMCST, ONLY : RG 7 USE phys_local_var_mod, ONLY : SO2_lifetime, sulf_convert 7 8 8 9 IMPLICIT NONE … … 17 18 REAL,DIMENSION(klon,klev+1),INTENT(IN) :: paprs ! pression pour chaque inter-couche (en Pa) 18 19 REAL,DIMENSION(klon,klev),INTENT(IN) :: sh ! humidite specifique 19 REAL,DIMENSION(klon,klev),INTENT(IN) :: SO2_lifetime 20 LOGICAL,DIMENSION(klon,klev),INTENT(IN) :: is_strato 20 LOGICAL,DIMENSION(klon,klev),INTENT(IN) :: is_strato ! stratospheric flag 21 21 22 22 ! local variables in coagulation routine 23 23 INTEGER :: i,j,k,nb,ilon,ilev 24 REAL,DIMENSION(klon,klev) :: sulf_convert_sub ! SO2 converted to H2SO4 [kg(S)/m2/layer/s]25 24 26 25 !--convert SO2 to H2SO4 26 sulf_convert(:,:)=0.0 27 27 DO ilon=1, klon 28 28 DO ilev=1, klev … … 30 30 IF (is_strato(ilon,ilev)) THEN 31 31 IF (SO2_lifetime(ilon,ilev).GT.0.0) THEN 32 sulf_convert_sub(ilon,ilev)=tr_seri(ilon,ilev,id_SO2_strat)*(1-exp(-pdtphys/SO2_lifetime(ilon,ilev))) 33 ELSE 34 sulf_convert_sub(ilon,ilev)=0.0 32 sulf_convert(ilon,ilev)=tr_seri(ilon,ilev,id_SO2_strat)*(1.0-exp(-pdtphys/SO2_lifetime(ilon,ilev))) 35 33 ENDIF 36 tr_seri(ilon,ilev,id_SO2_strat)=tr_seri(ilon,ilev,id_SO2_strat) - sulf_convert _sub(ilon,ilev)37 tr_seri(ilon,ilev,id_H2SO4_strat)=tr_seri(ilon,ilev,id_H2SO4_strat) + mH2SO4mol/mSO2mol*sulf_convert _sub(ilon,ilev)34 tr_seri(ilon,ilev,id_SO2_strat)=tr_seri(ilon,ilev,id_SO2_strat) - sulf_convert(ilon,ilev) 35 tr_seri(ilon,ilev,id_H2SO4_strat)=tr_seri(ilon,ilev,id_H2SO4_strat) + mH2SO4mol/mSO2mol*sulf_convert(ilon,ilev) 38 36 !convert sulf_convert from kg(SO2)/kgA to kg(S)/m2/layer/s for saving as diagnostic 39 sulf_convert _sub(ilon,ilev)=sulf_convert_sub(ilon,ilev)*mSatom/mSO2mol*(paprs(ilon,ilev)-paprs(ilon,ilev+1))/RG/pdtphys37 sulf_convert(ilon,ilev)=sulf_convert(ilon,ilev)*mSatom/mSO2mol*(paprs(ilon,ilev)-paprs(ilon,ilev+1))/RG/pdtphys 40 38 ENDIF 41 39 ENDDO
Note: See TracChangeset
for help on using the changeset viewer.