I have made sensitivity studies with ECHAM5. Now I would like to calculate the NAO-Index. First I have calculated the sea level pressure at two stations.
; surface pressure
aps = in->aps_ave(:,:,:)
tm1 = in->tm1_ave(:,:,:,:)
press = in->press_ave(:,:,:,:)
phis = geofile->GEOSP(:,:)
;**********************************************************************
; claculate sea level pressure at two points
;**********************************************************************
; sea level pressure at Akureyri
slp1 = pslec(tm1(:,30,:,:),phis,aps,press(:,30,:,:))
copy_VarMeta(aps,slp1)
slp1@long_name = "Sea Level Pressure"
printVarSummary(slp1)
slp1 = lonFlip(slp1)
printVarSummary(slp1)
slp1_ave = slp1(:,{65.7},{-18.1})
slp1_seas = month_to_season(slp1_ave,season)
nyrs = dimsizes(slp1_seas&time)
printVarSummary(slp1_seas)
; sea level pressure at Ponta Delgada
slp2 = pslec(tm1(:,30,:,:),phis,aps,press(:,30,:,:))
copy_VarMeta(aps,slp2)
slp2@long_name = "Sea Level Pressure"
printVarSummary(slp2)
slp2 = lonFlip(slp2)
printVarSummary(slp2)
slp2_ave = slp2(:,{65.7},{-18.1})
slp2_seas = month_to_season(slp2_ave,season)
nyrs = dimsizes(slp2_seas&time)
printVarSummary(slp2_seas)
How can I calculate the NAO-Index with ncl?
Thanks.
Lena