Re: ncl script problem

From: Dennis Shea <shea_at_nyahnyahspammersnyahnyah>
Date: Tue, 11 Jul 2006 16:38:48 -0600 (MDT)

>
>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