<html>
<head>
<meta content="text/html; charset=windows-1252"
http-equiv="Content-Type">
</head>
<body bgcolor="#FFFFFF" text="#000000">
Hi Mike,<br>
<br>
From where are you getting these files? The x,y coordinates are in
meters and not degrees. Without specifying the<br>
lat/lon it looks like you are getting the full disk image. The
netcdf files at JMA appear to be in lat/lon<br>
<a class="moz-txt-link-freetext" href="http://www.data.jma.go.jp/mscweb/en/himawari89/space_segment/sample_netcdf.html">http://www.data.jma.go.jp/mscweb/en/himawari89/space_segment/sample_netcdf.html</a>,
but as yet they only<br>
seem to have the Japan-centric images.<br>
<br>
Jim<br>
<br>
<div class="moz-cite-prefix">On 7/23/15 8:00 PM, Mike Manning wrote:<br>
</div>
<blockquote cite="mid:00b401d0c5d6$12747680$375d6380$@bsch.com.au"
type="cite">
<meta http-equiv="Content-Type" content="text/html;
charset=windows-1252">
<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]-->
<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>
<br>
<fieldset class="mimeAttachmentHeader"></fieldset>
<br>
<pre wrap="">_______________________________________________
gradsusr mailing list
<a class="moz-txt-link-abbreviated" href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a>
<a class="moz-txt-link-freetext" href="http://gradsusr.org/mailman/listinfo/gradsusr">http://gradsusr.org/mailman/listinfo/gradsusr</a>
</pre>
</blockquote>
<br>
</body>
</html>