[Thread Prev][Thread Next][Index]

Two conundrums





LAS is pretty slick once you get it running.
So far we've done a few test runs and most of it runs well. I'm having some 'issues', however.
But before we get into those, I've attached the xml I'm using. The GDS system is publicly accessible, so if you'd like to play around with the data we have, you're more than welcome to.
Ignore some of the...umm..oddness in the xml. I've been trying to track down errors this morning.

Anyway, when I attempt to plot "air temperature - 29 pressure levels [k]" I receive this error:
--- Calling: LAS::Server::Ferret::overlay
Sent command: set region/x="-124.0":"-52.0"/y="22.0":"63.0"/z="550":"550"/t="01-Jan-1979":"01-Jan-1979"
Got reply: yes?
Sent command: set win/size=.5
Got reply: yes?
Aspect: 0.854166666666667
Sent command: set win/asp=0.854166666666667
Got reply: yes?
Journal search path: custom:jnls:jnls/insitu:jnls/section

Searching for template/journal file: custom/std_gif.tmpl

Searching for template/journal file: custom/std_gif_xyref.jnl

Searching for template/journal file: custom/std_gif.tmpl

Searching for template/journal file: custom/std_gif_xyref.jnl

Searching for template/journal file: jnls/std_gif.tmpl

Running template file jnls/std_gif.tmpl
template: props is: HASH(0x93ef844) arg is: 'dataset_name' val is: ,
template: props is: HASH(0x93ef844) arg is: 'variable_name' val is: ,
template: props is: HASH(0x93ef844) arg is: 'fill_type' val is: fill
template: props is: HASH(0x93ef844) arg is: 'view' val is: xy
template: props is: HASH(0x93ef844) arg is: 'fill_levels' val is:
template: props is: HASH(0x93ef844) arg is: 'contour_levels' val is:
template: props is: HASH(0x93ef844) arg is: 'magnify' val is: 0
template: props is: HASH(0x93ef844) arg is: 'palette' val is:
template: props is: HASH(0x93ef844) arg is: 'draw_reference' val is: 1
template: props is: HASH(0x93ef844) arg is: 'do_shade' val is:
template: props is: HASH(0x93ef844) arg is: 'do_contour' val is: 1
template: props is: HASH(0x93ef844) arg is: 'do_overlay' val is:
template: props is: HASH(0x93ef844) arg is: 'rank' val is: 2
template: props is: HASH(0x93ef844) arg is: 'title' val is: @AS@TRair temperature - 29 pressure levels [k] ()
template: props is: HASH(0x93ef844) arg is: 'dump_args' val is:
template: props is: HASH(0x93ef844) arg is: 'dump_file' val is:
template: props is: HASH(0x93ef844) arg is: 'line_color' val is:
template: props is: HASH(0x93ef844) arg is: 'line_thickness' val is:
Template search path: custom:.:jnls:jnls/insitu:jnls/section
Sent command: GO "/tmp/lasgo50566485126.jnl"
Got reply: **ERROR: illegal limits: "TMPPRS" is not in the range T=-6.2451E+10
Axis extremes are T=0.5:75977
DEFINE VARIABLE STRIDE_istart = `tmpprs,return=istart`
Command file, command group, or REPEAT execution aborted

Error: transcript is
Adding an acceptable error string: "*** NOTE: ".
NOAA/PMEL TMAP
FERRET v5.80 Linux(g77) 2.4.20 - 01/03/05
24-Mar-05 11:51
yes? cancel mode verify
yes? cancel mode interp
yes? canc data/all
yes? go std_initialize "http://nomads4.ncdc.noaa.gov:9091/dods/NCEP_NARR_VAR_TEST/narr-a_221_tmpprs"; "1" "1" "tmpprs"
*** NOTE: Coordinates out of order or missing on axis time at subscript 2
*** NOTE: A dummy axis of subscripts will be used
yes? set region/x="-124.0":"-52.0"/y="22.0":"63.0"/z="550":"550"/t="01-Jan-1979":"01-Jan-1979"
yes? set win/size=.5
yes? set win/asp=0.854166666666667
yes? GO "/tmp/lasgo50566485126.jnl"
**ERROR: illegal limits: "TMPPRS" is not in the range T=-6.2451E+10
Axis extremes are T=0.5:75977
DEFINE VARIABLE STRIDE_istart = `tmpprs,return=istart`
Command file, command group, or REPEAT execution aborted

killing...Waiting for program to exit...
program exited...
killing...Waiting for program to exit...
program exited...
*********************

Also, attempting to overlay '0-3 hr acc convective precipitation [kg/m^2]' and 'anl. geopotential height [gpm] - 29 pressure levels' I recieve a "The region you chose is too small" error message. The debug log is below.


