PRO drawcoast_e, mask, xf, yf, nx, ny, CONT_THICK = cont_thick, YSEUIL = yseuil, XSEUIL = xseuil, onemore = onemore, _extra = ex
@common
   tempsun = systime(1)         ; pour key_performance
;---------------------------------------------------------
   if keyword_set(key_periodique) AND nx EQ jpi then begin
      mask = [mask, mask[0, *]]
      xf = [xf, xf[0, *]]
      yf = [yf, yf[0, *]]
      nx = nx+1
   ENDIF
;
; we plot the borders of the diamond in this sense : \
;
   if NOT keyword_set(onemore) then onemore = 0
   if NOT keyword_set(xseuil) then xseuil = 5.
   distanceseuil = (!p.position[2]-!p.position[0])/xseuil
; liste: liste des points i pourlesquels on va tracer un segment 
   index = lindgen(nx, ny)
   index = index[0:nx-2, 1:ny-1]
   indexbis = index-nx+((index/nx+onemore) MOD 2)
   liste = where(mask[index+1]+mask[indexbis] EQ 1 $
                 AND (xf[index]-xf[indexbis])^2+(yf[index]-yf[indexbis])^2 LE distanceseuil)
   IF liste[0] NE -1 THEN BEGIN
      index = index[liste]
      indexbis = indexbis[liste]
      for pt = 0, n_elements(index)-1 do begin
         plots, [xf[index[pt]], xf[indexbis[pt]]], [yf[index[pt]], yf[indexbis[pt]]] $
          , color=c_cote,thick=cont_thick, /normal, _extra = ex
      endfor
   ENDIF
;
; we plot the borders of the diamond in this sense : /
;
   if NOT keyword_set(xseuil) then xseuil = 5.
   distanceseuil = (!p.position[2]-!p.position[0])/xseuil
; liste: liste des points i pourlesquels on va tracer un segment 
   index = lindgen(nx, ny-1)
   index = index[0:nx-2, *]
   indexbis = index+nx+((index/nx+onemore) MOD 2)
   liste = where(mask[index+1]+mask[indexbis] EQ 1 $
                 AND (xf[index]-xf[indexbis])^2+(yf[index]-yf[indexbis])^2 LE distanceseuil)
   IF liste[0] NE -1 THEN BEGIN
      index = index[liste]
      indexbis = indexbis[liste]
      for pt = 0, n_elements(index)-1 do begin
         plots, [xf[index[pt]], xf[indexbis[pt]]], [yf[index[pt]], yf[indexbis[pt]]] $
          , color=c_cote,thick=cont_thick, /normal, _extra = ex
      endfor
   ENDIF

;---------------------------------------------------------
   if keyword_set(key_performance) THEN print, 'temps drawcoast_e', systime(1)-tempsun
   return
end