<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=utf-8">
<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:12.0pt;
font-family:"Times New Roman",serif;}
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;}
p.MsoNoSpacing, li.MsoNoSpacing, div.MsoNoSpacing
{mso-style-priority:1;
margin:0cm;
margin-bottom:.0001pt;
font-size:11.0pt;
font-family:"Calibri",sans-serif;
mso-fareast-language:EN-US;}
p.m6018092740948470009msonospacing, li.m6018092740948470009msonospacing, div.m6018092740948470009msonospacing
{mso-style-name:m_6018092740948470009msonospacing;
mso-margin-top-alt:auto;
margin-right:0cm;
mso-margin-bottom-alt:auto;
margin-left:0cm;
font-size:12.0pt;
font-family:"Times New Roman",serif;}
p.m6018092740948470009m4923127676996268827gmail-m7531163896496914688msonospacing, li.m6018092740948470009m4923127676996268827gmail-m7531163896496914688msonospacing, div.m6018092740948470009m4923127676996268827gmail-m7531163896496914688msonospacing
{mso-style-name:m_6018092740948470009m4923127676996268827gmail-m7531163896496914688msonospacing;
mso-margin-top-alt:auto;
margin-right:0cm;
mso-margin-bottom-alt:auto;
margin-left:0cm;
font-size:12.0pt;
font-family:"Times New Roman",serif;}
p.m6018092740948470009m4923127676996268827m7531163896496914688msonospacing, li.m6018092740948470009m4923127676996268827m7531163896496914688msonospacing, div.m6018092740948470009m4923127676996268827m7531163896496914688msonospacing
{mso-style-name:m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing;
mso-margin-top-alt:auto;
margin-right:0cm;
mso-margin-bottom-alt:auto;
margin-left:0cm;
font-size:12.0pt;
font-family:"Times New Roman",serif;}
span.EmailStyle21
{mso-style-type:personal;
font-family:"Calibri",sans-serif;
color:#1F497D;}
span.EmailStyle22
{mso-style-type:personal-reply;
font-family:"Calibri",sans-serif;
color:#1F497D;}
.MsoChpDefault
{mso-style-type:export-only;
font-size:10.0pt;}
@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="blue" vlink="purple">
<div class="WordSection1">
<p class="MsoNormal">Dear Dennis,<span style="font-size:11.0pt"><o:p></o:p></span></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">Sorry for my confusion. You are right that 9.96921e+36 is a float in version 6.x and later. I expected that the table.png will contain same missing value as the script (-999) accordingly version 5.2.x and earlier. Meantime I found a solution
to print -999 by modifying the taylor_metrics_table.ncl with:<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"> if (ismissing(values(icase,iseas,im)) .or. \ ; ???<o:p></o:p></p>
<p class="MsoNormal"> ismissing(values(icase,iseas,im)))then<o:p></o:p></p>
<p class="MsoNormal">; print ("Missing values skipped")<o:p></o:p></p>
<p class="MsoNormal"> values(icase,iseas,im) = -999<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">It produces table where the values can be better to see.<o:p></o:p></p>
<p class="MsoNormal">Have a nice day,<o:p></o:p></p>
<p class="MsoNormal">Beata<span style="font-family:"Calibri",sans-serif"><o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D;mso-fareast-language:EN-US"><o:p> </o:p></span></p>
<div>
<div style="border:none;border-top:solid #E1E1E1 1.0pt;padding:3.0pt 0cm 0cm 0cm">
<p class="MsoNormal"><b><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri",sans-serif">From:</span></b><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri",sans-serif"> ncl-talk-bounces@ucar.edu [mailto:ncl-talk-bounces@ucar.edu]
<b>On Behalf Of </b>Beáta Szabó-Takács<br>
<b>Sent:</b> Wednesday, June 21, 2017 9:06 AM<br>
<b>To:</b> Dennis Shea <shea@ucar.edu><br>
<b>Cc:</b> ncl-talk (ncl-talk@ucar.edu) <ncl-talk@ucar.edu>; Mary Haley <haley@ucar.edu><br>
<b>Subject:</b> Re: [ncl-talk] convert the double missing_value to float one<o:p></o:p></span></p>
</div>
</div>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNoSpacing">Hi Dennis,<o:p></o:p></p>
<p class="MsoNoSpacing"><o:p> </o:p></p>
<p class="MsoNoSpacing">Thank you for your suggestion. I printed from MPI-ESM-LR-CLM4.8.17_dry_power_sm.nc type of pr values as you suggested:<o:p></o:p></p>
<p class="MsoNoSpacing"><o:p> </o:p></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;mso-fareast-language:EN-US">print(typeof(pr))
<o:p></o:p></span></p>
<p class="MsoNoSpacing">print( typeof(pr@_FillValue))<o:p></o:p></p>
<p class="MsoNoSpacing">and I got the following results:<o:p></o:p></p>
<p class="MsoNoSpacing"> <o:p></o:p></p>
<p class="MsoNoSpacing">[szabo.b@fenix data]$ ncl taylor_1.ncl<o:p></o:p></p>
<p class="MsoNoSpacing">Copyright (C) 1995-2017 - All Rights Reserved<o:p></o:p></p>
<p class="MsoNoSpacing">University Corporation for Atmospheric Research<o:p></o:p></p>
<p class="MsoNoSpacing">NCAR Command Language Version 6.4.0<o:p></o:p></p>
<p class="MsoNoSpacing">The use of this software is governed by a License Agreement.<o:p></o:p></p>
<p class="MsoNoSpacing">See <a href="http://www.ncl.ucar.edu/">http://www.ncl.ucar.edu/</a> for more details.<o:p></o:p></p>
<p class="MsoNoSpacing">(0) float<o:p></o:p></p>
<p class="MsoNoSpacing">(0) float<o:p></o:p></p>
<p class="MsoNoSpacing"><o:p> </o:p></p>
<p class="MsoNoSpacing"><o:p> </o:p></p>
<p class="MsoNoSpacing">I also printed CD_biasn values which resulted _FillValues :
<o:p></o:p></p>
<p class="MsoNoSpacing">Variable: CD_biasn<o:p></o:p></p>
<p class="MsoNoSpacing">Type: float<o:p></o:p></p>
<p class="MsoNoSpacing">Total Size: 20 bytes<o:p></o:p></p>
<p class="MsoNoSpacing"> 5 values<o:p></o:p></p>
<p class="MsoNoSpacing">Number of Dimensions: 1<o:p></o:p></p>
<p class="MsoNoSpacing">Dimensions and sizes: [5]<o:p></o:p></p>
<p class="MsoNoSpacing">Coordinates:<o:p></o:p></p>
<p class="MsoNoSpacing">Number Of Attributes: 1<o:p></o:p></p>
<p class="MsoNoSpacing"> _FillValue : -999<o:p></o:p></p>
<p class="MsoNoSpacing">(0) -999<o:p></o:p></p>
<p class="MsoNoSpacing">(1) -999<o:p></o:p></p>
<p class="MsoNoSpacing">(2) -999<o:p></o:p></p>
<p class="MsoNoSpacing">(3) -7.358609<o:p></o:p></p>
<p class="MsoNoSpacing">(4) -3.67009<o:p></o:p></p>
<p class="MsoNoSpacing"><o:p> </o:p></p>
<p class="MsoNoSpacing">In the resulted taylor_bias_wet.png table by taylor_metric_table.ncl these _FillValues are double. I attched the table. Moreover I used taylor_diagram_lh.ncl for taylor diagram which printed the X, Y coordinate values of diagram:
<o:p></o:p></p>
<p class="MsoNoSpacing">(0) X(n,i) = 0.91817 Y(n,i) = 0.338616<o:p></o:p></p>
<p class="MsoNoSpacing">(0) X(n,i) = 0.728601 Y(n,i) = 0.418846<o:p></o:p></p>
<p class="MsoNoSpacing">(0) X(n,i) = 0.618642 Y(n,i) = 0.295433<o:p></o:p></p>
<p class="MsoNoSpacing">(0) X(n,i) = 1.09463 Y(n,i) = 0.428762<o:p></o:p></p>
<p class="MsoNoSpacing">(0) X(n,i) = 0.7988 Y(n,i) = 0.328584<o:p></o:p></p>
<p class="MsoNoSpacing">(0) X(n,i) = 1.27678 Y(n,i) = 0.44035<o:p></o:p></p>
<p class="MsoNoSpacing">(0) X(n,i) = 1.06292 Y(n,i) = 0.606676<o:p></o:p></p>
<p class="MsoNoSpacing">(0) X(n,i) = 0.920748 Y(n,i) = 0.41256<o:p></o:p></p>
<p class="MsoNoSpacing">(0) X(n,i) = 1.39532 Y(n,i) = 0.583074<o:p></o:p></p>
<p class="MsoNoSpacing">(0) X(n,i) = 0.842161 Y(n,i) = 0.348939<o:p></o:p></p>
<p class="MsoNoSpacing">(0) X(n,i) = 0.856069 Y(n,i) = 0.117605<o:p></o:p></p>
<p class="MsoNoSpacing">(0) X(n,i) = 0.85604 Y(n,i) = 0.117578<o:p></o:p></p>
<p class="MsoNoSpacing">(0) X(n,i) = 0.856069 Y(n,i) = 0.117605<o:p></o:p></p>
<p class="MsoNoSpacing">(0) X(n,i) = 0.856069 Y(n,i) = 0.117605<o:p></o:p></p>
<p class="MsoNoSpacing">(0) X(n,i) = 0.856069 Y(n,i) = 0.117605<o:p></o:p></p>
<p class="MsoNoSpacing">(0) X(n,i) = 0.602354 Y(n,i) = 0.316136<o:p></o:p></p>
<p class="MsoNoSpacing">(0) X(n,i) = 1.03864 Y(n,i) = 0.355518<o:p></o:p></p>
<p class="MsoNoSpacing">(0) X(n,i) = 0.902709 Y(n,i) = 0.284563<o:p></o:p></p>
<p class="MsoNoSpacing">(0) X(n,i) = 0.754752 Y(n,i) = 0.122997<o:p></o:p></p>
<p class="MsoNoSpacing">(0) X(n,i) = 0.791836 Y(n,i) = 0.275488<o:p></o:p></p>
<p class="MsoNoSpacing">(0) X(n,i) = 0.561078 Y(n,i) = 0.542241<o:p></o:p></p>
<p class="MsoNoSpacing">(0) X(n,i) = 0.55108 Y(n,i) = 0.558689<o:p></o:p></p>
<p class="MsoNoSpacing">(0) X(n,i) = 0.93554 Y(n,i) = 0.691576<o:p></o:p></p>
<p class="MsoNoSpacing">(0) X(n,i) = 0.857712 Y(n,i) = 0.495699<o:p></o:p></p>
<p class="MsoNoSpacing">(0) X(n,i) = 0.835235 Y(n,i) = 0.531927<o:p></o:p></p>
<p class="MsoNoSpacing">(0) X(n,i) = 0.929597 Y(n,i) = 0.663454<o:p></o:p></p>
<p class="MsoNoSpacing">(0) X(n,i) = 0.864055 Y(n,i) = 0.867048<o:p></o:p></p>
<p class="MsoNoSpacing">(0) X(n,i) = 1.27496 Y(n,i) = 0.861098<o:p></o:p></p>
<p class="MsoNoSpacing">(0) X(n,i) = 1.19209 Y(n,i) = 0.731529<o:p></o:p></p>
<p class="MsoNoSpacing">(0) X(n,i) = 0.903928 Y(n,i) = 0.5948<o:p></o:p></p>
<p class="MsoNoSpacing">(0) X(n,i) = 0.94527 Y(n,i) = 0.370124<o:p></o:p></p>
<p class="MsoNoSpacing">(0) X(n,i) = 0.94527 Y(n,i) = 0.370124<o:p></o:p></p>
<p class="MsoNoSpacing">(0) X(n,i) = 0.945275 Y(n,i) = 0.37013<o:p></o:p></p>
<p class="MsoNoSpacing">(0) X(n,i) = 0.94527 Y(n,i) = 0.370124<o:p></o:p></p>
<p class="MsoNoSpacing">(0) X(n,i) = 0.94527 Y(n,i) = 0.370124<o:p></o:p></p>
<p class="MsoNoSpacing">(0) X(n,i) = 9.96921e+36 Y(n,i) = 9.96921e+36<o:p></o:p></p>
<p class="MsoNoSpacing">(0) X(n,i) = 9.96921e+36 Y(n,i) = 9.96921e+36<o:p></o:p></p>
<p class="MsoNoSpacing">(0) X(n,i) = 9.96921e+36 Y(n,i) = 9.96921e+36<o:p></o:p></p>
<p class="MsoNoSpacing">(0) X(n,i) = 0.878889 Y(n,i) = 0.492841<o:p></o:p></p>
<p class="MsoNoSpacing">(0) X(n,i) = 0.953258 Y(n,i) = 0.563082<o:p></o:p></p>
<p class="MsoNoSpacing"><o:p> </o:p></p>
<p class="MsoNoSpacing">The _FillValues are also double. I attached my script, taylor_diagram_lh.ncl and taylor_metrics_table.ncl. I have sent the MPI-ESM-LR-CLM4.8.17_dry_power_sm.nc by ftp as you suggested.<o:p></o:p></p>
<p class="MsoNoSpacing">Thank you for your help in advance!<o:p></o:p></p>
<p class="MsoNoSpacing">Have a nice day,<o:p></o:p></p>
<p class="MsoNoSpacing">Beata<o:p></o:p></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D;mso-fareast-language:EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D;mso-fareast-language:EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><b><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri",sans-serif">From:</span></b><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri",sans-serif"> Dennis Shea [</span><a href="mailto:shea@ucar.edu"><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri",sans-serif">mailto:shea@ucar.edu</span></a><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri",sans-serif">]
<br>
<b>Sent:</b> Tuesday, June 20, 2017 3:35 PM<br>
<b>To:</b> Beáta Szabó-Takács <</span><a href="mailto:szabo.b@czechglobe.cz"><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri",sans-serif">szabo.b@czechglobe.cz</span></a><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri",sans-serif">><br>
<b>Cc:</b> Mary Haley <</span><a href="mailto:haley@ucar.edu"><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri",sans-serif">haley@ucar.edu</span></a><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri",sans-serif">>; ncl-talk (</span><a href="mailto:ncl-talk@ucar.edu"><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri",sans-serif">ncl-talk@ucar.edu</span></a><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri",sans-serif">)
<</span><a href="mailto:ncl-talk@ucar.edu"><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri",sans-serif">ncl-talk@ucar.edu</span></a><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri",sans-serif">><br>
<b>Subject:</b> Re: [ncl-talk] convert the double missing_value to float one<o:p></o:p></span></p>
<p class="MsoNormal"><o:p> </o:p></p>
<div>
<div>
<div>
<div>
<div>
<div>
<div>
<div>
<div>
<div>
<div>
<div>
<div>
<div>
<div>
<div>
<div>
<p class="MsoNormal" style="margin-bottom:12.0pt">I am Matlab ignorant so I can't answer what Matlab does/does-not do..<br>
<br>
NCL has some 'netCDF rules' built into the language:<br>
<br>
ncl 0> x := fspan(1,10,10) ; float<br>
ncl 1> x@_FillValue = -999d ; double<br>
fatal:Type Mismatch: The type of missing value could not be converted to type of variable (x)<br>
fatal:["Execute.c":8573]:Execute: Error occurred at or near line 1<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">NCL will 'help' in some cases<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><br>
ncl 5> x@_FillValue = -999 ; integer ... here NCL silently 'promotes' the integer to a float<br>
ncl 6> print(x)<br>
ncl 7> print(typeof(x@_FillValue)) ; float<br>
<br>
====<o:p></o:p></p>
</div>
<p class="MsoNormal" style="margin-bottom:12.0pt">ncdump -h <a href="http://MPI-ESM-LR-CLM4.8.17_dry_power.nc" target="_blank">
MPI-ESM-LR-CLM4.8.17_dry_power.nc</a><o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">[snip]<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">float pr(time, latitude, longitude) ; <br>
pr:_FillValue = -999.f ;<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">[snip]<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<p class="MsoNormal" style="margin-bottom:12.0pt">Don't be misled by the 'f' ... extension (-999.f).
<o:p></o:p></p>
</div>
<p class="MsoNormal"> f = addfile("MPI-ESM-LR-CLM4.8.17_dry_power","r")<o:p></o:p></p>
</div>
<p class="MsoNormal"> pr = f->pr<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> print(typeof(pr)) ; float<o:p></o:p></p>
</div>
<p class="MsoNormal"> print( typeof(pr@_FillValue)) ; I speculate this will be type double.<br>
<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> ; Note the := syntax<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> pr@_FillValue := totype( pr@_FillValue, typeof(pr) ) ; MH's suggestion<br>
print( typeof(pr@_FillValue)) <o:p></o:p></p>
</div>
<p class="MsoNormal" style="margin-bottom:12.0pt">[snip]<o:p></o:p></p>
</div>
<p class="MsoNormal"><br>
<br>
====<o:p></o:p></p>
</div>
<p class="MsoNormal" style="margin-bottom:12.0pt">You can send the file to<o:p></o:p></p>
</div>
<p class="MsoNormal">ftp <a href="http://ftp.cgd.ucar.edu" target="_blank">ftp.cgd.ucar.edu</a><o:p></o:p></p>
</div>
<p class="MsoNormal">anonymous<o:p></o:p></p>
</div>
<p class="MsoNormal">your_email<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">cd incoming<o:p></o:p></p>
</div>
<p class="MsoNormal">put .....<o:p></o:p></p>
</div>
<p class="MsoNormal" style="margin-bottom:12.0pt">quit<o:p></o:p></p>
</div>
<p class="MsoNormal" style="margin-bottom:12.0pt">Then notify ncl-talk after a successful file transfer.<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">You can attach a clean version of your script to the email.<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<p class="MsoNormal">Cheers<o:p></o:p></p>
</div>
<p class="MsoNormal">D<o:p></o:p></p>
<div>
<div>
<div>
<div>
<div>
<div>
<div>
<div>
<div>
<p class="MsoNormal" style="margin-bottom:12.0pt"><o:p> </o:p></p>
<div>
<div>
<p class="MsoNormal" style="margin-bottom:12.0pt"><o:p> </o:p></p>
<div>
<div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
<div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
<div>
<p class="MsoNormal">On Tue, Jun 20, 2017 at 1:32 AM, Beáta Szabó-Takács <<a href="mailto:szabo.b@czechglobe.cz" target="_blank">szabo.b@czechglobe.cz</a>> wrote:<o:p></o:p></p>
<blockquote style="border:none;border-left:solid #CCCCCC 1.0pt;padding:0cm 0cm 0cm 6.0pt;margin-left:4.8pt;margin-top:5.0pt;margin-right:0cm;margin-bottom:5.0pt">
<div>
<div>
<p class="m6018092740948470009msonospacing">Dear Dennis,<o:p></o:p></p>
<p class="m6018092740948470009msonospacing"> <o:p></o:p></p>
<p class="m6018092740948470009msonospacing">You are right I agree with you. To be honest I created the netCDF files in Matlab. The missing values were originally NaN values but I converted them to single -999 by:<o:p></o:p></p>
<p class="m6018092740948470009msonospacing">data(isnan(data))=single(-999);<o:p></o:p></p>
<p class="m6018092740948470009msonospacing">pr = permute(data,[2,3,1]);<o:p></o:p></p>
<p class="m6018092740948470009msonospacing">pr = single(pr);<o:p></o:p></p>
<p class="m6018092740948470009msonospacing"> <o:p></o:p></p>
<p class="m6018092740948470009msonospacing">and I defined -999 as _FillValue by:<o:p></o:p></p>
<p class="m6018092740948470009msonospacing">netcdf.putAtt(ncid1,3,'_FillValue', cast(-999, class(pr)));<o:p></o:p></p>
<p class="m6018092740948470009msonospacing"> <o:p></o:p></p>
<p class="m6018092740948470009msonospacing">I checked the netCDF files with ncdump -h which denoted that _FillValue and pr value are both float:<o:p></o:p></p>
<p class="m6018092740948470009msonospacing">$ ncdump -h <a href="http://MPI-ESM-LR-CLM4.8.17_dry_power.nc" target="_blank">
MPI-ESM-LR-CLM4.8.17_dry_power.nc</a><o:p></o:p></p>
<p class="m6018092740948470009msonospacing">netcdf MPI-ESM-LR-CLM4.8.17_dry_power {<o:p></o:p></p>
<p class="m6018092740948470009msonospacing">dimensions:<o:p></o:p></p>
<p class="m6018092740948470009msonospacing"> longitude = 464 ;<o:p></o:p></p>
<p class="m6018092740948470009msonospacing"> latitude = 201 ;<o:p></o:p></p>
<p class="m6018092740948470009msonospacing"> time = 5479 ;<o:p></o:p></p>
<p class="m6018092740948470009msonospacing">variables:<o:p></o:p></p>
<p class="m6018092740948470009msonospacing"> float longitude(longitude) ;<o:p></o:p></p>
<p class="m6018092740948470009msonospacing"> longitude:long_name = "longitude" ;<o:p></o:p></p>
<p class="m6018092740948470009msonospacing"> longitude:standard_name = "lon" ;<o:p></o:p></p>
<p class="m6018092740948470009msonospacing"> longitude:units = "degrees_east" ;<o:p></o:p></p>
<p class="m6018092740948470009msonospacing"> float latitude(latitude) ;<o:p></o:p></p>
<p class="m6018092740948470009msonospacing"> latitude:long_name = "latitude" ;<o:p></o:p></p>
<p class="m6018092740948470009msonospacing"> latitude:standard_name = "lat" ;<o:p></o:p></p>
<p class="m6018092740948470009msonospacing"> latitude:units = "degrees_north" ;<o:p></o:p></p>
<p class="m6018092740948470009msonospacing"> double time(time) ;<o:p></o:p></p>
<p class="m6018092740948470009msonospacing"> time:long_name = "time" ;<o:p></o:p></p>
<p class="m6018092740948470009msonospacing"> time:standard_name = "time" ;<o:p></o:p></p>
<p class="m6018092740948470009msonospacing"> time:units = "days since 1949-12-1 00:00:00" ;<o:p></o:p></p>
<p class="m6018092740948470009msonospacing"> float pr(time, latitude, longitude) ;<o:p></o:p></p>
<p class="m6018092740948470009msonospacing"> pr:_FillValue = -999.f ;<o:p></o:p></p>
<p class="m6018092740948470009msonospacing"> pr:long_name = "precipitation" ;<o:p></o:p></p>
<p class="m6018092740948470009msonospacing"> pr:standard_name = "precipitation" ;<o:p></o:p></p>
<p class="m6018092740948470009msonospacing"> pr:units = "mm" ;<o:p></o:p></p>
<p class="m6018092740948470009msonospacing"> pr:cell_methods = "time: mean" ;<o:p></o:p></p>
<p class="m6018092740948470009msonospacing">}<o:p></o:p></p>
<p class="m6018092740948470009msonospacing">When I tried to calculate the multi year seasonal mean values with CDO I got the following error message:<o:p></o:p></p>
<p class="m6018092740948470009msonospacing">cdo yseasmean <a href="http://MPI-ESM-LR-CLM4.8.17_dry_power.nc" target="_blank">
MPI-ESM-LR-CLM4.8.17_dry_power.nc</a> <a href="http://MPI-ESM-LR-CLM4.8.17_dry_power_sm.nc" target="_blank">
MPI-ESM-LR-CLM4.8.17_dry_power_sm.nc</a><o:p></o:p></p>
<p class="m6018092740948470009msonospacing">cdf_put_vara_double : ncid = 131072 varid = 3 val0 = -999.000000<o:p></o:p></p>
<p class="m6018092740948470009msonospacing">cdf_put_vara_double : varname = pr<o:p></o:p></p>
<p class="m6018092740948470009msonospacing"> <o:p></o:p></p>
<p class="m6018092740948470009msonospacing">Error (cdf_put_vara_double) : NetCDF: Numeric conversion not representable<o:p></o:p></p>
<p class="m6018092740948470009msonospacing">I could carry out this operator without error message if I use –b 64 option but due to this option the values become double precision. I do not understand why CDO considers _Fill Value to be double?<o:p></o:p></p>
<p class="m6018092740948470009msonospacing"> <o:p></o:p></p>
<p class="m6018092740948470009msonospacing">Beata<o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D"> </span><o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D"> </span><o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><b><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri",sans-serif">From:</span></b><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri",sans-serif"> Dennis
Shea [mailto:</span><a href="mailto:shea@ucar.edu" target="_blank"><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri",sans-serif">shea@ucar.edu</span></a><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri",sans-serif">]
<br>
<b>Sent:</b> Monday, June 19, 2017 6:39 PM<br>
<b>To:</b> Mary Haley <</span><a href="mailto:haley@ucar.edu" target="_blank"><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri",sans-serif">haley@ucar.edu</span></a><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri",sans-serif">><br>
<b>Cc:</b> Beáta Szabó-Takács <</span><a href="mailto:szabo.b@czechglobe.cz" target="_blank"><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri",sans-serif">szabo.b@czechglobe.cz</span></a><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri",sans-serif">>;
ncl-talk (</span><a href="mailto:ncl-talk@ucar.edu" target="_blank"><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri",sans-serif">ncl-talk@ucar.edu</span></a><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri",sans-serif">) <</span><a href="mailto:ncl-talk@ucar.edu" target="_blank"><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri",sans-serif">ncl-talk@ucar.edu</span></a><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri",sans-serif">><br>
<b>Subject:</b> Re: [ncl-talk] convert the double missing_value to float one</span><o:p></o:p></p>
<div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> <o:p></o:p></p>
<div>
<div>
<div>
<div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;margin-bottom:12.0pt">A comment about netCDF variables. Since NCL variables follow the netCDF variable model, it holds for NCL variables also.<o:p></o:p></p>
</div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;margin-bottom:12.0pt">By rule, the type of the _FillValue should match the type of the variable with which it is associated.
<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">short variable must have a short _Fillvalue<o:p></o:p></p>
</div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">integer variable ... integer _FillValue<o:p></o:p></p>
</div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">float variable ... float _FillValue<o:p></o:p></p>
</div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;margin-bottom:12.0pt">double variable ... double _FillValue<br>
<br>
---<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> <o:p></o:p></p>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">On Mon, Jun 19, 2017 at 9:38 AM, Mary Haley <<a href="mailto:haley@ucar.edu" target="_blank">haley@ucar.edu</a>> wrote:<o:p></o:p></p>
<blockquote style="border:none;border-left:solid #CCCCCC 1.0pt;padding:0cm 0cm 0cm 6.0pt;margin-left:4.8pt;margin-top:5.0pt;margin-right:0cm;margin-bottom:5.0pt">
<div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">Beata,<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">When going from a 'higher' type to a 'lower' type (double to float, for example), you need to use the reassignment operator (':=') to force the conversion:<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> <o:p></o:p></p>
</div>
<div>
<p class="m6018092740948470009m4923127676996268827gmail-m7531163896496914688msonospacing" style="margin:0cm;margin-bottom:.0001pt">
<span style="font-size:11.0pt;font-family:"Calibri",sans-serif"> p8@missing_value := doubletofloat(p8@missing_value)</span><o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827gmail-m7531163896496914688msonospacing" style="margin:0cm;margin-bottom:.0001pt">
<span style="font-size:11.0pt;font-family:"Calibri",sans-serif"> p8@_FillValue := p8@missing_value</span><o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827gmail-m7531163896496914688msonospacing" style="margin:0cm;margin-bottom:.0001pt">
<span style="font-size:11.0pt;font-family:"Calibri",sans-serif"> </span><o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827gmail-m7531163896496914688msonospacing" style="margin:0cm;margin-bottom:.0001pt">
<span style="font-size:11.0pt;font-family:"Calibri",sans-serif">It might be better to use the "totype" function:</span><o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827gmail-m7531163896496914688msonospacing" style="margin:0cm;margin-bottom:.0001pt">
<span style="font-size:11.0pt;font-family:"Calibri",sans-serif"> </span><o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827gmail-m7531163896496914688msonospacing" style="margin:0cm;margin-bottom:.0001pt">
<span style="font-size:11.0pt;font-family:"Calibri",sans-serif"> p8@_FillValue := totype(p8@_FillValue,typeof(p8))</span><o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827gmail-m7531163896496914688msonospacing" style="margin:0cm;margin-bottom:.0001pt">
<span style="font-size:11.0pt;font-family:"Calibri",sans-serif"> p8@missing_value := totype(p8@missing_value,typeof(p8))</span><o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827gmail-m7531163896496914688msonospacing" style="margin:0cm;margin-bottom:.0001pt">
<span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#888888"> </span><o:p></o:p></p>
</div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="color:#888888">--Mary</span><o:p></o:p></p>
</div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="color:#888888"> </span><o:p></o:p></p>
</div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="color:#888888"> </span><o:p></o:p></p>
</div>
</div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> <o:p></o:p></p>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">On Mon, Jun 19, 2017 at 1:04 AM, Beáta Szabó-Takács <<a href="mailto:szabo.b@czechglobe.cz" target="_blank">szabo.b@czechglobe.cz</a>> wrote:<o:p></o:p></p>
<blockquote style="border:none;border-left:solid #CCCCCC 1.0pt;padding:0cm 0cm 0cm 6.0pt;margin-left:4.8pt;margin-top:5.0pt;margin-right:0cm;margin-bottom:5.0pt">
<div>
<div>
<div>
<div>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
Dear NCL Users,<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
I would like to create a taylor_metrics_table where the missing values are denoted by -999. This part of script is:<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
season = (/ "DJF","JJA" /)<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
nSeason = dimsizes(season)<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
table = new ( (/nCase,nSeason,nSource/), typeof(ratio) )<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
table(0,:,:) = (/CA_bias, CA_biasn/)<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
table(1,:,:) = (/CB_bias, CB_biasn/)<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
table(2,:,:) = (/CC_bias, CC_biasn/)<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
table(3,:,:) = (/CD_bias, CD_biasn/)<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
tt_opt = True<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
tt_opt@tableTitle = "Bias(%)"<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
tt_opt@pltType= "png" ; "eps" [default], "pdf", "ps"<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
; "png", "gif" [if you have ImageMajik 'convert']<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
; tt_opt@color0 = "palegreen2"<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
; tt_opt@color1 = "tomato2"<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
tt_opt@color0 = "white"<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
tt_opt@color1 = "white"<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
taylor_metrics_table("taylor_bias_dry", source, case ,season, table, tt_opt)<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
I have some netcdf files where the precipitation (pr) values are float data type, but the _FillValue and missing_value are double. I tried convert these data to float with:<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
ncap2 -s 'pr=float(pr)' CM5A-MR_RCA4_dry_power_sm.nc CM5A-MR_RCA4_dry_power_sm2.nc<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
It converted the _FillValue to float but the missing_value are still double:<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
netcdf CM5A-MR_RCA4_dry_power_sm2 {<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
dimensions:<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
time = UNLIMITED ; // (4 currently)<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
latitude = 201 ;<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
longitude = 464 ;<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
variables:<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
float pr(time, latitude, longitude) ;<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
pr:_FillValue = -999.f ;<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
pr:cell_methods = "time: mean" ;<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
pr:long_name = "precipitation" ;<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
pr:missing_value = -999. ;<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
pr:standard_name = "precipitation_flux" ;<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
pr:units = "mm" ;<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
float longitude(longitude) ;<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
longitude:standard_name = "longitude" ;<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
longitude:long_name = "longitude" ;<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
longitude:units = "degrees_east" ;<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
longitude:axis = "X" ;<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
float latitude(latitude) ;<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
latitude:standard_name = "latitude" ;<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
latitude:long_name = "latitude" ;<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
latitude:units = "degrees_north" ;<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
latitude:axis = "Y" ;<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
double time(time) ;<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
time:standard_name = "time" ;<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
time:long_name = "time" ;<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
time:units = "days since 1949-12-1 00:00:00" ;<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
time:calendar = "standard" ;<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
time:axis = "T" ;<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
The resulted table is attached. In the table the missing values are denoted by double values in CD_biasn despite the fact that the printed CD_biasn contains _FillValues:<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
Variable: CD_biasn<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
Type: float<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
Total Size: 20 bytes<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
5 values<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
Number of Dimensions: 1<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
Dimensions and sizes: [5]<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
Coordinates:<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
Number Of Attributes: 1<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
_FillValue : -999<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
(0) -999<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
(1) -4.541727<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
(2) 12.45617<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
(3) 10.47692<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
(4) -999<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
I also tried to convert the missing values to float type by ncl doubletofloat function but it does not work:<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
p8@missing_value=doubletofloat(p8@missing_value)<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
p8@_FillValue = p8@missing_value<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
delete(p8@missing_value)<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
fatal:Type Mismatch: The type of missing value could not be converted to type of variable (p8)<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
p8_out = addfile("<a href="http://p8.nc" target="_blank">p8.nc</a>","c")<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
p8_out ->p8 = p8<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
ncdump -h <a href="http://p8.nc" target="_blank">p8.nc</a><o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
netcdf p8 {<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
dimensions:<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
time = 4 ;<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
latitude = 201 ;<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
longitude = 464 ;<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
variables:<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
float p8(time, latitude, longitude) ;<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
p8:_FillValue = -999.f ;<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
p8:cell_methods = "time: mean" ;<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
p8:long_name = "precipitation" ;<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
p8:missing_value = -999. ;<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
p8:standard_name = "precipitation_flux" ;<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
p8:units = "mm" ;<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
double time(time) ;<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
time:standard_name = "time" ;<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
time:long_name = "time" ;<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
time:units = "days since 1949-12-1 00:00:00" ;<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
time:calendar = "standard" ;<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
time:axis = "T" ;<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
float latitude(latitude) ;<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
latitude:standard_name = "latitude" ;<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
latitude:long_name = "latitude" ;<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
latitude:units = "degrees_north" ;<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
latitude:axis = "Y" ;<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
float longitude(longitude) ;<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
longitude:standard_name = "longitude" ;<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
longitude:long_name = "longitude" ;<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
longitude:units = "degrees_east" ;<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
longitude:axis = "X" ;<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
}<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
Can someone suggest me a solution?<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
Thank you for your help in advance!<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
Kind regards,<o:p></o:p></p>
<p class="m6018092740948470009m4923127676996268827m7531163896496914688msonospacing">
Beata<o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> <o:p></o:p></p>
</div>
</div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> <o:p></o:p></p>
</div>
</div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;margin-bottom:12.0pt">_______________________________________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><o:p></o:p></p>
</blockquote>
</div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> <o:p></o:p></p>
</div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;margin-bottom:12.0pt"><br>
_______________________________________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><o:p></o:p></p>
</blockquote>
</div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> <o:p></o:p></p>
</div>
</div>
</div>
</div>
</div>
</blockquote>
</div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
</div>
</body>
</html>