***************
--- Calling: LAS::Server::Ferret::accessDatabase
database_access prop: dsetname = file:dummy.nc

--- Returning from: LAS::Server::Ferret::accessDatabase
FER_GO=custom custom jnls jnls/insitu jnls/section . . /usr/local/ferret/go /usr/local/ferret/examples /usr/local/ferret/contrib
Starting:/usr/local/ferret/bin/ferret -gif -server
Got reply: NOAA/PMEL TMAP
FERRET v5.80 Linux(g77) 2.4.20 - 01/03/05
24-Mar-05 12:02
yes?
Sent command: cancel mode verify
Got reply: yes?
Sent command: cancel mode interp
Got reply: yes?
Sent command: canc data/all
Got reply: yes?
Sent command: go std_initialize "http://nomads4.ncdc.noaa.gov:9091/dods/NCEP_NARR_VAR_TEST/narr-a_221_pcp"; "1" "1" "acpcp"
Got reply: yes?
Sent command: go std_initialize "http://nomads4.ncdc.noaa.gov:9091/dods/NCEP_NARR_VAR_TEST/narr-a_221_hgtprs"; "2" "2" "hgtprs"
Got reply: yes?
LAS::Server::Handler::run(): 0.273446 secs to perform all 'preExecute' activities (database access and parsing)

--- Calling: LAS::Server::Ferret::overlay
names = "** 0-3 hr acc convective precipitation [kg/m^2],anl. geopotential height [gpm] - 29 pressure levels,North American Regional Reanalysis,North American Regional Reanalysis,,"
Compare var expr: 1:[y="22.0":"63.0",x="-124.0":"-52.0",d=1,t="01-Jan-1979"]
Compare var expr: 2:[y="22.0":"63.0",gx=var1_grid_,gy=var1_grid_,x="-124.0":"-52.0",d=2,t="01-Jan-1979"]
Sent command: let var1_ = acpcp[y="22.0":"63.0",x="-124.0":"-52.0",d=1,t="01-Jan-1979"]
Got reply: yes?
Sent command: let var2_ = hgtprs[y="22.0":"63.0",gx=var1_grid_,gy=var1_grid_,x="-124.0":"-52.0",d=2,t="01-Jan-1979"]
Got reply: yes?
Sent command: set var/units="`acpcp[y="22.0":"63.0",x="-124.0":"-52.0",d=1,t="01-Jan-1979"],return=units`" var1_
Got reply: yes?
Sent command: set var/units="`hgtprs[y="22.0":"63.0",gx=var1_grid_,gy=var1_grid_,x="-124.0":"-52.0",d=2,t="01-Jan-1979"],return=units`" var2_
Got reply: yes?
Sent command: let data1_ = var1_
Got reply: yes?
Sent command: let data2_ = var2_
Got reply: yes?
title part1,part2: Regional Reanalysis(1),<NL>@P2anl. geopotential height [gpm] - 29 pressure levels () from North American Regional ... Reanalysis(2)(contours)
title part1,part2: @AS@TR** 0-3 hr acc convective precipitation [kg/m^2] () from North American ... Regional Reanalysis(1),<NL>@P2anl. geopotential height [gpm] - 29 pressure levels () from North American Regional<NL>Reanalysis(2)(contours)
Sent command: set var/units="`var1_,return=units`" data1_
Got reply: yes?
Sent command: set var/units="`var2_,return=units`" data2_
Got reply: yes?
Sent command: set win/size=.5
Got reply: yes?
Aspect: 0.854166666666667
Sent command: set win/asp=0.854166666666667
Got reply: yes?
Journal search path: custom:jnls:jnls/insitu:jnls/section

Searching for template/journal file: custom/std_gif.tmpl

Searching for template/journal file: custom/std_gif_xyref.jnl

Searching for template/journal file: custom/std_gif.tmpl

Searching for template/journal file: custom/std_gif_xyref.jnl

Searching for template/journal file: jnls/std_gif.tmpl

