source: trunk/SRC/ToBeReviewed/PLOTS/DIVERS/axis4pltz.pro @ 72

;------------------------------------------------------------
;------------------------------------------------------------
;------------------------------------------------------------
;+
; NAME:axis4pltz
;
; PURPOSE:compute the mask and the axis for a vertical section
;
; CATEGORY:
;
; CALLING SEQUENCE:
; 
; INPUTS:mask: 3d mask
;        glam, gphi: 2d longitudes and latitudes
;        z:1d depth
;
; KEYWORD PARAMETERS:
;        XXAXIS to get the xaxis we need to use in pltbase
;        ZZAXIS to get the yaxis we need to use in pltbase
;
;        Others: see pltz 
;
; OUTPUTS:
;
; COMMON BLOCKS:
;
; SIDE EFFECTS:
;
; RESTRICTIONS:
;
; EXAMPLE:
;
; MODIFICATION HISTORY:Sebastien Masson (smasson@lodyc.jussieu.fr)
;                      June 24, 2002
;-
;------------------------------------------------------------
;------------------------------------------------------------
;------------------------------------------------------------

PRO axis4pltz, type, mask, glam, gphi, z, XXAXIS = xxaxis, ZZAXIS = zzaxis, SIN = sin, ZRATIO = zratio, ZOOM = zoom, PROFMAX = profmax, PROFMIN = profmin, _extra = ex
;
; include common
@cm_general
  IF NOT keyword_set(key_forgetold) THEN BEGIN
@updatekwd
  ENDIF
  
;--------------------------------------------------------------------
; 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
;--------------------------------------------------------------------
; define xxaxis and yyaxis the axis used for this section
;--------------------------------------------------------------------
  nx = (size(glam))[1]
  CASE (size(gphi))[0] OF
    1:ny = (size(gphi))[1]
    2:ny = (size(gphi))[2]
  ENDCASE
  CASE (size(z))[0] OF
    1:nz = (size(z))[1]
    2:nz = (size(z))[2]
  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 
;----------------------------------------------------------
; 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

; 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)

   return
end