<p><b>dwj07@fsu.edu</b> 2012-03-19 15:57:27 -0600 (Mon, 19 Mar 2012)</p><p><br>
        -- BRANCH COMMIT --<br>
<br>
        Adding a temporal convergence test case.<br>
</p><hr noshade><pre><font color="gray">Deleted: branches/ocean_projects/ocean_test_cases_staging/ocean/temporal_convergence/makeMeshes.sh
===================================================================
--- branches/ocean_projects/ocean_test_cases_staging/ocean/baroclinic_channel/makeMeshes.sh        2012-03-19 15:32:38 UTC (rev 1672)
+++ branches/ocean_projects/ocean_test_cases_staging/ocean/temporal_convergence/makeMeshes.sh        2012-03-19 21:57:27 UTC (rev 1676)
@@ -1,265 +0,0 @@
-#!/bin/bash
-
-########################################################################################
-## Change the horizontal spacings, vertical levels, and cells in the x direction here ##
-## To start, grids are periodic in the x direction. This may change later ##
-########################################################################################
-SPACINGS="1000 4000 10000"
-VERTLEVS="20"
-TCNAME="baroclinic_channel"
-
-###############################################################
-## Change reference spacing, time_step, and viscosities here ##
-###############################################################
-REF_TIME_STEP="60"
-REF_VISC_H="10.0"
-REF_VISC_V="0.0001"
-REF_SPACING="10000.0"
-TOTAL_TIME="17280000" # in seconds
-
-echo "Building ${TCNAME} meshes"
-
-CUR_DIR=`pwd`
-
-if [ -z "$3" ]; then
-        PROCS="2 4 8 16 32"
-else
-        PROCS="$3"
-fi
-
-X_EXTENT="160000"
-Y_EXTENT="500000"
-
-#######################################################
-## Setup Variables for different run initializations ##
-#######################################################
-TIME_STEPPERS["config_time_integration"]="rk4 se1 se20 use"
-TIME_STEPS["config_dt"]="1 5 10 20 40 80 160"
-
-############################################################
-## Setup Default namelist keys and values to update later ##
-############################################################
-i=0
-KEYS[$i]="config_run_duration"; VALUES[$i]="'0000_00:02:40'"; i=$i+1;
-KEYS[$i]="config_input_name"; VALUES[$i]="'grid.nc'"; i=$i+1;
-KEYS[$i]="config_output_name"; VALUES[$i]="'output.nc'"; i=$i+1;
-KEYS[$i]="config_restart_name"; VALUES[$i]="'restart.nc'"; i=$i+1;
-KEYS[$i]="config_output_interval"; VALUES[$i]="'00_00:02:40'"; i=$i+1;
-KEYS[$i]="config_do_restart"; VALUES[$i]=".false."; i=$i+1;
-KEYS[$i]="config_vert_grid_type"; VALUES[$i]="'zstar'"; i=$i+1;
-KEYS[$i]="config_rho0"; VALUES[$i]="1000"; i=$i+1;
-KEYS[$i]="config_bottom_drag_coeff"; VALUES[$i]="1.0e-2"; i=$i+1;
-KEYS[$i]="config_vert_visc_type"; VALUES[$i]="'const'"; i=$i+1;
-KEYS[$i]="config_vert_diff_type"; VALUES[$i]="'const'"; i=$i+1;
-KEYS[$i]="config_eos_type"; VALUES[$i]="'linear'"; i=$i+1;
-KEYS[$i]="config_monotonic"; VALUES[$i]=".true."; i=$i+1;
-KEYS[$i]="config_vert_tracer_adv_order"; VALUES[$i]="3"; i=$i+1;
-KEYS[$i]="config_horiz_tracer_adv_order"; VALUES[$i]="3"; i=$i+1;
-KEYS[$i]="config_thickness_adv_order"; VALUES[$i]="3"; i=$i+1;
-
-#################################################
-## Build perfect hex meshes using periodic_hex ##
-#################################################
-cd periodic_hex
-
-if [ -a Makefile.front ]; then
-        cat Makefile.front > Makefile
-        cat Makefile.end >> Makefile
-else
-        cp Makefile.bak Makefile
-fi
-
-
-make clean > /dev/null
-make > /dev/null
-
-for SPACING in $SPACINGS
-do
-        NAME=`echo "${SPACING}m"`
-        NY=`echo "(${Y_EXTENT} / ${SPACING}) + 2" | bc`
-        NX=`echo "(${X_EXTENT} / ${SPACING})" | bc`
-        DC=`echo "${SPACING}.0"`
-
-        echo " Creating ${NAME} base mesh"
-
-        cat namelist.input.template | sed "s/*NX/${NX}/g" | sed "s/*NY/${NY}/g" | sed "s/*DC/${DC}/g" > namelist.input
-        ./periodic_grid
-
-        mv grid.nc ${TCNAME}_${NAME}.grid.nc
-done
-
-mv ${TCNAME}*.grid.nc $CUR_DIR/.
-make clean > /dev/null
-
-rm Makefile
-
-##############################################################################
-## Clear run_paths file, which will be used to submit all test cases later. ##
-##############################################################################
-
-cd $CUR_DIR
-
-if [ -a run_paths ]; then
-        rm run_paths
-fi
-
-########################################
-## Setup Template Namelist input file ##
-########################################
-if [ -a MPAS-namelist.input.repo ]; then
-        cp MPAS-namelist.input.repo MPAS-namelist.input.temporary
-else
-        cp MPAS-namelist.input.default MPAS-namelist.input.temporary
-fi
-
-for ((i=0; i<${#KEYS[@]}; i++));
-do
-        cat MPAS-namelist.input.temporary | sed "s/${KEYS[$i]} .*/${KEYS[$i]} = ${VALUES[$i]}/g" > Temp
-        mv Temp MPAS-namelist.input.temporary
-done
-
-
-mv MPAS-namelist.input.temporary MPAS-namelist.input.template
-
-touch run_paths
-
-#################################################################
-## Generate full meshes, with initial conditions, using basin. ##
-#################################################################
-
-for VERTLEV in $VERTLEVS
-do
-        if [ -a grid.nc ]; then
-                unlink grid.nc
-        fi
-
-        if [ -a map ]; then
-                rm map
-        fi
-
-        echo " Bulding basin for ${VERTLEV} levels"
-
-        ## Build Basin
-        cd basin_src
-        cat basin-template.F | sed "s/*VERTLEVS/${VERTLEV}/g" > basin.F
-
-        if [ -a Makefile.front ]; then
-                cat Makefile.front > Makefile
-                cat Makefile.end >> Makefile
-        else
-                cp Makefile.bak Makefile
-        fi
-
-        make clean > /dev/null
-        make > /dev/null
-        cd ../
-        cp basin_src/map .
-
-        ## Call basin, for each perfect hex mesh.
-        for SPACING in $SPACINGS
-        do
-                NAME=`echo "${SPACING}m"`
-                NY=`echo "(${Y_EXTENT} / ${SPACING}) + 2" | bc`
-                NX=`echo "(${X_EXTENT} / ${SPACING})" | bc`
-                DC=`echo "${SPACING}.0"`
-
-                ## Compute scaled spacing, time_step, and viscosities
-                D_SPACING=`echo "scale=4; ${SPACING}/${REF_SPACING}" | bc`
-                TIME_STEP=`echo "scale=4; $D_SPACING * $REF_TIME_STEP" | bc`
-                VISC_H=`echo "scale=4; ($D_SPACING^2)*${REF_VISC_H}" | bc`
-                VISC_H=`echo ${REF_VISC_H}` # No Scaling in Horizontal
-                VISC_V=`echo ${REF_VISC_V}` # No Scaling in Vertocal
-                STATS=`echo " $TOTAL_TIME / $TIME_STEP / 10 " | bc`
-
-                echo " Converting ${NAME} mesh to have ${VERTLEV} levels"
-
-                ln -s ${TCNAME}_${NAME}.grid.nc grid.nc
-
-                mkdir -p dx
-                ./map > /dev/null
-
-                unlink grid.nc
-
-                mkdir -p ${NAME}_${VERTLEV}levs
-
-                mv ocean.nc ${TCNAME}_${NAME}_${VERTLEV}levs.grid.nc
-                mv graph.info ${TCNAME}_${NAME}_${VERTLEV}levs.graph.info
-
-                for TIME_STEPPER in ${TIME_STEPPERS}
-                do
-                        if [ ${TIME_STEPPER} == 'rk4' ]; then
-                                TIME_INTEGRATOR="'RK4'"
-                        elif [ ${TIME_STEPPER:0:2} == "se" ]; then
-                                TIME_INTEGRATOR="'split_explicit'"
-                                SUB_CYCLES=${TIME_STEPPER:2}
-                        elif [ ${TIME_STEPPER} == 'use' ]; then
-                                TIME_INTEGRATOR="'unsplit_explicit'"
-                        fi
-
-                        for TIME_STEP in ${TIME_STEPS}
-                        do
-
-                                for PROC in $PROCS
-                                do
-
-                                        RUN_NAME=${NAME}_${VERTLEV}levs/${TIME_STEPPER}_${TIME_STEP}/${PROC}procs
-
-                                        mkdir -p ${RUN_NAME}
-
-                                        pmetis ${TCNAME}_${NAME}_${VERTLEV}levs.graph.info $PROC > /dev/null
-
-                                        ln -f -s ${CUR_DIR}/${NAME}_${VERTLEV}levs/${TCNAME}_${NAME}_${VERTLEV}levs.graph.info ${RUN_NAME}/graph.info
-                                        ln -f -s ${CUR_DIR}/${NAME}_${VERTLEV}levs/${TCNAME}_${NAME}_${VERTLEV}levs.graph.info.part.${PROC} ${RUN_NAME}/graph.info.part.${PROC}
-                                        ln -f -s ${CUR_DIR}/${NAME}_${VERTLEV}levs/${TCNAME}_${NAME}_${VERTLEV}levs.grid.nc ${RUN_NAME}/grid.nc
-                                        ln -f -s ${CUR_DIR}/${NAME}_${VERTLEV}levs/dx ${RUN_NAME}/dx
-#                                        cp -R dx ${RUN_NAME}/dx
-
-                                        ## Copy executable to run directory
-                                        if [ $# -ge 1 ]; then
-                                                if [ -e ${RUN_NAME}/ocean_model.exe ]; then
-                                                        unlink ${RUN_NAME}/ocean_model.exe
-                                                fi
-
-                                                ln -s $1 ${RUN_NAME}/ocean_model.exe
-                                        fi
-
-                                        ## Copy run information to run directory
-                                        if [ $# -ge 2 ]; then
-                                                cp $2 ${RUN_NAME}/.
-                                        fi
-
-                                        echo "${CUR_DIR}/${RUN_NAME}" >> run_paths
-
-                                        ## Generate new namelist.input file for test case
-                                        cat MPAS-namelist.input.template \
-                                                | sed "s/config_time_integration .*/config_time_integration = ${TIME_INTEGRATOR}/g" \
-                                                | sed "s/config_n_btr_subcycles .*/config_n_btr_subcycles = ${SUB_CYCLES}/g" \
-                                                | sed "s/config_dt .*/config_dt = ${TIME_STEP}/g" \
-                                                | sed "s/config_h_mom_eddy_visc2 .*/config_h_mom_eddy_visc2 = ${VISC_H}/g" \
-                                                | sed "s/config_h_tracer_eddy_diff2 .*/config_h_tracer_eddy_diff2 = 0.0/g" \
-                                                | sed "s/config_vert_visc .*/config_vert_visc = ${VISC_V}/g" \
-                                                | sed "s/config_vert_diff .*/config_vert_diff = 0.0/g" \
-                                                | sed "s/config_dt .*/config_dt = ${TIME_STEP}/g" \
-                                                | sed "s/config_stats_interval .*/config_stats_interval = ${STATS}/g" \
-                                                > ${RUN_NAME}/namelist.input
-                                done
-                        done
-                done
-
-                mv ${TCNAME}_${NAME}_${VERTLEV}levs* ${NAME}_${VERTLEV}levs/.
-                mv dx ${NAME}_${VERTLEV}levs/.
-        done
-
-        rm map
-done
-
-rm -rf dx
-rm MPAS-namelist.input.template
-rm ${TCNAME}*
-rm fort.*
-
-cd basin_src
-make clean > /dev/null
-rm Makefile
-cd ${CUR_DIR}
-
Copied: branches/ocean_projects/ocean_test_cases_staging/ocean/temporal_convergence/makeMeshes.sh (from rev 1674, branches/ocean_projects/ocean_test_cases_staging/ocean/baroclinic_channel/makeMeshes.sh)
===================================================================
--- branches/ocean_projects/ocean_test_cases_staging/ocean/temporal_convergence/makeMeshes.sh         (rev 0)
+++ branches/ocean_projects/ocean_test_cases_staging/ocean/temporal_convergence/makeMeshes.sh        2012-03-19 21:57:27 UTC (rev 1676)
@@ -0,0 +1,294 @@
+#!/bin/bash
+
+########################################################################################
+## Change the horizontal spacings, vertical levels, and cells in the x direction here ##
+## To start, grids are periodic in the x direction. This may change later ##
+########################################################################################
+SPACINGS="1000 4000 10000"
+VERTLEVS="20"
+TCNAME="temporal_convergence"
+
+###############################################################
+## Change reference spacing, time_step, and viscosities here ##
+###############################################################
+REF_TIME_STEP="300"
+REF_VISC_H="10.0"
+REF_VISC_V="0.0001"
+REF_SPACING="10000.0"
+TOTAL_TIME="17280000" # in seconds
+
+echo "Building ${TCNAME} meshes"
+
+CUR_DIR=`pwd`
+
+if [ -z "$3" ]; then
+        PROCS="2 4 8 16 32"
+else
+        PROCS="$3"
+fi
+
+X_EXTENT="160000"
+Y_EXTENT="500000"
+
+#######################################################
+## Setup Variables for different run initializations ##
+#######################################################
+TIME_STEPPERS["config_time_integration"]="rk4 se1 se20 use"
+TIME_STEPS["config_dt"]="1 5 10 20 40 80 160"
+
+DEFAULT_TIME_INTEGRATOR="'split_explicit'"
+DEFAULT_SUBCYCLES=20
+DEFAULT_RUN_DURATION="'0200_00:00:00'"
+
+############################################################
+## Setup Default namelist keys and values to update later ##
+############################################################
+i=0
+KEYS[$i]="config_run_duration"; VALUES[$i]="'0000_00:02:40'"; i=$i+1;
+KEYS[$i]="config_input_name"; VALUES[$i]="'grid.nc'"; i=$i+1;
+KEYS[$i]="config_output_name"; VALUES[$i]="'output.nc'"; i=$i+1;
+KEYS[$i]="config_restart_name"; VALUES[$i]="'restart.nc'"; i=$i+1;
+KEYS[$i]="config_output_interval"; VALUES[$i]="'00_00:02:40'"; i=$i+1;
+KEYS[$i]="config_do_restart"; VALUES[$i]=".false."; i=$i+1;
+KEYS[$i]="config_vert_grid_type"; VALUES[$i]="'zstar'"; i=$i+1;
+KEYS[$i]="config_rho0"; VALUES[$i]="1000"; i=$i+1;
+KEYS[$i]="config_bottom_drag_coeff"; VALUES[$i]="1.0e-2"; i=$i+1;
+KEYS[$i]="config_vert_visc_type"; VALUES[$i]="'const'"; i=$i+1;
+KEYS[$i]="config_vert_diff_type"; VALUES[$i]="'const'"; i=$i+1;
+KEYS[$i]="config_eos_type"; VALUES[$i]="'linear'"; i=$i+1;
+KEYS[$i]="config_monotonic"; VALUES[$i]=".true."; i=$i+1;
+KEYS[$i]="config_vert_tracer_adv_order"; VALUES[$i]="3"; i=$i+1;
+KEYS[$i]="config_horiz_tracer_adv_order"; VALUES[$i]="3"; i=$i+1;
+KEYS[$i]="config_thickness_adv_order"; VALUES[$i]="3"; i=$i+1;
+
+#################################################
+## Build perfect hex meshes using periodic_hex ##
+#################################################
+cd periodic_hex
+
+if [ -a Makefile.front ]; then
+        cat Makefile.front > Makefile
+        cat Makefile.end >> Makefile
+else
+        cp Makefile.bak Makefile
+fi
+
+
+make clean > /dev/null
+make > /dev/null
+
+for SPACING in $SPACINGS
+do
+        NAME=`echo "${SPACING}m"`
+        NY=`echo "(${Y_EXTENT} / ${SPACING}) + 2" | bc`
+        NX=`echo "(${X_EXTENT} / ${SPACING})" | bc`
+        DC=`echo "${SPACING}.0"`
+
+        echo " Creating ${NAME} base mesh"
+
+        cat namelist.input.template | sed "s/*NX/${NX}/g" | sed "s/*NY/${NY}/g" | sed "s/*DC/${DC}/g" > namelist.input
+        ./periodic_grid
+
+        mv grid.nc ${TCNAME}_${NAME}.grid.nc
+done
+
+mv ${TCNAME}*.grid.nc $CUR_DIR/.
+make clean > /dev/null
+
+rm Makefile
+
+##############################################################################
+## Clear run_paths file, which will be used to submit all test cases later. ##
+##############################################################################
+
+cd $CUR_DIR
+
+if [ -a run_paths ]; then
+        rm run_paths
+fi
+
+########################################
+## Setup Template Namelist input file ##
+########################################
+if [ -a MPAS-namelist.input.repo ]; then
+        cp MPAS-namelist.input.repo MPAS-namelist.input.temporary
+else
+        cp MPAS-namelist.input.default MPAS-namelist.input.temporary
+fi
+
+for ((i=0; i<${#KEYS[@]}; i++));
+do
+        cat MPAS-namelist.input.temporary | sed "s/${KEYS[$i]} .*/${KEYS[$i]} = ${VALUES[$i]}/g" > Temp
+        mv Temp MPAS-namelist.input.temporary
+done
+
+
+mv MPAS-namelist.input.temporary MPAS-namelist.input.template
+
+touch run_paths
+
+#################################################################
+## Generate full meshes, with initial conditions, using basin. ##
+#################################################################
+
+for VERTLEV in $VERTLEVS
+do
+        if [ -a grid.nc ]; then
+                unlink grid.nc
+        fi
+
+        if [ -a map ]; then
+                rm map
+        fi
+
+        echo " Bulding basin for ${VERTLEV} levels"
+
+        ## Build Basin
+        cd basin_src
+        cat basin-template.F | sed "s/*VERTLEVS/${VERTLEV}/g" > basin.F
+
+        if [ -a Makefile.front ]; then
+                cat Makefile.front > Makefile
+                cat Makefile.end >> Makefile
+        else
+                cp Makefile.bak Makefile
+        fi
+
+        make clean > /dev/null
+        make > /dev/null
+        cd ../
+        cp basin_src/map .
+
+        ## Call basin, for each perfect hex mesh.
+        for SPACING in $SPACINGS
+        do
+                NAME=`echo "${SPACING}m"`
+                NY=`echo "(${Y_EXTENT} / ${SPACING}) + 2" | bc`
+                NX=`echo "(${X_EXTENT} / ${SPACING})" | bc`
+                DC=`echo "${SPACING}.0"`
+
+                ## Compute scaled spacing, time_step, and viscosities
+                D_SPACING=`echo "scale=4; ${SPACING}/${REF_SPACING}" | bc`
+                SCALED_TIME_STEP=`echo "scale=4; $D_SPACING * $REF_TIME_STEP" | bc`
+                VISC_H=`echo "scale=4; ($D_SPACING^2)*${REF_VISC_H}" | bc`
+                VISC_H=`echo ${REF_VISC_H}` # No Scaling in Horizontal
+                VISC_V=`echo ${REF_VISC_V}` # No Scaling in Vertocal
+
+                echo " Converting ${NAME} mesh to have ${VERTLEV} levels"
+
+                ln -s ${TCNAME}_${NAME}.grid.nc grid.nc
+
+                mkdir -p dx
+                ./map > /dev/null
+
+                unlink grid.nc
+
+                mkdir -p ${NAME}_${VERTLEV}levs
+                
+                for TIME_STEPPER in ${TIME_STEPPERS}
+                do
+                        if [ ${TIME_STEPPER} == 'rk4' ]; then
+                                TIME_INTEGRATOR="'RK4'"
+                                SUB_CYCLES=0
+                        elif [ ${TIME_STEPPER:0:2} == "se" ]; then
+                                TIME_INTEGRATOR="'split_explicit'"
+                                SUB_CYCLES=${TIME_STEPPER:2}
+                        elif [ ${TIME_STEPPER} == 'use' ]; then
+                                TIME_INTEGRATOR="'unsplit_explicit'"
+                                SUB_CYCLES=0
+                        fi
+
+                        for TIME_STEP in ${TIME_STEPS}
+                        do
+
+                                for PROC in $PROCS
+                                do
+
+                                        RUN_NAME=${NAME}_${VERTLEV}levs/.batch_runs/${TIME_STEPPER}_${TIME_STEP}/${PROC}procs
+
+                                        mkdir -p ${RUN_NAME}
+
+                                        ln -f -s ${CUR_DIR}/${NAME}_${VERTLEV}levs/graph.info ${RUN_NAME}/graph.info
+                                        ln -f -s ${CUR_DIR}/${NAME}_${VERTLEV}levs/graph.info.part.${PROC} ${RUN_NAME}/graph.info.part.${PROC}
+                                        ln -f -s ${CUR_DIR}/${NAME}_${VERTLEV}levs/grid.nc ${RUN_NAME}/grid.nc
+                                        ln -f -s ${CUR_DIR}/${NAME}_${VERTLEV}levs/dx ${RUN_NAME}/dx
+
+                                        ## Copy executable to run directory
+                                        if [ $# -ge 1 ]; then
+                                                if [ -e ${RUN_NAME}/ocean_model.exe ]; then
+                                                        unlink ${RUN_NAME}/ocean_model.exe
+                                                fi
+
+                                                ln -s $1 ${RUN_NAME}/ocean_model.exe
+                                        fi
+
+                                        ## Copy run information to run directory
+                                        if [ $# -ge 2 ]; then
+                                                cp $2 ${RUN_NAME}/.
+                                        fi
+
+                                        echo "${CUR_DIR}/${RUN_NAME}" >> run_paths
+
+                                        ## Generate new namelist.input file for test case
+                                        STATS=`echo " $TOTAL_TIME / $TIME_STEP / 10 " | bc`
+                                        cat MPAS-namelist.input.template \
+                                                | sed "s/config_time_integration .*/config_time_integration = ${TIME_INTEGRATOR}/g" \
+                                                | sed "s/config_n_btr_subcycles .*/config_n_btr_subcycles = ${SUB_CYCLES}/g" \
+                                                | sed "s/config_dt .*/config_dt = ${TIME_STEP}/g" \
+                                                | sed "s/config_h_mom_eddy_visc2 .*/config_h_mom_eddy_visc2 = ${VISC_H}/g" \
+                                                | sed "s/config_h_tracer_eddy_diff2 .*/config_h_tracer_eddy_diff2 = 0.0/g" \
+                                                | sed "s/config_vert_visc .*/config_vert_visc = ${VISC_V}/g" \
+                                                | sed "s/config_vert_diff .*/config_vert_diff = 0.0/g" \
+                                                | sed "s/config_stats_interval .*/config_stats_interval = ${STATS}/g" \
+                                                > ${RUN_NAME}/namelist.input
+                                done
+                        done
+                done
+
+                mv ocean.nc ${NAME}_${VERTLEV}levs/grid.nc
+                mv graph.info ${NAME}_${VERTLEV}levs/graph.info
+                mv dx ${NAME}_${VERTLEV}levs/.
+
+                cd ${NAME}_${VERTLEV}levs
+                for PROC in $PROCS
+                do
+                        pmetis graph.info $PROC > /dev/null
+                done
+                cd ${CUR_DIR}
+
+                ## Generate new namelist.input file for typical test case
+                STATS=`echo " $TOTAL_TIME / $SCALED_TIME_STEP / 10 " | bc`
+                cat MPAS-namelist.input.template \
+                        | sed "s/config_run_duration .*/config_run_duration = ${DEFAULT_RUN_DURATION}/g" \
+                        | sed "s/config_time_integration .*/config_time_integration = ${DEFAULT_TIME_INTEGRATOR}/g" \
+                        | sed "s/config_n_btr_subcycles .*/config_n_btr_subcycles = ${DEFAULT_SUBCYCLES}/g" \
+                        | sed "s/config_dt .*/config_dt = ${SCALED_TIME_STEP}/g" \
+                        | sed "s/config_h_mom_eddy_visc2 .*/config_h_mom_eddy_visc2 = ${VISC_H}/g" \
+                        | sed "s/config_h_tracer_eddy_diff2 .*/config_h_tracer_eddy_diff2 = ${VISC_H}/g" \
+                        | sed "s/config_vert_visc .*/config_vert_visc = ${VISC_V}/g" \
+                        | sed "s/config_vert_diff .*/config_vert_diff = ${VISC_V}/g" \
+                        | sed "s/config_stats_interval .*/config_stats_interval = ${STATS}/g" \
+                        > ${NAME}_${VERTLEV}levs/namelist.input
+                ## Copy executable to run directory
+                if [ $# -ge 1 ]; then
+                        if [ -e ${NAME}_${VERTLEV}levs/ocean_model.exe ]; then
+                                unlink ${NAME}_${VERTLEV}levs/ocean_model.exe
+                        fi
+
+                        ln -s $1 ${NAME}_${VERTLEV}levs/ocean_model.exe
+                fi
+        done
+
+        rm map
+done
+
+rm -rf dx
+rm MPAS-namelist.input.template
+rm ${TCNAME}*
+rm fort.*
+
+cd basin_src
+make clean > /dev/null
+rm Makefile
+cd ${CUR_DIR}
+
</font>
</pre>