Re: ncl script problem

total@_FillValue=1e+35
total=mask(total,total.le.0.,False)
fraction@_FillValue=1e+35
fraction=emission/total

On Jul 11, 2006, at 5:04 PM, Chao Luo wrote:

> Hi,
>
> I am trying to plot fraction of each dust bin emission to total
> emission. Since in many regions the total emissions are zero. So I
> need
> to set Fill value for zero numbers in these region in case of divided
> by zero. But the error message: Subscript out of range. Is there any
> better ways for this?
>
> Thanks,
>
> Chao
>
> part of my script is like:
>
> diri = "/data3/chaoluo/camdst/cam3019aersombf/"
> fl_mdl_3 = addfile(diri+"cam3019aersombf_annual.nc","r")
> sfmbl = fl_mdl_3->DSTSFMBL(0,:,:)
> sfmblx01 = fl_mdl_3->DSTX01SF(0,:,:)
>
> sfmbl@_FillValue=1.e+36
> nlon=128
> nlat=64
>
> frc1 = new((/nlat,nlon/),float)
> do ilon=0,nlon-1
> do ilat=0,nlat-1
> if (sfmbl(ilon,ilat).ne.0.) then
> frc1(ilon,ilat)=sfmblx01(ilon,ilat)/sfmbl(ilon,ilat)
> else
> frc1(ilon,ilat)=sfmbl@_FillValue
> end if
> end do
> end do
> print (frc1)
>
>
> _______________________________________________
> 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:11:18 MDT