Dear Erik,
I think you try to store the smaller data set in larger data variable in
your program. You have already assigned the dimension of v700 as [lat_3 |
181] x [lon_3 | 360].
Hence you need to store the smaller data in a new variable to get rid the
dimension size mismatch error.
So try
new_v700 = v700({-25:35},{-35:35})
Good luck
Debasish
_____
From: ncl-talk-bounces_at_ucar.edu [mailto:ncl-talk-bounces_at_ucar.edu] On Behalf
Of Erik Noble
Sent: Thursday, November 15, 2007 4:14 PM
To: ncl-talk_at_ucar.edu
Subject: 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:46:29 MST
This archive was generated by hypermail 2.2.0 : Sun Nov 25 2007 - 07:21:25 MST