Commit 8afe68e0 authored by Evans, Katherine's avatar Evans, Katherine
Browse files

small changes to script, now moving to a branch

parent 1de20877
......@@ -4,7 +4,7 @@
if [ "$postprocess" = "true" ]; then
#if [ "$version" == "e3sm" ]; then
if [ "$version" == "e3sm" ]; then
#ncwa -a -time $outpath/masks/mask_ne30np4_lnd.nc $outpath/masks/mask_ne30np4_lnd_notime.nc
#for fl in $path/lnd/hist/*h0*; do
#echo "first add mask in original grid to monthly lnd output in $fl"
......@@ -15,10 +15,11 @@ if [ "$postprocess" = "true" ]; then
#h0lnd_path=$outpath/lnd/hist_regrid
#mkdir -p "$h0lnd_path"
#ls $path/lnd/hist/*h0* | ncremap -m $map_file -O $h0lnd_path
#else
h0lnd_path=$path/archive/lnd/hist
#fi
else
h0lnd_path=$path/lnd/hist
fi
echo "make land h0 climos"
mkdir -p "$outpath/lnd/climos"
......
#!/bin/bash
#clear
# 1/18 to run this script on E3SM v1
# 8/18 to run this script on E3SM v1 on CORI at NERSC
#(1) qsub -X -I -A PROJ# -lnodes=1,partition=gpu,walltime=1:00:00
# alter parameters below to match desired settings:
# Key settings:
# job case, location (path), year start and end,
......@@ -12,30 +11,38 @@
# update specialty files used for masking and regridding
# total area of GrIS in model is 1.719e6
project=cli115
case=20180215.DECKv1b_H1.ne30_oEC.edison
year_start=1979
year_end=2005
base=/lustre/atlas1/$project/proj-shared/salil/archive
path=$base/$case/archive
obspath=/lustre/atlas1/$project/world-shared/$USER/obs_data
outpath=/lustre/atlas1/$project/world-shared/$USER/$case/postproc
area_file=/lustre/atlas1/$project/world-shared/$USER/surfdata_0.9x1.25_simyr1850_c110725.nc
map_file=/ccs/home/zender/data/maps/map_ne30np4_to_fv129x256_aave.20150901.nc
mask_file=/lustre/atlas1/$project/world-shared/4ue/$case/postproc/mask_ne30np4_notime.nc
# RHEA
#project=cli115
#case=20180215.DECKv1b_H1.ne30_oEC.edison
#year_start=1979
#year_end=2005
#base=/lustre/atlas1/$project/proj-shared/salil/archive
#path=$base/$case/archive
#obspath=/lustre/atlas1/$project/world-shared/$USER/obs_data
#outpath=/lustre/atlas1/$project/world-shared/$USER/$case/postproc
#area_file=/lustre/atlas1/$project/world-shared/$USER/surfdata_0.9x1.25_simyr1850_c110725.nc
#map_file=/ccs/home/zender/data/maps/map_ne30np4_to_fv129x256_aave.20150901.nc
#mask_file=/lustre/atlas1/$project/world-shared/4ue/$case/postproc/mask_ne30np4_notime.nc
# CORI
project=piscees
case=20180612.B_case.T62_oEC60to30v3wLI.modified_runoff_mapping.edison
year_start=0070
year_end=0075
base= /global/cscratch1/sd/dcomeau/acme_scratch/edison/archive
path=$base/$case/run
obspath=/project/projectdirs/piscees/postproc/obs_data
# outpath does not need to exist yet, it will be created below if its not already
outpath=/global/cscratch1/sd/$USER/$case/postproc
area_file=/project/projectdirs/piscees/postproc/surfdata_0.9x1.25_simyr1850_c110725.nc
map_file=/project/projectdirs/piscees/postproc/map_ne30np4_to_fv129x256_aave.20150901.nc
#mask_file=/project/projectdirs/piscees/postproc/mask_ne30np4_notime.nc
version=e3sm
if [[ $version = cesm1 ]]; then
atm_vars=("PSL" "CLDLOW" "CLDHGH" "CLDMED" "CLDTOT" "TREFHT" "TS" "TSMX" "TSMN" "PRECT" "PRECC" "PRECL")
lnd_h0_vars=("QICE" "TREFMNAV" "TREFMXAV" "TSA" "RAIN" "SNOW")
lnd_h1_vars=("QICE_FRZ" "QICE_MELT" "QSOIL_ICE" "QSNOMELT_ICE" "RAIN" "SNOW")
elif [[ $version = e3sm ]]; then
atm_vars=("PSL" "CLDLOW" "CLDHGH" "CLDMED" "CLDTOT" "TREFHT" "TS" "TSMX" "TSMN" "PRECC" "PRECL")
lnd_h0_vars=("FSDS" "QRUNOFF" "QSOIL" "TSA" "RAIN" "SNOW" "FIRA" "FLDS" "FSH" "QSNOMELT")
#glc_vars=("acab" "thk")
fi
postprocess=true
postprocess_obs=false
postprocess_obs=true
make_plots=true
mkdir -p $outpath
......@@ -85,6 +92,31 @@ source deactivate && source activate LIVVpy2
module list
elif [[ $(hostname -s) = cori* ]]; then
module unload PrgEnv-pgi PrgEnv-gnu PrgEnv-intel
module load PrgEnv-gnu/6.0.4
module unload nco ncl netcdf gcc
module load gcc
module load cray-netcdf/4.4.1.1.6
module load gsl
# module load hdf5-parallel/1.10.1 conflit with netcdf
module load nco/4.7.4
module load python/2.7-anaconda-4.4
module list
# load the conda environment
source /global/homes/k/kennedyj/LIVVkit/conda/etc/profile.d/conda.sh
# For python 3
conda activate lex_py3
# For python 2
# conda activate lex_py2
# . deactivate && . activate LIVVpy2
module list
elif [[ $(hostname -s) = titan* ]]; then
module load netcdf/cray
module load nco/4.6.9
......@@ -110,14 +142,14 @@ else
fi # if mask file needs to be created
fi # CESM vs E3SM
fi # CESM
if [ "$postprocess_obs" = "true" ]; then
echo "process and remap cloud obs"
. atm_obs.sh
../obs/atm_obs.sh
echo "process RACMO"
. lnd_obs.sh
../obs/lnd_obs.sh
if [ "$version" = "cesm1" ]; then
echo "remap RACMO to CESM"
......
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