config_ORCA2_L31_v36_triolith.sh 7.1 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211
  1. #!/bin/bash
  2. #==========================================================
  3. #
  4. # Configuration file for
  5. #
  6. # OCEAN MONITORING for NEMO v3.6 ORCA2 on 31 levels
  7. #
  8. # HPC: triolith
  9. #
  10. # L. Brodeau, 2015
  11. #
  12. #===========================================================
  13. export CONF=ORCA2 ; # horizontal global configuration
  14. export NBL=31 ; # number of levels
  15. # File system where NEMO config and output files are stored:
  16. export STORE_DIR="/proj/bolinc/users/x_laubr"
  17. # Is it an ec-earth run?
  18. export ece_run=0 ; # means that NEMO files in something like ${STORE_DIR}/<RUN>/output/nemo/<YYY>
  19. # # where YYY starts from '001' to
  20. export Y_INI_EC=1990 ; # initial year if ec-earth run...
  21. # List of suffixed of files that have been saved by NEMO and that are needed for the diags:
  22. export NEMO_SAVED_FILES="grid_T grid_U grid_V icemod SBC"
  23. # Directory structure in which to find NEMO output file (use <ORCA> and <RUN>):
  24. export NEMO_OUT_STRCT="${STORE_DIR}/<ORCA>/<ORCA>-<RUN>-S"
  25. export TSTAMP="1m" ; # output time-frequency stamp as in NEMO output files...
  26. # How does the nemo files prefix looks like
  27. # Everything before "<year_related_info>_grid_<X>" or "<year_related_info>_icemod"
  28. # use <ORCA>, <RUN> and <TSTAMP>=> Ex: export NEMO_FILE_PREFIX="<ORCA>-<RUN>_<TSTAMP>_"
  29. export NEMO_FILE_PREFIX="<ORCA>-<RUN>_<TSTAMP>_"
  30. # => should get rid of TSTAMP actually...
  31. # Temporary file system (scratch) on which to perform the job you can use <JOB_ID> if scracth depends on JOB ID:
  32. export SCRATCH="/scratch/local/<JOB_ID>"
  33. export PYTHON_HOME=${HOME}/opt/Canopy_64bit/User
  34. # If variables names in NEMO files are not the default...
  35. export NN_SST="tos"
  36. export NN_SSS="sos"
  37. export NN_SSH="zos"
  38. export NN_T="thetao"
  39. export NN_S="so"
  40. export NN_MLD="mldr10_1"
  41. export NN_U="uo"
  42. export NN_V="vo"
  43. #export NN_U_EIV="vozoeivu"
  44. #export NN_V_EIV="vomeeivv"
  45. export NN_U_EIV="0" ; # ignore
  46. export NN_V_EIV="0" ; # ignore
  47. export NN_TAUX="tauuo"
  48. export NN_TAUY="tauvo"
  49. export FILE_ICE_SUFFIX="icemod" ; # in what file to find ice fraction and volume?
  50. export NN_ICEF="siconc" ; # name of ice fraction in "FILE_ICE_SUFFIX" file...
  51. export NN_ICET="sivolu" ; # ice thickness or rather volume...
  52. export FILE_FLX_SUFFIX="SBC" ; # in what file to find surface fluxes ?
  53. export NN_FWF="wfo" ; # name of net freshwater flux (E-P-R) in "FILE_FLX_SUFFIX" file...
  54. export NN_EMP="emp_oce" ; # name of E-P in "FILE_FLX_SUFFIX" file...
  55. export NN_P="precip" ; # name of P in "FILE_FLX_SUFFIX" file...
  56. export NN_RNF="XXX" ; # name of continental runoffs in "FILE_FLX_SUFFIX" file...
  57. export NN_CLV="XXX" ; # calving from icebergs in "FILE_FLX_SUFFIX" file...
  58. export NN_E="XXX" ; # evaporation in "FILE_FLX_SUFFIX" file...
  59. export L_CONV2NC3=false ; # Set to true if your NEMO output is in Netcdf4 and your NCO does not support netcdf4!
  60. export L_RENAME=false ; # set to true if your ORCA output has old name convention (ex: votemper instead of thetao)
  61. export EXTRA_CONF="NEMO v3.6 ${CONF} (L${NBL}) - LIM3 / ocean-only experiment"
  62. # Land-sea mask and basins files:
  63. export MM_FILE="/proj/bolinc/users/x_laubr/${CONF}/mesh_mask_${CONF}_20151014.nc"
  64. export BM_FILE="/proj/bolinc/users/x_laubr/${CONF}/basin_mask_${CONF}.nc"
  65. # 3D monthly climatologies of potential temperature and salinity (can be those you used for the NEMO run):
  66. export F_T_CLIM_3D_12=${STORE_DIR}/${CONF}/${CONF}-I/data_1m_potential_temperature_nomask_${CONF}.nc
  67. export F_S_CLIM_3D_12=${STORE_DIR}/${CONF}/${CONF}-I/data_1m_salinity_nomask_${CONF}.nc
  68. export SST_CLIM_12=${STORE_DIR}/${CONF}/${CONF}-I/sst_data_${CONF}.nc
  69. export NN_T_CLIM="votemper"
  70. export NN_S_CLIM="vosaline"
  71. export NN_SST_CLIM="sst"
  72. export ICE_CLIM_12="${STORE_DIR}/${CONF}/${CONF}-I/ice_cover_180x360-${CONF}_Hurrell_monthly_mean1980-1999.nc4"
  73. export NN_ICEF_CLIM="ice_cover"
  74. # A text file where the vertical hydraugraphical sections of interest are defined :
  75. #export TRANSPORT_SECTION_FILE="${BARAKUDA_ROOT}/data/transportiz_${CONF}_light.dat"
  76. export TRANSPORT_SECTION_FILE="${BARAKUDA_ROOT}/data/transportiz_${CONF}_all.dat"
  77. # For transport by sigma-class:
  78. export DENSITY_SECTION_FILE="${BARAKUDA_ROOT}/data/dens_section_${CONF}.dat"
  79. # In what directory of the local machine to save the diagnostics:
  80. export DIAG_DIR="/proj/bolinc/users/x_laubr/tmp/barakuda/${CONF}_v36"
  81. # Files with the list of rectangular boxes to look at more closely:
  82. export FILE_DEF_BOXES="${BARAKUDA_ROOT}/data/def_boxes_convection_${CONF}.txt"
  83. # About remote HOST to install HTML pages to:
  84. ihttp=1 ; # do we export on a remote http server (1) or keep on the local machine (0)
  85. RHOST=misu228.misu.su.se ; # remote host to send diagnostic page to///
  86. RUSER=laurent ; # username associated to remote host (for file export)
  87. RWWWD=/data/www/barakuda/${CONF} ; # directory of the local or remote host to send the diagnostic page to
  88. #########################
  89. # Diags to be performed #
  90. #########################
  91. # In what format should figures be produced:
  92. export FIG_FORM="png"
  93. # Basic 3D and surface averages:
  94. i_do_mean=1
  95. # FreshWater fluxes at the surface spatially averaged over the ocean, E-P-R, E-P, R, P, ...
  96. #i_do_fwf=1
  97. # AMOC:
  98. i_do_amoc=1
  99. export LMOCLAT="20-23 30-33 40-43 45-48 50-53" ; # List of latitude bands to look in for max of AMOC
  100. # Transport of mass, heat and salt through specified sections (into TRANSPORT_SECTION_FILE):
  101. i_do_trsp=2 ; # transport of mass, heat and salt through specified sections
  102. # # i_do_trsp=2 => treat also different depths range!
  103. z1_trsp=100 ; # first depth: i_do_trsp must be set to 2
  104. z2_trsp=1000 ; # second depth: i_do_trsp must be set to 2
  105. # meridional heat/salt transport (advective)
  106. i_do_mht=1
  107. # Transport by sigma class
  108. i_do_sigt=1
  109. # sea-ice diags
  110. i_do_ice=1 ; # Sea-ice diags
  111. i_do_bb=1 ; # Budget and other stuffs on a given rectangular box!
  112. # # => needs file FILE_DEF_BOXES !!!
  113. # => produces time-series f(t) (mean of 2D fields)
  114. i_do_ssx_box=0 ; # zoom on given boxes (+spatially-averaged values) for surface properties
  115. # # boxes defined into barakuda_orca.py ...
  116. # Vertical profiles on of box-averaged as a function of time...
  117. i_do_box_TS_z=1 ; # do sigma vert. profiles on given boxes... # 1 => no figures, 2 => figures
  118. # # => needs file FILE_DEF_BOXES !!!
  119. # => produces time-series f(t,z)
  120. #
  121. # Deep Mixed volume in prescribed boxes:
  122. i_do_dmv=0
  123. export MLD_CRIT="1000,725,500"
  124. # Some nerdy stuffs about the critical depth in prescribed boxes:
  125. i_do_zcrit=0
  126. # BETA / TESTING:
  127. # Fresh-water transport associated to sea-ice transport
  128. # => must compile cdficeflux.x but depends on more recent CDFTOOLS module...
  129. i_do_icet=0 ; # treat sea-ice volume transport!
  130. export TRANSPORT_ICE_SECTION_FILE="${BARAKUDA_ROOT}/data/transportiz_${CONF}_ARCTIC.dat"
  131. i_do_flx=0 ; # surface fluxes diags
  132. i_do_amo=0 ; # buit a SST time serie usable to build Atlantic Multidecadal Oscilation index
  133. i_do_sect=0 ; # do sigma vert. profiles on given boxes...
  134. VSECT_NM=( "Indian_77p5_E" "Atlantic_21p5_W" )
  135. VSECT_JI=( "5,5" "266,266" ) ; # X range in C convention
  136. VSECT_JJ=( "25,170" "7,291" ) ; # Y range in C convention