Changeset 231 for trunk/SRC/ToBeReviewed/PLOTS/VECTEUR/ajoutvect.pro

Timestamp:

03/19/07 18:15:51 (17 years ago)

Author:

pinsard

Message:

improvements/corrections of some *.pro headers

File:

• trunk/SRC/ToBeReviewed/PLOTS/VECTEUR/ajoutvect.pro (modified) (13 diffs)

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

@@ -1 @@
-;------------------------------------------------------------
-;------------------------------------------------------------
-;------------------------------------------------------------
 ;+
 ;
@@ -11 +8 @@
 ;
 ; @param VECTEUR {in}{required}{type=vector}
-; It is a structure with 2 elements containing we 2 matrices U and V of 
-; values of the zonal and meridian component of the field of vectors to 
+; It is a structure with 2 elements containing we 2 matrices U and V of
+; values of the zonal and meridian component of the field of vectors to
 ; be traced.
 ;      For ex:
@@ -22 +19 @@
 ; It is a scalar n or an array with 2 elements [n1,n2].
 ; In the first case, we will trace a vector on n following x and y.
-; In the second case, we will trace a vector on n1 following x and a 
+; In the second case, we will trace a vector on n1 following x and a
 ; vector n2 following n2
-;   Comments: To trace all vectors following y and one vector on two 
+;   Comments: To trace all vectors following y and one vector on two
 ; following x, put unvectsur=[2,1]
 ;
 ; @keyword VECTMIN {in}{required}
-; Minimum norme of vectors to be traced 
+; Minimum norme of vectors to be traced
 ;
 ; @keyword VECTMAX {in}{required}
-; Maximum norme of vectors to be traced 
+; Maximum norme of vectors to be traced
 ;
 ; @keyword _EXTRA
-; Used to pass your keywords
-; 
+; Used to pass keywords
+;
 ; @uses
 ; common.pro            
@@ -46 +43 @@
 ;
 ; @version
-; $Id$ 
-; 
+; $Id$
+;
 ;-
-;------------------------------------------------------------
-;------------------------------------------------------------
-;------------------------------------------------------------
-pro ajoutvect,vecteur, vectlegende, UNVECTSUR=unvectsur,VECTMIN=vectmin, VECTMAX=vectmax, _EXTRA = ex
+;
+PRO ajoutvect,vecteur, vectlegende, UNVECTSUR=unvectsur,VECTMIN=vectmin, VECTMAX=vectmax, _EXTRA = ex
 ;
   compile_opt idl2, strictarrsubs
@@ -90 +85 @@
       indiceyv = (lindgen(jpj))[firstyv:firstyv+nyv-1]
       indicey = inter(indiceyu, indiceyv)
-      nx = n_elements(indicex) 
+      nx = n_elements(indicex)
       ny = n_elements(indicey)
       indice2d = lindgen(jpi, jpj)
@@ -100 +95 @@
          (size(u))[0] NE 2 OR (size(v))[0] NE 2: return
          (size(u))[1] EQ nxu AND (size(u))[2] EQ nyu AND $
-          (size(v))[1] EQ nxv AND (size(v))[2] EQ nyv:BEGIN 
+          (size(v))[1] EQ nxv AND (size(v))[2] EQ nyv:BEGIN
             if nxu NE nx then $
-             if indicex[0] EQ firstxu then u = u[0:nx-1, *] ELSE u = u[1: nx, *] 
+             if indicex[0] EQ firstxu then u = u[0:nx-1, *] ELSE u = u[1: nx, *]
             IF nxv NE nx THEN $
              if indicex[0] EQ firstxv then v = v[0:nx-1, *] ELSE v = v[1: nx, *]
             IF nyu NE ny THEN $
-             if indicey[0] EQ firstyu then u = u[*, 0:ny-1] ELSE u = u[*, 1: ny] 
+             if indicey[0] EQ firstyu then u = u[*, 0:ny-1] ELSE u = u[*, 1: ny]
             IF nyv NE ny THEN $
              if indicey[0] EQ firstyv then v = v[*, 0:ny-1] ELSE v = v[*, 1: ny]
          END
          (size(u))[1] EQ jpi AND (size(u))[2] EQ jpj AND $
-          (size(v))[1] EQ jpi AND (size(v))[2] EQ jpj:BEGIN 
+          (size(v))[1] EQ jpi AND (size(v))[2] EQ jpj:BEGIN
             u = u[indice2d]
             v = v[indice2d]
          END
-         ELSE:BEGIN 
+         ELSE:BEGIN
             ras = report('problemes d''adequation entre la taille du domaine et la taille des matrices necessaires a tracer des vecteurs')
             return
@@ -138 +133 @@
 ;----------------------------------------------------------------------------
 ; attribution of the mask and of longitude and latitude arrays.
-; We recuperate the complete grid to establish a big mask extensive 
-; in the four directions to cover points for which a land point has 
+; We recuperate the complete grid to establish a big mask extensive
+; in the four directions to cover points for which a land point has
 ; been considerated (do a small drawing)
 ;----------------------------------------------------------------------------
@@ -154 +149 @@
       endif
 ;-----------------------------------------------------------
-; We mask u and v et v the long of coasts (the place where we 
+; We mask u and v et v the long of coasts (the place where we
 ; can not calculate the average)
-;----------------------------------------------------------- 
+;-----------------------------------------------------------
 ; extention of the mask
       u = u*msku*shift(msku,1,0)
       v = v*mskv*shift(mskv,0,1)
-   ENDIF ELSE BEGIN 
+   ENDIF ELSE BEGIN
       u = u*tmask[firstxt:lastxt,firstyt:lastyt,firstzt]
       v = v*tmask[firstxt:lastxt,firstyt:lastyt,firstzt]
@@ -173 +168 @@
    nan = where(finite(v, /nan) EQ 1)
    if nan[0] NE -1 then v[nan] = 1e5
-   if keyword_set(vectmin) then BEGIN 
+   if keyword_set(vectmin) then BEGIN
       toosmall=where(tabnorme lt vectmin)
       if toosmall[0] NE -1 then begin
@@ -208 +203 @@
          indx = where((lindgen(nx) MOD unvectsur[0]) eq 0)
          indy = where((lindgen(ny) MOD unvectsur[0]) eq 0)
-      ENDIF ELSE BEGIN 
+      ENDIF ELSE BEGIN
          indx = where((lindgen(nx) MOD unvectsur[0]) eq 0)
          indy = where((lindgen(ny) MOD unvectsur[1]) eq 0)
@@ -222 +217 @@
    endif
 ;-----------------------------------------------------------
-; 
+;
 ;-----------------------------------------------------------
    if keyword_set(inverse) then begin
@@ -238 +233 @@
    if terre[0] ne -1 then mini = min(tabnorme[where(t2 eq 1)], max = maxi, /nan) $
    ELSE mini = min(tabnorme, max = maxi, /nan)
-   
+
    if litchamp(vecteur.(0), /u) NE '' then $
     vectlegende = {minmax:[mini, maxi], unite:litchamp(vecteur.(0), /u)} $
@@ -245 +240 @@
 
 sortie:
-   if keyword_set(key_performance) NE 0 THEN print, 'temps ajoutvect', systime(1)-tempsun 
+   if keyword_set(key_performance) NE 0 THEN print, 'temps ajoutvect', systime(1)-tempsun
    return
 end
-
-