#!/bin/bash

cmpset="2000_DATM%ICONNORESMhistorical_CLM51%BGC-CROP_SICE_SOCN_SROF_SGLC_SWAV"

tstmp="$(date +"%y%m%d-%H%M")"

fld=$1

# IF NO ARGUMENT PROVIDED,  set default case name:

if [[ -z ${fld} ]]; then
    
    fld="exice_cold_${cmpset}_$tstmp"
    echo fld
else
    fld=$1
fi

cfld="${fld}"

echo $cfld 

project=nn8002k

ctsmsrcfld=/cluster/home/irismuz/CTSM_EXICE/

#Remove dirs to rerun if previous run failed
apwd=$(pwd)
rm -r ./$cfld
rm -r $USERWORK/archive/$cfld
rm -r $USERWORK/noresm/$cfld

${ctsmsrcfld}/cime/scripts/create_newcase \
--case ./$cfld --res CLM_USRDAT --compset $cmpset \
        --run-unsupported --project nn8002k --machine betzy --handle-preexisting-dirs r --mpilib openmpi \
        --queue normal --pecount L 

cp $0 "$cfld/spinup_from_restart_not_accelerated_clean.sh"

cd $cfld

./xmlchange CLM_FORCE_COLDSTART=off
./xmlchange CLM_ACCELERATED_SPINUP=off
./xmlchange CALENDAR=NO_LEAP
./xmlchange RUN_STARTDATE=1985-01-01
./xmlchange REST_N=30,REST_OPTION=nyears
./xmlchange --subgroup case.run JOB_WALLCLOCK_TIME=55:00:00
./xmlchange SAVE_TIMING=TRUE

./xmlchange STOP_N=30,STOP_OPTION=nyears 
./xmlchange RESUBMIT=5 #180 years
./xmlchange DATM_YR_ALIGN=1985  
./xmlchange DATM_YR_START=1985  
./xmlchange DATM_YR_END=2014
./xmlchange DATM_SKIP_RESTART_READ=TRUE

#./xmlchange NTASKS=1024
./xmlchange NTASKS_CPL=2048,NTASKS_ATM=256,NTASKS_LND=1792,NTASKS_OCN=256,NTASKS_WAV=256,NTASKS_ESP=256,NTASKS_GLC=256,NTASKS_ROF=256,NTASKS_ICE=256
./xmlchange ROOTPE_LND=256

./xmlchange CLM_USRDAT_NAME=Arctic
./xmlchange CLM_USRDAT_DIR=/cluster/shared/noresm/inputdata
./xmlchange LND_DOMAIN_PATH=/cluster/shared/noresm/inputdata
./xmlchange ATM_DOMAIN_MESH=/cluster/shared/noresm/inputdata/wrf2clm_land_ESMFMesh_c290924.nc
./xmlchange PTS_DOMAINFILE=/cluster/shared/noresm/inputdata/domain.lnd.wrf2clm_lnd_noneg_wrf2clm_ocn_noneg.240929.nc
./xmlchange LND_DOMAIN_MESH=/cluster/shared/noresm/inputdata/wrf2clm_land_ESMFMesh_c290924.nc
./xmlchange LND_DOMAIN_FILE=/cluster/shared/noresm/inputdata/domain.lnd.wrf2clm_lnd_noneg_wrf2clm_ocn_noneg.240929.nc
./xmlchange MASK_MESH=/cluster/shared/noresm/inputdata/wrf2clm_land_ESMFMesh_c290924.nc

./xmlchange RUN_TYPE=hybrid
./xmlchange GET_REFCASE=TRUE
./xmlchange RUN_REFCASE="spinup_ICONNORESM_P_2" #180 years
./xmlchange RUN_REFDIR="/cluster/work/users/irismuz/archive/spinup_ICONNORESM_P_2/rest/2165-01-01-00000"
./xmlchange RUN_REFDATE=2165-01-01-00000


echo "done"

# Adjust namelists
echo "  
create_crop_landunit = .true.
enable_water_isotopes = .false.
use_excess_ice = .true.
!the line below should be true only for coldstart
use_excess_ice_streams = .false.
soil_layerstruct_predefined= '20SL_8.5m'
use_bedrock = .true.
hist_empty_htapes = .true.
hist_fincl1 = 'TWS', 'ALTMAX', 'TBOT', 'QBOT', 'TSOI', 'SOILICE', 'SOILLIQ', 'H2OSOI', 'GPP', 'TLAI', 'TOTECOSYSC', 'TOTVEGC', 'TOTSOMC'
hist_fincl1 += 'SUBSIDENCE', 'TOTEXICE_VOL', 'EXCESS_ICE'
hist_mfilt = 1
hist_nhtfrq = 0
fsurdat='/cluster/shared/noresm/inputdata/surfdata.nc'" >> user_nl_clm

./case.setup

./xmlchange ATM_NX=577
./xmlchange ATM_NY=582

./pelayout

./case.build

./case.submit