Shell Programming and Scripting

input files:
Monthly files: Jan.hdf ..Dec.hdf
location.dat containing the ff. information:
lon lat month
123.5 12.4 1
124.5 12.3 1
125.5 12.2 2
126.5 12.1 2
127.5 12.0 2... and so on

awk ' [awk script to give you your output] '  *.hdf > location.dat

cat *.hdf > location.dat

files ="$(ls *.hdf)"
for i in $files; do
    read lines of location.dat
    plot $i;    
    overlay location points corresponding to similar month
done

[root@hostname test]# ls -rt *.hdf
Jan.hdf  Feb.hdf

[root@hostname test]# cat location.dat
lon lat month
123.5 12.4 1
124.5 12.3 1
125.5 12.2 2
126.5 12.1 2
127.5 12.0 2

[root@hostname test]# ./test.sh
lon lat for the month of Jan.hdf
123.5 12.4
124.5 12.3
lon lat for the month of Feb.hdf
125.5 12.2
126.5 12.1
127.5 12.0

#!/bin/bash

files="$(ls *.hdf)"

for f in $files; do
   for line in location.dat; do
      grdimage $f $range -Csst.cpt -JM$scale -P -B10g5 -K -O >> "$f".ps (this will produce an image map for each month)
      psxy $line $range -JM$scale -B10g5 -Sc0.05 -Gblack -P -K -O "$f".ps (this will plot location points on the output image from earlier command-
this is where i need to tell the program that it will only plot location points from the list for the given month)
    done 
done

#! /bin/bash

mnths=( Jan Feb Mar Apr Jun Jul Aug Sep Oct Nov Dec )

for files in `ls -rt *.hdf`
do
    mt=`basename $files .hdf`
    echo "lon lat for the month of $files" >> $mt.out
    i=0
    for m in ${mnths[@]}
    do
        i=$(($i + 1))
        if [ "$mt" == "$m" ]
        then
            mt_num=$i
            break
        fi
    done
    while read LINE
    do
        [[ `echo $LINE | cut -d" " -f1` == "lon" ]] && continue
        [ `echo $LINE | cut -d" " -f3` -eq $mt_num ] && echo $LINE | cut -d" " -f1-2 >> $mt.out
    done < location.dat
done

lon lat for jan > month1
lon lat for feb > month2