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

From: Dennis Shea <shea_at_nyahnyahspammersnyahnyah>
Date: Mon May 06 2013 - 11:30:06 MDT

You can create your own function.

The attached is a modification of a function in shea_util

On 5/6/13 9:58 AM, 진현근 wrote:

> Unfortunately, I don't have any kind of available variables as "area"
> please, check the attached information about Sea Ice Concentration(SIC)
> Could you send me a mail to me again?
> sincerely yours
>
>
> dimensions:
> time = UNLIMITED ; // (12 currently)
> longitude = 360 ;
> latitude = 180 ;
> variables:
> double time(time) ;
> time:units = "days since 1979-1-1 0" ;
> double longitude(longitude) ;
> longitude:units = "degrees_east" ;
> longitude:long_name = "Longitude" ;
> double latitude(latitude) ;
> latitude:units = "degrees_north" ;
> latitude:long_name = "Latitude" ;
> float sic(time, latitude, longitude) ;
> sic:comment = " Observed Monthly Mean Sea Ice Concentration (%)" ;
> sic:long_name = "Sea-ice concentration" ;
> sic:units = "%" ;
> sic:grid_name = "linear" ;
> sic:grid_type = "uniform" ;
> sic:level_description = "Earth surface" ;
> sic:time_statistic = "average" ;
> sic:missing_value = 1.e+20f ;
>
>
>
> 2013/5/7 Dennis Shea <shea@ucar.edu>
>
>> 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<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 11:29:40 2013

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