<div dir="ltr"><div><div><div><div>Hi,<br><br></div>I'm having trouble reproducing the situation. However, you have these lines:<br><br><p class="gmail-m_7838753465502777532MsoNoSpacing"> p8@missing_value=<wbr>doubletofloat(p8@missing_value<wbr>)</p>
<p class="gmail-m_7838753465502777532MsoNoSpacing"> p8@_FillValue = p8@missing_value</p>
<p class="gmail-m_7838753465502777532MsoNoSpacing"> delete(p8@missing_value)</p>
<br></div>In my tests, it was possible to delete the missing_value attribute first, and then reassign it as a float value. Note that the line: <br><br> p8@missing_value=<wbr>doubletofloat(p8@missing_value<wbr>)<br><br></div>would not work, because as written, the missing_value is of type double, and even through doubletofloat() would return a float, the assignment would coerce it back to double. I think you need to delete the attribute and recreate it.<br><br></div>Rick<br><br></div><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Jun 19, 2017 at 1:04 AM, Beáta Szabó-Takács <span dir="ltr"><<a href="mailto:szabo.b@czechglobe.cz" target="_blank">szabo.b@czechglobe.cz</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div link="#0563C1" vlink="#954F72" lang="HU">
<div class="m_7838753465502777532WordSection1">
<p class="m_7838753465502777532MsoNoSpacing">Dear NCL Users,<u></u><u></u></p>
<p class="m_7838753465502777532MsoNoSpacing"><u></u> <u></u></p>
<p class="m_7838753465502777532MsoNoSpacing">I would like to create a taylor_metrics_table where the missing values are denoted by -999. This part of script is:<u></u><u></u></p>
<p class="m_7838753465502777532MsoNoSpacing"><u></u> <u></u></p>
<p class="m_7838753465502777532MsoNoSpacing"> season = (/ "DJF","JJA" /)<u></u><u></u></p>
<p class="m_7838753465502777532MsoNoSpacing"> nSeason = dimsizes(season)<u></u><u></u></p>
<p class="m_7838753465502777532MsoNoSpacing"><u></u> <u></u></p>
<p class="m_7838753465502777532MsoNoSpacing"> table = new ( (/nCase,nSeason,nSource/), typeof(ratio) )<u></u><u></u></p>
<p class="m_7838753465502777532MsoNoSpacing"> table(0,:,:) = (/CA_bias, CA_biasn/)<u></u><u></u></p>
<p class="m_7838753465502777532MsoNoSpacing"> table(1,:,:) = (/CB_bias, CB_biasn/)<u></u><u></u></p>
<p class="m_7838753465502777532MsoNoSpacing"> table(2,:,:) = (/CC_bias, CC_biasn/)<u></u><u></u></p>
<p class="m_7838753465502777532MsoNoSpacing"> table(3,:,:) = (/CD_bias, CD_biasn/)<u></u><u></u></p>
<p class="m_7838753465502777532MsoNoSpacing"><u></u> <u></u></p>
<p class="m_7838753465502777532MsoNoSpacing"><u></u> <u></u></p>
<p class="m_7838753465502777532MsoNoSpacing"> tt_opt = True<u></u><u></u></p>
<p class="m_7838753465502777532MsoNoSpacing"> tt_opt@tableTitle = "Bias(%)"<u></u><u></u></p>
<p class="m_7838753465502777532MsoNoSpacing"> tt_opt@pltType= "png" ; "eps" [default], "pdf", "ps"<u></u><u></u></p>
<p class="m_7838753465502777532MsoNoSpacing"> <wbr> ; "png", "gif" [if you have ImageMajik 'convert']<u></u><u></u></p>
<p class="m_7838753465502777532MsoNoSpacing">; tt_opt@color0 = "palegreen2"<u></u><u></u></p>
<p class="m_7838753465502777532MsoNoSpacing">; tt_opt@color1 = "tomato2"<u></u><u></u></p>
<p class="m_7838753465502777532MsoNoSpacing"> tt_opt@color0 = "white"<u></u><u></u></p>
<p class="m_7838753465502777532MsoNoSpacing"> tt_opt@color1 = "white"<u></u><u></u></p>
<p class="m_7838753465502777532MsoNoSpacing"><u></u> <u></u></p>
<p class="m_7838753465502777532MsoNoSpacing"><u></u> <u></u></p>
<p class="m_7838753465502777532MsoNoSpacing"> taylor_metrics_table("taylor_<wbr>bias_dry", source, case ,season, table, tt_opt)<u></u><u></u></p>
<p class="m_7838753465502777532MsoNoSpacing"><u></u> <u></u></p>
<p class="m_7838753465502777532MsoNoSpacing">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:<u></u><u></u></p>
<p class="m_7838753465502777532MsoNoSpacing"><u></u> <u></u></p>
<p class="m_7838753465502777532MsoNoSpacing">ncap2 -s 'pr=float(pr)' CM5A-MR_RCA4_dry_power_sm.nc CM5A-MR_RCA4_dry_power_sm2.nc<u></u><u></u></p>
<p class="m_7838753465502777532MsoNoSpacing"><u></u> <u></u></p>
<p class="m_7838753465502777532MsoNoSpacing">It converted the _FillValue to float but the missing_value are still double:<u></u><u></u></p>
<p class="m_7838753465502777532MsoNoSpacing">netcdf CM5A-MR_RCA4_dry_power_sm2 {<u></u><u></u></p>
<p class="m_7838753465502777532MsoNoSpacing">dimensions:<u></u><u></u></p>
<p class="m_7838753465502777532MsoNoSpacing"> time = UNLIMITED ; // (4 currently)<u></u><u></u></p>
<p class="m_7838753465502777532MsoNoSpacing"> latitude = 201 ;<u></u><u></u></p>
<p class="m_7838753465502777532MsoNoSpacing"> longitude = 464 ;<u></u><u></u></p>
<p class="m_7838753465502777532MsoNoSpacing">variables:<u></u><u></u></p>
<p class="m_7838753465502777532MsoNoSpacing"> float pr(time, latitude, longitude) ;<u></u><u></u></p>
<p class="m_7838753465502777532MsoNoSpacing"> pr:_FillValue = -999.f ;<u></u><u></u></p>
<p class="m_7838753465502777532MsoNoSpacing"> pr:cell_methods = "time: mean" ;<u></u><u></u></p>
<p class="m_7838753465502777532MsoNoSpacing"> pr:long_name = "precipitation" ;<u></u><u></u></p>
<p class="m_7838753465502777532MsoNoSpacing"> pr:missing_value = -999. ;<u></u><u></u></p>
<p class="m_7838753465502777532MsoNoSpacing"> pr:standard_name = "precipitation_flux" ;<u></u><u></u></p>
<p class="m_7838753465502777532MsoNoSpacing"> pr:units = "mm" ;<u></u><u></u></p>
<p class="m_7838753465502777532MsoNoSpacing"> float longitude(longitude) ;<u></u><u></u></p>
<p class="m_7838753465502777532MsoNoSpacing"> longitude:standard_name = "longitude" ;<u></u><u></u></p>
<p class="m_7838753465502777532MsoNoSpacing"> longitude:long_name = "longitude" ;<u></u><u></u></p>
<p class="m_7838753465502777532MsoNoSpacing"> longitude:units = "degrees_east" ;<u></u><u></u></p>
<p class="m_7838753465502777532MsoNoSpacing"> longitude:axis = "X" ;<u></u><u></u></p>
<p class="m_7838753465502777532MsoNoSpacing"> float latitude(latitude) ;<u></u><u></u></p>
<p class="m_7838753465502777532MsoNoSpacing"> latitude:standard_name = "latitude" ;<u></u><u></u></p>
<p class="m_7838753465502777532MsoNoSpacing"> latitude:long_name = "latitude" ;<u></u><u></u></p>
<p class="m_7838753465502777532MsoNoSpacing"> latitude:units = "degrees_north" ;<u></u><u></u></p>
<p class="m_7838753465502777532MsoNoSpacing"> latitude:axis = "Y" ;<u></u><u></u></p>
<p class="m_7838753465502777532MsoNoSpacing"> double time(time) ;<u></u><u></u></p>
<p class="m_7838753465502777532MsoNoSpacing"> time:standard_name = "time" ;<u></u><u></u></p>
<p class="m_7838753465502777532MsoNoSpacing"> time:long_name = "time" ;<u></u><u></u></p>
<p class="m_7838753465502777532MsoNoSpacing"> time:units = "days since 1949-12-1 00:00:00" ;<u></u><u></u></p>
<p class="m_7838753465502777532MsoNoSpacing"> time:calendar = "standard" ;<u></u><u></u></p>
<p class="m_7838753465502777532MsoNoSpacing"> time:axis = "T" ;<u></u><u></u></p>
<p class="m_7838753465502777532MsoNoSpacing"><u></u> <u></u></p>
<p class="m_7838753465502777532MsoNoSpacing">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:<u></u><u></u></p>
<p class="m_7838753465502777532MsoNoSpacing"><u></u> <u></u></p>
<p class="m_7838753465502777532MsoNoSpacing">Variable: CD_biasn<u></u><u></u></p>
<p class="m_7838753465502777532MsoNoSpacing">Type: float<u></u><u></u></p>
<p class="m_7838753465502777532MsoNoSpacing">Total Size: 20 bytes<u></u><u></u></p>
<p class="m_7838753465502777532MsoNoSpacing"> 5 values<u></u><u></u></p>
<p class="m_7838753465502777532MsoNoSpacing">Number of Dimensions: 1<u></u><u></u></p>
<p class="m_7838753465502777532MsoNoSpacing">Dimensions and sizes: [5]<u></u><u></u></p>
<p class="m_7838753465502777532MsoNoSpacing">Coordinates:<u></u><u></u></p>
<p class="m_7838753465502777532MsoNoSpacing">Number Of Attributes: 1<u></u><u></u></p>
<p class="m_7838753465502777532MsoNoSpacing"> _FillValue : -999<u></u><u></u></p>
<p class="m_7838753465502777532MsoNoSpacing">(0) -999<u></u><u></u></p>
<p class="m_7838753465502777532MsoNoSpacing">(1) -4.541727<u></u><u></u></p>
<p class="m_7838753465502777532MsoNoSpacing">(2) 12.45617<u></u><u></u></p>
<p class="m_7838753465502777532MsoNoSpacing">(3) 10.47692<u></u><u></u></p>
<p class="m_7838753465502777532MsoNoSpacing">(4) -999<u></u><u></u></p>
<p class="m_7838753465502777532MsoNoSpacing"><u></u> <u></u></p>
<p class="m_7838753465502777532MsoNoSpacing">I also tried to convert the missing values to float type by ncl doubletofloat function but it does not work:<u></u><u></u></p>
<p class="m_7838753465502777532MsoNoSpacing"><u></u> <u></u></p>
<p class="m_7838753465502777532MsoNoSpacing">p8@missing_value=<wbr>doubletofloat(p8@missing_value<wbr>)<u></u><u></u></p>
<p class="m_7838753465502777532MsoNoSpacing"> p8@_FillValue = p8@missing_value<u></u><u></u></p>
<p class="m_7838753465502777532MsoNoSpacing"> delete(p8@missing_value)<u></u><u></u></p>
<p class="m_7838753465502777532MsoNoSpacing"><u></u> <u></u></p>
<p class="m_7838753465502777532MsoNoSpacing">fatal:Type Mismatch: The type of missing value could not be converted to type of variable (p8)<u></u><u></u></p>
<p class="m_7838753465502777532MsoNoSpacing"><u></u> <u></u></p>
<p class="m_7838753465502777532MsoNoSpacing"> p8_out = addfile("<a href="http://p8.nc" target="_blank">p8.nc</a>","c")<u></u><u></u></p>
<p class="m_7838753465502777532MsoNoSpacing"> p8_out ->p8 = p8<u></u><u></u></p>
<p class="m_7838753465502777532MsoNoSpacing"><u></u> <u></u></p>
<p class="m_7838753465502777532MsoNoSpacing"><u></u> <u></u></p>
<p class="m_7838753465502777532MsoNoSpacing">ncdump -h <a href="http://p8.nc" target="_blank">p8.nc</a><u></u><u></u></p>
<p class="m_7838753465502777532MsoNoSpacing">netcdf p8 {<u></u><u></u></p>
<p class="m_7838753465502777532MsoNoSpacing">dimensions:<u></u><u></u></p>
<p class="m_7838753465502777532MsoNoSpacing"> time = 4 ;<u></u><u></u></p>
<p class="m_7838753465502777532MsoNoSpacing"> latitude = 201 ;<u></u><u></u></p>
<p class="m_7838753465502777532MsoNoSpacing"> longitude = 464 ;<u></u><u></u></p>
<p class="m_7838753465502777532MsoNoSpacing">variables:<u></u><u></u></p>
<p class="m_7838753465502777532MsoNoSpacing"> float p8(time, latitude, longitude) ;<u></u><u></u></p>
<p class="m_7838753465502777532MsoNoSpacing"> p8:_FillValue = -999.f ;<u></u><u></u></p>
<p class="m_7838753465502777532MsoNoSpacing"> p8:cell_methods = "time: mean" ;<u></u><u></u></p>
<p class="m_7838753465502777532MsoNoSpacing"> p8:long_name = "precipitation" ;<u></u><u></u></p>
<p class="m_7838753465502777532MsoNoSpacing"> p8:missing_value = -999. ;<u></u><u></u></p>
<p class="m_7838753465502777532MsoNoSpacing"> p8:standard_name = "precipitation_flux" ;<u></u><u></u></p>
<p class="m_7838753465502777532MsoNoSpacing"> p8:units = "mm" ;<u></u><u></u></p>
<p class="m_7838753465502777532MsoNoSpacing"> double time(time) ;<u></u><u></u></p>
<p class="m_7838753465502777532MsoNoSpacing"> time:standard_name = "time" ;<u></u><u></u></p>
<p class="m_7838753465502777532MsoNoSpacing"> time:long_name = "time" ;<u></u><u></u></p>
<p class="m_7838753465502777532MsoNoSpacing"> time:units = "days since 1949-12-1 00:00:00" ;<u></u><u></u></p>
<p class="m_7838753465502777532MsoNoSpacing"> time:calendar = "standard" ;<u></u><u></u></p>
<p class="m_7838753465502777532MsoNoSpacing"> time:axis = "T" ;<u></u><u></u></p>
<p class="m_7838753465502777532MsoNoSpacing"> float latitude(latitude) ;<u></u><u></u></p>
<p class="m_7838753465502777532MsoNoSpacing"> latitude:standard_name = "latitude" ;<u></u><u></u></p>
<p class="m_7838753465502777532MsoNoSpacing"> latitude:long_name = "latitude" ;<u></u><u></u></p>
<p class="m_7838753465502777532MsoNoSpacing"> latitude:units = "degrees_north" ;<u></u><u></u></p>
<p class="m_7838753465502777532MsoNoSpacing"> latitude:axis = "Y" ;<u></u><u></u></p>
<p class="m_7838753465502777532MsoNoSpacing"> float longitude(longitude) ;<u></u><u></u></p>
<p class="m_7838753465502777532MsoNoSpacing"> longitude:standard_name = "longitude" ;<u></u><u></u></p>
<p class="m_7838753465502777532MsoNoSpacing"> longitude:long_name = "longitude" ;<u></u><u></u></p>
<p class="m_7838753465502777532MsoNoSpacing"> longitude:units = "degrees_east" ;<u></u><u></u></p>
<p class="m_7838753465502777532MsoNoSpacing"> longitude:axis = "X" ;<u></u><u></u></p>
<p class="m_7838753465502777532MsoNoSpacing">}<u></u><u></u></p>
<p class="m_7838753465502777532MsoNoSpacing"><u></u> <u></u></p>
<p class="m_7838753465502777532MsoNoSpacing"><u></u> <u></u></p>
<p class="m_7838753465502777532MsoNoSpacing">Can someone suggest me a solution?<u></u><u></u></p>
<p class="m_7838753465502777532MsoNoSpacing">Thank you for your help in advance!<u></u><u></u></p>
<p class="m_7838753465502777532MsoNoSpacing">Kind regards,<u></u><u></u></p>
<p class="m_7838753465502777532MsoNoSpacing">Beata<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
</div>
<br>______________________________<wbr>_________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/<wbr>mailman/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div>