Index: trunk/LMDZ.MARS/libf/phymars/newstart.F
===================================================================
--- trunk/LMDZ.MARS/libf/phymars/newstart.F	(revision 1383)
+++ trunk/LMDZ.MARS/libf/phymars/newstart.F	(revision 1390)
@@ -35,4 +35,5 @@
       use comgeomfi_h, only: ini_fillgeom
 
+
       implicit none
 
@@ -168,4 +169,5 @@
       real :: Mair_old,Mair_new,vmr_old,vmr_new
       real,allocatable :: coefvmr(:)  ! Correction coefficient when changing composition
+      integer :: iloc(1), iqmax
 
 c sortie visu pour les champs dynamiques
@@ -1021,14 +1023,14 @@
 
               if (iq.eq.igcm_n2) then
-                write(*,*) "New vmr(n2)? (MSL: 1.89E-02 at Ls~186)"
+                write(*,*) "New vmr(n2)? (MSL: 2.03e-02 at Ls~184)"
               endif
               if (iq.eq.igcm_ar) then
-                write(*,*) "New vmr(ar)? (MSL: 1.93E-02 at Ls~186)"
+                write(*,*) "New vmr(ar)? (MSL: 2.07e-02 at Ls~184)"
               endif
               if (iq.eq.igcm_o2) then
-                write(*,*) "New vmr(o2)? (MSL: 1.46E-03) at Ls~186"
+                write(*,*) "New vmr(o2)? (MSL: 1.73e-03 at Ls~184)"
               endif
               if (iq.eq.igcm_co) then
-                write(*,*) "New vmr(co)? (~8.E-04)"
+                write(*,*) "New vmr(co)? (MSL: 7.49e-04 at Ls~184)"
               endif
  302          read(*,*,iostat=ierr) vmr_new
@@ -1055,5 +1057,6 @@
           end do
 
-        ! Recompute mass mixing ratios everywhere, and adjust mmr(CO2) to keep sum of mmr constant.
+        ! Recompute mass mixing ratios everywhere, and adjust mmr of most abundant species
+        ! to keep sum of mmr constant.
           do l=1,llm
            do j=1,jjp1
@@ -1069,5 +1072,7 @@
                 end if 
               enddo
-              q(i,j,l,igcm_co2) = q(i,j,l,igcm_co2) + Smmr_old-Smmr_new
+              iloc = maxloc(q(i,j,l,:))
+              iqmax = iloc(1)
+              q(i,j,l,iqmax) = q(i,j,l,iqmax) + Smmr_old - Smmr_new
             enddo
            enddo
@@ -1075,9 +1080,10 @@
 
           write(*,*)
-     &   "mmr(CO2) is modified everywhere to keep sum of mmr constant"
+     &   "The most abundant species is modified everywhere to keep "//
+     &   "sum of mmr constant"
           write(*,*) 'At reference site vmr(CO2)=', 
      &        q(iref,jref,1,igcm_co2)*Mair_new/mmol(igcm_co2)
-          write(*,*) "Compared to MSL observation: vmr(CO2)= 0.9597" 
-
+          write(*,*) "Compared to MSL observation: vmr(CO2)= 0.957 "//
+     &   "at Ls=184" 
 
 c      wetstart : wet atmosphere with a north to south gradient