Re: wrf_user_ll_to_ij problem

Hi Chris,

The wrf_user_ll_to_ij function calls an internal built-in function called "wrf_ll_to_ij". This is the one that appears to be undefined.

What version of NCL are you running? This internal function was added in version 5.1.0. You'll need to upgrade if you have an earlier version.

--Mary

On Jun 28, 2011, at 1:36 PM, Chris Klich wrote:

> I have been trying to start a script to create a cross section, and I am trying to use the function wrf_user_ll_to_ij to extract the i,j locations for specific lat and lon locations. However, I am getting the error:
>
> Variable: times
> Type: string
> Total Size: 4 bytes
> 1 values
> Number of Dimensions: 1
> Dimensions and sizes: [1]
> Coordinates:
> Number Of Attributes: 2
> description : times in file
> _FillValue : missing
> (0) 2008-04-21_00:00:00
> fatal:Undefined identifier: (wrf_ll_to_ij) is undefined, can't continue
> fatal:Execute: Error occurred at or near line 1456 in file ./WRFUserARW.ncl
>
> fatal:Execute: Error occurred at or near line 30
>
>
> I've checked the file and the function is there. Note that the WRFUserARW.ncl file is in my script directory as I've downloaded the newest version of it, and not the $NCARG_ROOT directory. Other WRF functions seem to work, but this one I am having problems with. Any help or suggestions would be appreciated.
>
> Thank you,
> Chris
>
> ------------------------------------------------------------------------------------------------------------------------------------------------------------------
>
> 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/wrf/WRF_contributed.ncl"
> load "./WRFUserARW.ncl"
>
> begin
>
>
> a = addfile("/sully/r0/cklich/WRF/DATA_FILES/WRFOUT_DATA/2_DOMS_45000_KM_35_LEVS_2008_0418/wrfout_d02_2008-04-21_00:00:00.nc","r")
>
> times = wrf_user_list_times(a) ; get time(s) in the file
> ntimes = dimsizes(times)
>
> lat=a->XLAT(0,:,:)
> lon=a->XLONG(0,:,:)
> hgt=a->HGT(0,:,:)
>
> dimll = dimsizes(lat) ; get size of dimensions
> nlat = dimll(0)
> mlon = dimll(1)
>
> start_lat = 27.1
> start_lon = 133.1
> end_lat = 38.1
> end_lon = 142.1
>
> resll = True
> resll@returnInt = True
>
> start_loc = wrf_user_ll_to_ij(a,start_lon,start_lat,resll)
> end_loc = wrf_user_ll_to_ij(a,end_lon,end_lat,resll)
>
> print("start_loc = "+start_loc)
> print("end_loc = "+end_loc)
>
> end
>
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Thu Jun 30 07:37:39 2011