Index: trunk/LMDZ.TITAN/libf/dynphy_lonlat/phytitan/newstart.F
===================================================================
--- trunk/LMDZ.TITAN/libf/dynphy_lonlat/phytitan/newstart.F	(revision 1848)
+++ trunk/LMDZ.TITAN/libf/dynphy_lonlat/phytitan/newstart.F	(revision 1849)
@@ -368,16 +368,20 @@
 c-----------------------------------------------------------------------
 
-      kappa = tab_cntrl(9) 
-      etot0 = tab_cntrl(12)
-      ptot0 = tab_cntrl(13)
-      ztot0 = tab_cntrl(14)
-      stot0 = tab_cntrl(15)
-      ang0 = tab_cntrl(16)
+      ! JVO 2017 :  For this part it seems to need a shift of 1 in the 
+      ! tab_cntrl indexes compared to dynetat0 ... I introduced it to have
+      ! correct outputs of newstart but I don't really understand why ...
+
+      kappa = tab_cntrl(10) 
+      etot0 = tab_cntrl(13)
+      ptot0 = tab_cntrl(14)
+      ztot0 = tab_cntrl(15)
+      stot0 = tab_cntrl(16)
+      ang0 = tab_cntrl(17)
       write(*,*) "Newstart: kappa,etot0,ptot0,ztot0,stot0,ang0"
       write(*,*) kappa,etot0,ptot0,ztot0,stot0,ang0
 
       ! for vertical coordinate
-      preff=tab_cntrl(18) ! reference surface pressure
-      pa=tab_cntrl(17)  ! reference pressure at which coord is purely pressure
+      preff=tab_cntrl(19) ! reference surface pressure
+      pa=tab_cntrl(18)  ! reference pressure at which coord is purely pressure
       !NB: in start_archive files tab_cntrl(17)=tab_cntrl(18)=0
       write(*,*) "Newstart: preff=",preff," pa=",pa
@@ -389,9 +393,9 @@
       if (yes.eq.'y') then
         write(*,*)"New value of reference surface pressure preff?"
-        write(*,*)"   (for Mars, typically preff=610)"
+        write(*,*)"   (for Titan, typically preff=140000)"
         read(*,*) preff
         write(*,*)"New value of reference pressure pa for purely"
         write(*,*)"pressure levels (for hybrid coordinates)?"
-        write(*,*)"   (for Mars, typically pa=20)"
+        write(*,*)"   (for Titan, typically pa=10000)"
         read(*,*) pa
       endif