#!/bin/bash for i in CHB do for id in `seq 0 30`; do echo "#PBS -N my_job"> job.sh; echo "#PBS -l walltime=96:00:00">> job.sh; echo "#PBS -l nodes=1:ppn=12">> job.sh; echo "#PBS -j oe">> job.sh; echo "#PBS -m abe">> job.sh; echo "perl /nfs/11/utl0343/1000Genome/GenomeDiff_population_v4.pl $i ALL.chr1.phase1_release_v3.20101123.snps_indels_svs.genotypes.vcf.gz 1 $id">> job.sh; qsub job.sh; done done