Running template file jnls/std_gif.tmpl
template: props is: HASH(0x8f632e0) arg is: 'dataset_name' val is: ,
template: props is: HASH(0x8f632e0) arg is: 'variable_name' val is: ,
template: props is: HASH(0x8f632e0) arg is: 'fill_type' val is: fill
template: props is: HASH(0x8f632e0) arg is: 'view' val is: xy
template: props is: HASH(0x8f632e0) arg is: 'fill_levels' val is:
template: props is: HASH(0x8f632e0) arg is: 'contour_levels' val is:
template: props is: HASH(0x8f632e0) arg is: 'magnify' val is: 0
template: props is: HASH(0x8f632e0) arg is: 'palette' val is:
template: props is: HASH(0x8f632e0) arg is: 'draw_reference' val is: 1
template: props is: HASH(0x8f632e0) arg is: 'do_shade' val is:
template: props is: HASH(0x8f632e0) arg is: 'do_contour' val is: 0
template: props is: HASH(0x8f632e0) arg is: 'do_overlay' val is: 1
template: props is: HASH(0x8f632e0) arg is: 'rank' val is: 2
template: props is: HASH(0x8f632e0) arg is: 'title' val is: @AS@TR** 0-3 hr acc convective precipitation [kg/m^2] () from North American<NL>Regional Reanalysis(1),<NL>@P2anl. geopotential height [gpm] - 29 pressure levels () from North American Regional<NL>Reanalysis(2)(contours)
template: props is: HASH(0x8f632e0) arg is: 'dump_args' val is:
template: props is: HASH(0x8f632e0) arg is: 'dump_file' val is:
template: props is: HASH(0x8f632e0) arg is: 'line_color' val is:
template: props is: HASH(0x8f632e0) arg is: 'line_thickness' val is:
Template search path: custom:.:jnls:jnls/insitu:jnls/section
Sent command: GO "/tmp/lasgo51741959833.jnl"
Got reply: **ERROR: dimensions improperly specified: must be a 2D region
contour/over/levels=20/nolabel data2_
contour/over/levels=20/nolabel data2_
Command file, command group, or REPEAT execution aborted

Error: transcript is
Adding an acceptable error string: "*** NOTE: ".
NOAA/PMEL TMAP
FERRET v5.80 Linux(g77) 2.4.20 - 01/03/05
24-Mar-05 12:02
yes? cancel mode verify
yes? cancel mode interp
yes? canc data/all
yes? go std_initialize "http://nomads4.ncdc.noaa.gov:9091/dods/NCEP_NARR_VAR_TEST/narr-a_221_pcp"; "1" "1" "acpcp"
yes? go std_initialize "http://nomads4.ncdc.noaa.gov:9091/dods/NCEP_NARR_VAR_TEST/narr-a_221_hgtprs"; "2" "2" "hgtprs"
yes? let var1_ = acpcp[y="22.0":"63.0",x="-124.0":"-52.0",d=1,t="01-Jan-1979"]
yes? let var2_ = hgtprs[y="22.0":"63.0",gx=var1_grid_,gy=var1_grid_,x="-124.0":"-52.0",d=2,t="01-Jan-1979"]
yes? set var/units="`acpcp[y="22.0":"63.0",x="-124.0":"-52.0",d=1,t="01-Jan-1979"],return=units`" var1_
yes? set var/units="`hgtprs[y="22.0":"63.0",gx=var1_grid_,gy=var1_grid_,x="-124.0":"-52.0",d=2,t="01-Jan-1979"],return=units`" var2_
yes? let data1_ = var1_
yes? let data2_ = var2_
yes? set var/units="`var1_,return=units`" data1_
yes? set var/units="`var2_,return=units`" data2_
yes? set win/size=.5
yes? set win/asp=0.854166666666667
yes? GO "/tmp/lasgo51741959833.jnl"
**ERROR: dimensions improperly specified: must be a 2D region
contour/over/levels=20/nolabel data2_
contour/over/levels=20/nolabel data2_
Command file, command group, or REPEAT execution aborted

killing...Waiting for program to exit...
program exited...
killing...Waiting for program to exit...
program exited...
**************


Any ideas on causes/fixes?

Adam
<datasets>
 <narr_a_221 name="North American Regional Reanalysis" url="file:dummy.nc" doc="">
  <variables>

   <hgtprs name="anl. geopotential height [gpm] - 29 pressure levels" url="http://nomads4.ncdc.noaa.gov:9091/dods/NCEP_NARR_VAR_TEST/narr-a_221_hgtprs";>
    <link match="/lasdata/grids/narr_a_221_time_lev_lat_lon_grid"/>
   </hgtprs>

   <apcpn name="** 0-3 hr acc total precipitation (nearest grid point)[kg/m^2] " url="http://nomads4.ncdc.noaa.gov:9091/dods/NCEP_NARR_VAR_TEST/narr-a_221_pcp";>
    <link match="/lasdata/grids/narr_a_221_time_lat_lon_grid"/>
   </apcpn>

   <acpcp name="** 0-3 hr acc convective precipitation [kg/m^2] " url="http://nomads4.ncdc.noaa.gov:9091/dods/NCEP_NARR_VAR_TEST/narr-a_221_pcp"; >
    <link match="/lasdata/grids/narr_a_221_time_lat_lon_grid"/>
   </acpcp>

   <pcp name="** 0-3 hr acc total precipitation [kg/m^2] " url="http://nomads4.ncdc.noaa.gov:9091/dods/NCEP_NARR_VAR_TEST/narr-a_221_pcp"; >
    <link match="/lasdata/grids/narr_a_221_time_lat_lon_grid"/>
   </pcp>

   <tmpprs name="air temperature - 29 pressure levels [k] " url="http://nomads4.ncdc.noaa.gov:9091/dods/NCEP_NARR_VAR_TEST/narr-a_221_tmpprs";>
    <link match="/lasdata/grids/narr_a_221_time_lev_lat_lon_grid"/>
   </tmpprs>

   <ugrdprs name="anl u wind [m/s] " url="http://nomads4.ncdc.noaa.gov:9091/dods/NCEP_NARR_VAR_TEST/narr-a_221_uvgrdprs";>
    <link match="/lasdata/grids/narr_a_221_wind_time_lev_lat_lon_grid"/>
   </ugrdprs>

   <vgrdprs name="anl v wind [m/s] " url="http://nomads4.ncdc.noaa.gov:9091/dods/NCEP_NARR_VAR_TEST/narr-a_221_uvgrdprs";>
    <link match="/lasdata/grids/narr_a_221_wind_time_lev_lat_lon_grid"/>
   </vgrdprs>

  </variables>
 </narr_a_221>
