#!/bin/bash


for NUM in $(seq 2011 2014)
do
infile="air.$NUM.nc"
coords="coords.txt"

while read line
do
   line=$(echo $line | sed -e 's/\r//g')
   station=$(echo $line | cut -d ' ' -f 1)
   #-- skip header line
   if [[ "$station" == "station" ]]; then continue; fi
   #-- select station coordinates
   lat=$(echo $line | cut -d ' ' -f 2)
   lon=$(echo $line | cut -d ' ' -f 3)
   station="${NUM}_${station}_${lat}_${lon}"
   #-- extract the station data
   cdo -remapnn,"lon=${lon}_lat=${lat}" ${infile} ${station}_out.nc
   cdo -s -outputtab,date,time,year,month,day,x,y,lon,lat,lev,name,value \
        ${station}_out.nc | \
        awk 'NR>1 {$8="'${lon}'"; $9="'${lat}'"}1' > ${station}_out.txt
done < $coords




done