Meeting notes, 2004-11-30

The modified ancil field

!p.multi=[0,2,1]
a=get_field('waterfix.ancil.32.pp',/nat)
b=get_field('waterfix.ancil.new.32.pp',/nat)

pp_plot,a,/blo,la=-45,ymarg=[2,0],xmarg=[0,0]    
pp_plot,b,/blo,la=-45,ymarg=[2,0],xmarg=[0,0]
gettwogifs,out='p1-new'

Sea ice conc and sfc T changes

!p.multi=[0,2,1]
pp_plot,pp_diff(gm('yacjf/64/x.03m',sc=146)),xmarg=[0,0]
pp_plot,pp_diff(gm('yacjf/64/x.03m',sc=24)),xmarg=[0,0]
gettwogifs,out='icec-and-t'

T, globally

!p.multi=[0,1,1]
pp_plot,pp_diff(gm('yacjf/64/x.03m',sc=24)),/glob
gettwogifs,out='t-global'

Depth-lat section at 180 E and 0 E

c=gm('yacjc/64/0.03',sc=101,nl=15,y='all',m=6)
c1=pp_avg1(reform(c,15,n_elements(c)/15),2)
f=gm('yacjf/64/0.03',sc=101,nl=15,y='all',m=6)
f1=pp_avg1(reform(f,15,n_elements(f)/15),2)

!p.multi=[0,1,2]
pp_plot,pp_xsect(lo=180,pp_diff(f1,c1)),lf=0.25
pp_plot,pp_xsect(lo=0,pp_diff(f1,c1)),lf=0.25
gettwogifs,out='t-sect-0-180'

Time-lat section (along 180E) at 447m depth

c=reform(c,15,n_elements(f)/15)
f=reform(f,15,n_elements(f)/15)

f2=xsect(f(10,*),lo=180,/y)
c2=xsect(c(10,*),lo=180,/y)
pp_plot,pp_diff(f2,c2),/no_extr
pp_plot,pp_diff(f2,c2),/no_extr,lf=0.25
gettwogifs,out='t-lat-180'

60day_xblfn

Ice conc (days 10, 30 , 60 minus day 1)

d=nc2pp('/users/icd/amjr/local/pp_fields/ocean/60day_xblfn_89b10_aice',lon='longitude',lat='latitude',data='iceconc')
!p.multi=[0,3,1]
pp_plot,pp_diff(d(1:3),d(0)),xmarg=[0,0]
gettwogifs,out='icec-5-10-30-60'

!p.multi=[0,3,1]
pp_plot,d(0:2),xmarg=[0,0]
gettwogifs,out='actual-ice-1-10-30'

Ice depth

d=nc2pp('/users/icd/amjr/local/pp_fields/ocean/60day_xblfn_89b10_hice',lon='longitude',lat='latitude',data='icedepth')
pp_plot,pp_diff(d(1:3),d(0)),xmarg=[0,0],lev=(indgen(11)-5)/10.,/noz
gettwogifs,out='iced-5-10-30-60'