#!/bin/csh
#PBS -l elapstim_req=10:00:00
#PBS -l cputim_job=10:00:00   # cputime per job
#PBS -l memsz_job=10gb      # memory per job
#PBS -l cpunum_job=1        # cpus per node
#PBS -q feque 
#PBS -j o
#PBS -l rlg0cap_job=0       # diese Zeile muss immer sein!
#PBS -p 0               #lowest/highest priority = -1024 => 1023
#PBS -N extr # job name


#some ini
 limit datasize unlimited
 limit stacksize unlimited
 module load nco
 setenv NC_BLKSZ 4M

#path to experiment
 set in_dir = /Volumes/work/mom4p1_pubrel_18dec2009/exp/iow05_603_588_47/out/iow05_603_588_47_longrun_I  
 set ou_dir = /Volumes/work/mom4p1_pubrel_18dec2009/exp/iow05_603_588_47/out/iow05_603_588_47_longrun_I 

#delete old output
 rm -f $ou_dir/s_by2.nc
 rm -f $ou_dir/t_by2.nc
 rm -f $ou_dir/s_by5.nc
 rm -f $ou_dir/t_by5.nc
 rm -f $ou_dir/s_by15.nc
 rm -f $ou_dir/t_by15.nc
 rm -f $ou_dir/s_by31.nc
 rm -f $ou_dir/t_by31.nc
 rm -f $ou_dir/s_ll7.nc
 rm -f $ou_dir/t_ll7.nc
 rm -f $ou_dir/s_ane.nc
 rm -f $ou_dir/t_ane.nc
 rm -f $ou_dir/s_f9.nc
 rm -f $ou_dir/t_f9.nc

#extraction
 echo extracing BY2
 ncrcat -d xt_ocean,14.0833 -d yt_ocean,55.     -v salt $in_dir/ocean_salt_19??_??_01.nc $ou_dir/s_by2.nc
 ncrcat -d xt_ocean,14.0833 -d yt_ocean,55.     -v temp $in_dir/ocean_temp_19??_??_01.nc $ou_dir/t_by2.nc

 echo extracting BY5
 ncrcat -d xt_ocean,15.9833 -d yt_ocean,55.2500 -v salt $in_dir/ocean_salt_19??_??_01.nc $ou_dir/s_by5.nc
 ncrcat -d xt_ocean,15.9833 -d yt_ocean,55.2500 -v temp $in_dir/ocean_temp_19??_??_01.nc $ou_dir/t_by5.nc

 echo extrcating BY15
 ncrcat -d xt_ocean,20.0500 -d yt_ocean,57.3333 -v salt $in_dir/ocean_salt_19??_??_01.nc $ou_dir/s_by15.nc
 ncrcat -d xt_ocean,20.0500 -d yt_ocean,57.3333 -v temp $in_dir/ocean_temp_19??_??_01.nc $ou_dir/t_by15.nc

 echo extracting BY31
 ncrcat -d xt_ocean,18.2333 -d yt_ocean,58.5833 -v salt $in_dir/ocean_salt_19??_??_01.nc $ou_dir/s_by31.nc
 ncrcat -d xt_ocean,18.2333 -d yt_ocean,58.5833 -v temp $in_dir/ocean_temp_19??_??_01.nc $ou_dir/t_by31.nc

 echo extracting LL7
 ncrcat -d xt_ocean,24.83   -d yt_ocean,59.85   -v salt $in_dir/ocean_salt_19??_??_01.nc $ou_dir/s_ll7.nc
 ncrcat -d xt_ocean,24.83   -d yt_ocean,59.85   -v temp $in_dir/ocean_temp_19??_??_01.nc $ou_dir/t_ll7.nc

 echo extracting ANHOLT E
 ncrcat -d xt_ocean,12.12   -d yt_ocean,56.67   -v salt $in_dir/ocean_salt_19??_??_01.nc $ou_dir/s_ane.nc
 ncrcat -d xt_ocean,12.12   -d yt_ocean,56.67   -v temp $in_dir/ocean_temp_19??_??_01.nc $ou_dir/t_ane.nc

 echo extracting F9
 ncrcat -d xt_ocean,22.04   -d yt_ocean,64.425  -v salt $in_dir/ocean_salt_19??_??_01.nc $ou_dir/s_f9.nc
 ncrcat -d xt_ocean,22.04   -d yt_ocean,64.425  -v temp $in_dir/ocean_temp_19??_??_01.nc $ou_dir/t_f9.nc

echo extrac_station.sh terminated gracefully
exit