runid = 'chan_M1e4'; Mmud = 1e-4; pmcrit = -1; 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 = 10'); fprintf(file, '%s\n', 'Thick = 15'); fprintf(file, '%s\n', 'Thick = 20'); fprintf(file, '%s\n', 'Thick = 30'); fprintf(file, '%s\n', 'Thick = 35'); 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', 'C02 = 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', '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