Changeset 41 for trunk

Timestamp:

05/02/06 17:12:07 (18 years ago)

Author:

pinsard

Message:

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

Location:

trunk

Files:

• PLOTS/VECTEUR (deleted)
• ToBeReviewed/PLOTS/VECTEUR/ajoutvect.pro (copied) (copied from trunk/PLOTS/VECTEUR/ajoutvect.pro) (9 diffs)
• ToBeReviewed/PLOTS/VECTEUR/vecteur.pro (copied) (copied from trunk/PLOTS/VECTEUR/vecteur.pro) (5 diffs)
• ToBeReviewed/PLOTS/VECTEUR/velovect.pro (copied) (copied from trunk/PLOTS/VECTEUR/velovect.pro)

trunk/ToBeReviewed/PLOTS/VECTEUR/ajoutvect.pro

@@ -62 +62 @@
    v = litchamp(vecteur.(1))
    v = checkfield(v, 'plt', TYPE = 'xy', /NOQUESTION)
-;         niv=''
-   IF (size(u))[0] EQ 3 THEN BEGIN 
-      niveau = xquestion('Le tableau d''entree est un tableau 3d, a quel niveau faut-il faire les vecteurs? (par defaut le premier niveau) ', strtrim(niveau, 1), /chkwidget)
-      niveau = 1 > fix(niveau) < jpk
-      u=u(*,*,niveau-1)
-      v=v(*,*,niveau-1)
-   endif
-
 ;-----------------------------------------------------------
 ;------------------------------------------------------------
@@ -78 +70 @@
    grillev = litchamp(vecteur.(1), /grid)
    if grillev EQ '' then grillev = 'V'
+
    IF grilleu EQ 'V' AND grillev EQ 'U' THEN inverse = 1
    IF grilleu EQ grillev THEN interpolle  = 0 ELSE interpolle = 1
@@ -89 +82 @@
 ;------------------------------------------------------------
    if interpolle then begin
-      indicexu = (lindgen(jpi))[premierxu:premierxu+nxu-1]
-      indicexv = (lindgen(jpi))[premierxv:premierxv+nxv-1]
+      indicexu = (lindgen(jpi))[firstxu:firstxu+nxu-1]
+      indicexv = (lindgen(jpi))[firstxv:firstxv+nxv-1]
       indicex = inter(indicexu, indicexv)
-      indiceyu = (lindgen(jpj))[premieryu:premieryu+nyu-1]
-      indiceyv = (lindgen(jpj))[premieryv:premieryv+nyv-1]
+      indiceyu = (lindgen(jpj))[firstyu:firstyu+nyu-1]
+      indiceyv = (lindgen(jpj))[firstyv:firstyv+nyv-1]
       indicey = inter(indiceyu, indiceyv)
       nx = n_elements(indicex) 
@@ -107 +100 @@
           (size(v))[1] EQ nxv AND (size(v))[2] EQ nyv:BEGIN 
             if nxu NE nx then $
-             if indicex[0] EQ premierxu 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 premierxv then v = v[0:nx-1, *] ELSE v = v[1: nx, *]
+             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 premieryu 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 premieryv then v = v[*, 0:ny-1] ELSE v = v[*, 1: ny]
+             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 $
@@ -137 +130 @@
       a=u(0,*)
       u=(u+shift(u,1,0))/2.
-      if NOT keyword_set(key_periodique) OR nx NE jpi then u(0,*)=a
+      if NOT keyword_set(key_periodic) OR nx NE jpi then u(0,*)=a
       a=v(*,0)
       v=(v+shift(v,0,1))/2.
-      if NOT keyword_set(key_periodique) OR nx NE jpi then v(*,0)=a
+      if NOT keyword_set(key_periodic) OR nx NE jpi then v(*,0)=a
 ;----------------------------------------------------------------------------
 ; attribution du mask et des tableau de longitude et latitude
@@ -148 +141 @@
 ;----------------------------------------------------------------------------
       vargrid='T'
-      msku = (umask())[indice2d+jpi*jpj*(niveau-1)]
-      mskv = (vmask())[indice2d+jpi*jpj*(niveau-1)]
+      msku = (umask())[indice2d+jpi*jpj*firstzt]
+      mskv = (vmask())[indice2d+jpi*jpj*firstzt]
       glam = glamt[indice2d]
       gphi = gphit[indice2d]
@@ -166 +159 @@
       v = v*mskv*shift(mskv,0,1)
    ENDIF ELSE BEGIN 
-      u = u*tmask[premierxt:dernierxt,premieryt:dernieryt,premierzt]
-      v = v*tmask[premierxt:dernierxt,premieryt:dernieryt,premierzt]
+      u = u*tmask[firstxt:lastxt,firstyt:lastyt,firstzt]
+      v = v*tmask[firstxt:lastxt,firstyt:lastyt,firstzt]
       indice2d = lindgen(jpi, jpj)
