[Thread Prev][Thread Next][Index]

Re: [ferret_users] regriding data from (lon lat) axis to (x_page y_page) axis (stereographic projection)



Dear Cécile,
As i understand you question you need to have your (rectilinear) data set in a curvilinear grid.
You can have it easily by the RECT_TO_CURV function. You only need two curvilinear grid variables (2D fields of lons and lats describing the curvilinear grid) from an other, curvilinear file or from yours:

go mp_stereographic_south
set grid your_old
!now you have variables like `mp_standard_parallel` and `mp_central_meridian` and x_page and y_page

!the `mp_standard_parallel and `mp_central_meridian` are constants (the middle of your map), the x_page and y_page are variables (they are just differences from the standard_parallel and central_meridian). so somehow you need to add them or sub them to get the correct curvilinear coordinates. i did not have time for further analyse.

let x_st=x_page+
`mp_standard_parallel`
let y_st=y_page+
`mp_central_meridian`

after this you can use:

let your_new=rect_to_curv(your_old,x_st,y_st,1)
save/file="..." your_new

hope this helps a bit,
Peter

[Thread Prev][Thread Next][Index]

Contact Us
Dept of Commerce / NOAA / OAR / PMEL / TMAP

Privacy Policy | Disclaimer | Accessibility Statement