Index: /trunk/LMDZ.PLUTO/libf/dynphy_lonlat/phypluto/newstart.F
===================================================================
--- /trunk/LMDZ.PLUTO/libf/dynphy_lonlat/phypluto/newstart.F	(revision 3969)
+++ /trunk/LMDZ.PLUTO/libf/dynphy_lonlat/phypluto/newstart.F	(revision 3970)
@@ -601,5 +601,5 @@
         ithfi(:,:)=inertiedat(:,:)
         ! copy topo
-        CALL gr_dyn_fi(1,iip1,jjp1,ngridmx,phisold_newgrid,phisfi)
+        ! CALL gr_dyn_fi(1,iip1,jjp1,ngridmx,phisold_newgrid,phisfi)
         ! Initialize n2frac
         n2fracfi(1:ngridmx)=1.
@@ -1458,5 +1458,5 @@
                 ! Loop where min
                 actualmin=1000.
-                do j=1,i
+                do j=1,i-1
                    if ( (array_angle(j).lt.angle+0.52).and.
      &                  (array_angle(j).gt.angle-0.52).and.
@@ -3150,16 +3150,8 @@
 
           ! get min max lon
-          print*, latfi*180/pi
-          print*, '***************************'
-          print*, '***************************'
-          print*, '***************************'
-          print*, '***************************'
-          print*, '***************************'
-          print*, lonfi*180/pi
-          print*, 'iip1=',iip1
-          do ig=2,ngridmx-1
+          ! do ig=2,ngridmx-1
+          do ig=2+iip1,ngridmx-1-iip1
              IF  (lonfi(ig)*180./pi.eq.-180.) THEN
-                 print*, lonfi(ig),lonfi(ig+iip1-2)
-                 phisfi(ig)=(phisfi(ig+1)+phisfi(ig+iip1))/2
+                 phisfi(ig)=(phisfi(ig-1)+phisfi(ig-iip1+2))/2.
              ENDIF
           enddo