Re: ncl script problem

From: Michael Notaro <mnotaro_at_nyahnyahspammersnyahnyah>
Date: Tue, 11 Jul 2006 17:49:28 -0500

When I sent the email, I hadn't realized that the NCL script
was at the bottom of Chao's email. I am in a proposal-writing fog.

Mike

On Jul 11, 2006, at 5:38 PM, Dennis Shea wrote:

>>
>> total@_FillValue=1e+35
>> total = mask(total,total.le.0.,False)
>> fraction@_FillValue=1e+35
>> fraction = emission/total
> ________________________________
>
> Mike's solution is a good way to do it.
>
> Minor modification to get meta data
>
> frac = sfmbl ; make array of same size with meta data
> frac@_FillValue=1e+20 ; **only *if sbmbl does not already
> ; have a _FillValue
>
> frac = mask(frac,frac.le.0.,False)
>
> frac = sfmblx01/frac
>
> frac_at_long_name = "ratio: sfmblx01/sfmbl"
> frac_at_units = ""
>
> printVarSummary(frac)
> printMinMax(frac, True ) ; contributed.ncl
>
> nmsg = num( ismissing( frac ) )
> print("frac: nmsg="+nmsg)
>
> _______________________________________________
> ncl-talk mailing list
> ncl-talk_at_ucar.edu
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk

_______________________________________________
ncl-talk mailing list
ncl-talk_at_ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Tue Jul 11 2006 - 16:49:28 MDT

This archive was generated by hypermail 2.2.0 : Thu Jul 13 2006 - 09:59:49 MDT