# Re: how to calculate the sea ice area with three-dimensional data

From: Dennis Shea <shea_at_nyahnyahspammersnyahnyah>
Date: Mon May 06 2013 - 09:07:52 MDT

Do the files you have available have an "area" variable?

eg:

float tarea(nlat, mlon) ;
tarea:units = "m^2" ;

If your variable is (say) "ice", then

ice = f->ice ; (time,nlat,mlon)
tarea = f->tarea ' (nlat,mlon)

TAREA = conform(ice, tarea, (/1,2/)) ; TAREA(ntim,nlat,mlon)
printVarSummary(TAREA)

iarea = where(ice.gt.0, TAREA, ice@_FillValue)
printVarSummary(iarea)

ice_area = dim_sum( iarea, (/1,2/) )
ice_area!0 = "time"
ice_area&time = ice&time
ice_area@long_name = "Total ice area"
ice_area@units = tarea@units
printVarSummary(ice_area)

On 5/6/13 6:21 AM, ì§„í˜„ê·¼ wrote:
> To whom it may concern
> I am trying to calculate total NH sea ice area with AMIP data are
> three-dimensional(time, lat, lon),
> but I don't know how can I match the array with "wgt_areasum2" function.
>
> I found some examples on the internet(but only used 4-dimensional data
> examples)
> but failed to calculate sea ice area
> because I don't have weight variables or weight files
>
>
>
>
> _______________________________________________
> 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 Mon May 6 09:07:28 2013

This archive was generated by hypermail 2.1.8 : Tue May 07 2013 - 08:38:33 MDT