Changeset 67 for trunk/ToBeReviewed/PLOTS/DESSINE/pltt.pro

Timestamp:

05/09/06 17:04:13 (18 years ago)

Author:

pinsard

Message:

miscellaneous modifications according to cerbere.lodyc.jussieu.fr: /usr/home/smasson/SAXO_RD/

File:

• trunk/ToBeReviewed/PLOTS/DESSINE/pltt.pro (modified) (26 diffs)

trunk/ToBeReviewed/PLOTS/DESSINE/pltt.pro

@@ -24 +24 @@
 ;            valmask) et type doit qd meme etre specifie pour qu''on
 ;            sache de quel trace il sagit. Pour avoir une legende
-;            correcte, respecifier la zone d''extraction via BOITE
+;            correcte, respecifier la zone d''extraction via BOXZOOM
 ;          * 1d: uniquement pour les traces de type 't'. Type doit qd
 ;            meme etre specifie pour qu''on sache de quel trace il
 ;            sagit. Pour avoir une legende correcte, respecifier la zone
-;            d''extraction via BOITE
+;            d''extraction via BOXZOOM
 ;       2) une structure repondant aux critaire specifies par
 ;       litchamp.pro. cf. IDL> xhelp,'litchamp'. Le tableau contennu ds
@@ -53 +53 @@
 ; KEYWORD PARAMETERS:
 ;
-;       BOITE:vecteur indiquant la zone geographique (3d) sur laquelle doit
+;       BOXZOOM:vecteur indiquant la zone geographique (3d) sur laquelle doit
 ; etre fait l''extraction du champ pour faire le hovmoeller.
-; Si BOITE a :
-; 1 element : l''extraction est faite sur [lon1, lon2, lat1, lat2, 0.,boite[0]]
-; 2 elements: l''extraction est faite sur [lon1, lon2, lat1, lat2, boite[0],boite[1]]
-; 4 elements: l''extraction est faite sur [Boite, prof1, prof2]
-; 5 elements: l''extraction est faite sur [Boite[0:3], 0, Boite[4]]
-; 6 elements: l''extraction est faite sur Boite
-;
-; Ou lon1, lon2,lat1,lat2,prof1,prof2 sont les variables globales
+; Si BOXZOOM a :
+; 1 element : l''extraction est faite sur [lon1, lon2, lat1, lat2, 0.,boxzoom[0]]
+; 2 elements: l''extraction est faite sur [lon1, lon2, lat1, lat2, boxzoom[0],boxzoom[1]]
+; 4 elements: l''extraction est faite sur [Boxzoom, vert1, vert2]
+; 5 elements: l''extraction est faite sur [Boxzoom[0:3], 0, Boxzoom[4]]
+; 6 elements: l''extraction est faite sur Boxzoom
+;
+; Ou lon1, lon2,lat1,lat2,vert1,vert2 sont les variables globales
 ; definies lors du dernier domdef!
 ;
@@ -71 +71 @@
 ;      CB_CHARSIZE: The character size of the color bar annotations
 ;
+;        COAST_COLOR: the color of the coastline.
+;                     defaut value is 0 => black
+;
+;        COAST_THICK: the thickness of the coastline.
+;                     defaut value is 1
+;
+;        /CONT_NOFILL: activer pour ne pas remplir les points masques
+;        pour les laisser en transparent! 
+;        Rq: on trace qd meme le contour du mask
+;
+;        CONT_COLOR: the color of the continent. defaut value is
+;        (!d.n_colors - 1) < 255 => white
+;
 ;       CONTINTERVALLE: lorsque CONTOUR est active, valeur d''un
-; intervalle entre deux isolignes traces par un trait ds la couleur
-; c_lab. Il peut ainsi etre different de celui specifie par INTERVALLE
+; intervalle entre deux isolignes traces par un trait. Il peut ainsi
+; etre different de celui specifie par INTERVALLE 
 ; qui, cas ce cas, ne controle que les isolignes coloriees en
 ; couleurs. Si aucun min n''est specifie, on choisit un contour min qui va
@@ -86 +99 @@
 ;       CONTMAX: lorsque CONTOUR est active, valeur maximum que l''on
 ; veut prendre en compte dans le trace des isolignes traces par un
