



REstar = [600 -600];

nR = length(REstar);

Bo = [4000];
nB = length(Bo);

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

for k=1:nday,
  switch prod_day(k),
    case 9,
      midstr = 'shoal_RE'
    case 18,
      midstr = 'RE';
    otherwise
      disp('*** error bad prod day')
      return
  end


  for i=1:nR,
      for j = 1:nB,
          prefixw = sprintf('cd /home/kghanson/basilisk/shoaling/runs/test_2022_%02d_%02d/%s%d_BO%d',prod_mon,prod_day(k),midstr,REstar(i),Bo(j));
          disp(prefixw)
          eval(prefixw)
          DD = dir('shoal_binary');
          if (isempty(DD))
              system('mkdir shoal_binary')
          end;
          %% Arguments: output_matrix  dump_num_start dump_num_end  Lx  NX  ymin ymax NY  scalar\n");          
          system('~/SurfRanch/Basilisk/basilisk_code/output_matrix 360 590 8 1600 -0.2 1.8 460 p focus') % delta = 0.005
          system('~/SurfRanch/Basilisk/basilisk_code/output_matrix 360 590 8 1600 -0.2 1.8 460 f focus')
          system('~/SurfRanch/Basilisk/basilisk_code/output_matrix 360 590 8 1600 -0.2 1.8 460 z focus')          
          %system('output_matrix 20 20 16 800 -0.37109 3 168 a broad') % delta = 0.02, t=1002
          %system('output_matrix 80 80 16 800 -0.37109 3 168 a broad') % delta = 0.02, t=1005
          %system('output_matrix 140 140 16 800 -0.37109 3 168 a broad') % delta = 0.02, t=1008                    
          %system('output_matrix 200 200 16 800 -0.37109 3 168 a broad') % delta = 0.02, t=1011          
          %system('output_matrix 280 280 16 800 -0.37109 3 168 a broad') % delta = 0.02, t=1015
          %system('output_matrix 360 360 16 800 -0.37109 3 168 a broad') % delta = 0.02, t=1019          
          %system('output_matrix 460 460 16 800 -0.37109 3 168 a broad') % delta = 0.02, t=1020
          %system('output_matrix 480 480 16 800 -0.37109 3 168 a broad') % delta = 0.02, t=1020.2          
          %system('output_matrix 500 500 16 800 -0.37109 3 168 a broad') % delta = 0.02, t=1020.4
          %system('output_matrix 520 520 16 800 -0.37109 3 168 a broad') % delta = 0.02, t=1020.6
          %system('output_matrix 540 540 16 800 -0.37109 3 168 a broad') % delta = 0.02, t=1020.8                
          %system('output_matrix 560 560 16 800 -0.37109 3 168 a broad') % delta = 0.02, t=1021
          %system('output_matrix 580 580 16 800 -0.37109 3 168 a broad') % delta = 0.02, t=1021.2
          %system('output_matrix 600 600 16 800 -0.37109 3 168 a broad') % delta = 0.02, t=1021.4                    
      end;
  end;

end;


cd ~/SurfRanch/Basilisk
