Index: LMDZ6/trunk/libf/phylmd/atke_exchange_coeff_mod.F90
===================================================================
--- LMDZ6/trunk/libf/phylmd/atke_exchange_coeff_mod.F90	(revision 4462)
+++ LMDZ6/trunk/libf/phylmd/atke_exchange_coeff_mod.F90	(revision 4463)
@@ -137,6 +137,6 @@
 ELSE ! TO DO
     
-    print*, 'traitement non-stationnaire de la tke pas encore prevu'
-    stop
+   call abort_physic("atke_compute_km_kh", &
+        'traitement non-stationnaire de la tke pas encore prevu', 1)
 
 END IF
Index: LMDZ6/trunk/libf/phylmd/phyaqua_mod.F90
===================================================================
--- LMDZ6/trunk/libf/phylmd/phyaqua_mod.F90	(revision 4462)
+++ LMDZ6/trunk/libf/phylmd/phyaqua_mod.F90	(revision 4463)
@@ -135,6 +135,5 @@
     if (year_len.ne.360) then
       write (*,*) year_len
-      write (*,*) 'iniaqua: 360 day calendar is required !'
-      stop
+      call abort_physic("iniaqua", 'iniaqua: 360 day calendar is required !', 1)
     endif
 
Index: LMDZ6/trunk/libf/phylmd/thermcell_down.F90
===================================================================
--- LMDZ6/trunk/libf/phylmd/thermcell_down.F90	(revision 4462)
+++ LMDZ6/trunk/libf/phylmd/thermcell_down.F90	(revision 4463)
@@ -58,5 +58,6 @@
 
    if ( iflag_thermals_down < 10 ) then
-        stop 'thermcell_down_dq = 0 or >= 10'
+      call abort_physic("thermcell_updown_dq", &
+           'thermcell_down_dq = 0 or >= 10', 1)
    else
         iflag_impl=iflag_thermals_down-10
@@ -128,6 +129,6 @@
          !!!! tentative de prise en compte d'un flux compensatoire montant  !!!!
          if (fup(ig,ilay)-fdn(ig,ilay) .lt. 0.) then
-            write(*,*) 'flux compensatoire montant, cas non traite par thermcell_updown_dq'
-            stop          
+            call abort_physic("thermcell_updown_dq", 'flux compensatoire '&
+                 // 'montant, cas non traite par thermcell_updown_dq', 1)
             !fthe(ig,ilay)=(fup(ig,ilay)-fdn(ig,ilay))*trac(ig,ilay-1)
          else
@@ -161,6 +162,6 @@
        do ig=1,ngrid
          if (fup(ig,ilay)-fdn(ig,ilay) .lt. 0.) then
-            write(*,*) 'flux compensatoire montant, cas non traite par thermcell_updown_dq'
-            stop          
+            call abort_physic("thermcell_updown_dq", 'flux compensatoire ' &
+                 // 'montant, cas non traite par thermcell_updown_dq', 1)
          else
             fthe(ig,ilay)=-(fup(ig,ilay)-fdn(ig,ilay))*trac(ig,ilay)
@@ -191,5 +192,5 @@
          if((fup(ig,ilay)-fdn(ig,ilay)) .lt. 0) then
             write(*,*) 'flux compensatoire montant, cas non traite par thermcell_updown_dq dans le cas d une resolution implicite, ilay : ', ilay
-            stop  
+            call abort_physic("thermcell_updown_dq", "", 1)
          else
            mstar_inv=ptimestep/masse(ig,ilay)
@@ -205,7 +206,6 @@
 
    else
-      write(*,*) 'valeur de iflag_impl non prevue'
-      stop
-
+      call abort_physic("thermcell_updown_dq", &
+           'valeur de iflag_impl non prevue', 1)
    endif