weighted area average of arctic

From: Mira Berdahl <mberdahl_at_nyahnyahspammersnyahnyah>
Date: Fri Jun 07 2013 - 17:14:03 MDT

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
Received on Fri Jun 7 17:14:11 2013

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