<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40"><head><META HTTP-EQUIV="Content-Type" CONTENT="text/html; charset=us-ascii"><meta name=Generator content="Microsoft Word 14 (filtered medium)"><style><!--
/* Font Definitions */
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0cm;
margin-bottom:.0001pt;
font-size:11.0pt;
font-family:"Calibri","sans-serif";
mso-fareast-language:EN-US;}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:blue;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-priority:99;
color:purple;
text-decoration:underline;}
span.EmailStyle17
{mso-style-type:personal-compose;
font-family:"Calibri","sans-serif";
color:windowtext;}
.MsoChpDefault
{mso-style-type:export-only;
font-family:"Calibri","sans-serif";
mso-fareast-language:EN-US;}
@page WordSection1
{size:612.0pt 792.0pt;
margin:72.0pt 72.0pt 72.0pt 72.0pt;}
div.WordSection1
{page:WordSection1;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]--></head><body lang=EN-AU link=blue vlink=purple><div class=WordSection1><p class=MsoNormal>Hi all,<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>I’m trying to display a Himawari-8 NetCDF file in Grads – here are the attributes of the nc file:<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>dimensions:<o:p></o:p></p><p class=MsoNormal> time = UNLIMITED ; // (1 currently)<o:p></o:p></p><p class=MsoNormal> y = 22000 ;<o:p></o:p></p><p class=MsoNormal> x = 22000 ;<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>variables:<o:p></o:p></p><p class=MsoNormal> double time(time) ;<o:p></o:p></p><p class=MsoNormal> time:coordinates = "time" ;<o:p></o:p></p><p class=MsoNormal> time:long_name = "time" ;<o:p></o:p></p><p class=MsoNormal> time:standard_name = "time" ;<o:p></o:p></p><p class=MsoNormal> time:units = "seconds since 1970-01-01 00:00:00" ;<o:p></o:p></p><p class=MsoNormal> time:calendar = "standard" ;<o:p></o:p></p><p class=MsoNormal> time:axis = "T" ;<o:p></o:p></p><p class=MsoNormal> time:unlimited = "True" ;<o:p></o:p></p><p class=MsoNormal> double y(y) ;<o:p></o:p></p><p class=MsoNormal> y:coordinates = "y" ;<o:p></o:p></p><p class=MsoNormal> y:long_name = "projection_y_coordinate" ;<o:p></o:p></p><p class=MsoNormal> y:standard_name = "projection_y_coordinate" ;<o:p></o:p></p><p class=MsoNormal> y:units = "m" ;<o:p></o:p></p><p class=MsoNormal> y:axis = "Y" ;<o:p></o:p></p><p class=MsoNormal> double x(x) ;<o:p></o:p></p><p class=MsoNormal> x:coordinates = "x" ;<o:p></o:p></p><p class=MsoNormal> x:long_name = "projection_x_coordinate" ;<o:p></o:p></p><p class=MsoNormal> x:standard_name = "projection_x_coordinate" ;<o:p></o:p></p><p class=MsoNormal> x:units = "m" ;<o:p></o:p></p><p class=MsoNormal> x:axis = "X" ;<o:p></o:p></p><p class=MsoNormal> float geostationary ;<o:p></o:p></p><p class=MsoNormal> geostationary:proj4 = "+proj=geos +lon_0=140.7 +h=35785863 +x_0=0 +y_0=0 +a=6378137 +b=6356752.3 +units=m +no_defs " ;<o:p></o:p></p><p class=MsoNormal> geostationary:proj_name = "GEOS141" ;<o:p></o:p></p><p class=MsoNormal> geostationary:grid_mapping_name = "geostationary" ;<o:p></o:p></p><p class=MsoNormal> geostationary:longitude_of_projection_origin = 140.7 ;<o:p></o:p></p><p class=MsoNormal> geostationary:perspective_point_height = 35785863. ;<o:p></o:p></p><p class=MsoNormal> geostationary:satellite_height = 35785863. ;<o:p></o:p></p><p class=MsoNormal> geostationary:semi_major_axis = 6378137. ;<o:p></o:p></p><p class=MsoNormal> geostationary:semi_minor_axis = 6356752.3 ;<o:p></o:p></p><p class=MsoNormal> geostationary:sweep_angle_axis = "y" ;<o:p></o:p></p><p class=MsoNormal> geostationary:GeoTransform = -5500000., 500., 0., 5500000., 0., -500. ;<o:p></o:p></p><p class=MsoNormal> geostationary:spatial_ref = "PROJCS[\"unnamed\",GEOGCS[\"unnamed ellipse\",DATUM[\"unknown\",SPHEROID[\"unnamed\",6378137,298.2570248822722]],PRIMEM[\"Greenwich\",0],UNIT[\"degree\",0.0174532925199433]],PROJECTION[\"Geostationary_Satellite\"],PARAMETER[\"central_meridian\",140.7],PARAMETER[\"satellite_height\",35785863],PARAMETER[\"false_easting\",0],PARAMETER[\"false_northing\",0]]" ;<o:p></o:p></p><p class=MsoNormal> double scan_line_time(y) ;<o:p></o:p></p><p class=MsoNormal> scan_line_time:units = "seconds since 2015-07-14 03:00:00" ;<o:p></o:p></p><p class=MsoNormal> scan_line_time:long_name = "estimated scan line time" ;<o:p></o:p></p><p class=MsoNormal> scan_line_time:calendar = "standard" ;<o:p></o:p></p><p class=MsoNormal> scan_line_time:coordinates = "y" ;<o:p></o:p></p><p class=MsoNormal> float channel_0003_brf(time, y, x) ;<o:p></o:p></p><p class=MsoNormal> channel_0003_brf:_FillValue = 1.e+20f ;<o:p></o:p></p><p class=MsoNormal> channel_0003_brf:least_significant_digit = 2s ;<o:p></o:p></p><p class=MsoNormal> channel_0003_brf:algorithm_name = "channel_3_brf" ;<o:p></o:p></p><p class=MsoNormal> channel_0003_brf:central_wavelength = 6.3915e-07 ;<o:p></o:p></p><p class=MsoNormal> channel_0003_brf:comment = "bidirectional reflectance factor observed by the Advanced Himawari Imager at 0.64 um (Band 3)" ;<o:p></o:p></p><p class=MsoNormal> channel_0003_brf:long_name = "bidirectional reflectance factor for channel 3 at 0.64 um" ;<o:p></o:p></p><p class=MsoNormal> channel_0003_brf:units = "1" ;<o:p></o:p></p><p class=MsoNormal> channel_0003_brf:valid_max = 1.2 ;<o:p></o:p></p><p class=MsoNormal> channel_0003_brf:valid_min = 0. ;<o:p></o:p></p><p class=MsoNormal> channel_0003_brf:coordinates = "time y x" ;<o:p></o:p></p><p class=MsoNormal> channel_0003_brf:grid_mapping = "geostationary" ;<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>I made a sample script as follows:<o:p></o:p></p><p class=MsoNormal>'sdfopen testfile.nc'<o:p></o:p></p><p class=MsoNormal>'d channel_0003_br'<o:p></o:p></p><p class=MsoNormal>'printim satimage.png'<o:p></o:p></p><p class=MsoNormal>'quit'<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>Ideally I want to be able to set a lat/lon region, but when I specify any lat/lon region I just get a chart with a red output. When I don’t specify any lat/lon output at all, I get a bizzare chart (I’m guessing because it’s using x and y coordinates instead of lat/lon?) Is there any way to correctly get this to display in grads at all?<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>Cheers, Mike<o:p></o:p></p></div></body></html>