Re: coordinate subscripting

Hi,
Thank you. I see my mistake.
-Erik
On Nov 15, 2007, at 8:46 PM, Debasish Pai Mazumder wrote:

> 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:23 MST