Re: How to calculate the area integrated evaporation in Atlantic

Well assuming there is a regional maskfor the region of interest

rmask[*][*] ; global mask that contains flags
                 ; to indicate different regions
wgt[*][*] ; global weights

evap(time,lats,lons)

      evap_area = new ( ntim, typeof(evap), getFillValue(evap) )
      do nt=0,ntim-1
         evap_area(nt) = sum( where(rmask.eq.?? , wgt*evap(nt,:,:), 0.) )
      end do

Wei Liu wrote:
> Hi, there
>
> I want to calculate the area integrated evaporation in Atlantic.
> The data EVAP_F is on POP grid which is not regular.
> So how to calculate the area sum of EVAP using
> two-dimensional weights? Thanks.
>
>
> V
>
>
> ------------------------------------------------------------------------
>
> _______________________________________________
> 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 Thu Jan 14 14:42:18 2010