[ncl-talk] reodering data

[Agnes Lim Huei Ni]

Hi Alan and Dennis,

This works. Thanks
Agnes

On 3/18/15 10:44 AM, Dennis Shea wrote:
> variable: nr2pgbq_T
> Type: float
> Total Size: 107976960 bytes
>              26994240 values
> Number of Dimensions: 3
> Dimensions and sizes:    [lv_ISBL3 | 26] x [lat | 721] x [lon | 1440]
> Coordinates:
>              lat: [-90..90]
>              lon: [-179.75..180]   <=====
> +++++
> Variable: anal_pgbq_T
> Type: float
> Total Size: 107976960 bytes
>              26994240 values
> Number of Dimensions: 3
> Dimensions and sizes:    [lv_ISBL3 | 26] x [g0_lat_0 | 721] x [g0_lon_1
> | 1440]
> Coordinates:
>              lv_ISBL3: [10..1000]
>              g0_lat_0: [-90..90]
>              g0_lon_1: [ 0..359.75] <=========
> +++++
> Assuming the levels are the same ... (no level coord variable for 
> anal_pgbq_T)
>
> anal_pgbq_T = lonFlip(anal_pgbq_T)
> printVarSummary(anal_pgbq_T)    ; look at lon coord variable
>
> diff = nr2pgbq_T - anal_pgbq_T
> diff at long_name = "Difference T"
> dif at units = anal_pgbq_T at units
> copy_VarCoords(anal_pgbq_T, diff)
> printVarSummary(diff)
>
>
>
>
> On Wed, Mar 18, 2015 at 9:29 AM, Agnes Lim Huei Ni <alim at ssec.wisc.edu 
> <mailto:alim at ssec.wisc.edu>> wrote:
>
>     Hi,
>
>
>     I have two arrays which I would like to take the difference. But looks
>     like i need to do some reordering before i can subtract them.
>
>     Variable: nr2pgbq_T
>     Type: float
>     Total Size: 107976960 bytes
>                  26994240 values
>     Number of Dimensions: 3
>     Dimensions and sizes:    [lv_ISBL3 | 26] x [lat | 721] x [lon | 1440]
>     Coordinates:
>                  lat: [-90..90]
>                  lon: [-179.75..180]
>     Number Of Attributes: 1
>        _FillValue :    1e+15
>
>     Variable: anal_pgbq_T
>     Type: float
>     Total Size: 107976960 bytes
>                  26994240 values
>     Number of Dimensions: 3
>     Dimensions and sizes:    [lv_ISBL3 | 26] x [g0_lat_0 | 721] x
>     [g0_lon_1
>     | 1440]
>     Coordinates:
>                  lv_ISBL3: [10..1000]
>                  g0_lat_0: [-90..90]
>                  g0_lon_1: [ 0..359.75]
>     Number Of Attributes: 12
>        center :    US National Weather Service - NCEP (WMC)
>        long_name :    Temperature
>        units :    K
>        _FillValue :    1e+20
>        level_indicator :    100
>        gds_grid_type :    0
>        parameter_table_version :    2
>        parameter_number :    11
>        model :    Spectral Statistical Interpolation (SSI) analysis from
>     "Final" run.
>        forecast_time :    0
>        forecast_time_units :    hours
>        initial_time :    04/12/2006 (18:00)
>
>     The following is a code segment that I have used.
>
>       anal_pgbq_T=anal_pgbq->TMP_GDS0_ISBL_10
>              lat=guess_pgbq->g0_lat_0
>              lon=guess_pgbq->g0_lon_1
>
>              nr2pgbq_T=nr2pgbq->regrid_int_nr_T
>              anal_pgbq_T=anal_pgbq_T(:,::-1,::-1)
>              anal_pgbq_T!1="lat"
>              anal_pgbq_T!2="lon"
>
>              dim=dimsizes(lon)
>              do i=1, dim(0)-1
>                 if(lon(i) .gt. 180.0) then
>                   lon(i)=lon(i)-360
>                 end if
>              end do
>              lon_1=new((/dim(0)/), "float")
>              lon_1(0:718)=lon(721:1439)
>              lon_1(719:1439)=lon(0:720)
>              delete(dim)
>
>              anal_pgbq_T&lon=lon_1
>
>              printVarSummary(nr2pgbq_T)
>              printVarSummary(anal_pgbq_T)
>
>     I am seeing regions of much larger difference than I expected, does
>     checking if I reorder my data incorrectly before subtracting.
>
>     Thanks
>     Agnes
>
>
>     _______________________________________________
>     ncl-talk mailing list
>     List instructions, subscriber options, unsubscribe:
>     http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>

-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20150318/f97b669f/attachment.html