runid = 'chan_M1e4';
Mmud  = 1e-4;
for i=0:10;
    filename = [runid,'_',num2str(i),'.sed'];
    file = fopen(filename,'w');
    pm = i/10;
    ps = 1-pm;
    fprintf(file, '%s\n', '[SedimentFileInformation]');
    fprintf(file, '%s\n', '   FileCreatedBy    = Delft3D FLOW-GUI, Version: 3.41.06.10981'); 
    fprintf(file, '%s\n', '   FileCreationDate = Mon May 09 2011, 09:50:10');
    fprintf(file, '%s\n', '   FileVersion      = 02.00');
    fprintf(file, '%s\n', '[SedimentOverall]');
    fprintf(file, '%s\n', '   Cref             = 1600       [kg/m3]  CSoil Reference density for hindered settling calculations');
    fprintf(file, '%s\n', '   IopSus           = 0                             If Iopsus = 1: susp. sediment size depends on local flow and wave condit');
    fprintf(file, '%s\n', '[Sediment]');
    fprintf(file, '%s\n', '   Name             = #Sediment sand#     Name of sediment fraction');
    fprintf(file, '%s\n', '   SedTyp           = sand                Must be "sand", "mud" or "bedload"');
    fprintf(file, '%s\n', '   RhoSol           =  2650      [kg/m3]  Specific density');
    fprintf(file, '%s\n', '   SedDia           =  1.0e-4    [m]      Median sediment diameter (D50)');
    fprintf(file, '%s\n', '   CDryB            =  1600      [kg/m3]  Dry bed density');
    fprintf(file, '%s%4.2f%s\n', '   IniSedThick      =  ',ps','      [m]      Initial sediment layer thickness at bed (uniform value or filename)');
    fprintf(file, '%s\n', '   FacDSS           =  1.0       [-]      FacDss * SedDia = Initial suspended sediment diameter. Range [0.6 - 1.0]');
    fprintf(file, '%s\n', '[Sediment]');
    fprintf(file, '%s\n', '   Name             = #Sediment mud#      Name of sediment fraction');
    fprintf(file, '%s\n', '   SedTyp           = mud                 Must be "sand", "mud" or "bedload"');
    fprintf(file, '%s\n', '   RhoSol           =  2650      [kg/m3]  Specific density');
    fprintf(file, '%s\n', '   SalMax           =  0.0       [ppt]    Salinity for saline settling velocity');
    fprintf(file, '%s\n', '   WS0              =  5.0e-4    [m/s]    Settling velocity fresh water');
    fprintf(file, '%s\n', '   WSM              =  5.0e-4    [m/s]    Settling velocity saline water');
    fprintf(file, '%s\n', '   TcrSed           =  1000      [N/m2]   Critical bed shear stress for sedimentation (uniform value or filename)');
    fprintf(file, '%s\n', '   TcrEro           =  1.0       [N/m2]   Critical bed shear stress for erosion       (uniform value or filename)');
    fprintf(file, '%s%6.1e%s\n', '   EroPar           =  ',Mmud,'  [kg/m2/s] Erosion parameter                           (uniform value or filename)');
    fprintf(file, '%s\n', '   CDryB            =  1600      [kg/m3]  Dry bed density');
    fprintf(file, '%s%4.2f%s\n', '   IniSedThick      =  ',pm','      [m]      Initial sediment layer thickness at bed (uniform value or filename)');
    fprintf(file, '%s\n', '   FacDSS           =  1.0       [-]      FacDss * SedDia = Initial suspended sediment diameter. Range [0.6 - 1.0]');
   
   fclose(file);
end
