Changeset 74 for trunk/SRC/ToBeReviewed/PLOTS/DIVERS/axis4pltz.pro

Timestamp:

05/23/06 15:45:47 (18 years ago)

Author:

smasson

Message:

debug xxx and cie + clean data file + rm perldoc_idl

File:

• trunk/SRC/ToBeReviewed/PLOTS/DIVERS/axis4pltz.pro (modified) (3 diffs)

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

@@ -41 +41 @@
 ;
 ; include common
-@cm_general
+@cm_4mesh
   IF NOT keyword_set(key_forgetold) THEN BEGIN
 @updatekwd
@@ -49 +49 @@
 ; define the mask used for this section
 ;--------------------------------------------------------------------
-   if mask[0] NE -1 AND (size(mask))[0] NE 2 then begin
-      if type EQ 'xz' then mask = total(mask, 2) < 1 $
-      ELSE mask = total(mask, 1) < 1
-   endif
+  if mask[0] NE -1 AND (size(mask))[0] NE 2 then begin
+    if type EQ 'xz' then mask = total(mask, 2) < 1 $
+    ELSE mask = total(mask, 1) < 1
+  endif
 ;--------------------------------------------------------------------
 ; define xxaxis and yyaxis the axis used for this section
@@ -66 +66 @@
   ENDCASE
 ;
-   if type eq 'yz' then BEGIN 
-      if (size(gphi))[0] EQ 2 then xxaxis=reform(gphi[0, *]) ELSE xxaxis=gphi 
-      if keyword_set(sin) then xxaxis = sin(!pi/180.*xxaxis)
-      if (size(z))[0] EQ 1 THEN zzaxis = z ELSE zzaxis = z
-   ENDIF ELSE BEGIN
-      xxaxis = glam[*, 0]
-      if (size(z))[0] EQ 1 then zzaxis = z ELSE zzaxis = z
-   ENDELSE 
+  if type eq 'yz' then BEGIN 
+    IF (size(gphi))[0] EQ 1 then xxaxis = gphi ELSE BEGIN 
+      IF keyword_set(key_irregular) THEN BEGIN 
+        cln = (where(gphi EQ max(gphi)))[0]
+        xxaxis = reform(gphi[cln MOD nx, *])
+      ENDIF ELSE xxaxis = reform(gphi[0, *])
+    ENDELSE
+    if keyword_set(sin) then xxaxis = sin(!pi/180.*xxaxis)
+    if (size(z))[0] EQ 1 THEN zzaxis = z ELSE zzaxis = z
+  ENDIF ELSE BEGIN
+    xxaxis = glam[*, 0]
+    if (size(z))[0] EQ 1 then zzaxis = z ELSE zzaxis = z
+  ENDELSE 
 ;----------------------------------------------------------
 ; on projette l''axe z dans [0,1]
 ;----------------------------------------------------------
-   if not keyword_set(zratio) then zratio = 2./3
-   if zoom ge profmax then zratio=1.
-   if zoom LT profmax then begin
-      mp = projsegment([profmin, zoom], [0, zratio], /mp)
-      zzaxis[where(zzaxis LE zoom)] = mp[0]*zzaxis[where(zzaxis LE zoom)]+mp[1]
-      mp = projsegment([zoom, profmax], [zratio, 1], /mp)
-      zzaxis[where(zzaxis GE zoom)] = mp[0]*zzaxis[where(zzaxis GE zoom)]+mp[1]
-   ENDIF ELSE BEGIN
-      mp = projsegment([profmin, profmax], [0, 1], /mp)
-      zzaxis = mp[0]*zzaxis+mp[1]
-   ENDELSE
+  if not keyword_set(zratio) then zratio = 2./3
+  if zoom ge profmax then zratio = 1.
+  if zoom LT profmax then begin
+    mp = projsegment([profmin, zoom], [0, zratio], /mp)
+    zzaxis[where(zzaxis LE zoom)] = mp[0]*zzaxis[where(zzaxis LE zoom)]+mp[1]
+    mp = projsegment([zoom, profmax], [zratio, 1], /mp)
+    zzaxis[where(zzaxis GE zoom)] = mp[0]*zzaxis[where(zzaxis GE zoom)]+mp[1]
+  ENDIF ELSE BEGIN
+    mp = projsegment([profmin, profmax], [0, 1], /mp)
+    zzaxis = mp[0]*zzaxis+mp[1]
+  ENDELSE
 
 ; to draw from bottom to top (avoid using cell_fill)
-   CASE size(zzaxis, /n_dimensions) OF
-     1:zzaxis=reverse(-zzaxis) 
-     2:zzaxis=reverse(-zzaxis,2) 
-   ENDCASE
-   if mask[0] NE -1 then mask=reverse(mask,2)
+  CASE size(zzaxis, /n_dimensions) OF
+    1:zzaxis = reverse(-zzaxis) 
+    2:zzaxis = reverse(-zzaxis, 2) 
+  ENDCASE
+  if mask[0] NE -1 then mask = reverse(mask, 2)
 
-   return
+  return
 end