<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;}
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;}
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="MsoNoSpacing">Dear NCL Users,<o:p></o:p></p>
<p class="MsoNoSpacing"><o:p> </o:p></p>
<p class="MsoNoSpacing">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="MsoNoSpacing"><o:p> </o:p></p>
<p class="MsoNoSpacing"> season = (/ "DJF","JJA" /)<o:p></o:p></p>
<p class="MsoNoSpacing"> nSeason = dimsizes(season)<o:p></o:p></p>
<p class="MsoNoSpacing"><o:p> </o:p></p>
<p class="MsoNoSpacing"> table = new ( (/nCase,nSeason,nSource/), typeof(ratio) )<o:p></o:p></p>
<p class="MsoNoSpacing"> table(0,:,:) = (/CA_bias, CA_biasn/)<o:p></o:p></p>
<p class="MsoNoSpacing"> table(1,:,:) = (/CB_bias, CB_biasn/)<o:p></o:p></p>
<p class="MsoNoSpacing"> table(2,:,:) = (/CC_bias, CC_biasn/)<o:p></o:p></p>
<p class="MsoNoSpacing"> table(3,:,:) = (/CD_bias, CD_biasn/)<o:p></o:p></p>
<p class="MsoNoSpacing"><o:p> </o:p></p>
<p class="MsoNoSpacing"><o:p> </o:p></p>
<p class="MsoNoSpacing"> tt_opt = True<o:p></o:p></p>
<p class="MsoNoSpacing"> tt_opt@tableTitle = "Bias(%)"<o:p></o:p></p>
<p class="MsoNoSpacing"> tt_opt@pltType= "png" ; "eps" [default], "pdf", "ps"<o:p></o:p></p>
<p class="MsoNoSpacing"> ; "png", "gif" [if you have ImageMajik 'convert']<o:p></o:p></p>
<p class="MsoNoSpacing">; tt_opt@color0 = "palegreen2"<o:p></o:p></p>
<p class="MsoNoSpacing">; tt_opt@color1 = "tomato2"<o:p></o:p></p>
<p class="MsoNoSpacing"> tt_opt@color0 = "white"<o:p></o:p></p>
<p class="MsoNoSpacing"> tt_opt@color1 = "white"<o:p></o:p></p>
<p class="MsoNoSpacing"><o:p> </o:p></p>
<p class="MsoNoSpacing"><o:p> </o:p></p>
<p class="MsoNoSpacing"> taylor_metrics_table("taylor_bias_dry", source, case ,season, table, tt_opt)<o:p></o:p></p>
<p class="MsoNoSpacing"><o:p> </o:p></p>
<p class="MsoNoSpacing">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="MsoNoSpacing"><o:p> </o:p></p>
<p class="MsoNoSpacing">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="MsoNoSpacing"><o:p> </o:p></p>
<p class="MsoNoSpacing">It converted the _FillValue to float but the missing_value are still double:<o:p></o:p></p>
<p class="MsoNoSpacing">netcdf CM5A-MR_RCA4_dry_power_sm2 {<o:p></o:p></p>
<p class="MsoNoSpacing">dimensions:<o:p></o:p></p>
<p class="MsoNoSpacing"> time = UNLIMITED ; // (4 currently)<o:p></o:p></p>
<p class="MsoNoSpacing"> latitude = 201 ;<o:p></o:p></p>
<p class="MsoNoSpacing"> longitude = 464 ;<o:p></o:p></p>
<p class="MsoNoSpacing">variables:<o:p></o:p></p>
<p class="MsoNoSpacing"> float pr(time, latitude, longitude) ;<o:p></o:p></p>
<p class="MsoNoSpacing"> pr:_FillValue = -999.f ;<o:p></o:p></p>
<p class="MsoNoSpacing"> pr:cell_methods = "time: mean" ;<o:p></o:p></p>
<p class="MsoNoSpacing"> pr:long_name = "precipitation" ;<o:p></o:p></p>
<p class="MsoNoSpacing"> pr:missing_value = -999. ;<o:p></o:p></p>
<p class="MsoNoSpacing"> pr:standard_name = "precipitation_flux" ;<o:p></o:p></p>
<p class="MsoNoSpacing"> pr:units = "mm" ;<o:p></o:p></p>
<p class="MsoNoSpacing"> float longitude(longitude) ;<o:p></o:p></p>
<p class="MsoNoSpacing"> longitude:standard_name = "longitude" ;<o:p></o:p></p>
<p class="MsoNoSpacing"> longitude:long_name = "longitude" ;<o:p></o:p></p>
<p class="MsoNoSpacing"> longitude:units = "degrees_east" ;<o:p></o:p></p>
<p class="MsoNoSpacing"> longitude:axis = "X" ;<o:p></o:p></p>
<p class="MsoNoSpacing"> float latitude(latitude) ;<o:p></o:p></p>
<p class="MsoNoSpacing"> latitude:standard_name = "latitude" ;<o:p></o:p></p>
<p class="MsoNoSpacing"> latitude:long_name = "latitude" ;<o:p></o:p></p>
<p class="MsoNoSpacing"> latitude:units = "degrees_north" ;<o:p></o:p></p>
<p class="MsoNoSpacing"> latitude:axis = "Y" ;<o:p></o:p></p>
<p class="MsoNoSpacing"> double time(time) ;<o:p></o:p></p>
<p class="MsoNoSpacing"> time:standard_name = "time" ;<o:p></o:p></p>
<p class="MsoNoSpacing"> time:long_name = "time" ;<o:p></o:p></p>
<p class="MsoNoSpacing"> time:units = "days since 1949-12-1 00:00:00" ;<o:p></o:p></p>
<p class="MsoNoSpacing"> time:calendar = "standard" ;<o:p></o:p></p>
<p class="MsoNoSpacing"> time:axis = "T" ;<o:p></o:p></p>
<p class="MsoNoSpacing"><o:p> </o:p></p>
<p class="MsoNoSpacing">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="MsoNoSpacing"><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) -4.541727<o:p></o:p></p>
<p class="MsoNoSpacing">(2) 12.45617<o:p></o:p></p>
<p class="MsoNoSpacing">(3) 10.47692<o:p></o:p></p>
<p class="MsoNoSpacing">(4) -999<o:p></o:p></p>
<p class="MsoNoSpacing"><o:p> </o:p></p>
<p class="MsoNoSpacing">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="MsoNoSpacing"><o:p> </o:p></p>
<p class="MsoNoSpacing">p8@missing_value=doubletofloat(p8@missing_value)<o:p></o:p></p>
<p class="MsoNoSpacing"> p8@_FillValue = p8@missing_value<o:p></o:p></p>
<p class="MsoNoSpacing"> delete(p8@missing_value)<o:p></o:p></p>
<p class="MsoNoSpacing"><o:p> </o:p></p>
<p class="MsoNoSpacing">fatal:Type Mismatch: The type of missing value could not be converted to type of variable (p8)<o:p></o:p></p>
<p class="MsoNoSpacing"><o:p> </o:p></p>
<p class="MsoNoSpacing"> p8_out = addfile("p8.nc","c")<o:p></o:p></p>
<p class="MsoNoSpacing"> p8_out ->p8 = p8<o:p></o:p></p>
<p class="MsoNoSpacing"><o:p> </o:p></p>
<p class="MsoNoSpacing"><o:p> </o:p></p>
<p class="MsoNoSpacing">ncdump -h p8.nc<o:p></o:p></p>
<p class="MsoNoSpacing">netcdf p8 {<o:p></o:p></p>
<p class="MsoNoSpacing">dimensions:<o:p></o:p></p>
<p class="MsoNoSpacing"> time = 4 ;<o:p></o:p></p>
<p class="MsoNoSpacing"> latitude = 201 ;<o:p></o:p></p>
<p class="MsoNoSpacing"> longitude = 464 ;<o:p></o:p></p>
<p class="MsoNoSpacing">variables:<o:p></o:p></p>
<p class="MsoNoSpacing"> float p8(time, latitude, longitude) ;<o:p></o:p></p>
<p class="MsoNoSpacing"> p8:_FillValue = -999.f ;<o:p></o:p></p>
<p class="MsoNoSpacing"> p8:cell_methods = "time: mean" ;<o:p></o:p></p>
<p class="MsoNoSpacing"> p8:long_name = "precipitation" ;<o:p></o:p></p>
<p class="MsoNoSpacing"> p8:missing_value = -999. ;<o:p></o:p></p>
<p class="MsoNoSpacing"> p8:standard_name = "precipitation_flux" ;<o:p></o:p></p>
<p class="MsoNoSpacing"> p8:units = "mm" ;<o:p></o:p></p>
<p class="MsoNoSpacing"> double time(time) ;<o:p></o:p></p>
<p class="MsoNoSpacing"> time:standard_name = "time" ;<o:p></o:p></p>
<p class="MsoNoSpacing"> time:long_name = "time" ;<o:p></o:p></p>
<p class="MsoNoSpacing"> time:units = "days since 1949-12-1 00:00:00" ;<o:p></o:p></p>
<p class="MsoNoSpacing"> time:calendar = "standard" ;<o:p></o:p></p>
<p class="MsoNoSpacing"> time:axis = "T" ;<o:p></o:p></p>
<p class="MsoNoSpacing"> float latitude(latitude) ;<o:p></o:p></p>
<p class="MsoNoSpacing"> latitude:standard_name = "latitude" ;<o:p></o:p></p>
<p class="MsoNoSpacing"> latitude:long_name = "latitude" ;<o:p></o:p></p>
<p class="MsoNoSpacing"> latitude:units = "degrees_north" ;<o:p></o:p></p>
<p class="MsoNoSpacing"> latitude:axis = "Y" ;<o:p></o:p></p>
<p class="MsoNoSpacing"> float longitude(longitude) ;<o:p></o:p></p>
<p class="MsoNoSpacing"> longitude:standard_name = "longitude" ;<o:p></o:p></p>
<p class="MsoNoSpacing"> longitude:long_name = "longitude" ;<o:p></o:p></p>
<p class="MsoNoSpacing"> longitude:units = "degrees_east" ;<o:p></o:p></p>
<p class="MsoNoSpacing"> longitude:axis = "X" ;<o:p></o:p></p>
<p class="MsoNoSpacing">}<o:p></o:p></p>
<p class="MsoNoSpacing"><o:p> </o:p></p>
<p class="MsoNoSpacing"><o:p> </o:p></p>
<p class="MsoNoSpacing">Can someone suggest me a solution?<o:p></o:p></p>
<p class="MsoNoSpacing">Thank you for your help in advance!<o:p></o:p></p>
<p class="MsoNoSpacing">Kind regards,<o:p></o:p></p>
<p class="MsoNoSpacing">Beata<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
</body>
</html>