[Thread Prev][Thread Next][Index]
[ferret_users] Reading NetcDF attributes with higher precision with ferret
Hello Ferret users,
I have a precision problem to plot a particular value of my variable
"chlorophyll_concentration" with ferret :
I read this variable from a NetCDF file, which contains a lat-lon grid
of 84*log(chlorophyll_concentration) stored on a byte.
This is the result of "show data" :
name title I
J K L
CHLOROPHYLL_CONCENTRATION
log(chlorophyll-a in mg/m3) 1:2627 1:2100 ... ...
I have set _FillValue=127 in the NetCDF file
I have also set two additional "out of range" values :126 to describe a
land mask, and 125 to describe the cloud mask
A scale factor = 1/84 is set in the NetCDF file, so ferret displays the
log of chlorophyll.
! reading the scale factor
LET sc = chlorophyll_concentration.scale_factor
list sc outputs at the screen 0.0190 (instead of 0.011904762 which is
the more exact value of 1/84)
shade/levels=(`sc*126`,`sc*126`,1)/title="land mask"
"chlorophyll_concentration" gives a correct result : (maybe because
sc*126 = 1.50000000)
The land mask value appear in pink, with only one pink colour on the
palette, with lower boundary = 1.5 and upper boundary = 1.5
shade/levels=(`sc*125`,`sc*125`,1)/title="cloud mask"
"chlorophyll_concentration" does not work (no cloud values displayed),
maybe because sc*125 = 1.4888095238 ?
The palette on the right is with only one pink colour on the palette,
with lower boundary = 1.49 and upper boundary = 1.49
Of course, I have checked that values = 125 are present in my variable !
Is there a way to improve the precision of the scale factor read with
ferret, or doing with other means ?
Thank you for your help
Jacques
--
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
JACQUES STUM ~~ E-mail : jstum@xxxxxx
Collecte Localisation Satellite,
8-10 rue Hermès, ~~ Phone : 33 5 61 39 48 12
Parc Technologique du Canal,
31526 Ramonville Saint-Agne, France ~~ Fax : 33 5 61 39 37 82
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
[Thread Prev][Thread Next][Index]
Contact Us
Dept of Commerce /
NOAA /
OAR /
PMEL /
TMAP
Privacy Policy | Disclaimer | Accessibility Statement