Index: trunk/LMDZ.MARS/libf/dyn3d/newstart.F
===================================================================
--- trunk/LMDZ.MARS/libf/dyn3d/newstart.F	(revision 165)
+++ trunk/LMDZ.MARS/libf/dyn3d/newstart.F	(revision 171)
@@ -734,4 +734,9 @@
              write(*,*) '(choose between 1 and ',nqmx,')'
              read(*,*) iq 
+             if ((iq.lt.1).or.(iq.gt.nqmx)) then
+               ! wrong value for iq, go back to menu
+               write(*,*) "wrong input value:",iq
+               cycle
+             endif
              write(*,*)'mixing ratio of tracer ',trim(tnom(iq)),
      &                 ' ? (kg/kg)'
@@ -1416,4 +1421,6 @@
      *rage est differente de la valeur parametree llm =',i4//)
 
+      write(*,*) "newstart: All is well that ends well."
+
       end
 
Index: trunk/LMDZ.MARS/libf/phymars/initracer.F
===================================================================
--- trunk/LMDZ.MARS/libf/phymars/initracer.F	(revision 165)
+++ trunk/LMDZ.MARS/libf/phymars/initracer.F	(revision 171)
@@ -231,4 +231,19 @@
       igcm_ar=0
       igcm_ar_n2=0
+      igcm_n=0
+      igcm_no=0
+      igcm_no2=0
+      igcm_n2d=0
+      igcm_co2plus=0
+      igcm_oplus=0
+      igcm_o2plus=0
+      igcm_coplus=0
+      igcm_cplus=0
+      igcm_nplus=0
+      igcm_noplus=0
+      igcm_n2plus=0
+      igcm_hplus=0
+      igcm_elec=0
+
 
       ! 1. find dust tracers
@@ -333,4 +348,74 @@
           count=count+1
         endif
+	if (noms(iq).eq."n") then
+          igcm_n=iq
+          mmol(igcm_n)=14.
+          count=count+1
+        endif
+        if (noms(iq).eq."no") then
+          igcm_no=iq
+          mmol(igcm_no)=30.
+          count=count+1
+        endif
+        if (noms(iq).eq."no2") then
+          igcm_no2=iq
+          mmol(igcm_no2)=46.
+          count=count+1
+        endif
+        if (noms(iq).eq."n2d") then
+          igcm_n2d=iq
+          mmol(igcm_n2d)=28.
+          count=count+1
+        endif
+        if (noms(iq).eq."co2plus") then
+          igcm_co2plus=iq
+          mmol(igcm_co2plus)=44.
+          count=count+1
+        endif
+        if (noms(iq).eq."oplus") then
+          igcm_oplus=iq
+          mmol(igcm_oplus)=16.
+          count=count+1
+        endif
+        if (noms(iq).eq."o2plus") then
+          igcm_o2plus=iq
+          mmol(igcm_o2plus)=32.
+          count=count+1
+        endif
+        if (noms(iq).eq."coplus") then
+          igcm_coplus=iq
+          mmol(igcm_coplus)=28.
+          count=count+1
+        endif
+        if (noms(iq).eq."cplus") then
+          igcm_cplus=iq
+          mmol(igcm_cplus)=12.
+          count=count+1
+        endif
+        if (noms(iq).eq."nplus") then
+          igcm_nplus=iq
+          mmol(igcm_nplus)=14.
+          count=count+1
+        endif
+        if (noms(iq).eq."noplus") then
+          igcm_noplus=iq
+          mmol(igcm_noplus)=30.
+          count=count+1
+        endif
+	if (noms(iq).eq."n2plus") then
+          igcm_n2plus=iq
+          mmol(igcm_n2plus)=28.
+          count=count+1
+        endif
+        if (noms(iq).eq."hplus") then
+          igcm_hplus=iq
+          mmol(igcm_hplus)=1.
+          count=count+1
+        endif
+        if (noms(iq).eq."elec") then
+          igcm_elec=iq
+          mmol(igcm_elec)=1./1822.89
+          count=count+1
+        endif
         if (noms(iq).eq."h2o_vap") then
           igcm_h2o_vap=iq
@@ -349,4 +434,7 @@
           count=count+1
         endif
+	
+
+
       enddo ! of do iq=1,nqmx
 !      count=count+nbqchem
Index: trunk/LMDZ.MARS/libf/phymars/tracer.h
===================================================================
--- trunk/LMDZ.MARS/libf/phymars/tracer.h	(revision 165)
+++ trunk/LMDZ.MARS/libf/phymars/tracer.h	(revision 171)
@@ -46,5 +46,19 @@
       integer :: igcm_n2
       integer :: igcm_ar
-      
+      integer :: igcm_n
+      integer :: igcm_no
+      integer :: igcm_no2
+      integer :: igcm_n2d
+      ! Ions
+      integer :: igcm_co2plus
+      integer :: igcm_oplus
+      integer :: igcm_o2plus
+      integer :: igcm_coplus
+      integer :: igcm_cplus
+      integer :: igcm_nplus
+      integer :: igcm_noplus
+      integer :: igcm_n2plus 
+      integer :: igcm_hplus
+      integer :: igcm_elec
       ! other tracers
       integer :: igcm_ar_n2 ! for simulations using co2 +neutral gas
@@ -59,5 +73,8 @@
      & igcm_h2o_vap,igcm_h2o_ice,igcm_co2,igcm_co,igcm_o,igcm_o1d,      &
      & igcm_o2,igcm_o3,igcm_h,igcm_h2,igcm_oh,igcm_ho2,igcm_h2o2,       &
-     & igcm_n2,igcm_ar,igcm_ar_n2!,nbqchem,niqchem
+     & igcm_n2,igcm_ar,igcm_n,igcm_no,igcm_no2,igcm_n2d,                &
+     & igcm_co2plus,igcm_oplus,igcm_o2plus,igcm_coplus,igcm_cplus,      &
+     & igcm_nplus,igcm_noplus,igcm_n2plus,igcm_hplus,igcm_elec,         &
+     & igcm_ar_n2!,nbqchem,niqchem
       COMMON/tracer3/noms
 !-----------------------------------------------------------------------