#!/bin/bash #========================================================== # # Configuration file for # # OCEAN MONITORING for NEMO v3.6 of EC-Earth 3.2.2 on 75 levels # # Machine: zenobe@cenaero # # F.Massonnet, 2017 # #=========================================================== export CONF=ORCA1.L75 ; # horizontal global ORCA configuration export NBL=75 ; # number of levels export HOST=ZENOBE ; # this has no importance at all, it will just become an "info" on the web-page! export MASTERMIND="UCL / $USER" ; # same here, who's the person who designed/ran this simulation? export EXTRA_CONF="NEMO 3.6 + LIM 3 (EC-Earth 3.2)" ; # // same here ... # Path / directory structure in which to find NEMO output file (you can use # and as substitute to your ORCA grid and experiment (EXP) name): export NEMO_OUT_STRCT="/SCRATCH/acad/ecearth/${USER}/run//output/nemo" # Path to root directory where to save the diagnostics (diagnostics for this "CONF"): export DIAG_DIR="/SCRATCH/acad/ecearth/${USER}/barakuda/${CONF}_ece32" # Path to directory containing some 2D and 3D climatologies on the relevant ORCA grid: export CONF_INI_DIR="/home/acad/ucl-elic/fmasson/ec-earth3.2_config_files/barakuda_clim/" # Temporary file system (scratch) on which to perform the job you can use if scracth depends on JOB ID: export SCRATCH="/SCRATCH/acad/ecearth/${USER}/TMP/" export PYTHON_HOME="/projects/acad/ecearth/opt/software/Python/2.7.12-foss-2016b/" ; # HOME to python distribution with matplotlib and basemap ! export DIR_NCVIEW_CMAP="${BARAKUDA_ROOT}/src/ncview_colormaps" # Is it an ec-earth experiment? export ece_exp=2 ; # 0 => not an EC-Earth experiment, it's a "pure" ocean-only NEMO experiment done from traditional NEMO setup # # 1 => it's an OCEAN-ONLY EC-Earth experiment done from a EC-Earth setup # # 2 => it's a COUPLED EC-Earth experiment # # Both 1 and 2 imply that NEMO files are stored in something like # # ${SOMEWHERE}//output/nemo/ # # where YYY starts from '001' to # # If you select '2', make sure 'cdo' is available and working!!! # # 10 => this experiment controled by AutoSubmit (so NEMO files are tared somerwhere?) # export Y_INI_EC=1990 ; # initial year if ece_exp /= 0 !!! export TRES_IFS=255 ; # spectral resolution for IFS, ex: T255 => TRES_IFS=255 ###--- end EC-Earth IFS relate section --- export ATMO_INFO="IFS T${TRES_IFS}" ; # Name of atmospheric model or forcing used (ex: COREv2, DFS5.2, IFS T255, ect...) # List of suffix of files that have been saved by NEMO and contain MONTHLY averages: export NEMO_SAVED_FILES="grid_T grid_U grid_V icemod SBC" export TSTAMP="1m" ; # output time-frequency stamp as in NEMO output files... # In case 3D fields have been saved on an annual mean basis rather than montly: export ANNUAL_3D="" ; # leave blanck "" if 3D fields are in monthly files... export NEMO_SAVED_FILES_3D="" ; # '' # How does the nemo files prefix looks like # Everything before "_grid_" or "_icemod" # use , and => Ex: export NEMO_FILE_PREFIX="-__" export NEMO_FILE_PREFIX="__" # => should get rid of TSTAMP actually... ####### NEMO => what fields in what files ??? ############ # ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ # => depends on the XIOS *.xml setup you used... # => always specify a string for the NN_* variables # USE "X" if the field is not present in your NEMO output # # State variables and others in grid_T files: export NN_SST="tos" export NN_SSS="sos" export NN_SSH="zos" export NN_T="thetao" export NN_S="so" export NN_MLD="mldr10_1" # # State variables and others in grid_U files: export NN_U="uo" export NN_TAUX="tauuo" export NN_U_EIV="0" ; # 0 => ignore # State variables and others in grid_V files: export NN_V="vo" export NN_TAUY="tauvo" export NN_V_EIV="0" ; # 0 => ignore # # Sea-ice fields: export FILE_ICE_SUFFIX="icemod" ; # in what file type extension to find ice fields export NN_ICEF="siconc" ; # name of ice fraction in "FILE_ICE_SUFFIX" file... export NN_ICET="sivolu" ; # ice thickness or rather volume... export NN_ICEU="sivelu" ; # ice U-velocity export NN_ICEV="sivelv" ; # ice V-velocity # # Surface fluxes: export FILE_FLX_SUFFIX="SBC" ; # in what file type extension to find surface fluxes # ++ Surface freswater fluxes: export NN_FWF="wfo" ; # *+/-* name of net freshwater flux (E-P-R) in "FILE_FLX_SUFFIX" file... export NN_EMP="evap_ao_cea+subl_ai_cea-precip" ; # *+/-* name of E-P in "FILE_FLX_SUFFIX" file... export NN_P="precip" ; # name of total precipitation (solid+liquid) in "FILE_FLX_SUFFIX" file... export NN_RNF="runoffs" ; # name of continental runoffs in "FILE_FLX_SUFFIX" file... export NN_CLV="calving" ; # calving from icebergs in "FILE_FLX_SUFFIX" file... export NN_E="evap_ao_cea+subl_ai_cea" ; # *+/-* name of total evaporation in "FILE_FLX_SUFFIX" file... # ++ Surface heat fluxes: export NN_QNET="qt_oce" ; # *+/-* name of total net surface heat flux in "FILE_FLX_SUFFIX" file... export NN_QSOL="rsntds" ; # name of net surface solar flux in "FILE_FLX_SUFFIX" file... # ++ Wind-stress module: export NN_TAUM="taum" ; # name of surface wind stress module in "FILE_FLX_SUFFIX" file... export NN_WNDM="windsp" ; # name of surface wind speed module in "FILE_FLX_SUFFIX" file... # ################################################################################################ # Land-sea mask and basins files: export MM_FILE=/home/acad/ucl-elic/fmasson/ec-earth3.2_config_files/mesh_mask.nc4 export BM_FILE=${BARAKUDA_ROOT}/data/basin_mask_ORCA1_ece3.2_2017.nc4 # OBSERVATIONS / REFERENCES # 3D monthly climatologies of potential temperature and salinity (can be those you used for the NEMO experiment): export NM_TS_OBS="EN4.2.0 [1990-2010]" export F_T_OBS_3D_12=${CONF_INI_DIR}/thetao_1degx1deg-ORCA1.L75_WOA2009_monthly_LB_20160223.nc4 export F_S_OBS_3D_12=${CONF_INI_DIR}/so_1degx1deg-ORCA1.L75_WOA2009_monthly_LB_20160223.nc4 export F_SST_OBS_12=${CONF_INI_DIR}/tos_180x360-ORCA1_Reynolds_monthly_mean1982-2005.nc4 export NN_T_OBS="thetao" export NN_S_OBS="so" export NN_SST_OBS="tos" export F_ICE_OBS_12=${CONF_INI_DIR}/ice_cover_180x360-ORCA1_Hurrell_monthly_mean1980-1999.nc4 export NN_ICEF_OBS="ice_cover" # A text file where the cross sections (to compute transports) are defined : export TRANSPORT_SECTION_FILE="${BARAKUDA_ROOT}/data/transportiz_ORCA1.dat" ; # set i_do_trsp=1 ! export TRANSPORT_SECTION_FILE_ICE="${BARAKUDA_ROOT}/data/transport_ice_ORCA1.dat" ; # set i_do_trsp_ice=1 ! # For transport by sigma-class: export DENSITY_SECTION_FILE="${BARAKUDA_ROOT}/data/dens_section_ORCA1.dat" # Files with the list of rectangular domains to "analyze" more closely: export FILE_DEF_BOXES="${BARAKUDA_ROOT}/data/def_boxes_convection_ORCA1.txt" export FILE_DMV_BOXES="${BARAKUDA_ROOT}/data/def_boxes_convection_ORCA1.txt" # In what format should figures be produced ('png' recommanded, but 'svg' supported!): export FIG_FORM="png" # About remote HOST to send/install HTML pages to: export ihttp=1 ; # do we export on a remote http server (1) or keep on the local machine (0) export RHOST=elixir%gwelic ; # remote host to send diagnostic page to/// export RUSER=$USER ; # username associated to remote host (for file export) export RWWWD=/elic/ftp/$USER/TMP/ ; # directory of the local or remote host to send the diagnostic page to ######################### # Diags to be performed # ######################### # Movies of SST and SSS compared to OBS: export i_do_movi=1 export iffmpeg_x264=0 ; # is, by chance, ffmpeg with support for x264 encoding available on your stystem? => 1 ! # Basic 3D and surface averages: export i_do_mean=1 # IFS surface fluxes of heat and freshwater export i_do_ifs_flx=1 ; # only relevant when ece_exp=2... # AMOC: export i_do_amoc=1 export LMOCLAT="20-23 30-33 40-43 45-48 50-53" ; # List of latitude bands to look in for max of AMOC # Sea-ice diags export i_do_ice=1 ; # Sea-ice diags # Transport of mass, heat and salt through specified sections (into TRANSPORT_SECTION_FILE): export i_do_trsp=1 ; # transport of mass, heat and salt through specified sections # # i_do_trsp=2 => treat also different depths range! z1_trsp=100 ; # first depth: i_do_trsp must be set to 2 z2_trsp=1000 ; # second depth: i_do_trsp must be set to 2 # Solid freshwater transport through sections due to sea-ice drift export i_do_trsp_ice=1 ; # must have i_do_ice=1 # Meridional heat/salt transport (advective) export i_do_mht=1 # Transport by sigma class export i_do_sigt=1 # Budget on pre-defined (FILE_DEF_BOXES) rectangular domains: export i_do_bb=0 ; # Budget and other stuffs on a given rectangular box! # # => needs file FILE_DEF_BOXES !!! # => produces time-series f(t) (mean of 2D fields) # Vertical profiles on of box-averaged as a function of time... export i_do_box_TS_z=0 ; # do sigma vert. profiles on given boxes... # 1 => no figures, 2 => figures # # => needs file FILE_DEF_BOXES !!! # => produces time-series f(t,z) # Deep Mixed volume in prescribed boxes: export i_do_dmv=0 export MLD_CRIT="1000,725,500" # User-defined meridional or zonal cross sections (for temperature and salinity) # => TS_SECTION_FILE must be defined! export i_do_sect=1 export TS_SECTION_FILE="${BARAKUDA_ROOT}/data/TS_sections.dat" # BETA / TESTING / NERDY (at your own risks...): # export i_do_ssx_box=0 ; # zoom on given boxes (+spatially-averaged values) for surface properties # # boxes defined into barakuda_orca.py ... # Some nerdy stuffs about the critical depth in prescribed boxes: export i_do_zcrit=0