Re: lat and lon grids

From: Jatin Kala <jatin.kala.jk_at_nyahnyahspammersnyahnyah>
Date: Mon Jan 27 2014 - 14:12:26 MST

Hi Anil,
see:
http://www.ncl.ucar.edu/Document/Language/reorder.shtml

This is not tested:
;since there is only 1 “lev”, you can probably get rid of this dim by:
array_no_lev = array(:,:,0,:) ; assuming array is (lon,lat,lev,time)
array_no_lev!0 = “lon”
array_no_lev!1 = “lat”
array_no_lev!2 = “time”
array_reoder = array_no_lev(time|:, lat|:, lon|:)

Cheers,
Jatin.

On 28 Jan 2014, at 4:29 am, Kumar, Anil (GSFC-617.0)[UNIV OF MARYLAND] <anil.kumar@nasa.gov> wrote:

> NCL folks,
> Please help me on this issue.
> I have single array of each like lat, lon, lev, and time. I want to create xlat_m and xlong_m with (time,lat,lon) as defined in the WRF like this XLAT_M(Time, south_north, west_east)
>
> Data in the file looks like this
>
> dimensions:
> lon = 1160 ;
> lat = 490 ;
> lev = 1 ;
> time = UNLIMITED ; // (46 currently)
> variables:
> double lon(lon) ;
> lon:standard_name = "longitude" ;
> lon:long_name = "longitude" ;
> lon:units = "degrees_east" ;
> lon:axis = "X" ;
> double lat(lat) ;
> lat:standard_name = "latitude" ;
> lat:long_name = "latitude" ;
> lat:units = "degrees_north" ;
> lat:axis = "Y" ;
> double lev(lev) ;
> lev:long_name = "generic" ;
> lev:units = "level" ;
> lev:axis = "Z" ;
> double time(time) ;
> time:standard_name = "time" ;
> time:units = "hours since 2001-01-01 00:00:00" ;
> time:calendar = "standard" ;
>
>
> Thanks,
> Anil
>
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Mon Jan 27 14:12:46 2014

This archive was generated by hypermail 2.1.8 : Fri Feb 07 2014 - 16:39:11 MST