[gradsusr] Help: Extent

[Jennifer Adams]

Please have a look at the atot() function  -- it was designed for this very purpose. Please see the second example in the doc page at http://iges.org/grads/gadoc/gradfuncatot.html. You have to use GrADS v2.0.2. 
--Jennifer


On Apr 23, 2013, at 10:20 AM, Fernanda Casagrande wrote:

> Dear All,
> 
> I have a problem with the calculation of the area (million km2).
> 
> I have a file of sea ice (ice concentration data from 0-100%) and I need to create a time series with a total extension of sea ice in the arctic and antartica (million km2). For this calculation, usually consider only areas with ice concentration equal to or above 15%, so I used the function maskout.
> 
> The. gs is running but the values ​​for the Arctic are wrong, can someone help me? Some other method to calculate the area (million km2)
> 
> these data are from satellite.
> 
> 
> 'sdfopen X172.16.9.158.108.11.23.58.nc'
> *calculate area (million km2)
> 'dx=2*3.1416*6371*cos(lat*3.1416/180)/360.0*cdiff(lon,x)/2 '
> 'dy=2*3.1416*6371/360.0*cdiff(lat,y)/2'
> 'area=(dx*dy)'  
> 'fracao=area*(icec/100)'
> 
> say 'calcula a area'
> 
> 'set lon 0 360'
> 'set lat -90 90'
> 'set t 1 last'
> 
> 'define ones = const(const(maskout(icec,icec-15),1),0,-u)' 
> 'set x 1'
> 'set y 1'
> 
> 'define xx=tloop(asum(ones*area,lon=0,lon=360,lat=-90,lat=0))'
> 'd xx'
> 
> 'ltrend xx a b c'
> 'set ccolor 2'
> 'd a';'d b'
> trend=subwrd(result,4)
>  say 'trend='trend
> ' d c'
> rmse=subwrd(result,4)
>  say 'rmse='rmse
> 
> 'draw title Antarctic Sea Ice obs incl.'trend' rmse. 'rmse' '
> 
> 
> regards
> 
> 
> 
> 
> 
> Fernanda
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr

--
Jennifer M. Adams
IGES/COLA
4041 Powder Mill Road, Suite 302
Calverton, MD 20705
jma at cola.iges.org



-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20130423/2d76d6cf/attachment-0003.html