[gradsusr] Ctl to open Satellite netcdf data

[andersen panjaitan]

Hi all,

I want to open an satellite data with format netcdf. This is header of this
file :

    netcdf file:/D:/Kerja/gradz/data/201008060032MTS.Indonesia--L1B.nc {
>
>      dimensions:
>        LATITUDE = 1082;
>        LONGITUDE = 1758;
>        MTSAT_VIS = 65536;
>        MTSAT_IR1 = 65536;
>        MTSAT_IR2 = 65536;
>        MTSAT_IR3 = 65536;
>        MTSAT_IR4 = 65536;
>
>      variables:
>        short MTSAT_VIS_DATA(LATITUDE=1082, LONGITUDE=1758);
>        float MTSAT_VIS_LUT(MTSAT_VIS=65536);
>        short MTSAT_IR1_DATA(LATITUDE=1082, LONGITUDE=1758);
>        float MTSAT_IR1_LUT(MTSAT_IR1=65536);
>        short MTSAT_IR2_DATA(LATITUDE=1082, LONGITUDE=1758);
>        float MTSAT_IR2_LUT(MTSAT_IR2=65536);
>        short MTSAT_IR3_DATA(LATITUDE=1082, LONGITUDE=1758);
>        float MTSAT_IR3_LUT(MTSAT_IR3=65536);
>        short MTSAT_IR4_DATA(LATITUDE=1082, LONGITUDE=1758);
>        float MTSAT_IR4_LUT(MTSAT_IR4=65536);
>
>      :MINLON = 85.0f; // float
>      :MAXLON = 150.0f; // float
>      :MINLAT = -20.0f; // float
>      :MAXLAT = 20.0f; // float
>      :DATA_MINLON = 85.0f; // float
>      :DATA_MAXLON = 150.0f; // float
>      :DATA_MINLAT = -20.0f; // float
>      :DATA_MAXLAT = 20.0f; // float
>      :COVER = 75.0f; // float
>      :RESOLUTION = 0.037f; // float
>      :NLON = 1758; // int
>      :NLAT = 1082; // int
>      :BPP = 2; // int
>      :NCHAN = 5; // int
>      :DATETIME = "201008060032";
>      :PRODUCT = "L1B";
>      :SETUP = "Indonesia-";
>      :SATELLITE = "MTS";
>      :CHANNELS = " MTSAT_VIS,MTSAT_IR1,MTSAT_IR2,MTSAT_IR3,MTSAT_IR4";
>      :SEPARATOR = "-";
>     }
>


I have create a ctl, like this :

    dset ^201008060032MTS.Indonesia--L1B.nc
>     dtype netcdf
>     undef -1e+15
>     title mtsat 201008020732
>     xdef 1758 linear  85.  0.037
>     ydef 1082 linear  -20.  0.037
>     zdef 1 levels 0
>     tdef 1 linear 00z06aug2010 1hr
>     vars 5
>     MTSAT_VIS_DATA=>vis 0 y,x visible
>     MTSAT_IR1_DATA=>ir1 0 y,x ir1
>     MTSAT_IR2_DATA=>ir2 0 y,x ir2
>     MTSAT_IR3_DATA=>ir3 0 y,x ir3
>     MTSAT_IR4_DATA=>ir4 0 y,x ir4
>     endvars
>

but it only can open count data variable. I hope i could open the
temperature variable from convert count data with LUT, but i don't know how
to write it in ctl. Could you help me find some suggestion to solve this??
Thanks

Meteorological Student

Andersen
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20100923/7f78d0fe/attachment-0003.html