Re: Library error?

Hey Mary
No, I'm running Ubuntu. I think I need to download a version that isn't
OpeNDap enabled.
I'm running version 5.1.1. I'll download 5.2.0 and see if that sorts it out.
Thanks
Brendan

On Mon, Jul 19, 2010 at 2:10 PM, Mary Haley <haley@ucar.edu> wrote:

> Hi Brendan,
>
> By "PC" do you mean a Windows system running Cygwin? If so, it's possible
> that your newer PC doesn't have the
> latest OPeNDAP libraries installed. Also, what version of NCL are you
> running? You can type:
>
> ncl -V
>
> to get this information. If you are running 5.1.x or earlier, I highly
> recommend upgrading to V5.2.0. Or, if you can wait
> a week or so, we'll be releasing V5.2.1.
>
> The reference to "libdap" makes me believe this is an older version of NCL.
> Version 5.2.0 shouldn't have a reference
> to this library.
>
> If you are already running V5.2.0, then I would recommend downloading the
> latest X/Cygwin "setup.exe" program,
> and making sure you have the latest versions of the external software
> needed.
>
> More information on getting "setup.exe" and installing the Cygwin software
> you need can be found at:
>
> http://www.ncl.ucar.edu/Download/cygwin.shtml#InstallCygwinX
>
> --Mary
>
>
> On Jul 19, 2010, at 2:58 AM, Brendan Argent wrote:
>
> Below is the script I'm stuck with.
> Keep getting the error:
>
> terminate called after throwing an instance of 'libdap::Error'
> Abort
>
> I've checked all the paths so don't thinks its that.
> It ran fine on another pc so it probably has something to do with the way
> NCL is set up on this PC, just not sure where to look?
>
>
> ;program to pair up observed variables from station files with wrf data.
> ;currently set for a short period from the george station with no header
> ;wrf data files are daily 48 hour projections which I split into first and
> second day forecasts and plot as time series
> ;this domain misses george so im pretending its in the domain since the
> bigger domain only has 3hourly data
> ;this programme finds the nearest gridpoint to george
> ;produces arrays for observed, 24hforecast and 48hforecast windspeed and a
> timestamp
>
>
> load "$NCARG_ROOT/share/ncarg/nclscripts/csm/contributed.ncl"
> load "$NCARG_ROOT/share/ncarg/nclscripts/csm/gsn_code.ncl"
> ;loads the necessary functions
> load "$NCARG_ROOT/share/ncarg/nclscripts/csm/gsn_csm.ncl"
> load "$NCARG_ROOT/share/ncarg/nclscripts/wrf/WRFUserARW.ncl"
>
> begin
> ;******Read in station data***********
> datas = asciiread("/media/My
> Passport/NCDC/688280",(/5012,1/),"string") ;open station
> data
> dataf = asciiread("/media/My Passport/NCDC/688280",(/5012,14/),"float")
> ;******Recode the time stamp**********
> timestamp =new((/5012,1/),"string")
> do d=0,5011
> time = stringtocharacter(datas(d,0))
> year = charactertostring(time(0:3))
> month = charactertostring(time(4:5))
> day = charactertostring(time(6:7))
> hour = charactertostring(time(8:9))
> name = year+"-"+month+"-"+day+"_"+hour
> ;print(name)
> timestamp(d,:) = name
> end do
>
> ;******Store the variable***********
> stationid = (dataf(0,1))
> latitude = (dataf(0,2))
> longitude = (dataf(0,3))
> elevation = (dataf(0,4))
>
> ppt = (dataf(:,5))
> P= (dataf(:,6))
> T= (dataf(:,7))
> Td= (dataf(:,8))
> wspd= (dataf(:,9))
> direction= (dataf(:,10))
> u= (dataf(:,11))
> v= (dataf(:,12))
>
> ;do n=0,10
> ;print (timestamp(n,0)+" "+wspd(n))
> ;print (timestamp(0:71,0))
>
> ;end do
> ;******Find the closest point in WRF**
> opt = True
> opt@MAP_PROJ = 1
> opt@TRUELAT1 = -30.0
> opt@TRUELAT2 = -60.0
> opt@STAND_LON = 20.0
> opt@REF_LAT = -36.7366
> opt@REF_LON = 14.477
> opt@KNOWNJ = 0
> opt@KNOWNI = 0
> opt@DX = 7500.
> opt@DY = 7500.
> loc = wrf_ll_to_ij( longitude, latitude, opt )
> ni =floattointeger(loc(0))
> nj =floattointeger(loc(1))
> ;print (ni + " "+ nj)
>
> ;******Loop through the WRF files*****
> File = "/media/My Passport/SAfrica/wrfout_d03_" ;get
> data
> periods=5012
> firstdayforecast =new((/periods,1/),"float")
> seconddayforecast =new((/periods,1/),"float")
> d = -1
>
> years = 2009
> do months = 9,10
>
> if (months .gt. 12) then
> months = months-12
> end if
>
> do days=1,5
> d = d+1
> hours = 12
>
> FName =
> File+sprinti("%0.4i",years)+"-"+sprinti("%0.2i",months)+"-"+sprinti("%0.2i",days)+"_"+sprinti("%0.2i",hours)
> ;check we are getting
> it right
> if (isfilepresent(FName)) then ;chceks if
> file exists
> ff = addfile(FName+".nc","r") ;ff is now
> the files
> Time = -1 ;takes all times
> ua = wrf_user_getvar(ff,"ua",Time) ; u on mass points
> va = wrf_user_getvar(ff,"va",Time) ; v on mass points
>
> Wwspd = sqrt(ua*ua + va*va)
> copy_VarMeta (ua,Wwspd) ;replaces
> metadata
> Wwspd@description = "wind speed"
> ;corrects the variable name
> ;printVarSummary(Wwspd)
> ;printMinMax(Wwspd,True) ;a useful
> data check
>
> if (days.eq.1) then ;if statement
> lat = ff->XLAT ;quick way to
> grab a variable
> lon = ff->XLONG
> Wwspd@lat2d = lat ;assigns
> coordinates
> Wwspd@lon2d = lon
> end if
>
> firstday = Wwspd(0:24,0,ni,nj)
> secondday = Wwspd(24:48,0,ni,nj)
>
> do t =0,24
> firstdayforecast(24*d+t,:) = firstday(t)
> seconddayforecast(24+24*d+t,:) = secondday(t)
> end do
>
> end if
>
> end do
> end do
>
> ;******Print Pairs to a file*******
> datamatrix =new((/5012,4/),"string")
> print ("Timestamp"+ " "+"obs wspd"+" "+ "24hf wspd"+" " +"48hf
> wsdp")
> do n=0,5011
> print(timestamp(n,0)+" "+wspd(n)+" "+firstdayforecast(n,0)+"
> "+seconddayforecast(n,0))
> datamatrix(n,0)= timestamp(n,0)
> datamatrix(n,1) = wspd(n)
> datamatrix(n,2) = flt2string(firstdayforecast(n,0))
> datamatrix(n,3) = flt2string(seconddayforecast(n,0))
> end do
>
> opt = True
> opt@fout = "Pairs.dat"
> write_matrix (datamatrix(1:3,:), "f3.2" , opt)
>
>
> end
>
> On Fri, Jul 16, 2010 at 9:37 PM, Rick Brownrigg <brownrig@ucar.edu> wrote:
>
>> Hi Brendan,
>>
>> It sounds like perhaps you are trying to open a remote file via OPeNDAP,
>> and the read of that file is failing? Its likely not the copy_VarMeta
>> function per se. It might be useful if you can post your script.
>>
>> By the way, not that it should make a difference in the libdap:Error
>> outcome, but in your attempt to replace
>>
>>
>> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
>>
>> with
>>
>>
>> load "./usr/share/ncarg/nclscripts/csm/contributed.ncl"
>>
>>>
>>> Note that the leading "." makes that a relative path, not the absolute
>> path I think you intended.
>>
>> FWIW...
>> Rick
>>
>>
>> On Jul 16, 2010, at 3:37 AM, Brendan Argent wrote:
>>
>> Hey all
>>> I'm quite a rookie with NCL so any help would be welcome.
>>> I've been working on a script on a computer overseas and brought it back
>>> with me only to find a bunch of error messages I wasn't getting before.
>>> Specifically:
>>>
>>> terminate called after throwing an instance of 'libdap::Error'
>>> Abort
>>>
>>> I think it has to do with the libraries being in a different place but
>>> I'm not sure how to tell NCL where to find the functions I need.
>>>
>>> One of the functions I am using is copy_VarMeta. This requires the line
>>>
>>> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
>>>
>>> If I comment this line out then I get the error:
>>>
>>> fatal:syntax error: line 95 in file
>>> Pairs3.ncl before or near ;replaces metadata
>>>
>>> copy_VarMeta (ua,Wwspd)
>>> ;replaces metadata
>>> ^
>>>
>>> fatal:syntax error: possibly an
>>> undefined procedure
>>> fatal:Syntax Error in block, block not
>>> executed
>>> fatal:error at line 136 in file
>>> Pairs3.ncl
>>>
>>> SO then the program actually runs which makes me think the problem is
>>> with loading the function.
>>>
>>> I've checked the path to the contributed.ncl function
>>> (/usr/share/ncarg/nclscripts/csm/contributed.ncl)
>>> but if I change the suggested line to :
>>> load "./usr/share/ncarg/nclscripts/csm/contributed.ncl"
>>>
>>> I get the error:
>>> fatal:Could not open
>>> (./usr/share/ncarg/nclscripts/csm/contributed.ncl)
>>> fatal:error at line 9 in file Pairs3.ncl
>>>
>>> As I said I'm quite new to this so maybe I'm missing something obvious.
>>> Any help would be appreciated.
>>> Thanks
>>> Brendan
>>> _______________________________________________
>>> 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 Mon Jul 19 06:17:32 2010