Re: How can I use equivalent latitudes?

Hi Theresa,

Attached is a sample script that will plot the data
using the conventional latitudes and the equivalent latitudes.
The key is the explicitly create the sfXArray and sfYArray
resources and convert the corresponding data array to one-dimensional.

       diri = "./"
       fili = "PV_Eqlat_600K_22503_2403.nc"
       f = addfile(diri+fili, "r")

       x = f->vortpot_0600K(nt,:,:) ; (48,64)
       elat = f->eqlat_0600K(nt,:,:) ; (48,64)
       elon = conform_dims(dimsizes(elat), x&lon, 1)

       res@sfXArray = ndtooned(elon)
       res@sfYArray = ndtooned(elat)
       x1d = ndtooned(x)

       res@gsnLeftString = "Equivalent Latitudes"
       plot = gsn_csm_contour_map_polar(wks,x1d,res)

Cheers
D

On 6/29/10 8:16 AM, Runde, Theresa wrote:
> Hello,
>
> I have a big problem with handling equivalent latitudes and using them
> on a regular grid for plotting. I calculated equivalent latitudes and
> now I have an equivalent latitude array of equlat=equlat(time,lat,lon)
> for each potential temperature level. That means: my new (equivalent)
> latitudes are 3 dimensional.
>
> Now I want to plot (for example) potential vorticity (PV) dependent on
> these new equivalent latitudes (to get a zonally symmetric polar vortex).
> At the moment my PV field is also 3 dimensional (time x lat x lon), but
> regular (latitudes and longitudes are always the same). How can I
> transform my PV field (time x lat(1d) x lon) to a PV field (time x
> equivalent lats(3d) x lon) and plot this? The time dependence is necessary.
>
> Thank you for your help,
> Theresa Runde
> _______________________________________________
> 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