Commit 0378d2ab authored by Stephan Lorenz's avatar Stephan Lorenz
Browse files

[paleo_quickplots]: add run_qp_driver_ape.sh (preliminary grid paths to be...

[paleo_quickplots]: add run_qp_driver_ape.sh (preliminary grid paths to be updated) including red_list and 2 options for APE ocean plots
parent 5248cff5
#!/bin/bash
#SBATCH --job-name=pyicon_qp_ape
#SBATCH --time=00:30:00
#SBATCH --output=log.o-%j.out
#SBATCH --error=log.o-%j.out
#SBATCH --ntasks=1
#SBATCH --partition=compute,compute2
#SBATCH --account=mh0033
path_pyicon=`(cd .. && pwd)`"/"
qp_driver="${path_pyicon}pyicon/quickplots/qp_driver.py"
source ./conda_act_mistral_pyicon_env.sh
which python
rand=$(cat /dev/urandom | tr -dc 'A-Z' | fold -w 3 | head -n 1)
config_file="./config_qp_${rand}.py"
qp_driver="${path_pyicon}pyicon/quickplots/qp_driver.py"
cat > ${config_file} << %eof%
# --- path to quickplots
path_quickplots = '../all_ape/'
# --- set this to True if the simulation is still running
omit_last_file = True
# --- do ocean and/or atmosphere plots
do_atmosphere_plots = False
do_ocean_plots = True
# --- grid information
gname = 'icon_grid_0036_R02B04_G'
lev = 'L40'
gname_atm = 'r2b4_atm_r0013'
lev_atm = 'L47'
# --- path to interpolation files
#path_grid = '/mnt/lustre01/work/mh0033/m300602/icon/grids/'+gname+'/'
path_grid = '/work/mh0469/m211032/Icon/PYICON/pyicon/Data/'+gname+'/'
path_grid_atm = '/mnt/lustre01/work/mh0033/m300602/icon/grids/'+gname_atm+'/'
path_ckdtree = path_grid+'/ckdtree/'
path_ckdtree_atm = path_grid_atm+'/ckdtree/'
# --- grid files and reference data
fpath_tgrid = path_grid + gname+'.nc'
fpath_tgrid_atm = '/pool/data/ICON/grids/public/mpim/0013/icon_grid_0013_R02B04_G.nc'
fpath_ref_data_oce = path_grid + 'ts_phc3.0_annual_icon_icon_grid_0036_R02B04_G_L40.nc'
fpath_ref_data_atm = '/mnt/lustre01/work/mh0033/m300602/icon/era/pyicon_prepare_era.nc'
#fpath_fx = path_grid + 'slo1400_oce_fx.nc'
fpath_fx = path_grid + 'slo1403APE_oce_fx.nc'
# --- nc file prefixes
oce_def = '_oce_def'
oce_moc = '_oce_moc'
oce_mon = '_oce_mon'
oce_ice = '_oce_ice'
oce_monthly = '_oce_dbg'
atm_2d = '_atm_2d_ml'
atm_3d = '_atm_3d_ml'
atm_mon = '_atm_mon'
# --- nc output
save_data = True
path_nc = '/work/mh0469/m211032/Icon/PYICON/pyicon/TMP/'
# --- parameter for aquaplanet
ocean_depth = 3000
land_facecolor = 'none'
coastlines_color = 'none'
red_list = [
'sst_bias', 'temp_bias_gzave', 'temp_bias_azave', 'temp_bias_ipzave',
'sss_bias', 'salt_bias_gzave', 'salt_bias_azave', 'salt_bias_ipzave',
'temp_azave', 'temp_ipzave',
'salt_azave', 'salt_ipzave',
'ice_concentration_nh', 'ice_thickness_nh', 'snow_thickness_nh',
'ice_concentration_sh', 'ice_thickness_sh', 'snow_thickness_sh',
'ts_ice_extent_nh', 'ts_ice_extent_sh',
'amoc', 'pmoc'
]
# --- time average information (can be overwritten by qp_driver call)
tave_ints = [
['1500-02-01', '1600-01-01'],
]
%eof%
# --- start qp_driver
startdate=`date +%Y-%m-%d\ %H:%M:%S`
path_data="/work/mh0287/m211032/Icon/Git_Icon/icon.nwp.cpldev/experiments/"
run="slo1403"; python -u ${qp_driver} --batch=True ${config_file} --path_data="$path_data/$run/" --run="$run" --tave_int='2100-02-01,2120-01-01'
#run="slo1404"; python -u ${qp_driver} --batch=True ${config_file} --path_data="$path_data/$run/" --run="$run" --tave_int='2030-02-01,2040-01-01'
enddate=`date +%Y-%m-%d\ %H:%M:%S`
rm ${config_file}
echo "--------------------------------------------------------------------------------"
echo "Started at ${startdate}"
echo "Ended at ${enddate}"
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment