$adidat diff1=1.0d0, !Diffusion constant for eq. 1 diff2=1.0d0, !Diffusion constant for eq. 2 par_a=1.5d0, !Parameter a in non linear term par_b=0.05d0, !Parameter b in non linear term par_c=5.0d0, !Parameter c in non linear term xdom=2.0d3, !Domain size in x direction ydom=2.0d3, !Domain size in y direction tstart=0.0d0, !Starting time step tfinal=4.0d3, !Ending time step ht0=0.5d0, !Size of time step htmax=5.0d0, !Maximum length of time step tslice=3.0d3, !Output time slice for solution file tdrop=1.0d2, !Output time for drop file tscale=1.0d0, !Min^{-1} time scale tss=1.0d9, !Steady state solution check filename1='eve_16.stat', !Stat file filename2='eve_16.soln', !Solution file filename3='eve_16.drop', !Drop file dynflg=.true., !sert to true to start ellflg=.false., !Set to true to use elliptic solver nointflg=.false., !Set to true to skip integrator outflg=.false., !Set to false to start restart=.false., !Set to true if starting from drop file varstp=.false., !Set to true for variable time step perturb=.true., !Set to true to perturb travelling wave ic=2, !Initial conditions: px=1500, ! 1 - cylinder px=py=my=mx mx=1, ! 2 - 1D wave px=py=1, mx=nx my=? py=100, ! my=1, ic1=0.0d0, !Predator density in front of wave ic2=0.8d0, !Prey density in front of wave per1=0.4d0, !Predator density behind wave per2=0.6d0, !Prey density behind wave per3=10, !Dimension of perturbation v_samp=3, !Subsampling for gridview title='Pred-Prey eq3; A=1.5,B=.05,C=5.0; Domain=2000 Eve-R16' $end