Index: LMDZ5/trunk/libf/phylmd/physiq.F90
===================================================================
--- LMDZ5/trunk/libf/phylmd/physiq.F90	(revision 2368)
+++ LMDZ5/trunk/libf/phylmd/physiq.F90	(revision 2369)
@@ -1100,5 +1100,12 @@
      ! pour obtenir le meme resultat.
      dtime=pdtphys
-     radpas = NINT( 86400./dtime/nbapp_rad)
+     IF (MOD(INT(86400./dtime),nbapp_rad).EQ.0) THEN
+       radpas = NINT( 86400./dtime/nbapp_rad)
+     ELSE 
+       WRITE(lunout,*) 'le nombre de pas de temps physique doit etre un multiple de nbapp_rad'
+       WRITE(lunout,*) 'changer nbapp_rad ou alors commenter ce test mais 1+1<>2'
+       abort_message='nbre de pas de temps physique n est pas multiple de nbapp_rad'
+       call abort_physic(modname,abort_message,1)
+     ENDIF
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
 
@@ -1848,6 +1855,6 @@
      !  radpas-1 pas de temps suivant. Petite ruse avec MOD pour prendre en
      !  compte le premier pas de temps de la physique pendant lequel itaprad=0 
-        zdtime1=dtime*REAL(-MOD(itaprad,4)-1)      
-        zdtime2=dtime*REAL(radpas-MOD(itaprad,4)-1) 
+        zdtime1=dtime*REAL(-MOD(itaprad,radpas)-1)      
+        zdtime2=dtime*REAL(radpas-MOD(itaprad,radpas)-1) 
         CALL zenang(zlongi,jH_cur,zdtime1,zdtime2,rlat,rlon,rmu0,fract)
      !