Re: track plots

No, that was not the problem. The error didn't occur when loading the
track file but when loading the map, i.e. in this section:

  wks = gsn_open_wks("pdf","traj") ; open workstation
    res = True ; map resources

    res_at_gsnDraw = False ; don't draw
    res_at_gsnFrame = False ; don't advance frame
    res_at_vpWidthF = 0.80 ; make map bigger
    res_at_vpHeightF = 0.80
    res_at_mpMaxLatF = -20 ; select subregion
    res_at_mpMinLatF = -60
    res_at_mpMinLonF = -75
    res_at_mpMaxLonF = -25

    res_at_tiMainString = "Example of a trajectory plot" ; title
    res_at_gsnCenterString = "markers every 4th time step" ; center string

    map = gsn_csm_map_ce(wks,res) ; create map

>
> lat = addfile("$PATH/TRACK_lat_200704.nc","r")
> lon = addfile("$PATH/TRACK_lon_200704.nc","r")
>
>
> As noted by Andrea, the above are incorrect. In fact, they are
> incorrect in two ways: (a) the path and (b) the way
> the file is opened.
>
> path = getenv("PATH")
>
> f = addfile (path+"TRACK_lat_200704.nc","r")
> lat = f->lat
> f = addfile (path+"TRACK_lon_200704.nc","r")
> lon = f->lon
>
> Please note that both "lat" and "lon" should have the missing value
> attribute _FillValue
>
> http://www.ncl.ucar.edu/Document/Functions/Built-in/getenv.shtml
>
> I suggest that you download/read the mini-language manual at:
> http://www.ncl.ucar.edu/Document/Manuals/
>
> Good luck
>
>
> Hahmann, Andrea N. wrote:
>> Christine
>>
>> The correct syntax for what you want to do is:
>>
>> ff = addfile(<lat file name>,"r")
>> lat = ff->lat (or name of field in netcdf file)
>> ff = addfile(<lon file name>,"r")
>> lon = ff->lon (or name of field in netcdf file)
>>
>> You need to open the files first, then read from them.
>> Hope this helps,
>>
>> Andrea
>>
>> Andrea Hahmann, PhD.
>> Senior Scientist
>> Phone direct +45 4677 5471
>> andrea.n.hahmann_at_risoe.dk
>>
>> Wind Energy Department
>> Risø National Laboratory for Sustainable Energy
>> Technical University of Denmark
>> Building 118, P.O. Box 49
>> DK-4000 Roskilde, Denmark
>> Tel +45 4677 4677
>> Fax +45 4677 5083
>>
>> -----Original Message-----
>> From: ncl-talk-bounces_at_ucar.edu [mailto:ncl-talk-bounces_at_ucar.edu]
>> On Behalf Of Christine.Radermacher_at_gfi.uib.no
>> Sent: Wednesday, April 22, 2009 4:05 PM
>> To: ncl-talk_at_ucar.edu
>> Subject: track plots
>>
>>
>> Hello,
>>
>> I would like to plot storm tracks with ncl.
>> The tracks are stored in NetCDF and the latitude and longitude
>> variables both have the form
>>
>> 1: lat1 lat2 lat3 missing missing ......
>> 2: lat1 lat2 lat3 lat4 missing missing .....
>> .
>> .
>> .
>>
>> so that each track is filled with missing values to a certain length.
>>
>> I tried to modify the traj_1.ncl script for that problem and modified
>> the beginning like that:
>>
>> ;*************************************************
>> ; traj_1.ncl
>> ;*************************************************
>> load "/usr/share/ncarg/nclscripts/csm/gsn_code.ncl"
>> load "/usr/share/ncarg/nclscripts/csm/gsn_csm.ncl"
>> ;**************************************************
>> begin
>> ;*************************************
>> ; read in data
>> ;************************************
>>
>> lat = addfile("$PATH/TRACK_lat_200704.nc","r")
>> lon = addfile("$PATH/TRACK_lon_200704.nc","r")
>> ntime=105
>>
>> Then I get the following error:
>> warning:TextItem:libncarg Error:PLCHHQ - DATASET NOT LOADED CORRECTLY
>> Segmentation fault
>>
>> I hope you could help me.
>>
>> Regards
>> Christine
>>
>>
>>
>>
>>
>> ----- Ende der weitergeleiteten Nachricht -----
>>
>>
>>
>> ----- Ende der weitergeleiteten Nachricht -----
>>
>>
>>
>> ----- Ende der weitergeleiteten Nachricht -----
>>
>> _______________________________________________
>> 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
>>
>
>
> --
> ======================================================
> Dennis J. Shea tel: 303-497-1361 |
> P.O. Box 3000 fax: 303-497-1333 |
> Climate Analysis Section |
> Climate & Global Dynamics Div. |
> National Center for Atmospheric Research |
> Boulder, CO 80307 |
> USA email: shea 'at' ucar.edu |
> ======================================================
>
> _______________________________________________
> 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 Apr 23 2009 - 01:48:00 MDT