grads指定经纬度读取nc数据进行输出

loving

路过的大神们，给点建议或帮助，小妹在这谢谢了，下面是nc文件结构，直接用sdfopen会报如下错：ga-> sdfopen E:\Wave_SMW_20141113.000000.nc
Scanning self-describing file:  E:\Wave_SMW_20141113.000000.nc
gadsdf: SDF file has no discernable X coordinate.

急需解决 ，请路过的大神们给点指点，我实在脑子转不过来了

nc文件结构：
netcdf E:\Wave_SMW_20141113.000000 {dimensions:
        time = UNLIMITED ; // (73 currently)
        points = 48560 ;
variables:
        int time(time) ;
                time:units = "seconds since 1970-01-01" ;
                time:calendar = "gregorian" ;
                time:standard_name = "time" ;
                time:long_name = "time" ;
        float longitude(points) ;
                longitude:units = "degrees_east" ;
                longitude:long_name = "longitude" ;
                longitude:standard_name = "longitude" ;
        float latitude(points) ;
                latitude:units = "degrees_north" ;
                latitude:long_name = "latitude" ;
                latitude:standard_name = "latitude" ;
        short hs(time, points) ;
                hs:units = "m" ;
                hs:standard_name = "sea_surface_wave_significant_height" ;
                hs:long_name = "hs" ;
                hs:_FillValue = -32768s ;
                hs:scale_factor = 0.0007629627f ;
                hs:add_offset = 25.f ;
        short theta0(time, points) ;
                theta0:units = "degrees" ;
                theta0:standard_name = "sea_surface_wave_to_direction" ;
                theta0:long_name = "theta0" ;
                theta0:_FillValue = -32768s ;
                theta0:scale_factor = 0.005493332f ;
                theta0:add_offset = 180.f ;
        short tm01(time, points) ;
                tm01:units = "s" ;
                tm01:standard_name = "sea_surface_wave_mean_period_from_variance
_spectral_density_first_frequency_moment" ;
                tm01:long_name = "tm01" ;
                tm01:_FillValue = -32768s ;
                tm01:scale_factor = 0.0007629627f ;
                tm01:add_offset = 25.f ;

// global attributes:
                :Conventions = "CF-1.5" ;
                :History = "Created with agioncmd version 1.2" ;
                :Directional convention = "cartesian" ;
                :project = "SMW" ;
                :run = "01" ;