Index: trunk/LMDZ.VENUS/libf/phyvenus/new_photochemistry_venus.F90
===================================================================
--- trunk/LMDZ.VENUS/libf/phyvenus/new_photochemistry_venus.F90	(revision 2135)
+++ trunk/LMDZ.VENUS/libf/phyvenus/new_photochemistry_venus.F90	(revision 2186)
@@ -29,4 +29,5 @@
 
 integer :: ij, i, iesp, jesp, iz, iiz
+real :: fact_SF        !  Increasin/Decreasing Solar Flux
 
 ! matrix
@@ -436,4 +437,13 @@
          read(30,'(7e11.4)') (jphot(iso2,isza,iz,ij), ij = 1,nj)
          do ij = 1,nj
+            !!!    
+            ! GG TESTS increasing/decreasing Solar Flux 
+            !fact_SF = 5.
+            ! jphot(iso2,isza,iz,ij)=jphot(iso2,isza,iz,ij)*0.2
+
+            !   print*,'!!! WARNING !!!: You are Modifying Solar flux by a factor:'
+            !   print*, ' Please check new_chemistry* and sw_venus_rh ' 
+            !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
+
             if (jphot(iso2,isza,iz,ij) == 1.E-30) then 
                jphot(iso2,isza,iz,ij) = 0.
@@ -2738,6 +2748,8 @@
 
 !     jpl 2003
-
-      b003(:) = 1.1E-10
+!      b003(:) = 1.1E-10
+
+!      jpl 2016      
+      b003(:) = 1.2E-10
 
       nb_reaction_4 = nb_reaction_4 + 1
@@ -2787,6 +2799,8 @@
 
 !     jpl 2003
-
-      c002(:) = 2.2E-11*exp(120./t(:))
+!      c002(:) = 2.2E-11*exp(120./t(:))
+
+!     jpl 2016  
+      c002(:) = 1.8E-11*exp(180./t(:))
 
       nb_reaction_4 = nb_reaction_4 + 1
@@ -2845,6 +2859,8 @@
 
 !     christensen et al., grl, 13, 2002
-
-      c008(:) = 1.5E-12*exp(19./t(:))
+!      c008(:) = 1.5E-12*exp(19./t(:))
+
+!     jpl 2016
+      c008(:) = 3.0E-13*exp(460./t(:))
 
       nb_reaction_3 = nb_reaction_3 + 1
@@ -2924,6 +2940,8 @@
 
 !     jpl 2003
-
-      c016(:) = 2.5*1.7E-33*exp(1000./t(:))*conc(:)
+!      c016(:) = 2.5*1.7E-33*exp(1000./t(:))*conc(:)
+
+!     jpl 2016
+      c016(:) = 2.5*2.1E-33*exp(920./t(:))*conc(:)
 
       nb_reaction_3 = nb_reaction_3 + 1
@@ -3072,6 +3090,8 @@
 
 !     jpl 2006
-
-      f004(:) = 2.6E-12*exp(-350./t(:))
+!      f004(:) = 2.6E-12*exp(-350./t(:))
+
+!     jpl 2016
+      f004(:) = 1.8E-12*exp(-250./t(:))
 
       nb_reaction_4 = nb_reaction_4 + 1
@@ -3126,7 +3146,9 @@
 
 !     jpl 2006
-
-      f010(:) = 4.1E-11*exp(-450./t(:))
-
+!      f010(:) = 4.1E-11*exp(-450./t(:))
+
+!     jpl 2016
+      f010(:) = 3.6E-11*exp(-375./t(:)) 
+      
       nb_reaction_4 = nb_reaction_4 + 1
       v_4(:,nb_reaction_4) = f010(:)
@@ -3135,6 +3157,8 @@
 
 !     jpl 2006
-
-      f011(:) = 1.8E-11*exp(170./t(:))
+!      f011(:) = 1.8E-11*exp(170./t(:))
+
+!     jpl 2016
+      f011(:) = 1.4E-11*exp(270./t(:))
 
       nb_reaction_4 = nb_reaction_4 + 1
@@ -3423,5 +3447,8 @@
 !---  g001: s + o2 -> so + o
 
-      g001(:) = 2.3E-12
+!      g001(:) = 2.3E-12
+
+!     jpl 2016
+      g001(:) = 1.6E-12*exp(100./t(:))
 
       nb_reaction_4 = nb_reaction_4 + 1
@@ -3437,5 +3464,8 @@
 !---  g003: so + o2 -> so2 + o
 
-      g003(:) = 1.25E-13*exp(-2190./t(:))
+!      g003(:) = 1.25E-13*exp(-2190./t(:))
+
+!     jpl 2016
+       g003(:) = 1.6E-13*exp(-2280./t(:))
 
       nb_reaction_4 = nb_reaction_4 + 1