[Thread Prev][Thread Next][Index]

Re: polar stereographic projection



Hi, Ansley,

Thanks. When I ran the script with your suggestions, I got the following
error messages:

 **ERROR: invalid command: grave accent doesnt evaluate to scalar
          x_page[i=@max,j=@max]
DEFINE VARIABLE/quiet x_axis_range = `x_page[i=@max,j=@max]` - `x_page[i=@min,j=@min]`
Command file, command group, or REPEAT execution aborted

Could you please to figure out what is wrong? Look foreward to your more 
suggestions.  Thanks.

Cheers, Xiangdong


> 
> Hi Xiangdong,
> You can do this kind of thing, if your variables are v1, lon, lat:
> 
> 
> ! Run the mp_stereographic_north script to define projection variables
> set region/y=60:90N
> go mp_stereographic_north
> 
> ! Reset the map projection coordinates to use the lon, lat grid
> 
> let/quiet mp_lambda = lon * deg2rad
> let/quiet mp_phi = lon* deg2rad
> set grid lon
> go mp_aspect
> 
> shade/noaxes v1, x_page, y_page
> 
> 
> Xiangdong Zhang wrote:
> 
> >Thanks. In this case, both plots are on the rectangular map. I hope
> >the plot will be on the polar stereographic projection map. More
> >suggestions? Thanks again -- Xiangdong
> >
> >  
> >
> >>You can use Ferret's curvilinear mapping procedures for this.  Let's  
> >>suppose you want to contour variable v1 which is on a rectangular i,j  
> >>grid that is not necessarily a global lat,lon grid, but you also have  
> >>lat and lon on that same rectangular grid.  Then the following will do  
> >>the job
> >>
> >>contour v1,lon,lat
> >>
> >>Your other data set with variable v2, already on a lon/lat grid, is  
> >>contoured by
> >>
> >>contour v2
> >>------------------------------------------------------------------------ 
> >>------------------
> >>Edward D. (Ned) Cokelet, Ph.D.	Edward.D.Cokelet@noaa.gov
> >>NOAA/PMEL					off:  (206) 526-6820
> >>7600 Sand Point Way NE			fax: (206) 526-6485
> >>Seattle, WA 98115-6439
> >>
> >>
> >>    
> >>
> >
> >  
> >
> 
> --------------020800080405000107010404
> Content-Type: text/html; charset=us-ascii
> Content-Transfer-Encoding: 7bit
> 
> <!DOCTYPE html PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">
> <html>
> <head>
>   <meta http-equiv="Content-Type" content="text/html;charset=ISO-8859-1">
>   <title></title>
> </head>
> <body text="#000000" bgcolor="#ffffff">
> Hi Xiangdong,<br>
> You can do this kind of thing, if your variables are v1, lon, lat:<br>
> <br>
> <br>
> ! Run the mp_stereographic_north script to define projection variables<br>
> set region/y=60:90N<br>
> go mp_stereographic_north<br>
> <br>
> ! Reset the map projection coordinates to use the lon, lat grid<br>
> <br>
> let/quiet mp_lambda = lon * deg2rad<br>
> let/quiet mp_phi = lon* deg2rad<br>
> set grid lon<br>
> go mp_aspect<br>
> <br>
> shade/noaxes v1, x_page, y_page<br>
> <br>
> <br>
> Xiangdong Zhang wrote:<br>
> <blockquote type="cite"
>  cite="mid20050212060428.A9E5140D28D@salix.frontier.iarc.uaf.edu">
>   <pre wrap="">Thanks. In this case, both plots are on the rectangular map. I hope
> the plot will be on the polar stereographic projection map. More
> suggestions? Thanks again -- Xiangdong
> 
>   </pre>
>   <blockquote type="cite">
>     <pre wrap="">You can use Ferret's curvilinear mapping procedures for this.  Let's  
> suppose you want to contour variable v1 which is on a rectangular i,j  
> grid that is not necessarily a global lat,lon grid, but you also have  
> lat and lon on that same rectangular grid.  Then the following will do  
> the job
> 
> contour v1,lon,lat
> 
> Your other data set with variable v2, already on a lon/lat grid, is  
> contoured by
> 
> contour v2
> ------------------------------------------------------------------------ 
> ------------------
> Edward D. (Ned) Cokelet, Ph.D.	<a class="moz-txt-link-abbreviated" href="mailto:Edward.D.Cokelet@noaa.gov";>Edward.D.Cokelet@noaa.gov</a>
> NOAA/PMEL					off:  (206) 526-6820
> 7600 Sand Point Way NE			fax: (206) 526-6485
> Seattle, WA 98115-6439
> 
> 
>     </pre>
>   </blockquote>
>   <pre wrap=""><!---->
>   </pre>
> </blockquote>
> </body>
> </html>
> 
> --------------020800080405000107010404--
> 
> 


[Thread Prev][Thread Next][Index]

Dept of Commerce / NOAA / OAR / PMEL / TMAP

Contact Us | Privacy Policy | Disclaimer | Accessibility Statement