



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

Bo = [1000];
%Bo = [4000 8000];

nB = length(Bo);

prod_mon = 11;
prod_day = [9 18];
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,
        prefix = sprintf('/home/kghanson/basilisk/shoaling/runs/test_2022_%02d_%02d/%s%d_BO%d/shoal_dump',prod_mon,prod_day(k),midstr,REstar(i),Bo(j));

        for dump_num=1:560;
          syscmd = sprintf('cd %s; /home/ffeddersen/SurfRanch/Basilisk/basilisk_code/get_offset_pressure %03d',prefix, dump_num );
          disp(syscmd)
          system(syscmd);
        end;
    end;
  end;
end;



