[ncl-talk] plotting problem

[Alan Brammer]

NCL will plot data on a map correctly regardless of the direction of the coordinates.  For some routines though ascending latitude coordinates are needed. 

The plot you attached looks correctly positioned to me. Why do you think it’s wrong?




> On 23 Oct 2015, at 08:14, mireiyue <mireiyue at gmail.com> wrote:
> 
> Dear NCL users
> 
> I am trying to plot JRA mass flux data. I can plot but it looks like upside down. So I checked latitude and it was 
> 
> Variable: lat
> Type: float
> Total Size: 580 bytes
>            145 values
> Number of Dimensions: 1
> Dimensions and sizes:	[g0_lat_2 | 145]
> Coordinates:
>            g0_lat_2: [90..-90]
> Number Of Attributes: 9
>  long_name :	latitude
>  GridType :	Cylindrical Equidistant Projection Grid
>  units :	degrees_north
>  Dj :	1.25
>  Di :	1.25
>  Lo2 :	-1.25
>  La2 :	-90
>  Lo1 :	 0
>  La1 :	90
> 
> Therefore, I used "::-1” to flip latitude. Then it shows 
> 
> Coordinates:
>            time: [0..83]
>            g0_lat_2: [-90..90]
>            g0_lon_3: [ 0..358.75]
> however, the figure is still same. how can I fix this problem?
> 
> Thank you for any helps.
> Sun-
> 
>  f = addfile("/Users/spark/vertical/00_DATA/20_JRA-55/jra/jramf20072013.nc/","r")
>   cmf = f->mf
>   lat = f->g0_lat_2
>   cmf at _FillValue = getFillValue(cmf)
>   cmf!0 = "time"
>   time = ispan(0,83,1)
>   cmf&time = time
> printVarSummary(cmf)
> 
>   mc = dim_avg_n_Wrap(cmf,1)
> printVarSummary(mc)
> ;print(lat)
> 
>   mfl= mc(:,::-1,:)
> printVarSummary(mfl)
> print(mfl&g0_lat_2(:))
> return
> 
> 
> 
> <jra_20072013.eps>_______________________________________________
> ncl-talk mailing list
> ncl-talk at ucar.edu
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk