runid = 'chan_M1e4';
Mmud  = 1e-4;
pmcrit = 0.3;
for i=0:10;
    if pmcrit>0
        filename = [runid,'_pm03_',num2str(i),'.mdf'];
    else
        filename = [runid,'_',num2str(i),'.mdf'];
    end
    sedname = [runid,'_',num2str(i),'.sed'];
    file = fopen(filename,'w');
    pm = i/10;
    ps = 1-pm;    
    fprintf(file, '%s\n', 'Ident = #Delft3D-FLOW  .03.02 3.41.06.10981#');
    fprintf(file, '%s\n', 'Commnt=                              ');
    fprintf(file, '%s\n', 'Runtxt= #Testing sand-mud interaction  #');
    fprintf(file, '%s\n', 'Filcco= #chan.grd#');
    fprintf(file, '%s\n', 'Fmtcco= #FR#');
    fprintf(file, '%s\n', 'Anglat=  0.0');
    fprintf(file, '%s\n', 'Grdang=  0.0');
    fprintf(file, '%s\n', 'Filgrd= #chan.enc#');
    fprintf(file, '%s\n', 'Fmtgrd= #FR#');
    fprintf(file, '%s\n', 'MNKmax= 152 12 5');
    fprintf(file, '%s\n', 'Thick =  25');
    fprintf(file, '%s\n', '         25');
    fprintf(file, '%s\n', '         25');
    fprintf(file, '%s\n', '         15');
    fprintf(file, '%s\n', '         10');
    fprintf(file, '%s\n', 'Commnt=                              ');
    fprintf(file, '%s\n', 'Fildep= #chan.dep#');
    fprintf(file, '%s\n', 'Fmtdep= #FR#');
    fprintf(file, '%s\n', 'Commnt=                              ');
    fprintf(file, '%s\n', 'Commnt=                              no. dry points: 0');
    fprintf(file, '%s\n', 'Commnt=                              no. thin dams: 0');
    fprintf(file, '%s\n', 'Commnt=                              ');
    fprintf(file, '%s\n', 'Itdate= #2011-05-09#');
    fprintf(file, '%s\n', 'Tunit = #M#');
    fprintf(file, '%s\n', 'Tstart= 0');
    fprintf(file, '%s\n', 'Tstop = 50');
    fprintf(file, '%s\n', 'Dt    = 0.5');
    fprintf(file, '%s\n', 'Tzone = 0');
    fprintf(file, '%s\n', 'Commnt=                              ');
    fprintf(file, '%s\n', 'Sub1  = #    #');
    fprintf(file, '%s\n', 'Sub2  = # C #');
    fprintf(file, '%s\n', 'Namc1 = #Sediment sand       #');
    fprintf(file, '%s\n', 'Namc2 = #Sediment mud        #');
    fprintf(file, '%s\n', 'Commnt=                              ');
    fprintf(file, '%s\n', 'Wnsvwp= #N#');
    fprintf(file, '%s\n', 'Wndint= #Y#');
    fprintf(file, '%s\n', 'Commnt=                              ');
    fprintf(file, '%s\n', 'Zeta0 =  0.0');
    fprintf(file, '%s\n', 'U0    = [.]');
    fprintf(file, '%s\n', 'V0    = [.]');
    fprintf(file, '%s\n', 'S0    = [.]');
    fprintf(file, '%s\n', 'C01   =  0.0');
    fprintf(file, '%s\n', '         0.0');
    fprintf(file, '%s\n', '         0.0');
    fprintf(file, '%s\n', '         0.0');
    fprintf(file, '%s\n', '         0.0');
    fprintf(file, '%s\n', 'C02   =  0.0');
    fprintf(file, '%s\n', '         0.0');
    fprintf(file, '%s\n', '         0.0');
    fprintf(file, '%s\n', '         0.0');
    fprintf(file, '%s\n', '         0.0');
    fprintf(file, '%s\n', 'Commnt=                              ');
    fprintf(file, '%s\n', 'Commnt=                              no. open boundaries: 2');
    fprintf(file, '%s\n', 'Filbnd= #chan.bnd#');
    fprintf(file, '%s\n', 'Fmtbnd= #FR#');
    fprintf(file, '%s\n', 'FilbcT= #chan.bct#');
    fprintf(file, '%s\n', 'FmtbcT= #FR#');
    fprintf(file, '%s\n', 'FilbcC= #chan.bcc#');
    fprintf(file, '%s\n', 'FmtbcC= #FR#');
    fprintf(file, '%s\n', 'Rettis=  0.0');
    fprintf(file, '%s\n', '         0.0');
    fprintf(file, '%s\n', 'Rettib=  0.0');
    fprintf(file, '%s\n', '         0.0');
    fprintf(file, '%s\n', 'Commnt=                              ');
    fprintf(file, '%s\n', 'Ag    =  9.81');
    fprintf(file, '%s\n', 'Rhow  =  1000');
    fprintf(file, '%s\n', 'Alph0 = [.]');
    fprintf(file, '%s\n', 'Tempw =  15');
    fprintf(file, '%s\n', 'Salw  =  31');
    fprintf(file, '%s\n', 'Rouwav= #    #');
    fprintf(file, '%s\n', 'Wstres=  6.3000000e-004  0.0000000e+000  7.2300000e-003  1.0000000e+002');
    fprintf(file, '%s\n', 'Rhoa  =  1.0000000e+000');
    fprintf(file, '%s\n', 'Betac =  5.0000000e-001');
    fprintf(file, '%s\n', 'Equili= #N#');
    fprintf(file, '%s\n', 'Tkemod= #            #');
    fprintf(file, '%s\n', 'Ktemp = 0');
    fprintf(file, '%s\n', 'Fclou =  0.0');
    fprintf(file, '%s\n', 'Sarea =  0.0');
    fprintf(file, '%s\n', 'Temint= #Y#');
    fprintf(file, '%s\n', 'Commnt=                  ');            
    fprintf(file, '%s\n', 'Roumet= #C#');
    fprintf(file, '%s\n', 'Ccofu =  50');
    fprintf(file, '%s\n', 'Ccofv =  50');
    fprintf(file, '%s\n', 'Xlo   =  0.0');
    fprintf(file, '%s\n', 'Vicouv=  1.0');
    fprintf(file, '%s\n', 'Dicouv=  10.0');
    fprintf(file, '%s\n', 'Htur2d= #N#');
    fprintf(file, '%s\n', 'Vicoww=  0.0');
    fprintf(file, '%s\n', 'Dicoww=  0.0');
    fprintf(file, '%s\n', 'Irov  = 0');
    fprintf(file, '%s%s%s\n', 'Filsed= #',sedname,'#');
    fprintf(file, '%s\n', 'Fmtsed= #FR#');
    fprintf(file, '%s\n', 'Filmor= #chan.mor#');
    fprintf(file, '%s\n', 'Fmtmor= #FR#');
    fprintf(file, '%s\n', 'Commnt=                              ');
    fprintf(file, '%s\n', 'Iter  =      2');
    fprintf(file, '%s\n', 'Dryflp= #YES#');
    fprintf(file, '%s\n', 'Dpsopt= #MAX#');
    fprintf(file, '%s\n', 'Dpuopt= #MOR#');
    fprintf(file, '%s\n', 'Dryflc=  1.0000000e-001');
    fprintf(file, '%s\n', 'Dco   = -9.9900000e+002');
    fprintf(file, '%s\n', 'Tlfsmo=  60');
    fprintf(file, '%s\n', 'ThetQH=  0.0');
    fprintf(file, '%s\n', 'Forfuv= #Y#');
    fprintf(file, '%s\n', 'Forfww= #N#');
    fprintf(file, '%s\n', 'Sigcor= #N#');
    fprintf(file, '%s\n', 'Trasol= #Cyclic-method#');
    fprintf(file, '%s\n', 'Momsol= #Cyclic#');
    fprintf(file, '%s\n', 'Commnt=                              ');
    fprintf(file, '%s\n', 'Commnt=                              no. discharges: 0');
    fprintf(file, '%s\n', 'Commnt=                              no. observation points: 14');
    fprintf(file, '%s\n', 'Filsta= #chan.obs#');
    fprintf(file, '%s\n', 'Fmtsta= #FR#');
    fprintf(file, '%s\n', 'Commnt=                              no. drogues: 0');
    fprintf(file, '%s\n', 'Commnt=                              ');
    fprintf(file, '%s\n', 'Commnt=                              ');
    fprintf(file, '%s\n', 'Commnt=                              no. cross sections: 0');
    fprintf(file, '%s\n', 'Commnt=                              ');
    fprintf(file, '%s\n', 'SMhydr= #YYYYY#     ');
    fprintf(file, '%s\n', 'SMderv= #YYYYYY#    ');
    fprintf(file, '%s\n', 'SMproc= #YYYYYYYYYY#');
    fprintf(file, '%s\n', 'PMhydr= #YYYYYY#    ');
    fprintf(file, '%s\n', 'PMderv= #YYY#       ');
    fprintf(file, '%s\n', 'PMproc= #YYYYYYYYYY#');
    fprintf(file, '%s\n', 'SHhydr= #YYYY#      ');
    fprintf(file, '%s\n', 'SHderv= #YYYYY#     ');
    fprintf(file, '%s\n', 'SHproc= #YYYYYYYYYY#');
    fprintf(file, '%s\n', 'SHflux= #YYYY#      ');
    fprintf(file, '%s\n', 'PHhydr= #YYYYYY#    ');
    fprintf(file, '%s\n', 'PHderv= #YYY#       ');
    fprintf(file, '%s\n', 'PHproc= #YYYYYYYYYY#');
    fprintf(file, '%s\n', 'PHflux= #YYYY#      ');
    fprintf(file, '%s\n', 'Online= #N#');
    fprintf(file, '%s\n', 'Flmap =  0.0 1.0  2.8800000e+003');
    fprintf(file, '%s\n', 'Flhis =  0.0 1.0  2.8800000e+003');
    fprintf(file, '%s\n', 'Commnt=                              ');
    fprintf(file, '%s\n', 'Commnt=');
    if pmcrit>0
        fprintf(file, '%s%f\n', 'PmCrit=  ',pmcrit);
    end
    fprintf(file, '%s\n', 'TraFrm= #chan.tra#');
    
   fclose(file);
end
