Re: missing values from grib variable

From: Dennis Shea <shea_at_nyahnyahspammersnyahnyah>
Date: Mon, 04 May 2009 17:26:01 -0600

re [2]
http://www.ncl.ucar.edu/Document/Manuals/Ref_Manual/NclFormatSupport.shtml

'ncl_scalar' is a reserved dimension name

To write a pure scalar

    fnc = addfile ("foo.nc", "r")

--
    r   = 6.37122e06
    r!0 = "ncl_scalar"
    fnc->R = r
or, better
    re = 6.37122e06
    re_at_long_name = "radius of earth"
    re_at_units = "m"
    filevardef(fnc, "REARTH", typeof(re), "ncl_scalar")
    filevarattdef(fnc, "REARTH", re)
    fnc->REARTH = (/ re /)           ; value only
Good luck
Wolfgang Langhans wrote:
> Hi,
> I have two questions:
> 1) I read in a variable from a grib file:
> LAI_MX=f->LAI_MX_GDS10_SFC(:,:)
> where f points at the grib file. The command printVarSummary(LAI_MX)
> shows the attribute _FillValue :  -999. There should actually be plenty
> of missing value entries in this variable, but these have the value zero
> instead, which is not good. How does it come that wrong missing values
> are assigned here?
> 
> 2) Following option 2 for netcdf output on the webpage I export to
> netcdf format. One variable is simply a dimensionless character and I
> run into troubles, if I use getvardims:
> 
> filevardef(fout,
> "rotated_pole",typeof(rotated_pole),getvardims(rotated_pole))
>  
> How should I proceed here, when no dimensions need to be assigned to
> "rotated_pole"?
> 
> Many thanks,
> Wolfgang
> 
> ************************************************
> Wolfgang Langhans
> 
> Institute for Atmospheric and Climate Science
> ETH Zürich
> Universitätstrasse 16
> CH-8092 Zürich, Switzerland
> 
> e-mail: wolfgang.langhans_at_env.ethz.ch
> phone : +41 44 632 78 18
> fax   : +41 44 632 13 11
> room  : CHN L16.3
> web   : http://www.iac.ethz.ch/people/wolangha
> ************************************************
> 
> 
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Mon May 04 2009 - 17:26:01 MDT

This archive was generated by hypermail 2.2.0 : Wed May 06 2009 - 16:48:59 MDT