Re: rotated grid

From: Dennis Shea <shea_at_nyahnyahspammersnyahnyah>
Date: Tue, 12 Feb 2008 15:13:08 -0700

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