Wolfgang Langhans wrote:
> Hi,
> I want to plot temperature fields from the Remo model, which has a rotated grid.
> When I use the rotaded lat/lon (which is both a 1-dimensional array) I can
> easily plot the output, but of course the overlayed map is not identical to the
> grid location. Thus I tried to plot with the real lat/lon arrays (which are
> 2-dimensional). Here I have the problem that the map is not placed correctly.
> Of course, I could plot it succesfully with the map-projections. But here the
> output shows areas with missing values, as the viewing area is not identical
> with the grid.
> Thanks for your help!!!
> Wolfgang
>
> P.S.:Here is my script:
>
>
[SNIP the script you sent]
====================================
The script you sent had little to do with the file.
It looks like you extracted some other script and
changed it.
[0] Always look at your data before proceding.
ncl_filedump -v time,rlat,rlon e001500mean1961-2000_167.nc
[1] NCL uses a netCDF variable model. When it imports data
it reads the array and attributes and, if applicable, netCDF
coordinate arrays.
see: http://www.cgd.ucar.edu/~shea/netCDF-varimodel.png
[2] By netCDF definition, rlat(rlat) and rlon(rlon) are coordinate
variables. Hence, NCL assigns these to the variable.
[3] The units associated with rlat and rlon are not standard.
"degrees_east" and "degrees_north" are commonly used.
[4] If you look at the file, there is only one level and one time step.
Your code has : T_surface = f->var167(100,0,:,:)
Where does the 100 come from? It should be:
T_surface = f->var167(0,0,:,:) ; => T_surface(rlat,rlon)
or
T_surface = f->var167 ;=>
T_surface(time,level,rlat,rlon)
[5] Unless you know what you are doing do not use:
resm_at_tfDoNDCOverlay = True
===============================================================
Try the attached script.
If new to NCL please download and read:
http://www.ncl.ucar.edu/Document/Manuals/
The mini-Language and Graphics manuals.
Good luck
_______________________________________________
ncl-talk mailing list
ncl-talk_at_ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Tue Feb 12 2008 - 15:13:08 MST
This archive was generated by hypermail 2.2.0 : Fri Feb 15 2008 - 17:17:57 MST