Re: Coordinate scripting for

From: Erik Noble <enoble_at_nyahnyahspammersnyahnyah>
Date: Tue, 11 Dec 2007 14:26:28 -0500

Hi Soyoung Ha.
Thank you.
This helps.

Now I have a different error (Below):
fatal:Coordinate variables must have one dimension only
fatal:No coordinate variable exists for dimension (lat) in variable (v_plane)
Bus error

This is all WRF data. My inquiry is related to the earlier inquiry
addressing interpolating WRF pressure-level data.
I interpolated pressure level data using the WRF_user_interp3d
command. The result of using the WRF_NCL code presents the
pressure-level data in 109 rows and 139 columns; It used to have
coordinates before using the ncl WRF_user_interp3d command.
Variable: v_plane
Type: float
Total Size: 60604 bytes
Number of Dimensions: 2
Dimensions and sizes: [109] x [139]
Coordinates:

I am simply trying to extract a smaller subset from a larger set of
data. Since ncl WRF_user_interp3d command took away coordinates, I
just want data in between rows 37-55 and columns 11-89.

My code keeps getting longer and longer. Yet with these errors, I am
having to resort to converting the large [109] x [139] files to ascii
and creating the the [37-55] x [11-89] subsets through Microsoft
Excel, which takes a while.
Any suggestions here?
-Erik

