<div dir="ltr"><div><div><div><div><div><div><div><div><div><div><div><div><div><div><div><div>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]:<wbr>Execute: Error occurred at or near line 1<br><br></div><div>NCL will 'help' in some cases<br></div><div><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>====<br></div>ncdump -h <a href="http://MPI-ESM-LR-CLM4.8.17_dry_power.nc" target="_blank">MPI-ESM-LR-CLM4.8.17_dry_<wbr>power.nc</a><br><br></div><div>[snip]<br></div><div>float pr(time, latitude, longitude) ;
<br>pr:_FillValue = -999.f ;<br></div><div>[snip]<br></div><div><br></div>Don't be misled by the 'f' ... extension (-999.f). <br><br></div> f = addfile("MPI-ESM-LR-CLM4.8.17_<wbr>dry_power","r")<br></div> pr = f->pr<br></div><div> print(typeof(pr)) <wbr> ; float<br></div> print( typeof(pr@_FillValue)) ; I speculate this will be type double.<br> <br></div><div> <wbr> ; Note the := syntax<br></div><div> pr@_FillValue := totype( pr@_FillValue, typeof(pr) ) ; MH's suggestion<br> print( typeof(pr@_FillValue)) <br></div>[snip]<br><br></div><br><br>====<br></div>You can send the file to<br><br></div>ftp <a href="http://ftp.cgd.ucar.edu" target="_blank">ftp.cgd.ucar.edu</a><br></div>anonymous<br></div>your_email<br></div><div>cd incoming<br></div>put .....<br></div>quit<br><br></div>Then notify ncl-talk after a successful file transfer.<br><br></div><div>You can attach a clean version of your script to the email.<br></div><div><br></div>Cheers<br></div>D<br><div><div><div><div><div><div><div><div><div><br><br><div><div><br><br><div><div><div><br><div><div><br></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Jun 20, 2017 at 1:32 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="blue" vlink="purple" lang="HU">
<div class="m_6018092740948470009WordSection1">
<p class="m_6018092740948470009MsoNoSpacing">Dear Dennis,<u></u><u></u></p>
<p class="m_6018092740948470009MsoNoSpacing"><u></u> <u></u></p>
<p class="m_6018092740948470009MsoNoSpacing">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:<u></u><u></u></p>
<p class="m_6018092740948470009MsoNoSpacing"><span>data(isnan(data))=single(-999)<wbr>;<u></u><u></u></span></p>
<p class="m_6018092740948470009MsoNoSpacing">pr = permute(data,[2,3,1]);<u></u><u></u></p>
<p class="m_6018092740948470009MsoNoSpacing">pr = single(pr);<u></u><u></u></p>
<p class="m_6018092740948470009MsoNoSpacing"><u></u> <u></u></p>
<p class="m_6018092740948470009MsoNoSpacing">and I defined -999 as _FillValue by:<u></u><u></u></p>
<p class="m_6018092740948470009MsoNoSpacing">netcdf.putAtt(ncid1,3,'_<wbr>FillValue', cast(-999, class(pr)));<u></u><u></u></p>
<p class="m_6018092740948470009MsoNoSpacing"><u></u> <u></u></p>
<p class="m_6018092740948470009MsoNoSpacing">I checked the netCDF files with ncdump -h which denoted that _FillValue and pr value are both float:<u></u><u></u></p>
<p class="m_6018092740948470009MsoNoSpacing">$ ncdump -h <a href="http://MPI-ESM-LR-CLM4.8.17_dry_power.nc" target="_blank">MPI-ESM-LR-CLM4.8.17_dry_<wbr>power.nc</a><u></u><u></u></p>
<p class="m_6018092740948470009MsoNoSpacing">netcdf MPI-ESM-LR-CLM4.8.17_dry_power {<u></u><u></u></p>
<p class="m_6018092740948470009MsoNoSpacing">dimensions:<u></u><u></u></p>
<p class="m_6018092740948470009MsoNoSpacing"> longitude = 464 ;<u></u><u></u></p>
<p class="m_6018092740948470009MsoNoSpacing"> latitude = 201 ;<u></u><u></u></p>
<p class="m_6018092740948470009MsoNoSpacing"> time = 5479 ;<u></u><u></u></p>
<p class="m_6018092740948470009MsoNoSpacing">variables:<u></u><u></u></p>
<p class="m_6018092740948470009MsoNoSpacing"> float longitude(longitude) ;<u></u><u></u></p>
<p class="m_6018092740948470009MsoNoSpacing"> longitude:long_name = "longitude" ;<u></u><u></u></p>
<p class="m_6018092740948470009MsoNoSpacing"> longitude:standard_name = "lon" ;<u></u><u></u></p>
<p class="m_6018092740948470009MsoNoSpacing"> longitude:units = "degrees_east" ;<u></u><u></u></p>
<p class="m_6018092740948470009MsoNoSpacing"> float latitude(latitude) ;<u></u><u></u></p>
<p class="m_6018092740948470009MsoNoSpacing"> latitude:long_name = "latitude" ;<u></u><u></u></p>
<p class="m_6018092740948470009MsoNoSpacing"> latitude:standard_name = "lat" ;<u></u><u></u></p>
<p class="m_6018092740948470009MsoNoSpacing"> latitude:units = "degrees_north" ;<u></u><u></u></p>
<p class="m_6018092740948470009MsoNoSpacing"> double time(time) ;<u></u><u></u></p>
<p class="m_6018092740948470009MsoNoSpacing"> time:long_name = "time" ;<u></u><u></u></p>
<p class="m_6018092740948470009MsoNoSpacing"> time:standard_name = "time" ;<u></u><u></u></p><span class="">
<p class="m_6018092740948470009MsoNoSpacing"> time:units = "days since 1949-12-1 00:00:00" ;<u></u><u></u></p>
</span><p class="m_6018092740948470009MsoNoSpacing"> float pr(time, latitude, longitude) ;<u></u><u></u></p>
<p class="m_6018092740948470009MsoNoSpacing"> pr:_FillValue = -999.f ;<u></u><u></u></p>
<p class="m_6018092740948470009MsoNoSpacing"> pr:long_name = "precipitation" ;<u></u><u></u></p>
<p class="m_6018092740948470009MsoNoSpacing"> pr:standard_name = "precipitation" ;<u></u><u></u></p>
<p class="m_6018092740948470009MsoNoSpacing"> pr:units = "mm" ;<u></u><u></u></p><span class="">
<p class="m_6018092740948470009MsoNoSpacing"> pr:cell_methods = "time: mean" ;<u></u><u></u></p>
<p class="m_6018092740948470009MsoNoSpacing">}<u></u><u></u></p>
</span><p class="m_6018092740948470009MsoNoSpacing">When I tried to calculate the multi year seasonal mean values with CDO I got the following error message:<u></u><u></u></p>
<p class="m_6018092740948470009MsoNoSpacing">cdo yseasmean <a href="http://MPI-ESM-LR-CLM4.8.17_dry_power.nc" target="_blank">MPI-ESM-LR-CLM4.8.17_dry_<wbr>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_<wbr>power_sm.nc</a><u></u><u></u></p>
<p class="m_6018092740948470009MsoNoSpacing">cdf_put_vara_double : ncid = 131072 varid = 3 val0 = -999.000000<u></u><u></u></p>
<p class="m_6018092740948470009MsoNoSpacing">cdf_put_vara_double : varname = pr<u></u><u></u></p>
<p class="m_6018092740948470009MsoNoSpacing"><u></u> <u></u></p>
<p class="m_6018092740948470009MsoNoSpacing">Error (cdf_put_vara_double) : NetCDF: Numeric conversion not representable<u></u><u></u></p>
<p class="m_6018092740948470009MsoNoSpacing">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?<u></u><u></u></p>
<p class="m_6018092740948470009MsoNoSpacing"><u></u> <u></u></p>
<p class="m_6018092740948470009MsoNoSpacing">Beata<u></u><u></u></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d"><u></u> <u></u></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d"><u></u> <u></u></span></p>
<p class="MsoNormal"><b><span style="font-size:11.0pt;font-family:"Calibri",sans-serif" lang="EN-US">From:</span></b><span style="font-size:11.0pt;font-family:"Calibri",sans-serif" lang="EN-US"> Dennis Shea [mailto:<a href="mailto:shea@ucar.edu" target="_blank">shea@ucar.edu</a>]
<br>
<b>Sent:</b> Monday, June 19, 2017 6:39 PM<br>
<b>To:</b> Mary Haley <<a href="mailto:haley@ucar.edu" target="_blank">haley@ucar.edu</a>><br>
<b>Cc:</b> Beáta Szabó-Takács <<a href="mailto:szabo.b@czechglobe.cz" target="_blank">szabo.b@czechglobe.cz</a>>; ncl-talk (<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>) <<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>><br>
<b>Subject:</b> Re: [ncl-talk] convert the double missing_value to float one<u></u><u></u></span></p><div><div class="h5">
<p class="MsoNormal"><u></u> <u></u></p>
<div>
<div>
<div>
<div>
<div>
<p class="MsoNormal" style="margin-bottom:12.0pt">A comment about netCDF variables. Since NCL variables follow the netCDF variable model, it holds for NCL variables also.<u></u><u></u></p>
</div>
<p class="MsoNormal" style="margin-bottom:12.0pt">By rule, the type of the _FillValue should match the type of the variable with which it is associated.
<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal">short variable must have a short _Fillvalue<u></u><u></u></p>
</div>
<p class="MsoNormal">integer variable ... integer _FillValue<u></u><u></u></p>
</div>
<p class="MsoNormal">float variable ... float _FillValue<u></u><u></u></p>
</div>
<p class="MsoNormal" style="margin-bottom:12.0pt">double variable ... double _FillValue<br>
<br>
---<br>
<br>
<br>
<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
<div>
<p class="MsoNormal">On Mon, Jun 19, 2017 at 9:38 AM, Mary Haley <<a href="mailto:haley@ucar.edu" target="_blank">haley@ucar.edu</a>> wrote:<u></u><u></u></p>
<blockquote style="border:none;border-left:solid #cccccc 1.0pt;padding:0cm 0cm 0cm 6.0pt;margin-left:4.8pt;margin-right:0cm">
<div>
<div>
<p class="MsoNormal">Beata,<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">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:<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="m_6018092740948470009m4923127676996268827gmail-m7531163896496914688msonospacing" style="margin:0cm;margin-bottom:.0001pt">
<span style="font-size:11.0pt;font-family:"Calibri",sans-serif"> p8@missing_value := doubletofloat(p8@missing_<wbr>value)<u></u><u></u></span></p>
<p class="m_6018092740948470009m4923127676996268827gmail-m7531163896496914688msonospacing" style="margin:0cm;margin-bottom:.0001pt">
<span style="font-size:11.0pt;font-family:"Calibri",sans-serif"> p8@_FillValue := p8@missing_value<u></u><u></u></span></p>
<p class="m_6018092740948470009m4923127676996268827gmail-m7531163896496914688msonospacing" style="margin:0cm;margin-bottom:.0001pt">
<span style="font-size:11.0pt;font-family:"Calibri",sans-serif"><u></u> <u></u></span></p>
<p class="m_6018092740948470009m4923127676996268827gmail-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:<u></u><u></u></span></p>
<p class="m_6018092740948470009m4923127676996268827gmail-m7531163896496914688msonospacing" style="margin:0cm;margin-bottom:.0001pt">
<span style="font-size:11.0pt;font-family:"Calibri",sans-serif"><u></u> <u></u></span></p>
<p class="m_6018092740948470009m4923127676996268827gmail-m7531163896496914688msonospacing" style="margin:0cm;margin-bottom:.0001pt">
<span style="font-size:11.0pt;font-family:"Calibri",sans-serif"> p8@_FillValue := totype(p8@_FillValue,<wbr>typeof(p8))<u></u><u></u></span></p>
<p class="m_6018092740948470009m4923127676996268827gmail-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,<wbr>typeof(p8))<u></u><u></u></span></p>
<p class="m_6018092740948470009m4923127676996268827gmail-m7531163896496914688msonospacing" style="margin:0cm;margin-bottom:.0001pt">
<span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#888888"><u></u> <u></u></span></p>
</div>
<div>
<p class="MsoNormal"><span style="color:#888888">--Mary<u></u><u></u></span></p>
</div>
<div>
<p class="MsoNormal"><span style="color:#888888"><u></u> <u></u></span></p>
</div>
<div>
<p class="MsoNormal"><span style="color:#888888"><u></u> <u></u></span></p>
</div>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
<div>
<p class="MsoNormal">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:<u></u><u></u></p>
<blockquote style="border:none;border-left:solid #cccccc 1.0pt;padding:0cm 0cm 0cm 6.0pt;margin-left:4.8pt;margin-right:0cm">
<div>
<div>
<div>
<div>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing">Dear NCL Users,<u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing"> <u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing">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_6018092740948470009m4923127676996268827m7531163896496914688msonospacing"> <u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing"> season = (/ "DJF","JJA" /)<u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing"> nSeason = dimsizes(season)<u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing"> <u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing"> table = new ( (/nCase,nSeason,nSource/), typeof(ratio) )<u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing"> table(0,:,:) = (/CA_bias, CA_biasn/)<u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing"> table(1,:,:) = (/CB_bias, CB_biasn/)<u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing"> table(2,:,:) = (/CC_bias, CC_biasn/)<u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing"> table(3,:,:) = (/CD_bias, CD_biasn/)<u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing"> <u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing"> <u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing"> tt_opt = True<u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing"> tt_opt@tableTitle = "Bias(%)"<u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing"> tt_opt@pltType= "png" ; "eps" [default], "pdf", "ps"<u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing"> <wbr> ; "png", "gif" [if you have ImageMajik 'convert']<u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing">; tt_opt@color0 = "palegreen2"<u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing">; tt_opt@color1 = "tomato2"<u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing"> tt_opt@color0 = "white"<u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing"> tt_opt@color1 = "white"<u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing"> <u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing"> <u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing"> taylor_metrics_table("taylor_<wbr>bias_dry", source, case ,season, table, tt_opt)<u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing"> <u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing">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_6018092740948470009m4923127676996268827m7531163896496914688msonospacing"> <u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing">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_6018092740948470009m4923127676996268827m7531163896496914688msonospacing"> <u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing">It converted the _FillValue to float but the missing_value are still double:<u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing">netcdf CM5A-MR_RCA4_dry_power_sm2 {<u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing">dimensions:<u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing"> time = UNLIMITED ; // (4 currently)<u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing"> latitude = 201 ;<u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing"> longitude = 464 ;<u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing">variables:<u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing"> float pr(time, latitude, longitude) ;<u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing"> pr:_FillValue = -999.f ;<u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing"> pr:cell_methods = "time: mean" ;<u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing"> pr:long_name = "precipitation" ;<u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing"> pr:missing_value = -999. ;<u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing"> pr:standard_name = "precipitation_flux" ;<u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing"> pr:units = "mm" ;<u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing"> float longitude(longitude) ;<u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing"> longitude:standard_name = "longitude" ;<u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing"> longitude:long_name = "longitude" ;<u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing"> longitude:units = "degrees_east" ;<u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing"> longitude:axis = "X" ;<u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing"> float latitude(latitude) ;<u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing"> latitude:standard_name = "latitude" ;<u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing"> latitude:long_name = "latitude" ;<u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing"> latitude:units = "degrees_north" ;<u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing"> latitude:axis = "Y" ;<u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing"> double time(time) ;<u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing"> time:standard_name = "time" ;<u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing"> time:long_name = "time" ;<u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing"> time:units = "days since 1949-12-1 00:00:00" ;<u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing"> time:calendar = "standard" ;<u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing"> time:axis = "T" ;<u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing"> <u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing">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_6018092740948470009m4923127676996268827m7531163896496914688msonospacing"> <u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing">Variable: CD_biasn<u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing">Type: float<u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing">Total Size: 20 bytes<u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing"> 5 values<u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing">Number of Dimensions: 1<u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing">Dimensions and sizes: [5]<u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing">Coordinates:<u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing">Number Of Attributes: 1<u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing"> _FillValue : -999<u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing">(0) -999<u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing">(1) -4.541727<u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing">(2) 12.45617<u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing">(3) 10.47692<u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing">(4) -999<u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing"> <u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing">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_6018092740948470009m4923127676996268827m7531163896496914688msonospacing"> <u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing">p8@missing_value=<wbr>doubletofloat(p8@missing_value<wbr>)<u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing"> p8@_FillValue = p8@missing_value<u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing"> delete(p8@missing_value)<u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing"> <u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing">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_6018092740948470009m4923127676996268827m7531163896496914688msonospacing"> <u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing"> p8_out = addfile("<a href="http://p8.nc" target="_blank">p8.nc</a>","c")<u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing"> p8_out ->p8 = p8<u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing"> <u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing"> <u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing">ncdump -h <a href="http://p8.nc" target="_blank">
p8.nc</a><u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing">netcdf p8 {<u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing">dimensions:<u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing"> time = 4 ;<u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing"> latitude = 201 ;<u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing"> longitude = 464 ;<u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing">variables:<u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing"> float p8(time, latitude, longitude) ;<u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing"> p8:_FillValue = -999.f ;<u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing"> p8:cell_methods = "time: mean" ;<u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing"> p8:long_name = "precipitation" ;<u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing"> p8:missing_value = -999. ;<u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing"> p8:standard_name = "precipitation_flux" ;<u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing"> p8:units = "mm" ;<u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing"> double time(time) ;<u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing"> time:standard_name = "time" ;<u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing"> time:long_name = "time" ;<u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing"> time:units = "days since 1949-12-1 00:00:00" ;<u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing"> time:calendar = "standard" ;<u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing"> time:axis = "T" ;<u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing"> float latitude(latitude) ;<u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing"> latitude:standard_name = "latitude" ;<u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing"> latitude:long_name = "latitude" ;<u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing"> latitude:units = "degrees_north" ;<u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing"> latitude:axis = "Y" ;<u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing"> float longitude(longitude) ;<u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing"> longitude:standard_name = "longitude" ;<u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing"> longitude:long_name = "longitude" ;<u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing"> longitude:units = "degrees_east" ;<u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing"> longitude:axis = "X" ;<u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing">}<u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing"> <u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing"> <u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing">Can someone suggest me a solution?<u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing">Thank you for your help in advance!<u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing">Kind regards,<u></u><u></u></p>
<p class="m_6018092740948470009m4923127676996268827m7531163896496914688msonospacing">Beata<u></u><u></u></p>
<p class="MsoNormal"> <u></u><u></u></p>
</div>
</div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
</div>
<p class="MsoNormal" style="margin-bottom:12.0pt">______________________________<wbr>_________________<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/<wbr>mailman/listinfo/ncl-talk</a><u></u><u></u></p>
</blockquote>
</div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<p class="MsoNormal" style="margin-bottom:12.0pt"><br>
______________________________<wbr>_________________<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/<wbr>mailman/listinfo/ncl-talk</a><u></u><u></u></p>
</blockquote>
</div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
</div></div></div>
</div>
</blockquote></div><br></div>