Changeset 41 for trunk/ToBeReviewed/PLOTS/VECTEUR/vecteur.pro

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

File:

• trunk/ToBeReviewed/PLOTS/VECTEUR/vecteur.pro (copied) (copied from trunk/PLOTS/VECTEUR/vecteur.pro) (5 diffs)

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