On Dec 11, 2007 1:57 PM, Soyoung Ha <syha_at_ucar.edu> wrote:
> Erik,
>
> You should define lat/lon info first.
> Add these lines at the 19th line, and try it again.
>
> rank = dimsizes(filevardimsizes(a,"XLAT")) ; # of dimensions
> if (rank.eq.3) then
> lat = a->XLAT(0,:,:) ; opt could bt "nt" f->XLAT(opt,:,:)
> lon = a->XLONG(0,:,:)
> else
> if (rank.eq.2) then
> lat = a->XLAT
> lon = a->XLONG
> else
> print ("Error: unexpected lat/lon rank: rank="+rank)
> exit
> end if
> end if
>
> Soyoung
>
>
>
> Erik Noble wrote:
> > Thank you Dr. Shea,
> > M variable has 109 rows and 139 columns.
> > Does that change things?
> >
> > i tried the suggestions from the two pages you gave below by adding
> >
> > v_plane = wrf_user_intrp3d( v,p,"h",pressure,0.,False)
> > printVarSummary(v_plane)
> > v_plane!0 = "lat"
> > v_plane!1 = "lon"
> > v_plane&lat = lat
> > v_plane&lon = lon
> > V = v_plane({5:15},{-30:10})
> > printVarSummary(V)
> >
> >
> > I am now getting a different error:
> >
> > Variable: v_plane
> > Type: float
> > Total Size: 60604 bytes
> > 15151 values
> > Number of Dimensions: 2
> > Dimensions and sizes: [109] x [139]
> > Coordinates:
> > Number Of Attributes: 4
> > description : v Velocity
> > units : m/s
> > _FillValue : -999999
> > PlotLevelID : 850 hPa
> > fatal:Variable (lat) is undefined
> > fatal:Execute: Error occurred at or near line 39 in file NCL_ascii_V700.ncl
> >
> > On Dec 11, 2007 1:13 PM, Dennis Shea <shea_at_ucar.edu> wrote:
> >
> >> [1] To use "coordinate subscripting" two conditions must be present:
> >> [a] the dimensions must be named. The arrays below do not have
> >> named dimensions.
> >> [b] The named dimensions must have 1D [not 2D] monotonically
> >> {in/de}creasing coordinates associated with each named
> >> dimension.
> >>
> >> Note: These are netCDF rules and NCL adheres to the netCDF
> >> variable model.
> >> [2] If you have variables that have 2d variable arrays that contain
> >> coordinates (like WRF XLAT and XLONG variables] you
> >> can not use coordinate subscripting [See 1b above ]
> >>
> >> [3] Pls read
> >> http://www.ncl.ucar.edu/Document/Language/cv.shtml
> >> http://www.ncl.ucar.edu/Document/Manuals/Ref_Manual/NclVariables.shtml
> >>
> >> Erik Noble wrote:
> >>
> >>> Hi.
> >>> I added a printVarSummary(v_plane) to the code to see the error more clearly.
> >>> The error still remains (below).
> >>>
> >>> any suggestions?
> >>> -erik
> >>>
> >>>
> >>> Variable: v_plane
> >>> Type: float
> >>> Total Size: 60604 bytes
> >>> 15151 values
> >>> Number of Dimensions: 2
> >>> Dimensions and sizes: [109] x [139]
> >>> Coordinates:
> >>> Number Of Attributes: 4
> >>> description : v Velocity
> >>> units : m/s
> >>> _FillValue : -999999
> >>> PlotLevelID : 850 hPa
> >>> fatal:Dimension (1) of (v_plane) is not named and therefore doesn't
> >>> have an associated coordinate variable
> >>> fatal:Execute: Error occurred at or near line 37 in file NCL_ascii_V700.ncl
> >>>
> >>>
> >>> On Dec 11, 2007 12:31 PM, Erik Noble <enoble_at_giss.nasa.gov> wrote:
> >>>
> >>>
> >>>> Dear NCL,
> >>>> I am trying to use coordinate subscripting to pull data from an array:
> >>>> V = v_plane({5:15},{-30:10})
> >>>> I am getting the error below.
> >>>> My code is attached.
> >>>> My I have some help?
> >>>> Thank you for your help.
> >>>> -Erik
> >>>>
> >>>>
> >>>> Erik_at_noble:/Volumes/Data_and_Models/ncl_scripts:ncl NCL_ascii_V700.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: times
> >>>> Type: string
> >>>> Total Size: 1948 bytes
> >>>> 487 values
> >>>> Number of Dimensions: 1
> >>>> Dimensions and sizes: [487]
> >>>> Coordinates:
> >>>> Number Of Attributes: 2
> >>>> description : times in file
> >>>> _FillValue : missing
> >>>>
> >>>> Variable: v
> >>>> Type: float
> >>>> Total Size: 1636308 bytes
> >>>> 409077 values
> >>>> Number of Dimensions: 3
> >>>> Dimensions and sizes: [27] x [109] x [139]
> >>>> Coordinates:
> >>>> Number Of Attributes: 2
> >>>> description : v Velocity
> >>>> units : m/s
> >>>>
> >>>>
> >>>> Variable: u
> >>>> Type: float
> >>>> Total Size: 1636308 bytes
> >>>> 409077 values
> >>>> Number of Dimensions: 3
> >>>> Dimensions and sizes: [27] x [109] x [139]
> >>>> Coordinates:
> >>>> Number Of Attributes: 2
> >>>> description : u Velocity
> >>>> units : m/s
> >>>>
> >>>>
> >>>> Variable: u
> >>>> Type: float
> >>>> Total Size: 1636308 bytes
> >>>> 409077 values
> >>>> Number of Dimensions: 3
> >>>> Dimensions and sizes: [27] x [109] x [139]
> >>>> Coordinates:
> >>>> Number Of Attributes: 2
> >>>> description : u Velocity
> >>>> units : m/s
> >>>> fatal:Dimension (1) of (v_plane) is not named and therefore doesn't have an
> >>>> associated coordinate variable
> >>>> fatal:Execute: Error occurred at or near line 37 in file NCL_ascii_V700.ncl
> >>>>
> >>>>
> >>>>
> >>>> Erik_at_noble:/Volumes/Data_and_Models/ncl_scripts:
> >>>>
> >>>>
> >>>>
> >>>>
> >>
> >>
> >>
> >> ------------------------------------------------------------------------
>
> >>
> >> _______________________________________________
> >> ncl-talk mailing list
> >> ncl-talk_at_ucar.edu
> >> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
> >>
>
>
>
_______________________________________________
ncl-talk mailing list
ncl-talk_at_ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Tue Dec 11 2007 - 12:26:28 MST

This archive was generated by hypermail 2.2.0 : Tue Dec 11 2007 - 14:19:21 MST