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

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

I attached the wrong function. Sorry

On 5/6/13 11:30 AM, Dennis Shea wrote:
> 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

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk

Received on Mon May 6 11:59:21 2013

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