Re: Data processing problem

From: Dennis Shea <shea_at_nyahnyahspammersnyahnyah>
Date: Tue Feb 12 2013 - 21:03:33 MST

The data on the *file* may be 1676 x 4320

BUT you sub-sampled the data (decimated)

t = short2flt(in->med_c0(::5,::10))

decimates the lat dimension to 1/5 the original size (1676/5)
and
decimates the lon dimension to 1/10 the original size (4320/10)
-----

If you want to retain the original dimensions

t = short2flt(in->med_c0) ; read all the data

Good Luck

On 2/12/13 8:51 PM, Jiang, Lifen wrote:
> Hello,
> I have a data set with two dimensions: lat=1676 and lon=4320; Variable type is short; Missing data are -32768 but many values of -9999 and -9998 should also be treated as missing values.
> I want to divide all values by 10 but the returned result has two dimensions: nc10=336 and nc11=432 instead of lat=1676 and lon=4320. And I canít plot it at all.
> I greatly appreciate if anyone can help check whatís wrong with my scripts as bellow.
> Regards,
> Lifen
> **************************************************************************
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
>
> in = addfile("med_c0.nc", "r")
>
>
> t = short2flt(in->med_c0(::5,::10)) ; convert to float
> t@_FillValue = -32768
> t@missing_value = t@_FillValue
> t = where(t .lt. -9997, t@_FillValue, t)
>
> y = t/10
> printVarSummary(y)
>
>
> system("/bin/rm -f med_c0_kgm-2.nc")
> fout = addfile("med_c0_kgm-2.nc", "c") ; new netCDF file
> fout@title = "kg/m2"
> fout->y = y
>
>
>
>
>
>
> _______________________________________________
> 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 Tue Feb 12 21:03:39 2013

This archive was generated by hypermail 2.1.8 : Wed Feb 13 2013 - 09:25:58 MST