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

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
>
> Thank you for your help in advance.
>
>
>
> _______________________________________________
> 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