<html>
<head>
<meta http-equiv="Content-Type" content="text/html;
charset=iso-8859-2">
</head>
<body text="#000000" bgcolor="#FFFFFF">
Hi Beata,<br>
<br>
I would expect that the variable for the rotated pole information
has to look like<br>
<br>
int rotated_pole ;<br>
rotated_pole:grid_mapping_name =
"rotated_latitude_longitude" ;<br>
rotated_pole:grid_north_pole_latitude = 39.25. ;<br>
rotated_pole:grid_north_pole_longitude = -162. ;<br>
<br>
-Karin<br>
<br>
<div class="moz-cite-prefix">Am 31.08.17 um 13:24 schrieb Beáta
Szabó-Takács:<br>
</div>
<blockquote type="cite"
cite="mid:VI1PR08MB26858F465815D125BE3E7F3EF09D0@VI1PR08MB2685.eurprd08.prod.outlook.com">
<meta http-equiv="Content-Type" content="text/html;
charset=iso-8859-2">
<meta name="Generator" content="Microsoft Word 15 (filtered
medium)">
<style><!--
/* Font Definitions */
@font-face
{font-family:"Cambria Math";
panose-1:2 4 5 3 5 4 6 3 2 4;}
@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:#0563C1;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-priority:99;
color:#954F72;
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:70.85pt 70.85pt 70.85pt 70.85pt;}
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">Dear NCL Users,<o:p></o:p></p>
<p class="MsoNormal">I have also a problem with netCDF writing.
I calculated the Hargreaves reference evapotranspiration in
each year from 1981 to 2000 and wrote out the results into a
netCDF file. I attached my script. The created netCDF file is
correct according to ncdump -h:<o:p></o:p></p>
<p class="MsoNormal">ncdump -h evtH_31.nc<o:p></o:p></p>
<p class="MsoNormal">netcdf evtH_31 {<o:p></o:p></p>
<p class="MsoNormal">dimensions:<o:p></o:p></p>
<p class="MsoNormal"> time = 364 ;<o:p></o:p></p>
<p class="MsoNormal"> rlat = 412 ;<o:p></o:p></p>
<p class="MsoNormal"> rlon = 424 ;<o:p></o:p></p>
<p class="MsoNormal"> vertices = 4 ;<o:p></o:p></p>
<p class="MsoNormal"> bnds = 2 ;<o:p></o:p></p>
<p class="MsoNormal">variables:<o:p></o:p></p>
<p class="MsoNormal"> float evtH(time, rlat, rlon) ;<o:p></o:p></p>
<p class="MsoNormal"> evtH:_FillValue = 1.e+20f ;<o:p></o:p></p>
<p class="MsoNormal"> evtH:long_name = "reference
evapotranspiration: Hargreaves " ;<o:p></o:p></p>
<p class="MsoNormal"> evtH:units = "mm/day" ;<o:p></o:p></p>
<p class="MsoNormal"> evtH:url =
<a class="moz-txt-link-rfc2396E" href="http://www.fao.org/docrep/X0490E/x0490e07.htm">"http://www.fao.org/docrep/X0490E/x0490e07.htm"</a> ;<o:p></o:p></p>
<p class="MsoNormal"> evtH:info = "FAO 56; EQN
52; refevt_hargreaves_fao56" ;<o:p></o:p></p>
<p class="MsoNormal"> float lat(rlat, rlon) ;<o:p></o:p></p>
<p class="MsoNormal"> lat:standard_name =
"latitude" ;<o:p></o:p></p>
<p class="MsoNormal"> lat:long_name = "latitude"
;<o:p></o:p></p>
<p class="MsoNormal"> lat:units = "degrees_north"
;<o:p></o:p></p>
<p class="MsoNormal"> lat:_CoordinateAxisType =
"Lat" ;<o:p></o:p></p>
<p class="MsoNormal"> lat:bounds = "lat_bnds" ;<o:p></o:p></p>
<p class="MsoNormal"> double rlat(rlat) ;<o:p></o:p></p>
<p class="MsoNormal"> rlat:standard_name =
"grid_latitude" ;<o:p></o:p></p>
<p class="MsoNormal"> rlat:long_name = "latitude
in rotated pole grid" ;<o:p></o:p></p>
<p class="MsoNormal"> rlat:units = "degrees" ;<o:p></o:p></p>
<p class="MsoNormal"> rlat:axis = "Y" ;<o:p></o:p></p>
<p class="MsoNormal"> double rlon(rlon) ;<o:p></o:p></p>
<p class="MsoNormal"> rlon:standard_name =
"grid_longitude" ;<o:p></o:p></p>
<p class="MsoNormal"> rlon:long_name = "longitude
in rotated pole grid" ;<o:p></o:p></p>
<p class="MsoNormal"> rlon:units = "degrees" ;<o:p></o:p></p>
<p class="MsoNormal"> rlon:axis = "X" ;<o:p></o:p></p>
<p class="MsoNormal"> float lat_bnds(rlat, rlon,
vertices) ;<o:p></o:p></p>
<p class="MsoNormal"> float lon(rlat, rlon) ;<o:p></o:p></p>
<p class="MsoNormal"> lon:standard_name =
"longitude" ;<o:p></o:p></p>
<p class="MsoNormal"> lon:long_name = "longitude"
;<o:p></o:p></p>
<p class="MsoNormal"> lon:units = "degrees_east"
;<o:p></o:p></p>
<p class="MsoNormal"> lon:_CoordinateAxisType =
"Lon" ;<o:p></o:p></p>
<p class="MsoNormal"> lon:bounds = "lon_bnds" ;<o:p></o:p></p>
<p class="MsoNormal"> float lon_bnds(rlat, rlon,
vertices) ;<o:p></o:p></p>
<p class="MsoNormal"> int rotated_latitude_longitude ;<o:p></o:p></p>
<p class="MsoNormal">
rotated_latitude_longitude:grid_mapping_name =
"rotated_latitude_longitude" ;<o:p></o:p></p>
<p class="MsoNormal">
rotated_latitude_longitude:grid_north_pole_latitude = 39.25 ;<o:p></o:p></p>
<p class="MsoNormal">
rotated_latitude_longitude:grid_north_pole_longitude = -162. ;<o:p></o:p></p>
<p class="MsoNormal">
rotated_latitude_longitude:north_pole_grid_longitude = 0. ;<o:p></o:p></p>
<p class="MsoNormal"> double time(time) ;<o:p></o:p></p>
<p class="MsoNormal"> time:standard_name = "time"
;<o:p></o:p></p>
<p class="MsoNormal"> time:long_name = "time" ;<o:p></o:p></p>
<p class="MsoNormal"> time:bounds = "time_bnds" ;<o:p></o:p></p>
<p class="MsoNormal"> time:units = "days since
1949-12-1 00:00:00" ;<o:p></o:p></p>
<p class="MsoNormal"> time:calendar =
"proleptic_gregorian" ;<o:p></o:p></p>
<p class="MsoNormal"> time:axis = "T" ;<o:p></o:p></p>
<p class="MsoNormal"> double time_bnds(time, bnds) ;<o:p></o:p></p>
<p class="MsoNormal">}<o:p></o:p></p>
<p class="MsoNormal">But when I use this netCDF file in Panoply
I get the following error:
<o:p></o:p></p>
<p class="MsoNormal">There was an error preparing the data:
Could not find pole attributes.
<o:p></o:p></p>
<p class="MsoNormal">Moreover, when I merged the created evtH.nc
files in CDO with mergetime operator and calculated the yearly
sum the longitude and latitude dimeansiones disappeared.
<o:p></o:p></p>
<p class="MsoNormal">ncdump -h PET_ys.nc<o:p></o:p></p>
<p class="MsoNormal">netcdf PET_ys {<o:p></o:p></p>
<p class="MsoNormal">dimensions:<o:p></o:p></p>
<p class="MsoNormal"> rlon = 424 ;<o:p></o:p></p>
<p class="MsoNormal"> rlat = 412 ;<o:p></o:p></p>
<p class="MsoNormal"> time = UNLIMITED ; // (20
currently)<o:p></o:p></p>
<p class="MsoNormal"> bnds = 2 ;<o:p></o:p></p>
<p class="MsoNormal">variables:<o:p></o:p></p>
<p class="MsoNormal"> double rlon(rlon) ;<o:p></o:p></p>
<p class="MsoNormal"> rlon:standard_name =
"longitude" ;<o:p></o:p></p>
<p class="MsoNormal"> rlon:long_name =
"longitude" ;<o:p></o:p></p>
<p class="MsoNormal"> rlon:units = "degrees_east"
;<o:p></o:p></p>
<p class="MsoNormal"> rlon:axis = "X" ;<o:p></o:p></p>
<p class="MsoNormal"> double rlat(rlat) ;<o:p></o:p></p>
<p class="MsoNormal"> rlat:standard_name =
"latitude" ;<o:p></o:p></p>
<p class="MsoNormal"> rlat:long_name = "latitude"
;<o:p></o:p></p>
<p class="MsoNormal"> rlat:units =
"degrees_north" ;<o:p></o:p></p>
<p class="MsoNormal"> rlat:axis = "Y" ;<o:p></o:p></p>
<p class="MsoNormal"> double time(time) ;<o:p></o:p></p>
<p class="MsoNormal"> time:standard_name = "time"
;<o:p></o:p></p>
<p class="MsoNormal"> time:long_name = "time" ;<o:p></o:p></p>
<p class="MsoNormal"> time:bounds = "time_bnds" ;<o:p></o:p></p>
<p class="MsoNormal"> time:units = "days since
1949-12-1 00:00:00" ;<o:p></o:p></p>
<p class="MsoNormal"> time:calendar =
"proleptic_gregorian" ;<o:p></o:p></p>
<p class="MsoNormal"> time:axis = "T" ;<o:p></o:p></p>
<p class="MsoNormal"> double time_bnds(time, bnds) ;<o:p></o:p></p>
<p class="MsoNormal"> float evtH(time, rlat, rlon) ;<o:p></o:p></p>
<p class="MsoNormal"> evtH:long_name = "reference
evapotranspiration: Hargreaves " ;<o:p></o:p></p>
<p class="MsoNormal"> evtH:units = "mm/day" ;<o:p></o:p></p>
<p class="MsoNormal"> evtH:_FillValue = 1.e+20f ;<o:p></o:p></p>
<p class="MsoNormal"> evtH:missing_value =
1.e+20f ;<o:p></o:p></p>
<p class="MsoNormal"> evtH:url =
<a class="moz-txt-link-rfc2396E" href="http://www.fao.org/docrep/X0490E/x0490e07.htm">"http://www.fao.org/docrep/X0490E/x0490e07.htm"</a> ;<o:p></o:p></p>
<p class="MsoNormal"> evtH:info = "FAO 56; EQN
52; refevt_hargreaves_fao56" ;<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">// global attributes:<o:p></o:p></p>
<p class="MsoNormal"> :CDI = "Climate Data
Interface version 1.8.2 (<a class="moz-txt-link-freetext" href="http://mpimet.mpg.de/cdi">http://mpimet.mpg.de/cdi</a>)" ;<o:p></o:p></p>
<p class="MsoNormal"> :Conventions = "CF-1.6" ;<o:p></o:p></p>
<p class="MsoNormal"> :history = "Wed Aug 30
08:00:03 2017: cdo yearsum PET_80-00.nc PET_ys.nc\n",<o:p></o:p></p>
<p class="MsoNormal"> "Wed Aug 30
07:56:15 2017: cdo mergetime evtH_11.nc evtH_12.nc evtH_13.nc
evtH_14.nc evtH_15.nc evtH_16.nc evtH_17.nc evtH_18.nc
evtH_19.nc evtH_20.nc evtH_21.nc evtH_22.nc evtH_23.nc
evtH_24.nc evtH_25.nc evtH_26.nc evtH_27.nc evtH_28.nc
evtH_29.nc evtH_30.nc PET_80-00.nc" ;<o:p></o:p></p>
<p class="MsoNormal"> :frequency = "year" ;<o:p></o:p></p>
<p class="MsoNormal"> :CDO = "Climate Data
Operators version 1.8.2 (<a class="moz-txt-link-freetext" href="http://mpimet.mpg.de/cdo">http://mpimet.mpg.de/cdo</a>)" ;<o:p></o:p></p>
<p class="MsoNormal">}<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">I do not understand why. <o:p></o:p></p>
<p class="MsoNormal">Can someone write me what I did wrong?<o:p></o:p></p>
<p class="MsoNormal">Thank you for your help in advance!<o:p></o:p></p>
<p class="MsoNormal">Best regards,<o:p></o:p></p>
<p class="MsoNormal">Beáta<o:p></o:p></p>
</div>
<br>
<fieldset class="mimeAttachmentHeader"></fieldset>
<br>
<pre wrap="">_______________________________________________
ncl-talk mailing list
<a class="moz-txt-link-abbreviated" href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a>
List instructions, subscriber options, unsubscribe:
<a class="moz-txt-link-freetext" href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a>
</pre>
</blockquote>
<br>
<pre class="moz-signature" cols="72">--
Dipl. Geophys. Karin Meier-Fleischer
Visualization, NCL, CDO
Application Support
Deutsches Klimarechenzentrum GmbH (DKRZ)
Bundesstrasse 45a - D20146 Hamburg - Germany
Phone: +49 (0)40 460094 126
Fax: +49 (0)40 460094 270
E-Mail: <a class="moz-txt-link-abbreviated" href="mailto:meier-fleischer@dkrz.de">meier-fleischer@dkrz.de</a>
URL: <a class="moz-txt-link-abbreviated" href="http://www.dkrz.de">www.dkrz.de</a>
Geschäftsführer: Prof. Dr. Thomas Ludwig
Sitz der Gesellschaft: Hamburg
Amtsgericht Hamburg HRB 39784
</pre>
</body>
</html>