[2690] | 1 | SUBROUTINE OCS_TO_SO2(pdtphys,tr_seri,t_seri,pplay,paprs,sh,OCS_lifetime,is_strato,ocs_convert_sub) |
---|
| 2 | |
---|
| 3 | USE dimphy, ONLY : klon,klev |
---|
| 4 | USE aerophys |
---|
| 5 | USE infotrac |
---|
| 6 | USE YOMCST, ONLY : RG |
---|
| 7 | |
---|
| 8 | IMPLICIT NONE |
---|
| 9 | |
---|
| 10 | !-------------------------------------------------------- |
---|
| 11 | ! transfer variables when calling this routine |
---|
| 12 | REAL,INTENT(IN) :: pdtphys ! Pas d'integration pour la physique (seconde) |
---|
| 13 | REAL,DIMENSION(klon,klev,nbtr),INTENT(INOUT) :: tr_seri ! Concentration Traceur [U/KgA] |
---|
| 14 | REAL,DIMENSION(klon,klev),INTENT(IN) :: t_seri ! Temperature |
---|
| 15 | REAL,DIMENSION(klon,klev),INTENT(IN) :: pplay ! pression pour le mileu de chaque couche (en Pa) |
---|
| 16 | REAL,DIMENSION(klon,klev+1),INTENT(IN) :: paprs ! pression pour chaque inter-couche (en Pa) |
---|
| 17 | REAL,DIMENSION(klon,klev),INTENT(IN) :: sh ! humidite specifique |
---|
| 18 | REAL,DIMENSION(klon,klev),INTENT(IN) :: OCS_lifetime |
---|
| 19 | LOGICAL,DIMENSION(klon,klev),INTENT(IN) :: is_strato |
---|
| 20 | |
---|
| 21 | ! local variables |
---|
| 22 | INTEGER :: i,j,k,nb,ilon,ilev |
---|
| 23 | REAL,DIMENSION(klon,klev) :: ocs_convert_sub ! OCS converted to SO2 [kg(S)/m2/layer/s] |
---|
| 24 | |
---|
| 25 | !--convert OCS to SO2 |
---|
| 26 | DO ilon=1, klon |
---|
| 27 | DO ilev=1, klev |
---|
| 28 | !only in the stratosphere |
---|
| 29 | IF (is_strato(ilon,ilev)) THEN |
---|
| 30 | IF (OCS_lifetime(ilon,ilev).GT.0.0) THEN |
---|
| 31 | ocs_convert_sub(ilon,ilev)=tr_seri(ilon,ilev,id_OCS_strat)*(1-exp(-pdtphys/OCS_lifetime(ilon,ilev))) |
---|
| 32 | ELSE |
---|
| 33 | ocs_convert_sub(ilon,ilev)=0.0 |
---|
| 34 | ENDIF |
---|
| 35 | tr_seri(ilon,ilev,id_OCS_strat)=tr_seri(ilon,ilev,id_OCS_strat) - ocs_convert_sub(ilon,ilev) |
---|
| 36 | tr_seri(ilon,ilev,id_SO2_strat)=tr_seri(ilon,ilev,id_SO2_strat) + mSO2mol/mOCSmol*ocs_convert_sub(ilon,ilev) |
---|
| 37 | !convert ocs_convert from kg(OCS)/kgA to kg(S)/m2/layer/s for saving as diagnostic |
---|
| 38 | ocs_convert_sub(ilon,ilev)=ocs_convert_sub(ilon,ilev)*mSatom/mOCSmol* & |
---|
| 39 | & (paprs(ilon,ilev)-paprs(ilon,ilev+1))/RG/pdtphys |
---|
| 40 | ENDIF |
---|
| 41 | ENDDO |
---|
| 42 | ENDDO |
---|
| 43 | |
---|
| 44 | END SUBROUTINE OCS_TO_SO2 |
---|