Index: trunk/LMDZ.MARS/libf/dyn3d/newstart.F
===================================================================
--- trunk/LMDZ.MARS/libf/dyn3d/newstart.F	(revision 311)
+++ trunk/LMDZ.MARS/libf/dyn3d/newstart.F	(revision 321)
@@ -347,4 +347,10 @@
       CALL gr_dyn_fi(1,iip1,jjp1,ngridmx,aire,airefi)
 
+! also initialize various physics flags/settings which might be needed
+!    (for instance initracer needs to know about some flags, and/or
+!      'datafile' path may be changed by user)
+      call inifis(ngridmx,llm,day_ini,daysec,dtphys,
+     &                latfi,lonfi,airefi,rad,g,r,cpp)
+
 c=======================================================================
 c   lecture topographie, albedo, inertie thermique, relief sous-maille
@@ -361,8 +367,4 @@
         relief="mola"
 c     enddo
-
-! before using datareadnc, "datafile" must be set (normaly done in inifis)
-      datafile="/u/forget/WWW/datagcm/datafile" ! default value
-      call getin("datadir",datafile) ! in case user specified another path
 
       CALL datareadnc(relief,phis,alb,surfith,z0S,
@@ -457,4 +459,7 @@
         tnom(1:nqmx)=noms(1:nqmx)
         write(*,*) 'Newstart: updated tracer names'
+      else
+       ! initialize tracer names and indexes (igcm_co2, igcm_h2o_vap, ...)
+        call initracer(qsurf,co2ice)
       endif
 
@@ -883,7 +888,9 @@
           DO l=1,llm
             DO j=1,jjp1
-              DO i=1,iip1
+              DO i=1,iip1-1
                 q(i,j,l,igcm_h2o_vap)=150.e-6 * (rlatu(j)+pi/2.) / pi
               ENDDO
+              ! We want to have the very same value at lon -180 and lon 180
+              q(iip1,j,l,igcm_h2o_vap) = q(1,j,l,igcm_h2o_vap)
             ENDDO
           ENDDO