-; trait de couleur c_lab. Par defaut on prend le max (sur les pts mer)
+; trait. Par defaut on prend le max (sur les pts mer)
 ; du tableau passe ds le mot cle CONTOUR. 
 ;
 ;       CONTMIN: lorsque CONTOUR est active, valeur minimum que l''on
 ; veut prendre en compte dans le trace des isolignes traces par un
-; trait de couleur c_lab. Par defaut on prend le max (sur les pts mer)
+; trait. Par defaut on prend le max (sur les pts mer)
 ; du tableau passe ds le mot cle CONTOUR. 
 ;
 ;       CONTNLEVEL: lorsque CONTOUR est active, nombre de contours
-; trace par un trait de couleur c_lab a dessiner (actif si
+; trace par un trait a dessiner (actif si
 ; CONTLABEL=0) par defaut =20.
 ;
@@ -118 +131 @@
 ; INTERVALLE ne specifie que intervalle entre 2 isolignes coloriees en
 ; couleur. Pour specifier l''intervalle entre 2 isolignes traces par un trait
-; de couleur c_lab utiliser CONTINTERVALLE.
+; utiliser CONTINTERVALLE.
 ;
 ;       /INV: inverse le vecteur couleur utilise pour colorier le graphe
@@ -126 +139 @@
 ; label correspondant aux cas n. cf label.pro
 ; Rq: Qd CONTOUR est active, ne specifie le type de label que pour les
-; isolignes coloriees en couleur. Pour celles tracees par un trait de
-; couleur c_lab utiliser CONTLABEL.
+; isolignes coloriees en couleur. Pour celles tracees par un trait
+; utiliser CONTLABEL. 
 ;
 ;       /LANDSCAPE: oblige la feuille ou le fenetre a l''ecran a etre
@@ -162 +175 @@
 ;       les couleurs
 ;
-;       /NOCOULEUR: activer si on veut juste les contours en noir et
+;       /NOFILL: activer si on veut juste les contours en noir et
 ;       blanc sur fond blanc
 ;
@@ -246 +259 @@
 ; d'ouvrir plusieurs fenetres). par defaut on ouvre 'IDL0'
 ;
-;       YSURX: rapport d'echelle entre y et x. Par defaut=1
+;       YXASPECT: rapport d'echelle entre y et x. Par defaut=1
 ; Par defaut on adapte pour occupe une grande partie de la feuille en sortie
 ; postScript ou pour cree une fenetre pas trop etiree
+;
+;       /ZEROMIDDLE: fo force the middle of the colorbar to be equal
+;       to 0 (force max=max(abs([min,max])) and min=-max)
 ;
 ; COMMON BLOCKS:
@@ -277 +293 @@
 ;------------------------------------------------------------
 ;------------------------------------------------------------
-pro pltt,tab,giventype,givenmin,givenmax,datmin,datmax,BOITE = boite, CONTOUR=contour $
+pro pltt,tab,giventype,givenmin,givenmax,datmin,datmax,BOXZOOM = boxzoom, CONTOUR=contour $
          ,ENDPOINTS=endpoints,INTERVALLE=intervalle,INV=inv  $
          ,CONTINTERVALLE=contintervalle,LABEL=label,CONTLABEL=contlabel $
-         ,STYLE=style, CONTMAX=contmax, USETRI = usetri $
+         ,STYLE=style, CONTMAX=contmax $
          ,CONTMIN=contmin,NLEVEL=nlevel,CONTNLEVEL=contnlevel $
          ,COL1D=col1d,STY1D=sty1d, MININ = minin, MAXIN = maxin $
          ,OV1D=ov1d, FILTER = filter, TREND_TYPE = trend_type $
          ,REPEAT_C = repeat_c, TYPEIN = typein, XT = XT, YT = YT, ZT = zt $
-         , TT = tt, STRICTFILL = strictfill, OVERPLOT = overplot, EXCHANGE_XY = exchange_xy $
+         , TT = tt, STRICTFILL = strictfill, OVERPLOT = overplot $
+         , EXCHANGE_XY = exchange_xy $
          , _extra = ex
-@common
+;---------------------------------------------------------
+; include common
+@cm_4mesh
+@cm_4data
+@cm_4cal
+  IF NOT keyword_set(key_forgetold) THEN BEGIN
+@updatenew
+@updatekwd
+  ENDIF
+;--------------------------------------------------------------
 ;------------------------------------------------------------
    tempsun = systime(1)         ; pour key_performance
@@ -299 +325 @@
 ; I1) lecture du champ
 ;--------------------------------------------------------------
