[ncid, status] = mexnc ( 'open', 'CPB_GRID.nc')
[lat] = mexnc ( 'get_var_double', ncid, 20 );
[lon] = mexnc ( 'get_var_double', ncid, 21 );
[h] = mexnc ( 'get_var_double', ncid, 6 );
h = h(1:81,:);

n = 0;
for i = 1:6:277
    if i<1000, blank = '0'; end
    if i<100,  blank = '00'; end
    if i<10,   blank = '000'; end
    fname = ['para1000',blank,num2str(i),'.csv'];
    data = load(fname);

    surf(lon,lat,-h), shading interp
    hold on
    plot3(data(:,3),data(:,4),data(:,1),'g.','MarkerSize',.8)
    axis([-76.7 -75.8 36.8 39.2 -30 0])
    hold off
    
    set(gca,'CameraPosition',[-76-.005*n 36-.1*n 30+.1*n])
    n = n+1;
    pause(.5)
end