lambda=2160 rm result$lambda\nm.dat for tau in `seq 1 30`; do for reff in `seq 1 25`; do sed 's/LAM/'$lambda'/' disort.tmplt | sed 's/TAU/'$tau'/' > disort.inp sed 's/REFF/'$reff'/' wc.tmplt > wc.dat uvspec15 < disort.inp | gawk '{print '$tau', '$reff', $1}' \ >> result$lambda\nm.dat done done