-   if keyword_set(boite) then oldboite = [lon1, lon2, lat1, lat2, prof1, prof2]
+   if (keyword_set(boxzoom) OR keyword_set(endpoints)) $
+     AND n_elements(contour) ne 4 THEN BEGIN
+     savedbox = 1b
+     saveboxparam, 'boxparam4pltt.dat'
+   ENDIF
    if n_elements(giventype) NE 0 then type = giventype
    if n_elements(givenmin) NE 0 then min = givenmin
@@ -311 +341 @@
       endif
       type = typein
-   endif
-   checktypeminmax, 'pltt', TYPE = type, MIN = min, MAX = max, XT = XT, YT = YT, ZT = zt, TT = tt, ENDPOINTS=endpoints, BOITE = boite, _extra = ex
+    ENDIF
+;
+   checktypeminmax, 'pltt', TYPE = type, MIN = min, MAX = max $
+     , XT = XT, YT = YT, ZT = zt, TT = tt, ENDPOINTS = endpoints, _extra = ex
+;
    if keyword_set(endpoints) then begin
-      section, tab, z2d, glam, gphi, ENDPOINTS = endpoints, TYPE = type, BOITE = boite, DIREC = direc
+      section, tab, z2d, glam, gphi, ENDPOINTS = endpoints, TYPE = type $
+        , BOXZOOM = boxzoom, DIREC = direc
       nx = n_elements(glam) 
       ny = nx
       if strupcase(vargrid) EQ 'W' then begin
-         z = gdepw[premierzw:dernierzw]
+         z = gdepw[firstzw:lastzw]
          nz = nzw
       ENDIF ELSE BEGIN
-         z = gdept[premierzt:dernierzt]
+         z = gdept[firstzt:lastzt]
          nz = nzt
       ENDELSE
    ENDIF ELSE BEGIN 
-      z2d = checkfield(tab, 'pltt', TYPE = type, BOITE = boite, direc = direc, _extra = ex)
+      z2d = checkfield(tab, 'pltt', TYPE = type, BOXZOOM = boxzoom $
+                       , direc = direc, _extra = ex)
       if z2d[0] EQ -1 then BEGIN 
-         if keyword_set(boite) AND n_elements(oldboite) NE 0 then domdef, oldboite,GRILLE=vargrid
-         return
+        IF keyword_set(savedbox) THEN restoreboxparam, 'boxparam4pltt.dat'
+        return
       endif
       grille, mask, glam, gphi, gdep, nx, ny,nz
@@ -371 +406 @@
    ENDIF 
 
-
 ;---------------------------------------------------------------
 ; selection du type de graphique. 
@@ -377 +411 @@
    taille=size(z2d)
    case taille[0] of
-      2 : begin
-         niveau=1
-         typdes='2d'
-      end
+      2 : typdes='2d'
       1 : begin
          z1d=z2d
@@ -401 +432 @@
 ; on construit le mask pour cela le tableau doit etre masque (fait
 ; automatiquement a la valeur valmask si on passe ds moyenne ou grossemoyenne)
-   mask = fltarr(taille[1], taille[2])
    nan = total(finite(z2d,/nan)) < 1
    if n_elements(valmask) EQ 0 then valmask = 1e20
@@ -407 +437 @@
       notanum = where(finite(z2d) EQ 0)
       z2d[notanum] = 0.
-      mask(where(z2d LT valmask/10.)) = 1
+      mask = z2d LT valmask/10.
       z2d[notanum] = !values.f_nan
-   ENDIF ELSE mask(where(z2d LT valmask/10)) = 1
+   ENDIF ELSE mask = z2d LT valmask/10.
 ;----------------------------------------------------------------------------
 ;   determination du mi:min et du ma:max de z2d ainsi que de max: max et
@@ -416 +446 @@
 ; faudra-t-il faire un autoscale ???
    autoscale = testvar(var = min) EQ testvar(var = max) AND NOT keyword_set(intervalle)
-   determineminmax, z2d, mask, mi, ma, MININ = min, MAXIN = max, nan = nan, INTERVALLE = intervalle
+   determineminmax, z2d, mask, mi, ma, MININ = min, MAXIN = max, nan = nan, INTERVALLE = intervalle, _extra = ex
    if z2d[0] EQ -1 THEN return
 ; on fait un autoscale si besoin
