Data processing problem

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
Received on Tue Feb 12 20:52:05 2013