>I would like to use the NCL-equivalent of grads "set gxout grid."
>For an XY plot (month versus year), I would like grid boxes to
>be plotted with the numerical value appearing at the center of
>each box. How do I do this?
>
>Example:
>
>y ------------
>e | 2 | 4 |
>a ------------
>r | 3 | 5 |
> ------------
> month
>_______________________________________________
So .... Assuming your grid is defined by lat[*], lon[*]
Don't know how Grads numbers grid boxes
A crude way which assumes lat [S->N]
I am assiming you have drawn a map [plot], possibly
with gridlines ..
tres = True ; text mods desired
tres_at_txFontHeightF = 0.015 ; make smaller than defailt
;tres_at_txJust = CenterCenter ; this is the default
nBox = 0
do nl=0,nlat-2
LAT = (lat(nl)+lat(nl+1))*0.5
do ml=0,mlon-2
nBox = nBox+1
LON = (lon(ml)+lon(ml+1))*0.5
gsn_text(wks,plot,sprinti("%0.3i", i) ,LON,LAT,tres)
end do
end do
there are fancier/faster ways of doing the above but
this will do the job.
--- If your grid is cyclic, you may want to create nLAT = nlat+1 mLON = mlon+1 latitude = new (nLAT, typeof(lat) ) longitude= mLON latitude(0:nlat-1) = (/ lat /) longitude(0:mlon-1)= (/ lon /) latitude(nlat) = lat(nlat-1) + dlat longitude(mlon)= lon(mlon-1) + dlon and use these in the above loop with max values nLAT-1 mLON-1 good luck D _______________________________________________ ncl-talk mailing list ncl-talk_at_ucar.edu http://mailman.ucar.edu/mailman/listinfo/ncl-talkReceived on Thu Jun 08 2006 - 09:16:53 MDT
This archive was generated by hypermail 2.2.0 : Thu Jun 08 2006 - 13:13:37 MDT