




REstar = [2400 1800 1200 600 0 -600  -1200 -1800];
%REstar = [2400]; %1800 1200 600 0 -600 -1200 -1800];
nR = length(REstar);

Bo=[4000];
%Bo = [200 1000 2000 4000 8000];
%Bo = [4000 8000];
nB = length(Bo);

prod_mon = 11;
prod_day = [9 18];
prod_day = [9 ];
nday = length(prod_day);

for k=1:nday,

  for i=1:nR,
      for j = 1:nB,
          disp(sprintf('** testing extract of Restar=%d, Bo=%d',REstar(i),Bo(j)))
          R = read_surface_extract(9, REstar(i),Bo(j), 1009.0);
        R = read_surface_extract(9, REstar(i),Bo(j), 1019.0);          
    end;
  end;
end;



