Index: trunk/LMDZ.VENUS/libf/phyvenus/chemparam_mod.F90 =================================================================== --- trunk/LMDZ.VENUS/libf/phyvenus/chemparam_mod.F90 (revision 3754) +++ trunk/LMDZ.VENUS/libf/phyvenus/chemparam_mod.F90 (revision 3755) @@ -21,6 +21,6 @@ i_osso_cis, i_osso_trans, i_s2o2_cyc, & i_ocs, i_hso3, i_h2so4, i_s2, & - i_clso2, i_oscl, i_n2, i_he, i_n, i_no, & - i_no2, i_n2d, & + i_clso2, i_cl2so2, i_oscl, i_n2, i_he, & + i_n, i_no, i_no2, i_n2d, & i_co2plus, i_coplus, i_oplus, i_o2plus, & i_n2plus, i_hplus, i_h2oplus, i_nplus, & @@ -1542,4 +1542,5 @@ i_s2 = 0 i_clso2 = 0 + i_cl2so2 = 0 i_oscl = 0 i_n2 = 0 @@ -1744,4 +1745,9 @@ m_tr(i_clso2) = 99.517 type_tr(i_clso2) = 1 + case('cl2so2') + i_cl2so2 = i + print*,'cl2so2',i_cl2so2 + m_tr(i_cl2so2) = 134.97 + type_tr(i_cl2so2) = 1 case('oscl') i_oscl = i Index: trunk/LMDZ.VENUS/libf/phyvenus/photochemistry_venus.F90 =================================================================== --- trunk/LMDZ.VENUS/libf/phyvenus/photochemistry_venus.F90 (revision 3754) +++ trunk/LMDZ.VENUS/libf/phyvenus/photochemistry_venus.F90 (revision 3755) @@ -212,5 +212,5 @@ i_cl, i_clo, i_cl2, i_hocl, i_so2, i_so, i_so3, i_s2, & i_osso_cis, i_osso_trans, i_s2o2_cyc, i_clso2, & - i_ocs, i_cocl2, i_h2so4, & + i_cl2so2, i_ocs, i_cocl2, i_h2so4, & i_no2, i_no, i_n2, i_n2d, & nesp, tr, sza_input, dist_sol, v_phot) @@ -652,4 +652,12 @@ !=========================================================== +! Cl2SO2 + hv -> Cl + ClSO2 +!=========================================================== + +nb_phot = nb_phot + 1 + +indice_phot(nb_phot) = z3spec(1.0, i_cl2so2, 1.0, i_cl, 1.0, i_clso2) + +!=========================================================== ! OCS + hv -> CO + S !=========================================================== @@ -1429,10 +1437,10 @@ !=========================================================== -! f030 : ClSO2 + ClSO2 -> Cl2 + SO2 + SO2 +! f030 : ClSO2 + ClSO2 -> Cl2SO2 + SO2 !=========================================================== nb_reaction_3 = nb_reaction_3 + 1 -indice_3(nb_reaction_3) = z3spec(2.0, i_clso2, 1.0, i_cl2, 2.0, i_so2) +indice_3(nb_reaction_3) = z3spec(2.0, i_clso2, 1.0, i_cl2so2, 1.0, i_so2) !=========================================================== @@ -3540,9 +3548,11 @@ v_4(:,nb_reaction_4) = f029(:) -!--- f030: clso2 + clso2 -> cl2 + so2 + so2 +!--- f030: clso2 + clso2 -> cl2so2 + so2 ! moses et al. 2002 ! f030(:) = 5.0E-13 + +! croce and cobos, 2018 do iz = 1,nz Index: trunk/LMDZ.VENUS/libf/phyvenus/photolysis_mod.F90 =================================================================== --- trunk/LMDZ.VENUS/libf/phyvenus/photolysis_mod.F90 (revision 3754) +++ trunk/LMDZ.VENUS/libf/phyvenus/photolysis_mod.F90 (revision 3755) @@ -5,9 +5,9 @@ ! photolysis - integer, save :: nphot = 28 ! number of photolysis + integer, save :: nphot = 29 ! number of photolysis !$OMP THREADPRIVATE(nphot) - integer, parameter :: nabs = 25 ! number of absorbing gases + integer, parameter :: nabs = 26 ! number of absorbing gases ! spectral grid @@ -48,4 +48,5 @@ real, dimension(nw), save :: xss2o2_cyc ! s2o2_cyc absorption cross-section (cm2) real, dimension(nw), save :: xsclso2 ! clso2 absorption cross-section (cm2) + real, dimension(nw), save :: xscl2so2 ! cl2so2 absorption cross-section (cm2) real, dimension(nw), save :: xsocs ! cos absorption cross-section (cm2) real, dimension(nw), save :: xscocl2 ! cocl2 absorption cross-section (cm2) @@ -62,5 +63,5 @@ real, dimension(nw), save :: albedo ! surface albedo -!$OMP THREADPRIVATE(f, xsco2_195, xsco2_295, xsco2_370, yieldco2, xso2_150, xso2_200, xso2_250, xso2_300, yieldo2, xso3_218, xso3_298, xsh2o, xsh2o2, xsho2, xsh2, yieldh2, xsno2, xsno2_220, xsno2_294, yldno2_248, yldno2_298) +!$OMP THREADPRIVATE(f, xsco2_195, xsco2_295, xsco2_370, yieldco2, xso2_150, xso2_200, xso2_250, xso2_300, yieldo2, xso3_218, xso3_298, xsclso2, xscl2so2, xsh2o, xsh2o2, xsho2, xsh2, yieldh2, xsno2, xsno2_220, xsno2_294, yldno2_248, yldno2_298) !$OMP THREADPRIVATE(xsno, yieldno, xsn2, yieldn2, xshdo, albedo) @@ -148,4 +149,8 @@ call rdxsclso2(nw,wl,wc,xsclso2) + +! read and grid cl2so2 cross-sections + + call rdxscl2so2(nw,wl,wc,xscl2so2) ! read and grid ocs cross-sections @@ -2795,4 +2800,87 @@ !============================================================================== + subroutine rdxscl2so2(nw, wl, wc, yg) + +!-----------------------------------------------------------------------------* +!= PURPOSE: =* +!= Read cl2so2 cross-sections =* +!= From Uthman et al., J. Phys. Chem., 1978 =* +!-----------------------------------------------------------------------------* +!= PARAMETERS: =* +!= NW - INTEGER, number of specified intervals + 1 in working (I)=* +!= wavelength grid =* +!-----------------------------------------------------------------------------* + + USE mod_phys_lmdz_para, ONLY: is_master + USE mod_phys_lmdz_transfert_para, ONLY: bcast + + IMPLICIT NONE + +! input + + integer :: nw ! number of wavelength grid points + real, dimension(nw) :: wl, wc ! lower and central wavelength for each interval + +! output + + real, dimension(nw) :: yg ! clso2 cross-sections (cm2) + +! local + + real, parameter :: deltax = 1.e-4 + integer, parameter :: kdata = 1500 + + real, dimension(kdata) :: x1, y1 + real :: qy, lambda + integer :: i, iw, n, ierr + integer :: kin, kout ! input/output logical units + character*100 fil + + kin = 10 + +!*** cross sections from Uthman et al., J. Phys. Chem., 1978 + + fil = 'cross_sections/cl2so2_Uthman1978_295K_190_300nm.txt' + print*, 'section efficace Cl2SO2: ', fil + + if (is_master) then + + n = 56 + + OPEN(kin,FILE=fil,STATUS='OLD') + + do i = 1,7 + read(kin,*) + end do + do i = 1,n + read(kin,*) x1(i), y1(i) + end do + close(kin) + + call addpnt(x1,y1,kdata,n,x1(1)*(1.-deltax),0.) + call addpnt(x1,y1,kdata,n, 0.,0.) + call addpnt(x1,y1,kdata,n,x1(n)*(1.+deltax),0.) + call addpnt(x1,y1,kdata,n, 1e38,0.) + + call inter2(nw,wl,yg,n,x1,y1,ierr) + + if (ierr /= 0) then + write(*,*) ierr, fil + stop + end if + + do iw = 1,nw - 1 + lambda = wc(iw) +! write(55,'(f8.3,3e12.4)') lambda, yg(iw) + end do + + end if !is_master + + call bcast(yg) + + end subroutine rdxscl2so2 + +!============================================================================== + subroutine rdxsclo(nw, wl, yg) Index: trunk/LMDZ.VENUS/libf/phyvenus/photolysis_online.F =================================================================== --- trunk/LMDZ.VENUS/libf/phyvenus/photolysis_online.F (revision 3754) +++ trunk/LMDZ.VENUS/libf/phyvenus/photolysis_online.F (revision 3755) @@ -6,6 +6,6 @@ $ i_oh, i_ho2, i_h2o2, i_h2o, i_h, i_hcl, $ i_cl, i_clo, i_cl2, i_hocl, i_so2, i_so, i_so3, i_s2, - $ i_osso_cis, i_osso_trans, i_s2o2_cyc, i_clso2, - $ i_ocs, i_cocl2, i_h2so4, + $ i_osso_cis, i_osso_trans, i_s2o2_cyc, i_clso2, + $ i_cl2so2, i_ocs, i_cocl2, i_h2so4, $ i_no2, i_no, i_n2, i_n2d, & nesp, rm, sza, dist_sol, v_phot) @@ -23,7 +23,7 @@ $ i_oh, i_ho2, i_h2o2, i_h2o, i_h, i_hcl, $ i_cl, i_clo, i_cl2, i_hocl, i_so2, i_so, - $ i_so3, i_s2, - $ i_osso_cis, i_osso_trans, i_s2o2_cyc, - $ i_clso2, i_ocs, i_cocl2, i_h2so4, + $ i_so3, i_s2, i_osso_cis, i_osso_trans, + $ i_s2o2_cyc, + $ i_clso2, i_cl2so2, i_ocs, i_cocl2, i_h2so4, $ i_no2, i_no, i_n2, i_n2d @@ -72,10 +72,10 @@ $ j_h2o, j_h2o2, j_ho2, j_h, j_hcl, j_cl2, j_hocl, j_clo, $ j_so2, j_so, j_so3, j_s2, j_osso_cis, j_osso_trans, - $ j_s2o2_cyc, j_clso2, j_ocs, j_cocl2, j_h2so4, + $ j_s2o2_cyc, j_clso2, j_cl2so2, j_ocs, j_cocl2, j_h2so4, $ j_no2, j_no, j_n2, j_h2 integer :: a_o2, a_co2, a_o3, a_h2o, a_h2o2, a_ho2, a_hcl, a_cl2, $ a_hocl, a_clo, a_so2, a_so, a_so3, a_s2, a_osso_cis, - $ a_osso_trans, a_s2o2_cyc, a_clso2, a_ocs, + $ a_osso_trans, a_s2o2_cyc, a_clso2, a_cl2so2, a_ocs, $ a_cocl2, a_h2so4, a_no2, a_no, a_n2, a_h2 @@ -105,10 +105,11 @@ a_s2o2_cyc = 18 ! s2o2_cyc a_clso2 = 19 ! clso2 - a_ocs = 20 ! ocs - a_cocl2 = 21 ! cocl2 - a_h2so4 = 22 ! h2so4 - a_no2 = 23 ! no2 - a_no = 24 ! no - a_n2 = 25 ! n2 + a_cl2so2 = 20 ! cl2so2 + a_ocs = 21 ! ocs + a_cocl2 = 22 ! cocl2 + a_h2so4 = 23 ! h2so4 + a_no2 = 24 ! no2 + a_no = 25 ! no + a_n2 = 26 ! n2 ! photodissociation rates numbering. @@ -137,10 +138,11 @@ j_s2o2_cyc = 21 ! s2o2_cyc + hv -> so + so j_clso2 = 22 ! clso2 + hv -> cl + so2 - j_ocs = 23 ! ocs + hv -> co + s - j_cocl2 = 24 ! cocl2 + hv -> 2cl + co - j_h2so4 = 25 ! h2so4 + hv -> so3 + h2o - j_no2 = 26 ! no2 + hv -> no + o - j_no = 27 ! no + hv -> n + o - j_n2 = 28 ! n2 + hv -> n(2d) + n + j_cl2so2 = 23 ! cl2so2 + hv -> cl + clso2 + j_ocs = 24 ! ocs + hv -> co + s + j_cocl2 = 25 ! cocl2 + hv -> 2cl + co + j_h2so4 = 26 ! h2so4 + hv -> so3 + h2o + j_no2 = 27 ! no2 + hv -> no + o + j_no = 28 ! no + hv -> n + o + j_n2 = 29 ! n2 + hv -> n(2d) + n ! j_hdo_od = ! hdo + hv -> od + h @@ -240,4 +242,6 @@ dtgas(ilay,iw,a_clso2) = $ colinc(ilay)*rm(ilay,i_clso2)*xsclso2(iw) + dtgas(ilay,iw,a_cl2so2) = + $ colinc(ilay)*rm(ilay,i_cl2so2)*xscl2so2(iw) dtgas(ilay,iw,a_ocs) = colinc(ilay)*rm(ilay,i_ocs)*xsocs(iw) dtgas(ilay,iw,a_cocl2) = colinc(ilay)*rm(ilay,i_cocl2) @@ -279,4 +283,5 @@ sj(ilay,iw,j_s2o2_cyc) = xss2o2_cyc(iw) ! s2o2_cyc sj(ilay,iw,j_clso2) = xsclso2(iw) ! clso2 + sj(ilay,iw,j_cl2so2) = xscl2so2(iw) ! cl2so2 sj(ilay,iw,j_ocs) = xsocs(iw) ! ocs sj(ilay,iw,j_cocl2) = xscocl2(iw) ! cocl2