</datasets>
<grids>
 <narr_a_221_time_lev_lat_lon_grid>
  <link match="/lasdata/axes/narr_a_221_time"/>
  <link match="/lasdata/axes/narr_a_221_lev"/>
  <link match="/lasdata/axes/narr_a_221_lat"/>
  <link match="/lasdata/axes/narr_a_221_lon"/>
 </narr_a_221_time_lev_lat_lon_grid>
 <narr_a_221_wind_time_lev_lat_lon_grid>
  <link match="/lasdata/axes/narr_a_221_time_wind"/>
  <link match="/lasdata/axes/narr_a_221_lev"/>
  <link match="/lasdata/axes/narr_a_221_lat"/>
  <link match="/lasdata/axes/narr_a_221_lon"/>
 </narr_a_221_wind_time_lev_lat_lon_grid>
 <narr_a_221_time_temp_lev_lat_lon_grid>
  <link match="/lasdata/axes/narr_a_221_time_wind"/>
  <link match="/lasdata/axes/narr_a_221_lev"/>
  <link match="/lasdata/axes/narr_a_221_lat"/>
  <link match="/lasdata/axes/narr_a_221_lon"/>
 </narr_a_221_time_temp_lev_lat_lon_grid>
 <narr_a_221_time_lat_lon_grid>
  <link match="/lasdata/axes/narr_a_221_time"/>
  <link match="/lasdata/axes/narr_a_221_lat"/>
  <link match="/lasdata/axes/narr_a_221_lon"/>
 </narr_a_221_time_lat_lon_grid>
</grids>
<axes>
 <narr_a_221_time type="t" units="hour">
  <arange start="1979-01-01 00:00:00" step="3" size="75960"/>
 </narr_a_221_time>
 <narr_a_221_time_temp type="t" units="hour">
  <arange start="1979-01-01 00:00:00" step="3" size="73040"/>
 </narr_a_221_time_temp>
 <narr_a_221_time_wind type="t" units="hour">
  <arange start="1979-01-01 00:00:00" step="3" size="75976"/>
 </narr_a_221_time_wind>
 <narr_a_221_lev type="z" units="millibar">
  <v>1000
  </v>
  <v>975
  </v>
  <v>950
  </v>
  <v>925
  </v>
  <v>900
  </v>
  <v>875
  </v>
  <v>850
  </v>
  <v>825
  </v>
  <v>800
  </v>
  <v>775
  </v>
  <v>750
  </v>
  <v>725
  </v>
  <v>700
  </v>
  <v>650
  </v>
  <v>600
  </v>
  <v>550
  </v>
  <v>500
  </v>
  <v>450
  </v>
  <v>400
  </v>
  <v>350
  </v>
  <v>300
  </v>
  <v>275
  </v>
  <v>250
  </v>
  <v>225
  </v>
  <v>200
  </v>
  <v>175
  </v>
  <v>150
  </v>
  <v>125
  </v>
  <v>100
  </v>
 </narr_a_221_lev>
 <narr_a_221_lat type="y" units="degrees_north">
  <arange start="0" step="0.375" size="240"/>
 </narr_a_221_lat>
 <narr_a_221_lon type="x" units="degrees_east">
  <arange start="-220" step="0.375" size="586"/>
 </narr_a_221_lon>
</axes>



[Thread Prev][Thread Next][Index]

Dept of Commerce / NOAA / OAR / PMEL / TMAP
Contact Us | Privacy Policy | Disclaimer | Accessibility Statement