Index: /trunk/LMDZ.VENUS/libf/phyvenus/photochemistry_venus.F90
===================================================================
--- /trunk/LMDZ.VENUS/libf/phyvenus/photochemistry_venus.F90	(revision 2319)
+++ /trunk/LMDZ.VENUS/libf/phyvenus/photochemistry_venus.F90	(revision 2320)
@@ -1,3 +1,3 @@
-subroutine photochemistry_venus(nz, n_lon, ptimestep, p, t, tr, mumean, sza_input, nesp, iter)
+subroutine photochemistry_venus(nz, n_lon, ptimestep, p, t, tr, mumean, sza_input, lon, lat, nesp, iter)
 
 use chemparam_mod
@@ -17,4 +17,5 @@
 real :: ptimestep                 ! physics timestep (s)
 real :: sza_input                 ! solar zenith angle (degrees)
+real :: lon, lat                  ! longitude and latitude (degrees)
 
 integer :: n_lon                  ! for 1D test
@@ -176,5 +177,5 @@
 
    call define_dt(nesp, dt_corrected, dt_guess, ctimestep, cold(:), c(iz,:), &
-                  mat1, prod, loss, conc(iz))
+                  mat1, prod, loss, conc(iz), lon, lat)
 
    if (time + dt_corrected > ptimestep) then
@@ -2649,17 +2650,19 @@
 !      v_phot(:,nb_phot) = g023(:)*conc(:)
 
-!     Changement pour g023 ==> s2 + hv -> s + s
-!     Pas encore inclu dans la table jphot
-
-!     Pas de photodissociation sous le nuage photochimique 
-      g023(1:28)  = 0.0E+0
-
-!     Dependance en sza pour la photodissociation
-!     moins de W.m-2      
-      IF (sza_input.LT.90.0) THEN
-        g023(29:50) = 6.5E-3*COS(sza_input*pi/180.0)
-      ELSE
-        g023(29:50) = 0.0E+0
-      END IF
+!---  temporary treatment of s2 photolysis
+
+!     s2 + hv -> s + s
+
+!     under the clouds: no photolysis
+
+      g023(1:28)  = 0.
+
+!     above the clouds: value derived from Mills, 1998
+
+      if (sza_input < 90.) then
+        g023(29:nz) = 6.5e-3*cos(sza_input*pi/180.)
+      else
+        g023(29:nz) = 0.
+      end if
            
       nb_phot = nb_phot + 1
@@ -2980,5 +2983,5 @@
 
  subroutine define_dt(nesp, dtnew, dtold, ctimestep, cold, ccur, mat1, &
-                      prod, loss, dens)
+                      prod, loss, dens, lon, lat)
 
 !================================================================
@@ -3000,4 +3003,5 @@
 real, dimension(nesp)      :: prod, loss
 real                       :: dens
+real                       :: lon, lat
 
 ! output
Index: /trunk/LMDZ.VENUS/libf/phyvenus/phytrac_chimie.F
===================================================================
--- /trunk/LMDZ.VENUS/libf/phyvenus/phytrac_chimie.F	(revision 2319)
+++ /trunk/LMDZ.VENUS/libf/phyvenus/phytrac_chimie.F	(revision 2320)
@@ -83,5 +83,5 @@
                trac(:,1:22,i_co)  = 25.e-6
                trac(:,:,i_hcl)    = 0.4e-6
-               trac(:,1:22,i_so2) = 10.e-6
+               trac(:,1:22,i_so2) = 7.e-6
                trac(:,1:22,i_h2o) = 30.e-6
                trac(:,:,i_n2)     = 0.35e-1
@@ -257,5 +257,7 @@
      $                                ztrac(ilon,:,:),
      $                                mmean(ilon,:),
-     $                                sza_local, nqmax, iter(ilon,:))
+     $                                sza_local, 
+     $                                lon(ilon), lat(ilon),
+     $                                nqmax, iter(ilon,:))
 
          end do