Re: question on how to calculate the area sum with two-dimensional coordinates data

Hi Tingting

Can you please point me toward a dataset so I may look at it?
Ideally, the file would contain an area variable. For example,
the POP models have

         float TAREA(nlat, nlon) ;
                 TAREA:long_name = "area of T cells" ;
                 TAREA:units = "centimeter^2" ;
                 TAREA:coordinates = "TLONG TLAT" ;
                 TAREA:_FillValue = 9.96921e+36f ;
                 TAREA:missing_value = 9.96921e+36f ;

Consider:

     x(time,nlat,nlon) [ dimension #s (0,1,2) ]

A global weighted mean would be

     xAVG = dim_sum_n(x*conform(x,TAREA,(/1,2/)),0)/ sum(TAREA)

Then xAVG(ntim)

On 3/27/13 11:40 AM, tingting@ucar.edu wrote:
> Hello,
> I am doing the MRI-CGCM3 model output validation and need to calculate
> SIE(sea ice extent), the coordinates of MRI output are two-dimensional,
> while the NCL function(wgt_areasum2) seems to deal with data with regular
> coordinates, right?
> So,my question is that how to calculate the area sum? Or, do I need to
> regrid the MRI to regular-grid data?
> Thank you for any coming help in advance.
> Tingting
>
>
> _______________________________________________
> 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 Wed Mar 27 13:12:19 2013