@@ -505 +535 @@
       if keyword_set(intervalle) AND NOT keyword_set(label) then label=1
       if keyword_set(label) eq 0 then cas=0 else cas=label
-      label,cas,min,max,ncontour,level_z2d,NLEVEL=nlevel,INTERVALLE=intervalle, strictfill = strictfill
+      label, cas, min, max, ncontour, level_z2d, colnumb, NLEVEL = nlevel $
+             , INTERVALLE = intervalle, strictfill = strictfill
 ;--------------------------------------------------------------
 ; choix de style
@@ -511 +542 @@
       if not keyword_set(style) then style=0
       style,style,level_z2d,linestyle,thick
-      if keyword_set(inv) then couleur=reverse(couleur)
+      if keyword_set(inv) then colnumb=reverse(colnumb)
 ;----------------------------------------------------------------------
       nby = n_elements(yy)
       nbx = n_elements(xx)
-      xx = xx#replicate(1, nby) ; on passe les axes en tableaux 2d
-      yy = replicate(1, nbx)#yy
-;
-      if keyword_set(nan) then begin
-         notanumber = where(finite(z2d, /nan) EQ 1)
-         z2d[notanumber] = max
-         case type of
-            'xt':masknan = replicate(1, nx, jpt)
-            'yt':masknan = replicate(1, ny, jpt)
-            'zt':masknan = replicate(1, nz, jpt)
-         endcase
-         masknan[notanumber] = 0
-      ENDIF ELSE masknan = 1
-      z2d = remplit(z2d,nite=2+keyword_set(nan), mask = mask*masknan, /basique, _extra=ex)
+;
+      if keyword_set(nan) then BEGIN
+;         xx = xx#replicate(1, nby) ; on passe les axes en tableaux 2d
+;         yy = replicate(1, nbx)#yy
+         masknan = finite(z2d)
+         z2d[where(masknan EQ 0)] = max
+         mask = temporary(mask)*temporary(masknan)
+      ENDIF
+      z2d = remplit(z2d,nite=2+keyword_set(nan), mask = mask, /basique, _extra=ex)
       if NOT keyword_set(strictfill) then z2d = min > z2d <  max
-      if NOT keyword_set(usetri) then triangulation = -1 ELSE $
-       triangulation = triangule(mask*masknan,/regulier,coinmonte=coinmontemask,coindescend=coindescendmask)
-;
+      if keyword_set(nan) then $
+       triangulation = triangule(mask,/basic,coinmonte=coinmontemask $
+                                 ,coindescend=coindescendmask) $
+       ELSE triangulation = -1 
 ;----------------------------------------------------------------------
-      pltbase,z2d, xx, yy, mask,xx, yy, level_z2d,couleur, contour = contour,/noerase $
+      pltbase,z2d, xx, yy, mask,xx, yy, level_z2d,colnumb, contour = contour,/noerase $
        , c_linestyle=linestyle,c_labels=1-(indgen(n_elements(level_z2d)) MOD 2) $
        , trichamp = triangulation, trimsk = triangulation, overplot = overplot $
-       , c_thick=thick, performance = key_performance $
-       , coinmontemask=-1, coindescendmask=-1, masknan = masknan, _extra = ex
+       , c_thick=thick, performance = key_performance, usetri = keyword_set(nan) $
+       , coinmontemask=coinmontemask, coindescendmask=coindescendmask, _extra = ex
 ;------------------------------------------------------------
       if n_elements(contour) eq 4 then BEGIN ; c''est la 2eme fois que je passe ds pltt
@@ -548 +575 @@
          pourlegende = [1, 1, 1, 1]
          oldattributs = saveatt()
-         oldcouleur = couleur
+         oldcolnumb = colnumb
          pltt,contour,type, contmin,contmax,CONTOUR=pourlegende, /noerase, USETRI = usetri $
           ,INTERVALLE=contintervalle,LABEL=contlabel,STYLE=style, ENDPOINTS = endpoints $
-          ,NLEVEL=contnlevel,YSURX=ysurx, BOITE = boite, STRICTFILL = strictfill, _extra = ex
-         couleur = oldcouleur
+           , NLEVEL = contnlevel, BOXZOOM = boxzoom, STRICTFILL = strictfill, _extra = ex
+         colnumb = oldcolnumb
          restoreatt, oldattributs
       endif
