ARRAY(0x...) ?

To: las_users@xxxxxxxx
Subject: ARRAY(0x...) ?
From: Steve Cousins <cousins@xxxxxxxxxxxxxxxxxxxxxxx>
Date: Fri, 9 Jul 2004 19:01:52 -0400 (EDT)
Reply-to: cousins@xxxxxxxxxxxxxx
Sender: owner-las_users@xxxxxxxxxxxxx
Hi All,

In LAS 5 I was able to edit the draw2d.tmpl file to get LAS to show data
correctly that is on a curvilinear grid and with Sigma coordinate depth.
I'm now trying to get the Sigma depth part to work in LAS 6.3.  When I
try it though, it produces:

----------------------------------------------------------------------

LAS Error

 The following error message was received from LAS: 

 Adding an acceptable error string: "*** NOTE: ".
         NOAA/PMEL TMAP
         FERRET v5.53  
         Linux(g77) 2.4.x - 09/26/03
          9-Jul-04 18:03     

 yes? cancel mode verify
 yes? cancel mode interp
 yes? define symbol refmap_xhi = 290.0
 yes? define symbol refmap_yhi = -17.0
 yes? define symbol land_type = shade
 yes? define symbol refmap_xlo = 100.0
 yes? define symbol jnl_t = ARRAY(0x8cde4c4)
 yes? define symbol refmap_view = xline
 yes? define symbol refmap_ylo = -17.0
 yes? define symbol variable_name_0 = no3
 yes? define symbol fill_type = fill
 yes? define symbol jnl_x = ARRAY(0x8cde488)
 yes? define symbol variable_name = ARRAY(0x88605dc)
 yes? define symbol jnl_y = ARRAY(0x8cde4b8)
 yes? define symbol jnl_z = ARRAY(0x8cde2b4)
 yes? define symbol diag_direction = down
 yes? define symbol dataset_name_0 =
/juno/usr4/lshi/RESULT_NP_1990_2001_NCEP_daily/NP_avg.001.nc
 yes? define symbol dataset_name = ARRAY(0x8860600)
 yes? define symbol rank = 2
 yes? canc data/all
 yes? go std_initialize
"/juno/usr4/lshi/RESULT_NP_1990_2001_NCEP_daily/NP_avg.001.nc" "1" "1"
"no3"
 yes? set
region/x="100.0":"290.0"/y="-17.0"/z="1":"20"/t="01-Jan-1990":"01-Jan-1990"
 yes? set win/size=.5
 yes? GO "/tmp/lasgo223198511390.jnl"
  **ERROR: command syntax:
