Re: weighted area average of arctic

Not sure that ncl-talk is the place for this sort of question.

Is the grid rectilinear or curvilinear? Did you get any
errpr messages?

What do the folowing yield?

printVarSummary(tas_bnu_rcp45)
printVarSummary(cellarea_bnu)

On 6/7/13 5:14 PM, Mira Berdahl wrote:
> Hi everyone,
>
> I was just wondering if someone could confirm for me that I am properly
> taking the area average of the Arctic in the below script. Thank you very
> much in advance!
>
> Cheers,
> MB
>
>
> ;***************
> ; Start script
> ;***************
>
> begin
>
>
> in = addfile("tas_Amon_BNU-ESM_rcp45_r1i1p1_200601-210012.nc","r")
> tas_bnu_rcp45 = in->tas
>
> ;%%%%% read in the weights (area of grid cells %%%%%%%%%%%%%%%%
> in = addfile("areacella_fx_BNU-ESM_G1_r0i0p0_new.nc","r"); area of each
> grid cell
> cellarea_bnu = in->areacella
>
> wgty = cellarea_bnu(:,1)
> printVarSummary(wgty)
>
> ;%%% calculate the area weighted Arctic temperature --- based on area of
> grid cells.
> Arctic_tas_bnu_rcp45 =
> wgt_areaave(tas_bnu_rcp45(:,{60:90},:),wgty({60:90}),1.0,0)
> ;print("Arctic Average TAS BNU="+Arctic_tas_bnu)
>
> ; write out the Arctic average data for bcc into an ascii file
> asciiwrite("Arctic.tas.rcp45.BNU.asc",Arctic_tas_bnu_rcp45)
>
> end
> _______________________________________________
> 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 Fri Jun 7 17:35:14 2013