Changeset 29 for trunk/ToBeReviewed/TRIANGULATION/tracemask.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/tracemask.pro (copied) (copied from trunk/TRIANGULATION/tracemask.pro) (5 diffs)

trunk/ToBeReviewed/TRIANGULATION/tracemask.pro

@@ -17 +17 @@
 ; KEYWORD PARAMETERS:
 ;
-;        CONT_THICK: l''epaisseur du trait pour tracer les
+;        COAST_COLOR: the color of the coastline.
+;                     defaut value is 0 => black
+;
+;        COAST_THICK: l''epaisseur du trait pour tracer les
 ;        continents. par defaut c''est 1.
 ;
@@ -36 +39 @@
 ;------------------------------------------------------------
 ;------------------------------------------------------------
-PRO tracemask, maskentree, xin, yin, CONT_THICK = cont_thick, OVERPLOT = overplot, _extra = ex
-   xentree = xin
-   yentree = yin
+PRO tracemask, maskentree, xin, yin, COAST_COLOR = coast_color, COAST_THICK = coast_thick, OVERPLOT = overplot, _extra = ex
+;
    if keyword_set(overplot) then return
-@common
+;---------------------------------------------------------
+@cm_general
+  IF NOT keyword_set(key_forgetold) THEN BEGIN
+@updatekwd
+  ENDIF
+;--------------------------------------------------------- 
    tempsun = systime(1)         ; pour key_performance
 ; on s''afranchit des problemes de bord:
@@ -47 +54 @@
    nx = tailleentree[1]+1
    ny = tailleentree[2]+1
-; on agrandi le mask de une colonne a gauche et de une colonne a droite.
+; we check the input axis
+  IF n_elements(xin) EQ 0 THEN xentree = findgen(nx-1) ELSE xentree = xin
+  IF (size(xentree))[0] EQ 1 THEN xentree = xentree#replicate(1,ny-1)
+  IF n_elements(yin) EQ 0 THEN yentree = findgen(ny-1) ELSE yentree = yin
+  IF (size(yentree))[0] EQ 1 THEN yentree = replicate(1,nx-1)#yentree
+; on agrandi le mask de une colonne a gauche et de une colonne en bas
    mask = intarr(tailleentree[1]+1, tailleentree[2]+1)
    mask[1:tailleentree[1], 1:tailleentree[2]] = maskentree
@@ -105 +117 @@
 ; boucle sur les points concernes et trace du segment
 ; rq: on utilise plots au lieu de plot car plots est bcp plus rapide.
-         for pt = 0, n_elements(lx)-1 do BEGIN 
+         for pt = 0L, n_elements(lx)-1 do BEGIN 
             i = lx[pt] & j = ly[pt]
             plots, [xf[i, j-1], xf[i, j]], [yf[i, j-1], yf[i, j]] $
-             , color=c_cote,thick=cont_thick, _extra = ex
+              , color = coast_color, thick = coast_thick, _extra = ex
             if pt LT 5 then begin
             endif
@@ -129 +141 @@
        print, 'temps tracemask: liste traits horizontaux', systime(1)-tempdeux
       tempdeux = systime(1)     ; pour key_performance =2
-      for pt = 0, n_elements(lx)-1 do BEGIN 
+      for pt = 0L, n_elements(lx)-1 do BEGIN 
          i = lx[pt] & j = ly[pt]
          plots, [xf[i-1, j], xf[i, j]], [yf[i-1, j], yf[i, j]] $
-          , color=c_cote,thick=cont_thick, _extra = ex
+           , color = coast_color, thick = coast_thick, _extra = ex
       endfor
       IF testvar(var = key_performance) EQ 2 THEN $