Re: weighted area average of arctic

From: Dennis Shea <shea_at_nyahnyahspammersnyahnyah>
Date: Mon Jun 10 2013 - 15:14:26 MDT

Speculation: the answer is "no" this is not correct.
Likely you have a curvilinear grid. You are using {...} which
applies to rectilinear grids

   cellarea_bnu = in->areacella
   wgty = cellarea_bnu(:,1)

The areacella variable is rank 2 .. characterisic of curvilinear grid.

Howevr, you should include

cellarea_bnu = in->areacella
printVarSummary(cellarea_bnu)

tas_bnu_rcp45 = in->tas

printVarSummary(tas_bnu_rcp45)

On 06/07/2013 05: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 Mon Jun 10 15:14:34 2013

This archive was generated by hypermail 2.1.8 : Tue Jun 11 2013 - 12:03:57 MDT