Changeset 3462

Timestamp:

Oct 17, 2024, 5:48:37 PM (2 months ago)

Author:

flefevre

Message:

Mars PCM: reverting photochemistry to clean version 3288

File:

• trunk/LMDZ.MARS/libf/aeronomars/photochemistry.F90 (modified) (11 diffs)

trunk/LMDZ.MARS/libf/aeronomars/photochemistry.F90

@@ -161 +161 @@
 !integer,parameter :: i_hdo2    = 38
 
-integer, save :: i_last
-!$OMP THREADPRIVATE(i_last)
+integer :: i_last
+
 !Tracer indexes for photionization coeffs
 
@@ -180 +180 @@
 
 integer :: ilay
-integer, save :: ind_norec
-integer, save :: ind_orec
-!$OMP THREADPRIVATE(ind_norec, ind_orec)
+integer :: ind_norec
+integer :: ind_orec
+
 real :: ctimestep           ! standard timestep for the chemistry (s) 
 real :: dt_guess            ! first-guess timestep (s) 
@@ -503 +503 @@
 em_no(:)=c(:,i_o)*c(:,i_n)*v_4(:,ind_norec)   !2.8e-17*(300./temp(:)))**0.5
 em_o2(:)=0.75*c(:,i_o)*c(:,i_o)*c(:,i_co2)*v_3(:,ind_orec)   !2.5*9.46e-34*exp(485./temp(:))*dens(:)
-!print*, "ind_norec= ", ind_norec
-!print*, "ind_orec= ", ind_orec
+
 !===================================================================
 !     save chemical species for the gcm       
@@ -660 +659 @@
 
 use comcstfi_h
-!use photolysis_mod, only : nphot
+use photolysis_mod, only : nphot
 
 implicit none
@@ -676 +675 @@
 logical, intent(in)     :: deutchem
 integer                 :: lswitch           ! interface level between lower
-integer,parameter :: nphot0 =13
                                              ! atmosphere and thermosphere chemistries
 real, dimension(nlayer) :: dens              ! total number density (molecule.cm-3)
@@ -697 +695 @@
 real (kind = 8), dimension(nlayer,nb_reaction_3_max) :: v_3
 real (kind = 8), dimension(nlayer,nb_reaction_4_max) :: v_4
-integer, intent(out) :: ind_norec
-integer, intent(out) :: ind_orec
+integer :: ind_norec
+integer :: ind_orec
 
 !----------------------------------------------------------------------
@@ -705 +703 @@
 
 integer          :: ilev
-integer, save    :: nb_phot, nb_reaction_3, nb_reaction_4
-!$OMP THREADPRIVATE(nb_phot, nb_reaction_3, nb_reaction_4)
+integer          :: nb_phot, nb_reaction_3, nb_reaction_4
 real :: ak0, ak1, xpo, rate, rate1, rate2
 real :: k1a0, k1b0, k1ainf, k1a, k1b, fc, fx, x, y, gam
@@ -740 +737 @@
 !     initialisation
 !----------------------------------------------------------------------
-      nb_phot       = 0
-      nb_phot       = nphot0 + nphotion ! initialised to the number of photolysis + number of photoionization rates
+
+      nb_phot       = nphot + nphotion ! initialised to the number of photolysis + number of photoionization rates
       nb_reaction_3 = 0
       nb_reaction_4 = 0
@@ -779 +776 @@
       v_3(:,nb_reaction_3) = a002(:)
       ind_orec=nb_reaction_3
-!      print*,"ind_orec",ind_orec
+
 !---  a003: o + o3 -> o2 + o2
 
@@ -1203 +1200 @@
       v_4(:,nb_reaction_4) = d007(:)
       ind_norec = nb_reaction_4
-!      print*,"ind_norec=",ind_norec
+
 !---  d008: n + ho2 -> no + oh
 
@@ -2461 +2458 @@
 ! local
 
-integer, save :: nb_phot, nb_reaction_3, nb_reaction_4
-!$OMP THREADPRIVATE(nb_phot, nb_reaction_3, nb_reaction_4)
+integer :: nb_phot, nb_reaction_3, nb_reaction_4
 integer :: i_dummy