>
>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
Received on Tue Jul 11 2006 - 16:38:48 MDT
This archive was generated by hypermail 2.2.0 : Thu Jul 13 2006 - 09:59:49 MDT