Index: ICOSA_LMDZ/src/phylmd/icolmdz_param_gravity_wave.f90
===================================================================
--- ICOSA_LMDZ/src/phylmd/icolmdz_param_gravity_wave.f90	(revision 4897)
+++ ICOSA_LMDZ/src/phylmd/icolmdz_param_gravity_wave.f90	(revision 4966)
@@ -401,9 +401,9 @@
     INTEGER :: i,j,ij
 
-    DO j=jj_begin,jj_end
-      DO i=ii_begin,ii_end
+    DO j=jj_begin-1,jj_end+1
+      DO i=ii_begin-1,ii_end+1
         ij=(j-1)*iim+i
         zstd(ij) = ztz(ij) - zmea(ij)*zmea(ij)
-        IF (zstd(ij)<0) zstd(ij)=0
+        IF (zstd(ij)<=0) zstd(ij)=0
         zstd(ij)=SQRT(zstd(ij))
       ENDDO