-      indice2d = indice2d[premierxt:dernierxt, premieryt:dernieryt]
+      indice2d = indice2d[firstxt:lastxt, firstyt:lastyt]
       nx = nxt
       ny = nyt
@@ -179 +172 @@
    if nan[0] NE -1 then v[nan] = 1e5
    if keyword_set(vectmin) then BEGIN 
-      troppetit=where(tabnorme lt vectmin)
-      if troppetit[0] NE -1 then begin
-         u[troppetit] = 1e5
-         v[troppetit] = 1e5
+      toosmall=where(tabnorme lt vectmin)
+      if toosmall[0] NE -1 then begin
+         u[toosmall] = 1e5
+         v[toosmall] = 1e5
       ENDIF
    endif
    if keyword_set(vectmax) then BEGIN
-      tropgrand=where(tabnorme gt vectmax)
-      if tropgrand[0] NE -1 then begin
-         u[tropgrand] = 1e5
-         v[tropgrand] = 1e5
+      toobig=where(tabnorme gt vectmax)
+      if toobig[0] NE -1 then begin
+         u[toobig] = 1e5
+         v[toobig] = 1e5
       ENDIF
    ENDIF
@@ -197 +190 @@
 ;-----------------------------------------------------------
    if interpolle then t2 = msku*shift(msku,1,0)*mskv*shift(mskv,0,1) $
-   ELSE t2 = tmask[premierxt:dernierxt,premieryt:dernieryt,premierzt]
-   if NOT keyword_set(key_periodique) OR nx NE jpi then t2(0, *)=0.
+   ELSE t2 = tmask[firstxt:lastxt,firstyt:lastyt,firstzt]
+   if NOT keyword_set(key_periodic) OR nx NE jpi then t2(0, *)=0.
    t2(*,0)=0.
    terre=where(t2 eq 0)

trunk/ToBeReviewed/PLOTS/VECTEUR/vecteur.pro

@@ -84 +84 @@
 ; quelle est la longeur en coordonnees normales d''un trait qui fera 1
 ; cm sur le parier et qui est parallele a x?
-   sizexfeuille = petitfeuille*key_portrait+grandfeuille*(1-key_portrait)
-   sizeyfeuille = grandfeuille*key_portrait+petitfeuille*(1-key_portrait)
+  mipgsz = min(page_size, max = mapgsz)
+   sizexfeuille = mipgsz*key_portrait+mapgsz*(1-key_portrait)
+   sizeyfeuille = mapgsz*key_portrait+mipgsz*(1-key_portrait)
    cm_en_normal = 1./sizexfeuille
 ;
@@ -217 +218 @@
 ; definition de nu
    radius = replicate(1,nxgd*nygd)
-   coord_sphe = transpose([ [(glamu[indice2d])[*]], [(gphiu[indice2d])[*]], [radius(*)] ])
+   IF finite(glamu[0]*gphiu[0]) NE 0 THEN $
+     coord_sphe = transpose([ [(glamu[indice2d])[*]], [(gphiu[indice2d])[*]], [radius(*)] ]) $
+   ELSE coord_sphe = transpose([ [(glamf[indice2d])[*]], [(gphit[indice2d])[*]], [radius(*)] ])
    r = cv_coord(from_sphere=coord_sphe,/to_rect,/degrees)
 ; coordonnes de points de la grille u en cartesien
@@ -228 +231 @@
    nuz = uz-shift(uz, 1, 0)
 ; conditions aux limites
-   if NOT keyword_set(key_periodique) OR nxgd NE jpi then begin
+   if NOT keyword_set(key_periodic) OR nxgd NE jpi then begin
       nux[0, *] = nux[1, *]
       nuy[0, *] = nuy[1, *]
@@ -238 +241 @@
    nuz = nuz[reduitindice2d]
 ; definition de nv
-   coord_sphe = transpose([ [(glamv[indice2d])[*]], [(gphiv[indice2d])[*]], [radius(*)] ])
+   IF finite(glamv[0]*gphiv[0]) NE 0 THEN $
+   coord_sphe = transpose([ [(glamv[indice2d])[*]], [(gphiv[indice2d])[*]], [radius(*)] ]) $
+   ELSE coord_sphe = transpose([ [(glamt[indice2d])[*]], [(gphif[indice2d])[*]], [radius(*)] ])                
    r = cv_coord(from_sphere=coord_sphe,/to_rect,/degrees)
 ; coordonnes de points de la grille v en cartesien
@@ -353 +358 @@
 ;
    if NOT keyword_set(cmref) then BEGIN 
-      sizexfeuille = petitfeuille*key_portrait+grandfeuille*(1-key_portrait)
+      mipgsz = min(page_size, max = mapgsz)
+      sizexfeuille = mipgsz*key_portrait+mapgsz*(1-key_portrait)
       sizexfeuille = 10.*sizexfeuille
       cmref = 5 > floor(sizexfeuille/10.) < 15