Re: rcm2points

Hi Dennis: Thank you very much for your help!  It is quite odd. For the information shown below, I used Version 4.2.0.a033. If I set lat =  (/-22.3,-22.31/) and lon = (/166.4,166.41/), it works. If I set lat = (/-22.3/) and lon = (/166.4/), it has the following error infor.. Actually, I tested this with two different ncl scripts and data. The results are the same. Could this be a lower version problem? Thanks! Regards Ed  Copyright (C) 1995-2006 - All Rights Reserved  University Corporation for Atmospheric Research  NCAR Command Language Version 4.2.0.a033  The use of this software is governed by a License Agreement.  See http://www.ncl.ucar.edu/ for more details. fatal:rcm2points: The input lat/lon grids must have at least 2 elements, and the output lat/lon arrays 1 element fatal:Execute: Error occurred at or near line 64 in file yytest.ncl In my script, ;  lat =  (/-22.3,-22.31/)  ; works ;  lon = (/166.4,166.41/) ; works lat = (/-22.3/)                 ; does not work lon = (/166.4/)                ; does not work .......... ul  = rcm2points (tlat, tlon, u, lat, lon, 1)  ; line 64 ....................................... ----- Original Message ---- From: Dennis Shea <shea@ucar.edu> To: Yang Yang <yang123yang@yahoo.com> Cc: ncl-talk@ucar.edu Sent: Friday, June 13, 2008 2:33:13 AM Subject: Re: [ncl-talk] rcm2points Yang Yang wrote: > Hello: > I ran rcm2points by setting lat and lon with 2 or more elements without any problem. > However, if I set lat and lon with only one element, the ncl script failed with the following error information: > > fatal:rcm2points: The input lat/lon grids must have at least 2 elements, and the output lat/lon arrays 1 element > > How to solve this problem? Thanks! > Regards > Ed >  The following simple test script works fine. ================================= load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl" begin   nt    = 0   kl    = 20   a    = addfile("/project/cas/shea/WRF/wrfout_d01_2003-07-13_12:00:00.nc","r")   x    = a->T(nt,kl,:,:)   lat2d = a->XLAT(nt,:,:)      ; min=  16.789  max= 34.542   lon2d = a->XLONG(nt,:,:)    ; min=-105.153  max=-82.847   latPts = (/22., 29./)   lonPts = (/-100., -93./)   xNew  = rcm2points(lat2d, lon2d, x, latPts, lonPts, 1)   print(xNew) ;  Test Ed's ncl-talk question   yPts = (/22./)   xPts = (/-100./)   xTst  = rcm2points(lat2d, lon2d, x, yPts, xPts, 1)   print(xTst) end ================== Variable: xNew Dimensions and sizes:  [2] Coordinates: 0)    39.80119 (1)    39.21795 ================== Variable: xTst Dimensions and sizes:  [1] (0)    39.80119
_______________________________________________
ncl-talk mailing list
ncl-talk_at_ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Thu Jun 12 2008 - 15:34:49 MDT