Changeset 29 for trunk/ToBeReviewed/TRIANGULATION/drawcoast_e.pro

Timestamp:

05/02/06 15:32:01 (18 years ago)

Author:

pinsard

Message:

upgrade of TRIANGULATION according to cerbere.lodyc.jussieu.fr: /usr/home/smasson/SAXO_RD/ : files

File:

• trunk/ToBeReviewed/TRIANGULATION/drawcoast_e.pro (copied) (copied from trunk/TRIANGULATION/drawcoast_e.pro) (7 diffs)

trunk/ToBeReviewed/TRIANGULATION/drawcoast_e.pro

@@ -1 @@
-PRO drawcoast_e, mask, xf, yf, nx, ny, CONT_THICK = cont_thick, YSEUIL = yseuil, XSEUIL = xseuil, onemore = onemore, _extra = ex
-@common
+PRO drawcoast_e, mask, xf, yf, nx, ny, COAST_COLOR = coast_color, COAST_THICK = coast_thick, YSEUIL = yseuil, XSEUIL = xseuil, onemore = onemore, _extra = ex
+;---------------------------------------------------------
+@cm_4mesh
+  IF NOT keyword_set(key_forgetold) THEN BEGIN
+@updatenew
+@updatekwd
+  ENDIF
+;---------------------------------------------------------
    tempsun = systime(1)         ; pour key_performance
 ;---------------------------------------------------------
-   if keyword_set(key_periodique) AND nx EQ jpi then begin
+   if keyword_set(key_periodic) AND nx EQ jpi then begin
       mask = [mask, mask[0, *]]
       xf = [xf, xf[0, *]]
@@ -13 +19 @@
 ;
    if NOT keyword_set(onemore) then onemore = 0
-   if NOT keyword_set(xseuil) then xseuil = 5.
+   if NOT keyword_set(xseuil) then xseuil = 5. < (min([nx, ny])-2)
    distanceseuil = (!p.position[2]-!p.position[0])/xseuil
 ; liste: liste des points i pourlesquels on va tracer un segment 
@@ -20 +26 @@
    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)
+                 AND (xf[index]-xf[indexbis])^2+(yf[index]-yf[indexbis])^2 LE distanceseuil^2)
    IF liste[0] NE -1 THEN BEGIN
       index = index[liste]
@@ -26 +32 @@
       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
+           , color = coast_color, thick = coast_thick, /normal, _extra = ex
       endfor
    ENDIF
@@ -32 +38 @@
 ; we plot the borders of the diamond in this sense : /
 ;
-   if NOT keyword_set(xseuil) then xseuil = 5.
+   if NOT keyword_set(xseuil) then xseuil = 5. < (min([nx, ny])-2)
    distanceseuil = (!p.position[2]-!p.position[0])/xseuil
 ; liste: liste des points i pourlesquels on va tracer un segment 
@@ -39 +45 @@
    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)
+                 AND (xf[index]-xf[indexbis])^2+(yf[index]-yf[indexbis])^2 LE distanceseuil^2)
    IF liste[0] NE -1 THEN BEGIN
       index = index[liste]
@@ -45 +51 @@
       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
+           , color = coast_color, thick = coast_thick, /normal, _extra = ex
       endfor
    ENDIF