prep.sh

#!/usr/bin/env bash

Do all the preprocessing needed for each of the years for which I want to run

WRF.

set -eu

For metgrid.exe to run in parallel.

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