coordinate subscripting

I am trying to pull out a smaller data set from a larger one via
coordinates. I am not sure what the error means because I have used
the brackets correctly for coordinate subscripting.
How do I get around this? ( I eventually will write out to an ascii
file).
-Erik

ERROR: fatal:Dimension sizes of left hand side and right hand side of
assignment do not match
fatal:Execute: Error occurred at or near line 26 in file FNL-ascii.ncl

Terminal:

Erik_at_noble:/Volumes/Data_and_Models/ncl_scripts:ncl FNL-ascii.ncl
  Copyright (C) 1995-2007 - All Rights Reserved
  University Corporation for Atmospheric Research
  NCAR Command Language Version 5.0.0
  The use of this software is governed by a License Agreement.
  See http://www.ncl.ucar.edu/ for more details.

Variable: v700
Type: float
Total Size: 260640 bytes
             65160 values
Number of Dimensions: 2
Dimensions and sizes: [lat_3 | 181] x [lon_3 | 360]
Coordinates:
             lat_3: [90..-90]
             lon_3: [ 0..359]
Number Of Attributes: 13
   lv_ISBL3 : 700
   center : US National Weather Service - NCEP (WMC)
   long_name : v-component of wind
   units : m/s
   _FillValue : -999
   level_indicator : 100
   grid_number : 3
   parameter_table_version : 2
   parameter_number : 34
   model : Spectral Statistical Interpolation (SSI) analysis from
"Final" run.
   forecast_time : 0
   forecast_time_units : hours
   initial_time : 09/10/2006 (00:00)

Variable: v700
Type: float
Total Size: 260640 bytes
             65160 values
Number of Dimensions: 2
Dimensions and sizes: [lat_3 | 181] x [lon_3 | 360]
Coordinates:
             lat_3: [90..-90]
             lon_3: [-180..179]
Number Of Attributes: 14
   lonFlip : longitude coordinate variable has been reordered via
lonFlip
   initial_time : 09/10/2006 (00:00)
   forecast_time_units : hours
   forecast_time : 0
   model : Spectral Statistical Interpolation (SSI) analysis from
"Final" run.
   parameter_number : 34
   parameter_table_version : 2
   grid_number : 3
   level_indicator : 100
   units : m/s
   long_name : v-component of wind
   center : US National Weather Service - NCEP (WMC)
   lv_ISBL3 : 700
   _FillValue : -999
fatal:Dimension sizes of left hand side and right hand side of
assignment do not match
fatal:Execute: Error occurred at or near line 24 in file FNL-ascii.ncl

Erik_at_noble:/Volumes/Data_and_Models/ncl_scripts:

Code:
;*************************************************
; FNL_Plot.ncl
;************************************************
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"
;************************************************
begin
;************************************************
; read in netCDF file
;************************************************
  a = addfile("/Volumes/Data_and_Models/FNL_SOP3/FNL_09_grb/
fnl_060910_00_00.grb","r")
   title1 = "FNL Reanalysis - 10 AUG 2006 - 00 UTC"
   type = "x11" ; type = "pdf"
   wks = gsn_open_wks(type,"FNL-0910-12UTC_SLP_Hgt_Panel")
   ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
   v700 = a->V_GRD_3_ISBL_10(17,:,:) ; 700mb =level 17
         printVarSummary( v700 )
        v700 = lonFlip( v700 )
         lon2d_Flip=v700&lon_3
                printVarSummary( v700 )
v700 = v700({-25:35},{-35:35})

end

_______________________________________________
ncl-talk mailing list
ncl-talk_at_ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Thu Nov 15 2007 - 18:14:28 MST