



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

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

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

end;


