Changeset 2554 for trunk/LMDZ.MARS/libf/aeronomars

Timestamp:

Aug 18, 2021, 5:01:18 PM (3 years ago)

Author:

emillour

Message:

Mars GCM:
Make NO and O2 Nightglow emissions computations more flexible (reaction index

no longer hard coded in photochemistry.F90).

FGG

File:

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

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

@@ -175 +175 @@
 
 integer :: ilay
+integer :: ind_norec
+integer :: ind_orec
 
 real :: ctimestep           ! standard timestep for the chemistry (s) 
@@ -380 +382 @@
                    c(:,i_o2), c(:,i_o), c(:,i_n2), press, temp,           &
                    temp_elect, hetero_dust, hetero_ice,                   &
-                   surfdust1d, surfice1d, v_phot, v_3, v_4)
+                   surfdust1d, surfice1d, v_phot, v_3, v_4, ind_norec,    &
+                   ind_orec)
 
 !===================================================================
@@ -492 +495 @@
 
 !add calculation of NO and O2 nightglows
-em_no(:)=c(:,i_o)*c(:,i_n)*v_4(:,26)   !2.8e-17*(300./temp(:)))**0.5
-em_o2(:)=0.75*c(:,i_o)*c(:,i_o)*c(:,i_co2)*v_3(:,1)   !2.5*9.46e-34*exp(485./temp(:))*dens(:)
+
+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(:)
 
 !===================================================================
@@ -635 +639 @@
                           nb_phot_max, nphotion, lswitch, dens, co2, o2, o,      &
                           n2, press, t, t_elect, hetero_dust, hetero_ice,        &
-                          surfdust1d, surfice1d, v_phot, v_3, v_4)
+                          surfdust1d, surfice1d, v_phot, v_3, v_4,ind_norec,   &
+                          ind_orec)
  
 !================================================================
@@ -685 +690 @@
 real (kind = 8), dimension(nlayer,nb_reaction_3_max) :: v_3
 real (kind = 8), dimension(nlayer,nb_reaction_4_max) :: v_4
+integer :: ind_norec
+integer :: ind_orec
 
 !----------------------------------------------------------------------
@@ -758 +765 @@
       nb_reaction_3 = nb_reaction_3 + 1
       v_3(:,nb_reaction_3) = a002(:)
+      ind_orec=nb_reaction_3
 
 !---  a003: o + o3 -> o2 + o2
@@ -1136 +1144 @@
       nb_reaction_4 = nb_reaction_4 + 1
       v_4(:,nb_reaction_4) = d007(:)
+      ind_norec = nb_reaction_4
 
 !---  d008: n + ho2 -> no + oh