[Science]   [BAS home]   [Met home]	Antarctic Meteorology

Verification of the coupled model (freedrift version)

The runs here are from HADcm3 (aaxzc, for fans of exactness) which is ocean-drift and yaben, a coupled free-drift run done at BAS.

See yabeo for a verification of the ocean-drift code.

Note that when comparing 32- and 64- bit below, the 32-bit run has about twice as many years.

32-bit against HADcm3

Sea ice	JJA	DJF
SH
NH

64-bit against HADcm3

Sea ice	JJA	DJF
SH
NH

Comparison of time series of free drift and ocean drift

IDL stuff to generate the above...

en64_12=gm('yaben/64/0.03',y='all',m=12,var='icec')                
en64=gm('yaben/64/0.03',y='all',m=6,var='icec')    
en32_12=gm('yaben/32/0.03',y='all',m=12,var='icec')
en32=gm('yaben/32/0.03',y='all',m=6,var='icec')    

ctl3=gm('ctl3/x.03m',m=6,var='icec')
ctl3_12=gm('ctl3/x.03m',m=12,var='icec')
ctl3_series_12=gm('ctl3/0.03',y=2100+indgen(19),m=12,var='icec')                                   
ctl3_series_6=gm('ctl3/0.03',y=2100+indgen(19),m=6,var='icec')  

pp_plot,pp_diff(pp_avg(en32),ctl3,ctl=3),xm=[0,0]
gettwogifs,out='ice-jja-sh'
pp_plot,pp_diff(pp_avg(en32_12),ctl3_12,ctl=3),xm=[0,0]
gettwogifs,out='ice-djf-sh'
pp_plot,pp_diff(pp_avg(en64),ctl3,ctl=3),xm=[0,0]
gettwogifs,out='ice-jja-sh-64'
pp_plot,pp_diff(pp_avg(en64_12),ctl3_12,ctl=3),xm=[0,0]
gettwogifs,out='ice-jja-sh-64'
pp_plot,pp_diff(pp_avg(en32),ctl3,ctl=3),xm=[0,0],/nh
gettwogifs,out='ice-jja-nh'
pp_plot,pp_diff(pp_avg(en32_12),ctl3_12,ctl=3),xm=[0,0],/nh
gettwogifs,out='ice-djf-nh'
pp_plot,pp_diff(pp_avg(en64),ctl3,ctl=3),xm=[0,0],/nh
gettwogifs,out='ice-jja-nh-64'
pp_plot,pp_diff(pp_avg(en64_12),ctl3_12,ctl=3),xm=[0,0],/nh
gettwogifs,out='ice-jja-nh-64'

@time_series
gettwogifs,out='time_series'