Re: Unknown and unreasonable values when making

I just did the following and it worked fine ... including the plot

   diri = "./" ; input directory
   f = addfile (diri+"air.199001.nc", "r")

   xShort = f->air(10,:,:,:) ; (time=11,lev,lat,lon)
   if (isatt(xShort,"_FillValue") .and. isatt(xShort,"missing_value")) then
     if (xShort@_FillValue .ne. xShort_at_missing_value) then
       xShort@_FillValue = xShort_at_missing_value ; use alt. missing value
     end if ; only if needed
   end if

   x = short2flt(xShort) ; convert short=>float

   printVarSummary(x)

===
I speculate the problem with specific humidity, is related to
something you are doing. Use printVarSummary and printMinMax
to examine the variabe

http://www.ncl.ucar.edu/Document/Functions/Built-in/printVarSummary.shtml
http://www.ncl.ucar.edu/Document/Functions/Contributed/printMinMax.shtml

Yan Bao wrote:
> Hi,
> I have tried the way Dave had shown me , but the plots I made are
> same as before. One intresting thing is when I tried to use NCEP data to
> make the same figures, it works, but for specific humility, if you
> change the unit (from g/kg to kg/kg), you will get the totally different
> plots. I know that's the same problem both in NARR and NCEP, but how can
> I advoid the problem? Any ideas for that?
> Thanks,
>
> Yan
>
>
> ------------------------------------------------------------------------
>
> _______________________________________________
> 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 Wed Jul 01 2009 - 09:41:14 MDT