



%REstar = [2400 1800 1200 600 0 -600 -1200 -1800];

REstar = [2400  -1800];
nR = length(REstar);

%Bo = [200 1000]; %4000 8000];
Bo = [4001];
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_all 50 590 2 5 50 4500 -1 5 600 a') % 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')          
      
      end;
  end;

end;


cd ~/SurfRanch/Basilisk
