#!/bin/tcsh
#
# Purpose
# -------
# Runs the CSIRO Mk3L climate system model for one day, in coupled mode.
#
# Usage
# -----
# ./test_cpl
#
# History
# -------
# 2006 May 16	Steven Phipps	Original version
# 2007 Oct 24	Steven Phipps	Updated for change to directory structure
# 2008 Feb 5	Steven Phipps	Modified for the conversion of the coupled
#				model auxiliary files to netCDF
# 2008 Feb 6	Steven Phipps	Modified for the conversion of the atmosphere
#				model auxiliary files to netCDF
# 2008 Mar 8	Steven Phipps	Modified for the conversion of the ocean model
#				and coupled model restart files to netCDF
# 2008 Nov 10	Steven Phipps	Added the line "limit stacksize unlimited"
# 2008 Nov 22	Steven Phipps	Added the line "setenv KMP_STACKSIZE 16M"

# Set the stack sizes
limit stacksize unlimited
setenv KMP_STACKSIZE 16M

# Create a run directory if it doesn't already exist, and delete the contents
# of the directory if it does
set RUN_DIR = ../../run3/
if (-e $RUN_DIR) /bin/rm $RUN_DIR/*
if (! -e $RUN_DIR) mkdir $RUN_DIR

# Copy the model executable to the run directory
cp ../bin/model $RUN_DIR
# Copy code to calculate atmos co2 concentration
cp ../bin/atmCO2conc $RUN_DIR
# Copy code for computing radiation information from atm co2
cp ../bin/radint $RUN_DIR
cp ../bin/pset $RUN_DIR

# Copy the run script to the run directory
cp ./run_coal  $RUN_DIR

# Copy the control file to the run directory
cp ../control/input_cpl_1day $RUN_DIR/input
#cp ../control/input_atm_cable_1yr $RUN_DIR/input
cat ../control/input_obgc >> $RUN_DIR/input

# Copy the restart files to the run directory
cp ../data/atmosphere/restart/rest.start_default $RUN_DIR/rest.start
cp ../data/coupled/restart/oflux.nc_default $RUN_DIR/oflux.nc
cp ../data/ocean/restart/orest.nc_sync $RUN_DIR/orest.org.nc
cp ../data/ocean/restart/nuts.nc $RUN_DIR
cp ../data/ocean/restart/restart.sh $RUN_DIR


# Copy the basic data files to the run directory
cp ../data/atmosphere/basic/* $RUN_DIR
cp ../data/ocean/basic/* $RUN_DIR

# Copy the runoff relocation data to the run directory
cp ../data/atmosphere/runoff/landrun21 $RUN_DIR

# Copy the CO2 radiative data file for 280ppm to the run directory
cp ../data/atmosphere/co2/co2_data.280ppm.18l $RUN_DIR/co2_data.18l

# Copy the sea surface temperatures to the run directory (these are only used
# for diagnostic purposes)
cp ../data/atmosphere/sst/ssta.nc_default $RUN_DIR/ssta.nc

# Copy the flux adjustments to the run directory
cp ../data/coupled/flux_adjustments/hfcor.nc_default $RUN_DIR/hfcor.nc
cp ../data/coupled/flux_adjustments/sfcor.nc_default $RUN_DIR/sfcor.nc
cp ../data/coupled/flux_adjustments/txcor.nc_default $RUN_DIR/txcor.nc
cp ../data/coupled/flux_adjustments/tycor.nc_default $RUN_DIR/tycor.nc
cp ../data/coupled/flux_adjustments/sstcor.nc_default $RUN_DIR/sstcor.nc
cp ../data/coupled/flux_adjustments/ssscor.nc_default $RUN_DIR/ssscor.nc
cp ../data/coupled/flux_adjustments/dtm.nc_default $RUN_DIR/dtm.nc

# Copy the obgc input file
 cp ../data/ocean/obgc/gasx_ocmip2.nc $RUN_DIR
 cp ../data/ocean/obgc/co2_history.dat $RUN_DIR

# Copy the CABLE files
# new atmos restart file to be consistent with atmos only run
cp ../data/atmosphere/restart/csiro_rest.start $RUN_DIR/rest.start
# Copy the CABLE control file to the run directory
cp ../control/cable_mk3l.nml $RUN_DIR
# Copy the CABLE restart file to the run directory
cp ../data/CABLE/restart/csiro_restart_out.nc $RUN_DIR/restart_out.nc
# Copy the basic CABLE data files to the run directory
cp ../data/CABLE/surface_data/*  $RUN_DIR

# files for atmos only run
cp ../data/atmosphere/sst/ssta.nc_default $RUN_DIR/ssta.nc
cp ../data/atmosphere/currents/ocuv.nc_default $RUN_DIR/ocuv.nc

# files for an ocean only run
cp ../data/ocean/sst/sst.nc_monthly $RUN_DIR/sst.nc
cp ../data/ocean/sss/sss.nc_monthly $RUN_DIR/sss.nc
cp ../data/ocean/stress/stress.nc_monthly $RUN_DIR/stress.nc

# Copy files for the computing the evolving atm CO2 concentration
cp ../../data/grids/atmosphere/csiro_agcm_area.nc $RUN_DIR
cp ../../data/grids/ocean/grid_spec_mk3l_128_112_21_v2.nc $RUN_DIR
cp ../data/co2/*.txt $RUN_DIR

# Copy the files to compute the radiative forcing from co2 concentration
cp ../../pre/co2/*ppmv* $RUN_DIR


# Change to the run directory
cd $RUN_DIR

# make restart file for coupled obgc run
restart.sh 

# make the radiation file
# pset only needs to be run once
pset -t n -n 18
radint -c 280
mv co2_data co2_data.18l

# Set model year number to 00001
echo '00001' > year

exit

# Run the model
./model < input | tee output