BEGIN {
	k = 209.0;
	theta = 45.0;
	moistv = 0.2;
	moistg = 0.9057;
	height = 0.255;
	rad1 = 0.05;
	rad2 = 0.02;
	nd = 713;

	if(wave) k=wave;
	if(m) m=model;
	if(t) theta=t;
	if(mv) moistv=mv;
	if(mg) moistg=mg;
	if(h) height=h;
	if(r1) rad1=r1;
	if(r2) rad2=r2;
	if(n) nd=n;

	if(k==111) {model = "#ksp_g"}  		
	if(k==209) {model = "#ksp_g"}
	if(k==26.2) {model = "#spm_g"}
	if(k==9.0) {model = "#spm_g"}
	
	print "rt";
	print "";
	print " k ",k;
	print " theta",theta;
	print "";
 	print " alg0 std";
 	print "  loops 8.0e3";
 	print "  tol 1.0e-5";
 	print " alg1 std";
	print "  loops 1.0e1";

	print "";
 	print "ground";
 	print "model",model;
 	print " ls   0.15";
 	print " sigmah  0.01";
 	print " chi hallik";
	print "  moistv ",moistv;
 	print "  sand 41.0";
 	print "  clay 22.0";
	print "";
 	print "layer0";
	print "  height ",height; 
 	print "#disc_rg";
	print "  len  3.000e-4";
	print "  rad1 ",rad1;
	print "  rad2 ",rad2;
	print "  nd",nd;	
 	print "dista uniform";
 	print "disti spherical";
	print "chi leafy";
 	print "  temp 22.0";
 	print "  moistg",moistg;
}