Ignore:
Timestamp:
06/19/06 16:14:56 (18 years ago)
Author:
smasson
Message:

new compilation options (compile_opt idl2, strictarrsubs) in each routine

Location:
trunk/SRC/ToBeReviewed/PLOTS
Files:
41 edited

Legend:

Unmodified
Added
Removed

  • trunk/SRC/ToBeReviewed/PLOTS/DESSINE/bar_plot.pro

    r97 r114  
    131131          outline=outline,overplot=overplot,background=background, $ 
    132132          rotate=rotate, _EXTRA = ex 
     133; 
     134  compile_opt idl2, strictarrsubs 
     135; 
    133136if (n_params(d) eq 0) then begin  ;Print call & return if no parameters 
    134137  print,'bar_test,values,baselines=baselines,colors=colors,barnames=barnames,$' 

  • trunk/SRC/ToBeReviewed/PLOTS/DESSINE/plt.pro

    r67 r114  
    304304;--------------------------------------------------------- 
    305305; include common 
     306; 
     307  compile_opt idl2, strictarrsubs 
     308; 
    306309@cm_4mesh 
    307310@cm_4data 

  • trunk/SRC/ToBeReviewed/PLOTS/DESSINE/plt1d.pro

    r74 r114  
    160160;--------------------------------------------------------- 
    161161; include common 
     162; 
     163  compile_opt idl2, strictarrsubs 
     164; 
    162165@cm_4mesh 
    163166@cm_4data 
    … …  
    239242    notanum = where(finite(z1d) EQ 0) 
    240243    z1d[notanum] = 0 
    241     mask(where(z1d LT valmask/10)) = 1 
     244    mask[where(z1d LT valmask/10)] = 1 
    242245    z1d[notanum] = !values.f_nan 
    243   ENDIF ELSE mask(where(z1d LT valmask/10)) = 1 
     246  ENDIF ELSE mask[where(z1d LT valmask/10)] = 1 
    244247;----------------------------------------------------------------------------- 
    245248; determination du min et du max apres la moyenne 

  • trunk/SRC/ToBeReviewed/PLOTS/DESSINE/pltbase.pro

    r67 r114  
    135135             , _EXTRA = ex 
    136136;--------------------------------------------------------- 
     137; 
     138  compile_opt idl2, strictarrsubs 
     139; 
    137140@cm_4mesh 
    138141  IF NOT keyword_set(key_forgetold) THEN BEGIN 
    … …  
    215218; ds le cas on unsur2 est active on reduit levels 
    216219    if NOT keyword_set(nocontour) then begin 
    217       IF keyword_set(unsur2) THEN levels = levels(where(zeroun(n_elements(levels) ) eq 1)) 
     220      IF keyword_set(unsur2) THEN levels = levels[where(zeroun(n_elements(levels) ) eq 1)] 
    218221; unlabsur est active?  C_LABEL est passe via _EXTRA? 
    219222      if keyword_set(unlabsur) THEN IF chkstru(ex, 'C_LABELS') THEN $ 

  • trunk/SRC/ToBeReviewed/PLOTS/DESSINE/pltsc.pro

    r35 r114  
    33; 
    44; scatter plot (inspired from plt1d) 
     5; 
     6; 
     7  compile_opt idl2, strictarrsubs 
    58; 
    69 

  • trunk/SRC/ToBeReviewed/PLOTS/DESSINE/pltt.pro

    r74 r114  
    306306;--------------------------------------------------------- 
    307307; include common 
     308; 
     309  compile_opt idl2, strictarrsubs 
     310; 
    308311@cm_4mesh 
    309312@cm_4data 

  • trunk/SRC/ToBeReviewed/PLOTS/DESSINE/pltz.pro

    r35 r114  
    266266;--------------------------------------------------------- 
    267267; include common 
     268; 
     269  compile_opt idl2, strictarrsubs 
     270; 
    268271@cm_4mesh 
    269272@cm_4data 

  • trunk/SRC/ToBeReviewed/PLOTS/DESSINE/sbar_plot.pro

    r69 r114  
    5151 
    5252PRO sbar_plot, Values, COLORS = colors, NOREINITPLT = noreinitplt, _extra = ex  
     53; 
     54  compile_opt idl2, strictarrsubs 
     55; 
    5356@common 
    5457; 1) je reinitialise l''environnememt graphique (les variables !x, !y et !p):  

  • trunk/SRC/ToBeReviewed/PLOTS/DESSINE/scontour.pro

    r35 r114  
    4141;------------------------------------------------------------ 
    4242PRO scontour, x, y, z, NOREINITPLT = noreinitplt, _EXTRA = ex 
     43; 
     44  compile_opt idl2, strictarrsubs 
     45; 
    4346@common 
    4447; 1) je reinitialise l''environnememt graphique (les variables !x, !y et !p):  

  • trunk/SRC/ToBeReviewed/PLOTS/DESSINE/splot.pro

    r67 r114  
    4646;------------------------------------------------------------ 
    4747PRO splot,  x, y, NOREINITPLT = noreinitplt, _EXTRA = ex 
     48; 
     49  compile_opt idl2, strictarrsubs 
     50; 
    4851@common 
    4952; 1) je reinitialise l''environnememt graphique (les variables !x, !y et !p):  

  • trunk/SRC/ToBeReviewed/PLOTS/DESSINE/tvplus.pro

    r69 r114  
    8686            , MIN = min, MAX = max, MASK = mask, OFFSET = offset, NOUSEINFOS = NOUSEINFOS $ 
    8787            , NCOLORS = ncolors, NOINTERP = nointerp, _EXTRA = ex 
     88; 
     89; 
     90  compile_opt idl2, strictarrsubs 
    8891; 
    8992  IF n_elements(z2d) EQ 0 THEN return 
    … …  
    136139    if abs(mask) LT 1e6 then BEGIN  
    137140      masked = where(arr EQ mask) 
    138       if masked[0] NE -1 then arr(masked) = min(arr(where(arr NE mask))) 
     141      if masked[0] NE -1 then arr[masked] = min(arr[where(arr NE mask)]) 
    139142    ENDIF ELSE BEGIN  
    140143      masked = where(abs(arr) GE abs(mask)/10.) 
    141       if masked[0] NE -1 then arr(masked) = min(arr(where(abs(arr) LT abs(mask)/10.))) 
     144      if masked[0] NE -1 then arr[masked] = min(arr[where(abs(arr) LT abs(mask)/10.)]) 
    142145    ENDELSE 
    143146  ENDIF ELSE masked = -1 

  • trunk/SRC/ToBeReviewed/PLOTS/DIVERS/addaxe.pro

    r37 r114  
    3737;------------------------------------------------------------ 
    3838PRO addaxe, endpoints, type, posfenetre, _EXTRA = ex 
     39; 
     40  compile_opt idl2, strictarrsubs 
     41; 
    3942@common 
    4043;--------------------------------------- 

  • trunk/SRC/ToBeReviewed/PLOTS/DIVERS/autoscale.pro

    r37 r114  
    4242; ce CI est un multiple de l'unite en unite log de 10. 
    4343; 
     44; 
     45  compile_opt idl2, strictarrsubs 
     46; 
    4447    ci = (max-min)/20. 
    4548    ci = 10.^floor(alog10(ci)) 

  • trunk/SRC/ToBeReviewed/PLOTS/DIVERS/axis4pltz.pro

    r74 r114  
    4141; 
    4242; include common 
     43; 
     44  compile_opt idl2, strictarrsubs 
     45; 
    4346@cm_4mesh 
    4447  IF NOT keyword_set(key_forgetold) THEN BEGIN 

  • trunk/SRC/ToBeReviewed/PLOTS/DIVERS/barrecouleur.pro

    r37 r114  
    2525                  , CB_SUBTITLE = cb_subtitle, POST = post, _extra = ex 
    2626;------------------------------------------------------------ 
     27; 
     28  compile_opt idl2, strictarrsubs 
     29; 
    2730@cm_general 
    2831  IF NOT keyword_set(key_forgetold) THEN BEGIN 

  • trunk/SRC/ToBeReviewed/PLOTS/DIVERS/checkfield.pro

    r69 r114  
    4040;------------------------------------------------------------ 
    4141FUNCTION err_1d, type, n1, name, n2 
     42; 
     43  compile_opt idl2, strictarrsubs 
     44; 
    4245  return, report(['Error in "' + type + '" type plot with a 1D input array:' $ 
    4346                  , 'the number of elements of the input vector ('+strtrim(n1, 1)+') ' $ 
    … …  
    4649; 
    4750FUNCTION err_2d, type, sz, nx, ny, nz 
     51; 
     52  compile_opt idl2, strictarrsubs 
     53; 
    4854  @cm_4mesh 
    4955  @cm_4cal 
    … …  
    5965; 
    6066FUNCTION err_3d, type, sz, nx, ny, nz 
     67; 
     68  compile_opt idl2, strictarrsubs 
     69; 
    6170  @cm_4mesh 
    6271  @cm_4cal 
    … …  
    7483FUNCTION checkfield, field, procedure, TYPE = type, BOXZOOM = boxzoom, DIREC = direc, NOQUESTION = noquestion, VECTEUR = vecteur, WDEPTH = wdepth, _EXTRA = ex 
    7584;-------------------------------------------------------------- 
     85; 
     86  compile_opt idl2, strictarrsubs 
     87; 
    7688; include commons 
    7789@cm_4mesh 

  • trunk/SRC/ToBeReviewed/PLOTS/DIVERS/checktypeminmax.pro

    r37 r114  
    3636                     , XINDEX = xindex, YINDEX = yindex $ 
    3737                     , ENDPOINTS = endpoints, _extra = ex 
     38; 
     39; 
     40  compile_opt idl2, strictarrsubs 
    3841; 
    3942@common 

  • trunk/SRC/ToBeReviewed/PLOTS/DIVERS/determineminmax.pro

    r37 r114  
    4444;------------------------------------------------------------ 
    4545PRO determineminmax, tab, mask, vraimin, vraimax, glam, gphi, MAXIN = maxin, MININ = minin, INTERVALLE = intervalle, usetri = usetri, ZEROMIDDLE = zeromiddle, _extra = ex 
     46; 
     47  compile_opt idl2, strictarrsubs 
     48; 
    4649@common 
    4750;----------------------------------------------------------------------------- 
    … …  
    7174  endif 
    7275; ma et mi : max et min sur les points mer 
    73   vraimax = max(tab(mer), min = vraimin, _extra = ex) 
     76  vraimax = max(tab[mer], min = vraimin, _extra = ex) 
    7477  sameminmax = testvar(var = minin) EQ testvar(var = maxin)  
    7578  if n_elements(maxin) EQ 0 OR sameminmax then maxin = vraimax 

  • trunk/SRC/ToBeReviewed/PLOTS/DIVERS/givewindowsize.pro

    r37 r114  
    33;------------------------------------------------------------ 
    44; include commons 
     5; 
     6  compile_opt idl2, strictarrsubs 
     7; 
    58@cm_4ps 
    69IF NOT keyword_set(key_forgetold) THEN BEGIN 

  • trunk/SRC/ToBeReviewed/PLOTS/DIVERS/meridienparallele.pro

    r37 r114  
    2323;------------------------------------------------------------ 
    2424PRO meridienparallele, coupe 
     25; 
     26  compile_opt idl2, strictarrsubs 
     27; 
    2528@common 
    2629;------------------------------------------------------------ 

  • trunk/SRC/ToBeReviewed/PLOTS/DIVERS/placecolor.pro

    r37 r114  
    4242;------------------------------------------------------------ 
    4343; include commons 
     44; 
     45  compile_opt idl2, strictarrsubs 
     46; 
    4447@cm_4ps 
    4548   IF NOT keyword_set(key_forgetold) THEN BEGIN 

  • trunk/SRC/ToBeReviewed/PLOTS/DIVERS/placedessin.pro

    r69 r114  
    5656;--------------------------------------------------------- 
    5757; include common 
     58; 
     59  compile_opt idl2, strictarrsubs 
     60; 
    5861@cm_4ps 
    5962@cm_4mesh 

  • trunk/SRC/ToBeReviewed/PLOTS/DIVERS/projsegment.pro

    r37 r114  
    5050FUNCTION projsegment, vecteur, bornes, MP = mp 
    5151;-------------------------------------------------------------- 
     52; 
     53  compile_opt idl2, strictarrsubs 
     54; 
    5255   a1 = float(vecteur[0]) 
    5356   b1 = float(vecteur[n_elements(vecteur)-1]) 

  • trunk/SRC/ToBeReviewed/PLOTS/DIVERS/restoreatt.pro

    r37 r114  
    2929;------------------------------------------------------------ 
    3030PRO restoreatt, struct 
     31; 
     32  compile_opt idl2, strictarrsubs 
     33; 
    3134@common 
    3235;------------------------------------------------------------ 

  • trunk/SRC/ToBeReviewed/PLOTS/DIVERS/rotation.pro

    r37 r114  
    2828 
    2929PRO ROTATION,X,Y,DEG,NX,NY 
     30; 
     31  compile_opt idl2, strictarrsubs 
     32; 
    3033ang=deg*!dtor 
    3134 
    … …  
    3538;get angle in for loop so that zero radii will be left as zero angle 
    3639for i = 0,n_elements(r)-1 do $ 
    37 if r(i) ne 0 then theta(i) = atan(y(i),x(i))  ;range from -pi to +pi 
     40if r[i] ne 0 then theta[i] = atan(y[i],x[i])  ;range from -pi to +pi 
    3841; 
    3942;add rotation angle 

  • trunk/SRC/ToBeReviewed/PLOTS/DIVERS/saveatt.pro

    r37 r114  
    2626;------------------------------------------------------------ 
    2727FUNCTION saveatt 
     28; 
     29  compile_opt idl2, strictarrsubs 
     30; 
    2831@common 
    2932   return, {n:varname,g:vargrid,d:vardate,e:varexp,u:varunit,m:valmask} 

  • trunk/SRC/ToBeReviewed/PLOTS/DIVERS/terminedessin.pro

    r37 r114  
    3434PRO terminedessin, POST = post, SMALL = small, _extra = ex 
    3535;--------------------------------------------------------- 
     36; 
     37  compile_opt idl2, strictarrsubs 
     38; 
    3639@cm_4ps 
    3740  IF NOT keyword_set(key_forgetold) THEN BEGIN 

  • trunk/SRC/ToBeReviewed/PLOTS/LABEL/label.pro

    r39 r114  
    5050pro label, cas, min, max, ncontour, level_z2d, colnumb, NLEVEL = nlevel $ 
    5151          ,INTERVALLE=intervalle, STRICTFILL = strictfill 
     52; 
     53  compile_opt idl2, strictarrsubs 
     54; 
    5255@common 
    5356   if !d.name EQ 'PS' OR !d.name EQ 'Z' then BEGIN 
    … …  
    8386         level_z2d = min + intervalle*findgen(Ncontour)  
    8487         colnumb   = ncoul*(findgen(Ncontour))/Ncontour+ncoul/(2*ncontour) 
    85          max=level_z2d(Ncontour-1)+intervalle 
     88         max=level_z2d[Ncontour-1]+intervalle 
    8689      end 
    8790; label pour faire les memes sss que dessier 

  • trunk/SRC/ToBeReviewed/PLOTS/LABEL/label_date.pro

    r97 r114  
    9999FUNCTION LABEL_DATE, axis, index, x, DATE_FORMAT = format, MONTHS = months, $ 
    100100              OFFSET= offs, _EXTRA = ex 
     101; 
     102  compile_opt idl2, strictarrsubs 
     103; 
    101104COMMON label_date_com, fmt, month_chr, offset 
    102105 

  • trunk/SRC/ToBeReviewed/PLOTS/LABEL/label_gmt.pro

    r39 r114  
    33; 
    44PRO label_gmt, min, max, intervalle, ncoul, ncontour, level_z2d, coul 
     5; 
     6  compile_opt idl2, strictarrsubs 
     7; 
    58@common 
    69@com_eg 
    … …  
    2225      ncontour  = fix((max-min)/intervalle) 
    2326      level_z2d = min + intervalle*findgen(Ncontour)  
    24       max=level_z2d(Ncontour-1)+intervalle 
     27      max=level_z2d[Ncontour-1]+intervalle 
    2528 
    2629      print, '     Number of contour intervals, plotting min & max ', ncontour, min, max 
    … …  
    4043           ; difference plot : lighter below first negative interval 
    4144            red[51:98] = long((100.-float(grey_shade_2))/100.*255.) 
    42             red(1:48) = long((100.-float(grey_shade))/100.*255.) 
     45            red[1:48] = long((100.-float(grey_shade))/100.*255.) 
    4346            red[50] = 255 
    4447 

  • trunk/SRC/ToBeReviewed/PLOTS/LABEL/lataxe.pro

    r39 r114  
    3838;------------------------------------------------------------ 
    3939; on ramenne value ds le segment [0,180] 
     40; 
     41  compile_opt idl2, strictarrsubs 
     42; 
    4043   lat=value mod 360 
    4144   if lat lt 0 then lat=lat+360 

  • trunk/SRC/ToBeReviewed/PLOTS/LABEL/lonaxe.pro

    r39 r114  
    3838;------------------------------------------------------------ 
    3939; on ramenne value ds le segment [0,360[ 
     40; 
     41  compile_opt idl2, strictarrsubs 
     42; 
    4043   lon=value mod 360 
    4144   if lon lt 0 then lon=lon+360 

  • trunk/SRC/ToBeReviewed/PLOTS/VECTEUR/ajoutvect.pro

    r41 r114  
    5454;------------------------------------------------------------ 
    5555pro ajoutvect,vecteur, vectlegende, UNVECTSUR=unvectsur,VECTMIN=vectmin, VECTMAX=vectmax, _EXTRA = ex 
     56; 
     57  compile_opt idl2, strictarrsubs 
     58; 
    5659@common 
    5760   tempsun = systime(1)         ; pour key_performance 
    … …  
    128131; construction de u et v aux pts T 
    129132;----------------------------------------------------------- 
    130       a=u(0,*) 
     133      a=u[0,*] 
    131134      u=(u+shift(u,1,0))/2. 
    132       if NOT keyword_set(key_periodic) OR nx NE jpi then u(0,*)=a 
    133       a=v(*,0) 
     135      if NOT keyword_set(key_periodic) OR nx NE jpi then u[0,*]=a 
     136      a=v[*,0] 
    134137      v=(v+shift(v,0,1))/2. 
    135       if NOT keyword_set(key_periodic) OR nx NE jpi then v(*,0)=a 
     138      if NOT keyword_set(key_periodic) OR nx NE jpi then v[*,0]=a 
    136139;---------------------------------------------------------------------------- 
    137140; attribution du mask et des tableau de longitude et latitude 
    … …  
    191194   if interpolle then t2 = msku*shift(msku,1,0)*mskv*shift(mskv,0,1) $ 
    192195   ELSE t2 = tmask[firstxt:lastxt,firstyt:lastyt,firstzt] 
    193    if NOT keyword_set(key_periodic) OR nx NE jpi then t2(0, *)=0. 
    194    t2(*,0)=0. 
     196   if NOT keyword_set(key_periodic) OR nx NE jpi then t2[0, *]=0. 
     197   t2[*,0]=0. 
    195198   terre=where(t2 eq 0) 
    196199   if terre[0] ne -1 then begin 
    197       u(terre)=1e5 
    198       v(terre)=1e5 
     200      u[terre]=1e5 
     201      v[terre]=1e5 
    199202   ENDIF 
    200203;----------------------------------------------------------- 

  • trunk/SRC/ToBeReviewed/PLOTS/VECTEUR/vecteur.pro

    r41 r114  
    8181; angle peut etre un tableau. 
    8282; 
     83; 
     84  compile_opt idl2, strictarrsubs 
     85; 
    8386@common 
    8487; quelle est la longeur en coordonnees normales d''un trait qui fera 1 
    … …  
    103106; normalise le vecteur 
    104107; 
     108; 
     109  compile_opt idl2, strictarrsubs 
     110; 
    105111   IF n_elements(w) NE 0 THEN BEGIN  
    106112      norme = sqrt(u^2.+v^2.+w^2.) 
    107113      ind = where(norme NE 0) 
    108       u(ind) = u(ind)/norme[ind] 
    109       v(ind) = v(ind)/norme[ind] 
    110       w(ind) = w(ind)/norme[ind] 
     114      u[ind] = u[ind]/norme[ind] 
     115      v[ind] = v[ind]/norme[ind] 
     116      w[ind] = w[ind]/norme[ind] 
    111117   ENDIF ELSE BEGIN 
    112118      norme = sqrt(u^2.+v^2.) 
    113119      ind = where(norme NE 0) 
    114       u(ind) = u(ind)/norme[ind] 
    115       v(ind) = v(ind)/norme[ind] 
     120      u[ind] = u[ind]/norme[ind] 
     121      v[ind] = v[ind]/norme[ind] 
    116122   ENDELSE  
    117123END 
    … …  
    120126             , VECTCOLOR = vectcolor, VECTTHICK = vectthick, VECTREFPOS = vectrefpos $ 
    121127             , VECTREFFORMAT = vectrefformat, NOVECTREF = novectref, _extra = extra 
     128; 
     129  compile_opt idl2, strictarrsubs 
     130; 
    122131@common 
    123132   tempsun = systime(1)         ; pour key_performance 
    … …  
    202211   r = cv_coord(from_sphere=coord_sphe,/to_rect,/degrees) 
    203212; 
    204    x0 = reform(r(0, *), nx, ny) 
    205    y0 = reform(r(1, *), nx, ny) 
    206    z0 = reform(r(2, *), nx, ny) 
     213   x0 = reform(r[0, *], nx, ny) 
     214   y0 = reform(r[1, *], nx, ny) 
     215   z0 = reform(r[2, *], nx, ny) 
    207216; 
    208217; etape 1, b) 
    209218; 
    210219; Construction du vecteur nu (resp. nv), vecteur norme porte par 
    211 ; l''axe des points u(i,j) et u(i-1,j) (resp v(i,j) et v(i,j-1)) 
     220; l''axe des points u[i,j] et u[i-1,j] (resp v[i,j] et v[i,j-1]) 
    212221; qui definissent pour chaque point sur la shere les directions locales 
    213222; associee a u et v. ces vecteurs definissent un repere orthonorme 
    … …  
    219228   radius = replicate(1,nxgd*nygd) 
    220229   IF finite(glamu[0]*gphiu[0]) NE 0 THEN $ 
    221      coord_sphe = transpose([ [(glamu[indice2d])[*]], [(gphiu[indice2d])[*]], [radius(*)] ]) $ 
    222    ELSE coord_sphe = transpose([ [(glamf[indice2d])[*]], [(gphit[indice2d])[*]], [radius(*)] ]) 
     230     coord_sphe = transpose([ [(glamu[indice2d])[*]], [(gphiu[indice2d])[*]], [radius[*]] ]) $ 
     231   ELSE coord_sphe = transpose([ [(glamf[indice2d])[*]], [(gphit[indice2d])[*]], [radius[*]] ]) 
    223232   r = cv_coord(from_sphere=coord_sphe,/to_rect,/degrees) 
    224233; coordonnes de points de la grille u en cartesien 
    225    ux = reform(r(0, *), nxgd, nygd) 
    226    uy = reform(r(1, *), nxgd, nygd) 
    227    uz = reform(r(2, *), nxgd, nygd) 
     234   ux = reform(r[0, *], nxgd, nygd) 
     235   uy = reform(r[1, *], nxgd, nygd) 
     236   uz = reform(r[2, *], nxgd, nygd) 
    228237; calcul de nu  
    229238   nux = ux-shift(ux, 1, 0) 
    … …  
    242251; definition de nv 
    243252   IF finite(glamv[0]*gphiv[0]) NE 0 THEN $ 
    244    coord_sphe = transpose([ [(glamv[indice2d])[*]], [(gphiv[indice2d])[*]], [radius(*)] ]) $ 
    245    ELSE coord_sphe = transpose([ [(glamt[indice2d])[*]], [(gphif[indice2d])[*]], [radius(*)] ])                 
     253   coord_sphe = transpose([ [(glamv[indice2d])[*]], [(gphiv[indice2d])[*]], [radius[*]] ]) $ 
     254   ELSE coord_sphe = transpose([ [(glamt[indice2d])[*]], [(gphif[indice2d])[*]], [radius[*]] ])                 
    246255   r = cv_coord(from_sphere=coord_sphe,/to_rect,/degrees) 
    247256; coordonnes de points de la grille v en cartesien 
    248    vx = reform(r(0, *), nxgd, nygd) 
    249    vy = reform(r(1, *), nxgd, nygd) 
    250    vz = reform(r(2, *), nxgd, nygd) 
     257   vx = reform(r[0, *], nxgd, nygd) 
     258   vy = reform(r[1, *], nxgd, nygd) 
     259   vz = reform(r[2, *], nxgd, nygd) 
    251260; calcul de nv  
    252261   nvx = vx-shift(vx, 0, 1) 
    … …  
    290299; coordonnees de la pointe en spherique 
    291300 
    292    coord_rect = transpose([ [x1(*)], [y1(*)], [z1(*)] ]) 
     301   coord_rect = transpose([ [x1[*]], [y1[*]], [z1[*]] ]) 
    293302   r = cv_coord(from_rect=coord_rect,/to_sphere,/degrees) 
    294    glam1 = reform(r(0, *), nx, ny) 
    295    gphi1 = reform(r(1, *), nx, ny) 
     303   glam1 = reform(r[0, *], nx, ny) 
     304   gphi1 = reform(r[1, *], nx, ny) 
    296305 
    297306; 
    … …  
    303312; 
    304313   ind = where(glam1 LT !x.range[0] AND glam1+360. LE !x.range[1]) 
    305    if ind[0] NE -1 then glam1(ind) = glam1(ind)+360. 
     314   if ind[0] NE -1 then glam1[ind] = glam1[ind]+360. 
    306315   ind = where(glam1 GT !x.range[1] AND glam1-360. GE !x.range[0]) 
    307    if ind[0] NE -1 then glam1(ind) = glam1(ind)-360. 
     316   if ind[0] NE -1 then glam1[ind] = glam1[ind]-360. 
    308317 
    309318   ind = where(glam LT !x.range[0] AND glam+360. LE !x.range[1]) 
    310    if ind[0] NE -1 then glam(ind) = glam(ind)+360. 
     319   if ind[0] NE -1 then glam[ind] = glam[ind]+360. 
    311320   ind = where(glam  GT !x.range[1] AND glam-360. GE !x.range[0]) 
    312    if ind[0] NE -1 then glam(ind) = glam(ind)-360. 
     321   if ind[0] NE -1 then glam[ind] = glam[ind]-360. 
    313322; 
    314323; 
    … …  
    316325; 
    317326   r = convert_coord(glam,gphi,/data,/to_normal)  
    318    x0 = r(0, *)                 ; coordonnes normales du debut de la fleche 
    319    y0 = r(1, *)                 ;  
     327   x0 = r[0, *]                 ; coordonnes normales du debut de la fleche 
     328   y0 = r[1, *]                 ;  
    320329    
    321330   r = convert_coord(glam1,gphi1,/data,/to_normal)  
    322    x1 = r(0, *)                 ; coordonnes normales de la fin de la fleche (avant scaling) 
    323    y1 = r(1, *)                 ;  
     331   x1 = r[0, *]                 ; coordonnes normales de la fin de la fleche (avant scaling) 
     332   y1 = r[1, *]                 ;  
    324333; 
    325334; tests pour eviter que des fleches soient dessineees hors du domaine 
    … …  
    384393; 
    385394   r = cv_coord(from_polar = transpose([ [dirpol[*]], [norme[*]] ]), /to_rect) 
    386    composantex = r(0, *) 
    387    composantey = r(1, *) 
     395   composantex = r[0, *] 
     396   composantey = r[1, *] 
    388397; 
    389398   x1 = x0+composantex 
    … …  
    395404 
    396405   points = where(msk EQ 1) 
    397    IF points[0] NE -1 THEN arrow, x0(points), y0(points), x1(points), y1(points), /norm $ 
     406   IF points[0] NE -1 THEN arrow, x0[points], y0[points], x1[points], y1[points], /norm $ 
    398407    , hsize = -.2, COLOR = vectcolor, THICK = vectthick 
    399408; 

  • trunk/SRC/ToBeReviewed/PLOTS/VECTEUR/velovect.pro

    r97 r114  
    100100PRO VELOVECT,U,V,X,Y, Missing = Missing, Length = length, Dots = dots,  $ 
    101101        Color=color, CLIP=clip, NOCLIP=noclip, OVERPLOT=overplot, _EXTRA=extra 
     102; 
     103  compile_opt idl2, strictarrsubs 
     104; 
    102105        on_error,2                      ;Return to caller if an error occurs 
    103106        s = size(u) 

  • trunk/SRC/ToBeReviewed/PLOTS/axe.pro

    r42 r114  
    5050;------------------------------------------------------------ 
    5151PRO axe, coupe,tempsmin,tempsmax, REVERSE_X = reverse_x, REVERSE_Y = reverse_y, SIN = sin, SEPDATE = sepdate, DIGITSYEAR = digitsyear, _EXTRA = ex 
     52; 
     53  compile_opt idl2, strictarrsubs 
     54; 
    5255@common 
    5356   tempsun = systime(1)         ; pour key_performance 

  • trunk/SRC/ToBeReviewed/PLOTS/legende.pro

    r42 r114  
    4747            , INTERVALLE = intervalle, TYPE_YZ = type_yz, VARNAME2 = varname2 $ 
    4848            , NPTS = npts, _EXTRA = ex 
     49; 
     50  compile_opt idl2, strictarrsubs 
     51; 
    4952@common 
    5053   tempsun = systime(1)         ; pour key_performance 

  • trunk/SRC/ToBeReviewed/PLOTS/plotsym.pro

    r97 r114  
    3232                  angle=angle, box=box, line=line, scale=scale, $ 
    3333                  _extra=extra 
     34; 
     35  compile_opt idl2, strictarrsubs 
     36; 
    3437if not keyword_set(scale) then scale=1.0 
    3538if not keyword_set(angle) then angle=0.0 

  • trunk/SRC/ToBeReviewed/PLOTS/reinitplt.pro

    r42 r114  
    4141pro reinitplt, all=all,x=x,y=y,z=z,p=p, invert=invert 
    4242;------------------------------------------------------------ 
     43; 
     44  compile_opt idl2, strictarrsubs 
     45; 
    4346        clearx = 0 
    4447        cleary = 0 

  • trunk/SRC/ToBeReviewed/PLOTS/style.pro

    r42 r114  
    4444;------------------------------------------------------------ 
    4545pro style,labstyle,level_z2d,linestyle,thick 
     46; 
     47  compile_opt idl2, strictarrsubs 
     48; 
    4649   case labstyle of 
    4750      0: begin 
    … …  
    7578         a=replicate(0,n_elements(level_z2d)-n) 
    7679         c=replicate(2,n)          
    77          if seuil eq level_z2d(n) then begin 
     80         if seuil eq level_z2d[n] then begin 
    7881            thick=[replicate(1,n),2,replicate(1,n_elements(level_z2d)-1-n)] 
    7982            linestyle=[c,a] 
    … …  
    8689      3: begin 
    8790         n = n_elements(level_z2d) 
    88          seuil = level_z2d(1+n/2) 
     91         seuil = level_z2d[1+n/2] 
    8992 
    9093         thick = intarr(n) 
    91          thick(indgen((n)/4)*4) = 1 
    92          thick(indgen((n)/4)*4+1) = 1 
    93          thick(indgen((n)/4)*4+2) = 2 
    94          thick(indgen((n)/4)*4+3) = 1 
     94         thick[indgen(n/4)*4] = 1 
     95         thick[indgen(n/4)*4+1] = 1 
     96         thick[indgen(n/4)*4+2] = 2 
     97         thick[indgen(n/4)*4+3] = 1 
    9598 
    9699         linestyle = intarr(n) 
    97          linestyle(indgen((n)/4)*4) = 3 
    98          linestyle(indgen((n)/4)*4+1) = 0 
    99          linestyle(indgen((n)/4)*4+2) = 0 
    100          linestyle(indgen((n)/4)*4+3) = 0 
     100         linestyle[indgen(n/4)*4] = 3 
     101         linestyle[indgen(n/4)*4+1] = 0 
     102         linestyle[indgen(n/4)*4+2] = 0 
     103         linestyle[indgen(n/4)*4+3] = 0 
    101104 
    102105         labels = intarr(n) 
    103          labels(indgen((n)/2)*2) = 1 
    104          labels(n/2) = 0 
     106         labels[indgen(n/2)*2] = 1 
     107         labels[n/2] = 0 
    105108 
    106109         return 

  • trunk/SRC/ToBeReviewed/PLOTS/symbols.pro

    r42 r114  
    5050;- 
    5151pro symbols,nsym,scale,color=col 
     52; 
     53  compile_opt idl2, strictarrsubs 
     54; 
    5255on_error,2 
    5356fill = 0 
    … …  
    6568                      xarr = fltarr(5) 
    6669                      yarr = xarr 
    67                       xarr(1) = 10. 
    68                       xarr(2) = 6. 
    69                       yarr(2) = 2. 
     70                      xarr[1] = 10. 
     71                      xarr[2] = 6. 
     72                      yarr[2] = 2. 
    7073                      ;nsyms greater than 10 should be filled arrows 
    7174                      if nsym gt 10 then begin 
    72                          xarr(3) = 6.  
    73                          xarr(4) = 10. 
    74                          yarr(3) = -2. 
     75                         xarr[3] = 6.  
     76                         xarr[4] = 10. 
     77                         yarr[3] = -2. 
    7578                         fill = 1 
    7679                      endif else begin 
    77                          xarr(3) = 10. 
    78                          xarr(4) = 6. 
    79                          yarr(4) = -2. 
     80                         xarr[3] = 10. 
     81                         xarr[4] = 6. 
     82                         yarr[4] = -2. 
    8083                      endelse 
    8184                      case 1 of 
    … …  
    8588                            xarr = extrac(xarr,0,11) 
    8689                            yarr = extrac(yarr,0,11) 
    87                             yarr(6) = 0.5 
    88                             xarr(7) = 6 
    89                             yarr(7) = 0.5 
    90                             xarr(8) = 6 
    91                             yarr(8) = -0.5 
    92                             yarr(9) = -0.5 
     90                            yarr[6] = 0.5 
     91                            xarr[7] = 6 
     92                            yarr[7] = 0.5 
     93                            xarr[8] = 6 
     94                            yarr[8] = -0.5 
     95                            yarr[9] = -0.5 
    9396                            if nsym eq 12 then begin 
    9497                               rotation,xarr,yarr,180,nx,ny 
    … …  
    124127                      xarr = fltarr(5) + 3 
    125128                      yarr = xarr 
    126                       xarr(1) = -3. 
    127                       xarr(2) = -3. 
    128                       yarr(2) = -3. 
    129                       yarr(3) = -3. 
     129                      xarr[1] = -3. 
     130                      xarr[2] = -3. 
     131                      yarr[2] = -3. 
     132                      yarr[3] = -3. 
    130133                      if (nsym eq 21)+(nsym eq 31) then begin 
    131134                         rotation,xarr,yarr,45,nx,ny 
    … …  
    139142                      yarr = fltarr(4) - 6./4. 
    140143                      xarr = fltarr(4) - 6./2. 
    141                       xarr(1) = 6./2. 
    142                       xarr(2) = 0. 
    143                       yarr(2) = 6.*sqrt(3.)/2. - 6./4. 
     144                      xarr[1] = 6./2. 
     145                      xarr[2] = 0. 
     146                      yarr[2] = 6.*sqrt(3.)/2. - 6./4. 
    144147                      if nsym eq 32 then fill = 1 
    145148                                    end 
    … …  
    147150                      xarr = fltarr(2) + 1 
    148151                      yarr = xarr * 0. 
    149                       xarr(1) = -1. 
     152                      xarr[1] = -1. 
    150153                                    end 
    151154endcase 
Note: See TracChangeset for help on using the changeset viewer.