Re: Lambert coordinates

From: Dennis Shea <shea_at_nyahnyahspammersnyahnyah>
Date: Wed Apr 24 2013 - 06:56:38 MDT

As I noted before, please send WRF questions to wrfhelp@ucar.edu

===

Please look at your WRF created netCDF file(s). The variables
XLAT and XLONG contain the grid point locations of
whatever map projection is being used.

[snip]
         float XLAT(Time, south_north, west_east) ;
                 XLAT:FieldType = 104 ;
                 XLAT:MemoryOrder = "XY " ;
                 XLAT:description = "LATITUDE, SOUTH IS NEGATIVE" ;
                 XLAT:units = "degree" ;
                 XLAT:stagger = "" ;
         float XLONG(Time, south_north, west_east) ;
                 XLONG:FieldType = 104 ;
                 XLONG:MemoryOrder = "XY " ;
                 XLONG:description = "LONGITUDE, WEST IS NEGATIVE" ;
                 XLONG:units = "degree" ;
                 XLONG:stagger = "" ;
[snip]

Generally, the coordinates do not change with time so

    f = addfile("wrfout...nc , "r"
    xlat = f->XLAT(0,:,:)
    xlon = f->XLONG(0,:,:)

====
Did you look at any of the examples? If not, please look and read
the documentation.

http://www.ncl.ucar.edu/Applications/wrf.shtml

http://www.mmm.ucar.edu/wrf/OnLineTutorial/Graphics/NCL/
http://www.mmm.ucar.edu/wrf/OnLineTutorial/Graphics/NCL/NCL_examples.htm

On 4/24/13 3:27 AM, Virginie Hergault wrote:
> Thank you for your answer but I do not know how to get my data in the
> lambert coordinates.
>
> In namelist.wps, I put:
>
> map_proj = 'lambert',
> ref_lat = 42.5,
> ref_lon = 140.0,
> truelat1 = 25.0,
> truelat2 = 60.0,
> stand_lon = 140.0
>
> I just want to get the lambert coordinates for all the points of my
> grid. Is it possible easily with a function .ncl?
>
> Thanks again!
> virginie
>
> 2013/4/24, yagnesh@lowtem.hokudai.ac.jp <yagnesh@lowtem.hokudai.ac.jp>:
>>
>> Hello Virginie,
>>
>> On Apr 24 2013, Virginie Hergault <virginie.hergault@gmail.com> wrote:
>>
>>> Do you know which outputs correspond to Lambert coordinates?
>>
>> You can use something like this.
>>
>> ╭─────
>> │;; lifted straight from the WRFUserARW.ncl
>> │map_projection = nc_file@MAP_PROJ
>> │if( map_projection .eq. 1) then ; Lambert Conformal mapping
>> │ blah blah
>> │end if
>> ╰─────
>>
>> and also I would dare to say that, its good idea to glance through
>> "WRFUserARW.ncl" once at-least.
>>
>> Cheers.,
>> --
>> ఎందరో మహానుభావులు అందరికి వందనములు.
>> YYR
>>
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Wed Apr 24 06:56:47 2013

This archive was generated by hypermail 2.1.8 : Fri Apr 26 2013 - 17:10:25 MDT