need help (regional boundary)

[Mac MacLeod]

Hi Alessandro.


We have such a program as you require.  It is a "C" Program that takes the
lat and long of boundary points and encodes them in a binary file which you
then attach by set mpdset  and control through  mpt commands in the script.

Here is the program attached with a separate instruction file. The
makemaporig.c is the C code that must be compiled before use.  In the
program before compiling you set the : #define MPTID 24  in the program to
#define MPTID ??  where ?? = the number you want to use to specify the
background number in the mpt command.  Note do not use 1 or 2 as these are
used by GrADS for coasts and national boundaries.

Good luck.

Mac MacLeod

M.A.(Mac) MacLeod
President and General Manager,
Scotia Weather Services Inc.
192 Wyse Road, Suite 8,
Dartmouth, Nova Scotia, B3A 1M9
Tele: 902-468-3866
Cell: 902-225-5747
Fax: 902-461-1768
E-mail: mmacleod at scotiaweather.com
visit us at: www.scotiaweather.com

----- Original Message -----
From: "Alessandro Raolil" <raolil at LE.ISAC.CNR.IT>
To: <GRADSUSR at LIST.CINECA.IT>
Sent: Thursday, July 06, 2006 4:18 PM
Subject: Re: need help (regional boundary)


> Thanks a lot for your response, but i've already implemented by myself a
> script many similar. My problem is that the script is too slow to draw all
> italian regional boundary ; so i need a script that produce a binary file,
> like "hires", faster and optimized, like the one proposed by Francis in
> his previous mail.
>
> Regards,
> Alessandro
>
> > Hello Alessandro,
> >
> > Here is a simple script snippet that you can put in a
> > grads script that will draw a boundary using a lonlat
> > coordinate pairs from an ascii file. It is not elegant
> > code, but it works.
> > The lonlatfile format is one lon/lat coordinate pair
> > on each line. Maxpts is the total number of coordinate
> > points. If you  want a closed boundary, make the last
> > coordinate pair in the ascii file the same as the
> > first coordinate pair.
> > ------------------------------------------------
> > *Read lon lat coords of boundary points,calculate xy
> > coords and plot points
> >
> > count = 1
> > while (count<=maxpts)
> >   llpoint= sublin(read(lonlatfile),2)
> >
> >   ln = subwrd(llpoint,1)
> >
> >   lt = subwrd(llpoint,2)
> >
> >   'q w2xy 'ln' 'lt''
> >
> >   x1=subwrd(result,3)
> >   y1=subwrd(result,6)
> >
> > if count=1
> >   oldx=x1
> >   oldy=y1
> > else
> >   'draw line 'oldx' 'oldy' 'x1' 'y1''
> > endif
> >   oldx=x1
> >   oldy=y1
> >   count=count+1
> > endwhile
> >
> > --------------------------------------------------
> > This might work for your purposes.
> >
> > Regards,
> > Joe Marlow
> > Tohono O'odham Community College
> > jmarlow at tocc.cc.az.us
> > jemz57 at yahoo.com
> >
> >
> > --- Alessandro Raolil <raolil at LE.ISAC.CNR.IT> wrote:
> >
> >> Hi Francis,
> >> i have the same problem of Vikram, but for italian
> >> regions. Can you send
> >> me the correct code? i have all the lon/lat
> >> coordinates in ascii format.
> >> Thanks
> >> Alessandro
> >>
> >
> >
> > __________________________________________________
> > Do You Yahoo!?
> > Tired of spam?  Yahoo! Mail has the best spam protection around
> > http://mail.yahoo.com
> >
> >
>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: mpt.doc
Type: application/msword
Size: 2349 bytes
Desc: not available
Url : http://gradsusr.org/pipermail/gradsusr/attachments/20060706/c7e3f313/attachment.doc 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: map.doc
Type: application/msword
Size: 6133 bytes
Desc: not available
Url : http://gradsusr.org/pipermail/gradsusr/attachments/20060706/c7e3f313/attachment-0001.doc 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: makemaporig.c
Type: application/octet-stream
Size: 5320 bytes
Desc: not available
Url : http://gradsusr.org/pipermail/gradsusr/attachments/20060706/c7e3f313/attachment.obj