Due to the current gap in continued funding from the U.S. National Science Foundation (NSF), the NSF Unidata Program Center has temporarily paused most operations. See NSF Unidata Pause in Most Operations for details.
Due to the current gap in continued funding from the U.S. National Science Foundation (NSF), the NSF Unidata Program Center has temporarily paused most operations. See NSF Unidata Pause in Most Operations for details.
I'm trying to use ncwa to get an area average with latitude weighting (weights
are cos(lat)), but I keep getting infinite value as the result. Here is the
ncwa command I used:
ncwa -h -O -w weights -a lat,lon in2.nc out.nc
And here is my in2.nc dump (sorry for the long text):
netcdf in2 {
dimensions:
lat = 11 ;
lon = 11 ;
variables:
double weights(lat) ;
weights:grads_dim = "y" ;
weights:grads_mapping = "linear" ;
weights:grads_size = "180" ;
weights:long_name = "latitude" ;
weights:maximum = 90.f ;
weights:minimum = -90.f ;
weights:resolution = 1.f ;
weights:units = "degrees_north" ;
float lon(lon) ;
lon:grads_dim = "x" ;
lon:grads_mapping = "linear" ;
lon:grads_size = "360" ;
lon:units = "degrees_east" ;
lon:long_name = "longitude" ;
lon:minimum = -180.f ;
lon:maximum = 180.f ;
lon:resolution = 1.f ;
float lat(lat) ;
lat:grads_dim = "y" ;
lat:grads_mapping = "linear" ;
lat:grads_size = "180" ;
lat:units = "degrees_north" ;
lat:long_name = "latitude" ;
lat:minimum = -90.f ;
lat:maximum = 90.f ;
lat:resolution = 1.f ;
float
h5__HDFEOS_GRIDS_Aerosol_NearUV_Grid_Data_Fields_FinalAerosolAbsOpticalDepth388(lat,
lon) ;
h5__HDFEOS_GRIDS_Aerosol_NearUV_Grid_Data_Fields_FinalAerosolAbsOpticalDepth388:h5__HDFEOS_GRIDS_Aerosol_NearUV_Grid_Data_Fields_FinalAerosolAbsOpticalDepth388._FillValue
= -1.267651e+30f ;
h5__HDFEOS_GRIDS_Aerosol_NearUV_Grid_Data_Fields_FinalAerosolAbsOpticalDepth388:h5__HDFEOS_GRIDS_Aerosol_NearUV_Grid_Data_Fields_FinalAerosolAbsOpticalDepth388.units
= "NoUnits" ;
h5__HDFEOS_GRIDS_Aerosol_NearUV_Grid_Data_Fields_FinalAerosolAbsOpticalDepth388:h5__HDFEOS_GRIDS_Aerosol_NearUV_Grid_Data_Fields_FinalAerosolAbsOpticalDepth388.title
= "Final Aerosol Absorption Optical Depth at 388 nm" ;
h5__HDFEOS_GRIDS_Aerosol_NearUV_Grid_Data_Fields_FinalAerosolAbsOpticalDepth388:h5__HDFEOS_GRIDS_Aerosol_NearUV_Grid_Data_Fields_FinalAerosolAbsOpticalDepth388.UniqueFieldDefinition
= "OMI-Specific" ;
h5__HDFEOS_GRIDS_Aerosol_NearUV_Grid_Data_Fields_FinalAerosolAbsOpticalDepth388:h5__HDFEOS_GRIDS_Aerosol_NearUV_Grid_Data_Fields_FinalAerosolAbsOpticalDepth388.scale_factor
= 1. ;
h5__HDFEOS_GRIDS_Aerosol_NearUV_Grid_Data_Fields_FinalAerosolAbsOpticalDepth388:h5__HDFEOS_GRIDS_Aerosol_NearUV_Grid_Data_Fields_FinalAerosolAbsOpticalDepth388.add_offset
= 0. ;
h5__HDFEOS_GRIDS_Aerosol_NearUV_Grid_Data_Fields_FinalAerosolAbsOpticalDepth388:h5__HDFEOS_GRIDS_Aerosol_NearUV_Grid_Data_Fields_FinalAerosolAbsOpticalDepth388.missing_value
= -1.267651e+30f ;
h5__HDFEOS_GRIDS_Aerosol_NearUV_Grid_Data_Fields_FinalAerosolAbsOpticalDepth388:h5__FillValue
= -1.267651e+30f ;
h5__HDFEOS_GRIDS_Aerosol_NearUV_Grid_Data_Fields_FinalAerosolAbsOpticalDepth388:units
= "NoUnits" ;
h5__HDFEOS_GRIDS_Aerosol_NearUV_Grid_Data_Fields_FinalAerosolAbsOpticalDepth388:title
= "Final Aerosol Absorption Optical Depth at 388 nm" ;
h5__HDFEOS_GRIDS_Aerosol_NearUV_Grid_Data_Fields_FinalAerosolAbsOpticalDepth388:UniqueFieldDefinition
= "OMI-Specific" ;
h5__HDFEOS_GRIDS_Aerosol_NearUV_Grid_Data_Fields_FinalAerosolAbsOpticalDepth388:scale_factor
= 1. ;
h5__HDFEOS_GRIDS_Aerosol_NearUV_Grid_Data_Fields_FinalAerosolAbsOpticalDepth388:add_offset
= 0. ;
h5__HDFEOS_GRIDS_Aerosol_NearUV_Grid_Data_Fields_FinalAerosolAbsOpticalDepth388:missing_value
= -1.267651e+30f ;
h5__HDFEOS_GRIDS_Aerosol_NearUV_Grid_Data_Fields_FinalAerosolAbsOpticalDepth388:fonc_original_name
=
"_HDFEOS_GRIDS_Aerosol_NearUV_Grid_Data_Fields_FinalAerosolAbsOpticalDepth388" ;
// global attributes:
:NC_GLOBAL.title = "NASA HDFEOS5 Grid" ;
:NC_GLOBAL.Conventions = "CF-1.4" ;
:NC_GLOBAL.dataType = "Grid" ;
data:
weights = 0.995396469710366, 0.996917515468375, 0.998134908405148,
0.999048277712875, 0.999657345186828, 0.999961925310096,
0.999961925310096, 0.999657345186828, 0.999048277712875,
0.998134908405148, 0.996917515468375 ;
lon = -4.5, -3.5, -2.5, -1.5, -0.5, 0.5, 1.5, 2.5, 3.5, 4.5, 5.5 ;
lat = -5.5, -4.5, -3.5, -2.5, -1.5, -0.5, 0.5, 1.5, 2.5, 3.5, 4.5 ;
h5__HDFEOS_GRIDS_Aerosol_NearUV_Grid_Data_Fields_FinalAerosolAbsOpticalDepth388
=
-1.267651e+30, -1.267651e+30, -1.267651e+30, -1.267651e+30, -1.267651e+30,
-1.267651e+30, -1.267651e+30, -1.267651e+30, -1.267651e+30,
-1.267651e+30, -1.267651e+30,
-1.267651e+30, -1.267651e+30, -1.267651e+30, -1.267651e+30, -1.267651e+30,
-1.267651e+30, -1.267651e+30, -1.267651e+30, -1.267651e+30,
-1.267651e+30, -1.267651e+30,
-1.267651e+30, 0.0684, -1.267651e+30, -1.267651e+30, 0.1196, 0.0661,
0.0648, 0.072, 0.0896, 0.085, -1.267651e+30,
0.0672, 0.0772, 0.0745, 0.0676, 0.0717, 0.0701, 0.0757, 0.0879, 0.1147,
-1.267651e+30, -1.267651e+30,
0.0905, 0.0833, 0.0734, 0.0749, 0.0857, 0.0838, 0.0874, 0.0716, 0.0848,
0.1245, -1.267651e+30,
0.0771, 0.0859, 0.0727, 0.0672, 0.0836, -1.267651e+30, 0.0696, 0.0708,
0.0885, -1.267651e+30, -1.267651e+30,
0.073, 0.0859, 0.0801, 0.0686, 0.0676, -1.267651e+30, -1.267651e+30,
-1.267651e+30, -1.267651e+30, -1.267651e+30, -1.267651e+30,
0.1016, 0.1106, 0.0852, 0.0678, -1.267651e+30, -1.267651e+30,
-1.267651e+30, -1.267651e+30, 0.0898, -1.267651e+30, -1.267651e+30,
0.0693, 0.0755, 0.0971, 0.0733, 0.0845, -1.267651e+30, 0.146, 0.1097,
0.1307, -1.267651e+30, -1.267651e+30,
0.0817, 0.1014, 0.0964, 0.0961, 0.1, 0.0859, 0.1092, 0.1306, -1.267651e+30,
-1.267651e+30, -1.267651e+30,
0.1381, 0.1455, -1.267651e+30, 0.0711, 0.0705, 0.1011, 0.1402, 0.1481,
-1.267651e+30, -1.267651e+30, 0.1522 ;
}
Can anyone help? Thanks.
Jianfu
netcdfgroup archives: