Index: LMDZ6/trunk/libf/dynphy_lonlat/calfis.F
===================================================================
--- LMDZ6/trunk/libf/dynphy_lonlat/calfis.F	(revision 4463)
+++ LMDZ6/trunk/libf/dynphy_lonlat/calfis.F	(revision 4464)
@@ -186,5 +186,5 @@
          write(lunout,*) '  ngridmx  jjm   iim   '
          write(lunout,*) ngridmx,jjm,iim
-         STOP
+         call abort_gcm("calfis", "", 1)
         ENDIF
       ELSE
Index: LMDZ6/trunk/libf/dynphy_lonlat/calfis_loc.F
===================================================================
--- LMDZ6/trunk/libf/dynphy_lonlat/calfis_loc.F	(revision 4463)
+++ LMDZ6/trunk/libf/dynphy_lonlat/calfis_loc.F	(revision 4464)
@@ -262,5 +262,5 @@
           write(lunout,*) '  ngridmx  jjm   iim   ' 
           write(lunout,*) ngridmx,jjm,iim
-          STOP
+          call abort_gcm("calfis_loc", "", 1)
         ENDIF
 c$OMP MASTER
@@ -1218,7 +1218,6 @@
 
 #else 
-      write(lunout,*) 
-     & "calfis_p: for now can only work with parallel physics" 
-      stop 
+      call abort_gcm("calfis_loc",
+     & "calfis_p: for now can only work with parallel physics", 1)
 #endif 
 ! of #ifdef CPP_PHYS
Index: LMDZ6/trunk/libf/dynphy_lonlat/gr_dyn_fi.F
===================================================================
--- LMDZ6/trunk/libf/dynphy_lonlat/gr_dyn_fi.F	(revision 4463)
+++ LMDZ6/trunk/libf/dynphy_lonlat/gr_dyn_fi.F	(revision 4464)
@@ -23,6 +23,5 @@
 
       IF (ngrid.NE.2+(jm-2)*(im-1)) then
-         print *, 'probleme de dim'
-         stop 1
+         call abort_gcm("gr_dyn_fi", 'probleme de dim', 1)
       end if
 c   traitement des poles
Index: LMDZ6/trunk/libf/dynphy_lonlat/phylmd/limit_netcdf.F90
===================================================================
--- LMDZ6/trunk/libf/dynphy_lonlat/phylmd/limit_netcdf.F90	(revision 4463)
+++ LMDZ6/trunk/libf/dynphy_lonlat/phylmd/limit_netcdf.F90	(revision 4464)
@@ -599,5 +599,5 @@
          CALL pchfe_95(timeyear, champtime(i, j, :), yder, skip, &
               arth(0.5, real(ndays_in) / ndays, ndays), champan(i, j, :), ierr)
-         if (ierr < 0) stop 1
+         if (ierr < 0) call abort_physiq("get_2Dfield", "", 1)
          n_extrap = n_extrap + ierr
        END DO