prep.sh - UK-FVCOM-Usergroup/uk-fvcom GitHub Wiki
#!/usr/bin/env bash
set -eu
module load mpi/mpich-x86_64 np=$(grep -c physical\ id /proc/cpuinfo) if [ -z $np ]; then np=1 fi
namelist=templates/namelist.wps Vtable=templates/Vtable years=({2002..2004}) grids=(grids/geo_em.d??.nc)
for year in ${years[@]}; do if [ ! -d ./$year ]; then mkdir ./$year fi
(
cd $year
# Link to the model grids.
for grid in ${grids[@]}; do
if [ ! -h $(basename $grid) ]; then
ln -s ../$grid
fi
done
# Put the Vtable in the current directory.
if [ -h Vtable ]; then
rm -f Vtable
fi
if [ ! -f ./Vtable ]; then
ln -s ../$Vtable
else
echo "WARNING: Vtable found but not a symlink to the default one."
fi
# Do not clobber old namelists.
if [ -f namelist.wps ]; then
mv namelist.wps namelist.wps.$$
fi
cp ../$namelist namelist.wps
sed -i 's/2002-'/$year-'/g' namelist.wps
# Run the link script, ungrib and then finally metgrid. Do grib1 first,
# then grib2 for 2007, otherwise just do whatever we find.
if [ $year -eq 2007 ]; then
../bin/link_grib.csh ../../$year/*.grib1
../bin/ungrib.exe
../bin/link_grib.csh ../../$year/*.grib2
../bin/ungrib.exe
else
../bin/link_grib.csh ../../$year/*.grib?
../bin/ungrib.exe
fi
mpirun -n $np ../bin/metgrid.exe
# Tidy up
rm Vtable namelist.wps GRIBFILE.???
for grid in ${grids[@]}; do
rm $(basename $grid)
done
)
done