<html>
<head>
<meta http-equiv="Content-Type" content="text/html;
charset=iso-8859-2">
</head>
<body text="#000000" bgcolor="#FFFFFF">
Hi Beata,<br>
<br>
by the way, if your data variables would have the coordinates
attribute like<br>
<br>
<font color="#3333ff"><tt> t:coordinates = "lat lon" ;</tt></font><br>
<br>
which can be set for e.g. variable t with the CDO operator
setattribute<br>
<div class="page" title="Page 55">
<div class="layoutArea">
<div class="column">
<p><font color="#3333ff"><tt><span style="font-size: 10pt;
background-color: rgb(255, 236, 179);"> cdo
setattribute,t@coordinates="lat lon" infile outfile
</span></tt></font></p>
</div>
</div>
</div>
the lat and lon coordinates won't be deleted ('disappear'). ;)<br>
<br>
-Karin<br>
<br>
<br>
<div class="moz-cite-prefix">Am 31.08.17 um 13:38 schrieb Karin
Meier-Fleischer:<br>
</div>
<blockquote type="cite"
cite="mid:e5621e13-9789-8ea3-4d30-29e2c024f394@dkrz.de">
<meta http-equiv="Content-Type" content="text/html;
charset=iso-8859-2">
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"
moz-do-not-send="true">"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"
moz-do-not-send="true">"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" moz-do-not-send="true">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" moz-do-not-send="true">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" moz-do-not-send="true">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" moz-do-not-send="true">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" moz-do-not-send="true">meier-fleischer@dkrz.de</a>
URL: <a class="moz-txt-link-abbreviated" href="http://www.dkrz.de" moz-do-not-send="true">www.dkrz.de</a>
Geschäftsführer: Prof. Dr. Thomas Ludwig
Sitz der Gesellschaft: Hamburg
Amtsgericht Hamburg HRB 39784
</pre>
<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>
</body>
</html>