I have been trying to plot some data on a latitude longitude slice,
and I want to generate a mirror image of the results in longitude i.e.
so that the results at for example 30 W are plotted at 30 E and so on.
For some reason, whatever combination of subscripts I use, or however
I try to reorder the elements and put them into another variable, the
longitudinal ordering always remains exactly the same, and the plots I
generate always look identical. I think there must be something
obvious I am doing wrong, but I am not sure what it might be.
Does anyone have any ideas how I can reverse the direction of my
results in longitude on my plot ?. I have included below a section of
my script in which I am trying to reorder the results. (I also tried
the "plot" commands without including the -1, but the results were
always exactly the same).
Thanks,
Helen.
kl = 49
ureorder = U(lev | kl, lat | :, lon | :, time | :) ;
(lev,lat,lon,time)
utime = dim_avg_Wrap(ureorder) ; (lev,lat,lon)
utimelon = new((/192,288/),"float")
utimelon2 = new((/192,288/),"float")
do i = 0,287
utimelon(:,i) = utime(:,287-i)
utimelon2(:,i) = utime(:,i)
end do
plot = gsn_csm_contour(wks, utimelon(:,287:0:-1), res ) ;
(lev,lat)
plot = gsn_csm_contour(wks, utimelon(:,0:287:-1), res ) ;
(lev,lat)
plot = gsn_csm_contour(wks, utimelon2(:,287:0:-1), res ) ;
(lev,lat)
plot = gsn_csm_contour(wks, utimelon2(:,0:287:-1), res ) ;
(lev,lat)
plot = gsn_csm_contour(wks, utime(:,0:287:-1), res ) ;
(lev,lat)
plot = gsn_csm_contour(wks, utime(:,287:0:-1), res ) ;
(lev,lat)
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Fri Jan 27 18:13:04 2012
This archive was generated by hypermail 2.1.8 : Thu Feb 02 2012 - 03:10:31 MST