ZAXREPLACE(ARRAY(0x88605dc),depth_at_sigma,z[gz=zdepth])
           RRAY(0x88
 DEFINE VARIABLE varonsigma=ZAXREPLACE(ARRAY(0x88605dc),depth_at_sigma,z[gz=zdepth])
 Command file, command group, or REPEAT execution aborted

----------------------------------------------------------------------

The thing that sticks out is the ARRAY(0x88605dc).  I hadn't seen this
before. The /tmp/lasgo... file looks like:

----------------------------------------------------------------------

   can view

   define view/x=0,1/y=0.25,1 vt
   define view/x=0.029,.975/y=0,0.25 vb

   set view vt

   set MEMORY/SIZE=20.0

   DEFINE AXIS/Z=-5000.0:0.0:10.0/UNIT=meters/DEPTH zdepth 
   LET varonsigma=ZAXREPLACE(ARRAY(0x88605cc),depth_at_sigma,z[gz=zdepth])
   LET var_at_depth = varonsigma[  ]

   let var_kend = `var_at_depth,RETURN=KEND`

   go preplot_setup
     
   fill/levels=30/set no3,'X',depth_at_sigma
    
   go postplot_setup
   palette default
     
   ppl title,0.15,averaged nitrate (mMol N m-3)
   ppl fill
   ppl %range 'PPL$ZMIN' 'PPL$ZMAX' 30
   ppl set PPL$RANGE_INC2 'PPL$RANGE_INC' * 2
 
contour/over/levels=('PPL$RANGE_LOW','PPL$RANGE_HIGH','PPL$RANGE_INC2')/nolabel no3

   set view vb
   ppl origin ,.7
   ppl axlen ,2

----------------------------------------------------------------------

In order to compare it with my LAS 5 server, I put a typo in the
draw_2d.tmpl file so the /tmp/lasgo... file would not be deleted.  Here is
what I get with LAS 5:

----------------------------------------------------------------------

   can view

   define view/x=0,1/y=0.25,1 vt
   define view/x=0.029,.975/y=0,0.25 vb

   set view vt

   set MEMORY/SIZE=20.0
      
   DEFINE AXIS/Z=-5000.0:0.0:10.0/UNIT=meters/DEPTH zdepth 
   LET varonsigma= ZAXREPLACE(temp,depth_at_sigma,z[gz=zdepth])
   LET var_at_depth = varonsigma[  ]

   let var_kend = `var_at_depth,RETURN=KEND`

    go preplot_setup
     
    fill/levels=30/set temp,'X',depth_at_igma
    
    go postplot_setup
    palette default

    ppl title,0.2,temp
    
    ppl fill
    ppl %range 'PPL$ZMIN' 'PPL$ZMAX' 30
    ppl set PPL$RANGE_INC2 'PPL$RANGE_INC' * 2
	
    contour/over/levels=1/nolabel temp,'X',depth_at_sigma
    set view vb

    ppl origin ,.7
    ppl axlen ,2

----------------------------------------------------------------------

So you see that it doesn't use ARRAY(0x88605cc) in the ZAXREPLACE line.

Is this the problem?  Ferret doesn't know what ARRAY(0x88605cc) is?  Maybe
I'm going down the wrong path.  Just to make this email even longer, I'll
include the draw_2d.tmpl file:

----------------------------------------------------------------------

[% DEFAULT
    args.fill_type = 'fill'
    args.fill_levels = '30'
    args.palette = 'default'
    args.do_shade = 1
    args.do_contour = 0
    args.do_overlay = 0
%]


[% IF args.fill_type == 'default';
    args.fill_type = 'fill' ;
ELSIF args.fill_type == 'none';
   args.do_shade = 0 ;
END %]

! Support for <contour_levels>none</contour_levels>
! Must be overridden for overlay plots.
[% IF args.contour_levels == 'none';
    args.do_contour = 0 ;
    args.contour_levels = '' ;
END %]



! NetCDF strides should be used only for XY views
! NetCDF strides should not be used for comparison.
[% IF args.view == 'xy' && args.variable_name.0 != 'diff_' &&
args.variable_name.0 != 'data1_' && args.variable_name.0 != 'data2_' %]
    go stride_create_symbols [% args.variable_name.0 %] 300 200
[% END %]

! SDC: Dealing with Sigma (actually S-coordinates in this case) 

[% IF args.view == 'xy' %]
      [% needz = 'NO' %]
   [% ELSIF args.view == 'xz' %]
      [% needz = 'YES' %]
   [% ELSIF args.view == 'yz' %]
      [% needz = 'YES' %]
   [% ELSE %]
      [% needz = 'NO' %]
[% END %]
    
[% IF needz == 'YES' %]
   set MEMORY/SIZE=20.0

!   let sigma = z[gz=zpos]
!   let depth_x_sigma = sigma * depth
      
   DEFINE AXIS/Z=-5000.0:0.0:10.0/UNIT=meters/DEPTH zdepth 
!   LET varonsigma= ZAXREPLACE([% args.variable_name
%],depth_x_sigma,z[gz=zdepth])
   LET varonsigma= ZAXREPLACE([% args.variable_name
%],depth_at_sigma,z[gz=zdepth])
   LET var_at_depth = varonsigma[ [% args.z_lo %] ]

   let var_kend = `var_at_depth,RETURN=KEND`
[% END %]

[% IF args.do_shade %]
    go preplot_setup
    [% IF args.view == 'xy' && args.variable_name.0 != 'diff_' &&
args.variable_name.0 != 'data1_' && args.variable_name.0 != 'data2_' %]
        [% args.fill_type %]/levels=[% args.fill_levels %]/set [%
args.variable_name.0 %][i=($STRIDE_irange),j=($STRIDE_jrange)]
    [% ELSIF args.view == 'xz' %] 
        [% args.fill_type %]/levels=[% args.fill_levels %]/set [%
args.variable_name.0 %],'X',depth_at_sigma
    [% ELSIF args.view == 'yz' %]
        [% args.fill_type %]/levels=[% args.fill_levels %]/set [%
args.variable_name.0 %],'Y',depth_at_sigma    
    [% ELSE %]
        [% args.fill_type %]/levels=[% args.fill_levels %]/set [%
args.variable_name.0 %]
    [% END %]
    go postplot_setup
    palette [% args.palette %]
    [% IF args.title %] 
        ppl title,0.15,[% args.title %]
    [% END %]
    [% INCLUDE labels.tmpl args = args %]
    ppl [% args.fill_type %]
[% END %]

[% IF args.do_contour %]
    [% IF args.do_shade %]
	ppl %range 'PPL$ZMIN' 'PPL$ZMAX' 30
	ppl set PPL$RANGE_INC2 'PPL$RANGE_INC' * 2
	[% IF ! args.contour_levels %]
	    [% args.contour_levels =
"('PPL\$RANGE_LOW','PPL\$RANGE_HIGH','PPL\$RANGE_INC2')" %]
	[% END %]
    [% ELSE %]
	[% IF ! args.contour_levels %]
	    [% args.contour_levels = "20" %]
	[% END %]
    [% END %]
    [% IF args.view == 'xy' && args.variable_name.0 != 'diff_' &&
args.variable_name.0 != 'data1_' && args.variable_name.0 != 'data2_' %]
        contour/over/levels=[% args.contour_levels %]/nolabel [%
args.variable_name.0 %][i=($STRIDE_irange),j=($STRIDE_jrange)]
    [% ELSE %]
        contour/over/levels=[% args.contour_levels %]/nolabel [%
args.variable_name.0 %]
    [% END %]
[% END %]

[% IF args.do_overlay %]
    [% IF ! args.contour_levels %]
        [% args.contour_levels = '20' %]
    [% END %]
    [% IF args.view == 'xy' && args.variable_name.0 != 'diff_' &&
args.variable_name.0 != 'data1_' && args.variable_name.0 != 'data2_' %]
        contour/over/levels=[% args.contour_levels %]/nolabel [%
args.overlay_variable_name.0 %][i=($STRIDE_irange),j=($STRIDE_jrange)]
    [% ELSE %]
        contour/over/levels=[% args.contour_levels %]/nolabel [%
args.overlay_variable_name.0 %]
    [% END %]

[% END %]
----------------------------------------------------------------------

Unfortunately it doesn't wrap very well.  If you want me to just attach it
let me know.

Can anyone tell me what the problem is?

Thanks very much and have a good weekend.

Steve
______________________________________________________________________
 Steve Cousins, Ocean Modeling Group    Email: cousins@umit.maine.edu
 Marine Sciences, 208 Libby Hall        http://rocky.umeoce.maine.edu
 Univ. of Maine, Orono, ME 04469        Phone: (207) 581-4302
Follow-Ups:
- installing on sun-os 8
  - From: J-M Epitalon
- Re: ARRAY(0x...) ?
  - From: Joe McLean
Previous by thread: in-situ longitude problems.
Next by thread: Re: ARRAY(0x...) ?
[Thread Prev][Thread Next][Index]
Dept of Commerce / NOAA / OAR / PMEL / TMAP
Contact Us | Privacy Policy | Disclaimer | Accessibility Statement