@@ -559 +586 @@
       if keyword_set(overplot) then GOTO, fini
 ;------------------------------------------------------------
-; Trace de la ligne de changement de date
-;------------------------------------------------------------
-      if (type eq 'xt') then $ 
-       plot, [180,180],[tempsmin,tempsmax],/noerase,color=c_cote
+; Trace de la ligne de changement de date or the equator
+;------------------------------------------------------------
+      IF key_onearth THEN BEGIN
+        CASE type OF
+          'xt':oplot, [180, 180] - 360*(!x.range[1] LT 180), !y.range 
+          'yt':oplot, !x.range, [0, 0]
+          ELSE:
+        ENDCASE
+      ENDIF
 ;------------------------------------------------------------
 ; legendes + affichage de celles-ci
@@ -569 +601 @@
        , ENDPOINTS = endpoints, _extra = ex
 ;
-      plot,[0],[0],/noerase,color=c_cote, xstyle = 1+4*(keyword_set(endpoints) AND type EQ 'xt' AND lat1 NE lat2) $
-       , ystyle = 1+4*(keyword_set(endpoints) AND type EQ 'yt'), _extra = ex
+; we want to draw the time axis with tick length = 1. if we simply
+; use [xy]style=1 then each axis (up and down or left and right) have
+; majorticks with a length of 1 (-> not very nice when xgridstyle=2),
+; same if xticklen=0.5 (not very nice in the middle). => so we draw
+; the top (right) axis by hand using axis.
+;
+      if n_elements(ex) NE 0 then BEGIN 
+; pour avoir un cadre de la couleur noire
+        if (where(tag_names(ex) EQ 'COLOR'))[0] NE -1 then ex.COLOR = 0
+      ENDIF
+;
+      plot, [0], [0], /nodata, /noerase, _extra = ex $
+        , xstyle = 1+4*(keyword_set(endpoints) AND type EQ 'xt' AND lat1 NE lat2)+8*(type EQ 'yt' OR type EQ 'zt') $
+        , ystyle = 1+4*(keyword_set(endpoints) AND type EQ 'yt')+8*(type EQ 'xt')
+; call axis for the missing axis.
+      IF type EQ 'xt' AND NOT keyword_set(endpoints) THEN BEGIN
+        if n_elements(ex) NE 0 then $ 
+          if (where(tag_names(ex) EQ 'YTICKNAME'))[0] NE -1 then $
+          ex.YTICKNAME = replicate(' ', n_elements(ex.YTICKNAME))
+        axis, yaxis = 1, ystyle = 1, yticklen = 0 $
+          , ytickname = replicate(' ', !y.ticks+1), _extra = ex
+      ENDIF
+      IF (type EQ 'yt' OR type EQ 'zt') AND NOT keyword_set(endpoints) THEN BEGIN
+        if n_elements(ex) NE 0 then $ 
+          if (where(tag_names(ex) EQ 'XTICKNAME'))[0] NE -1 then $
+          ex.XTICKNAME = replicate(' ', n_elements(ex.XTICKNAME))
+       axis, xaxis = 1, xstyle = 1, xticklen = 0 $
+        , xtickname = replicate(' ', !x.ticks+1), _extra = ex
+      ENDIF
 ; ajout d'un axe ds le cas ou l'on utilise endpoints
       if keyword_set(endpoints) then addaxe, endpoints, type, posfenetre, _EXTRA = ex
@@ -576 +635 @@
 ; barre de couleur
 ;------------------------------------------------------------
-      couleur = couleur[0:ncontour-1-keyword_set(strictfill)]
-      barrecouleur, position=posbar,divisions=(ncontour-keyword_set(strictfill))/2,min=min $
-       , max=max,discret=couleur,_extra = ex
+      colnumb = colnumb[0:ncontour-1-keyword_set(strictfill)]
+      barrecouleur, colnumb, min,  max, (ncontour-keyword_set(strictfill))/2 $
+                    , position = posbar, _extra = ex
 ;------------------------------------------------------------
    endif
@@ -592 +651 @@
             bar_plot, yy, background = (!d.n_colors-1) < 255,  $
              baselines = replicate(!y.range[0], n_elements(yy)), barnames = ['', ''], $
