<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=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]-->
</head>
<body lang="HU" link="#0563C1" vlink="#954F72">
<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 = "http://www.fao.org/docrep/X0490E/x0490e07.htm" ;<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 = "http://www.fao.org/docrep/X0490E/x0490e07.htm" ;<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 (http://mpimet.mpg.de/cdi)" ;<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 (http://mpimet.mpg.de/cdo)" ;<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>
</body>
</html>