[Thread Prev][Thread Next][Index]
Re: Use of 'X' with Contour in LAS 6.3
Thanks Joe (and Cyndy). This has allowed me to continue but I now get a
different error:
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
12-Jul-04 11:46
yes? cancel mode verify
yes? cancel mode interp
yes? define symbol refmap_xhi = 290.071990966797
yes? define symbol refmap_yhi = 10.030584335327148
yes? define symbol land_type = shade
yes? define symbol refmap_xlo = 99.7431030273438
yes? define symbol jnl_t = ARRAY(0x8cde450)
yes? define symbol refmap_view = xline
yes? define symbol refmap_ylo = 10.030584335327148
yes? define symbol variable_name_0 = no3
yes? define symbol fill_type = fill
yes? define symbol jnl_x = ARRAY(0x8cde2c4)
yes? define symbol variable_name = ARRAY(0x886064c)
yes? define symbol jnl_y = ARRAY(0x8cde3fc)
yes? define symbol jnl_z = ARRAY(0x8cde414)
yes? define symbol diag_direction = down
yes? define symbol dataset_name_0 =
/usr2/fer_dsets/descr/pacific-roms-90-01.des
yes? define symbol dataset_name = ARRAY(0x8860670)
yes? define symbol rank = 2
yes? canc data/all
yes? go std_initialize "/usr2/fer_dsets/descr/pacific-roms-90-01.des" "1"
"1" "no3"
yes? set
region/x="99.7431030273438":"290.071990966797"/y="10.030584335327148"/z="1":"20"/l="1":"1"
yes? set win/size=.5
yes? GO "/tmp/lasgo31108960517.jnl"
**ERROR: variable unknown or not in data set: 'X'
CONTOUR/FILL/levels=30/set no3,'X',depth_at_sigma
Command file, command group, or REPEAT execution aborted
I'm just adding the 'X', and depth_at_sigma variables to the end of the
CONTOUR line. 'X' is supposed to refer to the X dimension variable as I
understand it. It works this way in LAS 5.
Any ideas?
Thanks,
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
On Sun, 11 Jul 2004, Joe McLean wrote:
> Hi Steve,
>
> As of LAS v6.0 template access to args.variable_name refers to an array
> rather than a string variable . Easy fix though.
> You'll notice in 'draw_2d.tmpl' which you sent that most references to
> the variable name are in the following syntax: args.variable_name.0
>
> There is still one place in your template which refers to 'args.variable_name' instead of 'args.variable_name.0'. Make this change and all should be well.
>
> Joe
> ---------------------------
> [...]
> ! SDC: Dealing with Sigma (actually S-coordinates in this case)
> [...]
> 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 %]
> --------------------------
>
>
> Steve Cousins wrote:
>
> >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
> >
> >
> >
> >
> >
>
[Thread Prev][Thread Next][Index]
Dept of Commerce /
NOAA /
OAR /
PMEL /
TMAP
Contact Us | Privacy Policy | Disclaimer | Accessibility Statement