-             colors    = replicate(col1d      , n_elements(yy)), /outline
-            plot, [0], [0],/noerase,color=c_cote, _extra = ex
+             colors    = replicate(col1d, n_elements(yy)), /outline
+            if n_elements(ex) NE 0 then BEGIN 
+; pour avoir un cadre de la couleur noire
+              if (where(tag_names(ex) EQ 'COLOR'))[0] NE -1 then ex.COLOR = 0
+            ENDIF
+            plot, [0], [0], /nodata, /noerase, _extra = ex
             GOTO, fini
          ENDIF 
@@ -609 +672 @@
             yy = rien
          endif
-         plot,xx,yy,color=col1d,linestyle=sty1d,thick=2, title = '', subtitle = '', _extra = ex
+         plot, xx, yy, color = col1d, linestyle = sty1d, thick = 2 $
+           , title = '', subtitle = '', _extra = ex
+; draw 0 line
+         if keyword_set(exchange_xy) then $
+           oplot, [0, 0], !y.range ELSE oplot, !x.range, [0, 0]
+; we want to draw the axis with the time axis length = 1. if we simply
+; use xstyle=1 then each axis (up and down) have majorticks with a
+; length of 1 (-> not very nice when xgridstyle=2), same if
+; xticklen=0.5 (not very nice in the middle). 
+; => so we draw the top axis by hand using axis.
          if n_elements(ex) NE 0 then BEGIN 
-; pour avoir un cadre de la couleur c_cote
-            if (where(tag_names(ex) EQ 'COLOR'))[0] NE -1 then ex.COLOR = c_cote
+; pour avoir un cadre de la couleur noire
+            if (where(tag_names(ex) EQ 'COLOR'))[0] NE -1 then ex.COLOR = 0
          ENDIF
-         if keyword_set(exchange_xy) then $
-          plot, [0, 0],!y.range ,/noerase,color=c_cote, xstyle = 1, ystyle = 1, _extra = ex $
-         ELSE plot, !x.range, [0, 0],/noerase,color=c_cote, xstyle = 1, ystyle = 1, _extra = ex
-      ENDIF ELSE oplot,xx,yy,color=col1d,linestyle=sty1d,thick=2, _extra = ex
+         plot, [0], [0], /nodata, /noerase $
+           , xstyle = 1+8*(1-keyword_set(exchange_xy)) $
+           , ystyle = 1+8*keyword_set(exchange_xy), _extra = ex
+; call axis for the missing axis.
+         if n_elements(ex) NE 0 then BEGIN 
+; force tickname to blank array
+            if (where(tag_names(ex) EQ 'YTICKNAME'))[0] NE -1 AND keyword_set(exchange_xy) then ex.YTICKNAME = replicate(' ', n_elements(ex.YTICKNAME))
+            if (where(tag_names(ex) EQ 'XTICKNAME'))[0] NE -1 AND NOT keyword_set(exchange_xy) then ex.XTICKNAME = replicate(' ', n_elements(ex.XTICKNAME))
+         ENDIF
+         if keyword_set(exchange_xy) then  axis, yaxis = 1, ystyle = 1 $
+           , yticklen = 0, ytickname = replicate(' ', !y.ticks+1), _extra = ex $
+         ELSE axis, xaxis = 1, xstyle = 1, xticklen = 0 $
+           , xtickname = replicate(' ', !x.ticks+1), _extra = ex
+       ENDIF ELSE oplot, xx, yy, color = col1d, linestyle = sty1d, thick = 2, _extra = ex
    endif
 fini:
@@ -633 +715 @@
 ; on fait un faut plot pour que ces valeurs soient prises en
 ; consideration
-   plot,[0], [0], /noerase,xstyle=5, ystyle = 5, title = '', subtitle = '', ytitle = '', xtitle = ''
+   plot, [0], [0], /nodata, /noerase, xstyle = 5, ystyle = 5 $
+     , title = '', subtitle = '', ytitle = '', xtitle = ''
 ;------------------------------------------------------------
 ;------------------------------------------------------------
@@ -640 +723 @@
 ;------------------------------------------------------------
    terminedessin, _extra = ex
-   if keyword_set(oldboite) then domdef, oldboite,GRILLE=[vargrid, 'T']
+   if keyword_set(savedbox) THEN restoreboxparam, 'boxparam4pltt.dat'
 ;------------------------------------------------------------
    if n_elements(key_performance) NE 0 then $