******************** eORCA025.L121-OPM026 ******************** Summary ======= Compare to eORCA025.L121-OPM021, we improve the model behavior at the tip of Antarctic Peninsula and sss restoring is applied as a corretion from 021 run. Input files =========== All the inputs files needed to run this configuration are here: https://ige-meom-opendap.univ-grenoble-alpes.fr/thredds/catalog/meomopendap/extract/eORCA025.L121/eORCA025.L121-OPM026/catalog.html?dataset=meomscanpublic/eORCA025.L121/eORCA025.L121-OPM026/eORCA025.L121-OPM026_input.tar A detailed description of each file is given in :ref:`eORCA025_input_files` Namelist ======== Only the change compare to the reference (eORCA025.L121-OPM021) are mentioned in the next sections. Full files are available here: https://github.com/pmathiot/NEMOCFG/tree/main/cfgs/eORCA025.L121-OPM026/cfgs/eORCA025.L121-OPM026/EXPREF namelist_oce ------------ - in order to have clean comparison later on, we enforce the restoring term to be a specified corrective term of emp (nn_ssr = 3, need code changes wrt NEMO distri). In this case sssr output from eORCA021.L121-OPM021 is used. .. code-block:: console !----------------------------------------------------------------------- &namsbc_ssr ! surface boundary condition : sea surface restoring (ln_ssr =T) !----------------------------------------------------------------------- nn_sstr = 0 ! add a retroaction term to the surface heat flux (=1) or not (=0) rn_dqdt = -40. ! magnitude of the retroaction on temperature [W/m2/K] nn_sssr = 3 ! add a damping term to the surface freshwater flux (=2) ! ! or to SSS only (=1) or no damping term (=0) or read from a file (=3) rn_deds = -166.67 ! magnitude of the damping on salinity [mm/day] ln_sssr_bnd = .true. ! flag to bound erp term (associated with nn_sssr=2) rn_sssr_bnd = 4.e0 ! ABS(Max/Min) value of the damping erp term [mm/day] nn_sssr_ice = 1 ! control of sea surface restoring under sea-ice ! 0 = no restoration under ice : * (1-icefrac) ! 1 = restoration everywhere ! >1 = enhanced restoration under ice : 1+(nn_icedmp-1)*icefrac cn_dir = './' ! root directory for the SST/SSS data location !___________!_________________________!___________________!___________!_____________!________!___________!___________!__________!_______________! ! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights e ! rotation ! land/sea mask ! ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! filename ! sn_sst = 'NOT_USED.nc' , 24. , 'sst' , .false. , .false., 'yearly' , '' , '' , '' sn_sss = 'eORCA025_sss_WOA2018_c3.0_v19812010.5.1' , -1. , 'sosaline' , .true. , .true. , 'yearly' , '' , '' , '' / !----------------------------------------------------------------------- &namsbc_ssr_drk ! surface boundary condition : sea surface restoring (ln_ssr =T) !----------------------------------------------------------------------- ln_sssr_flt = .false. ! use filtering of SSS model for sss restoring nn_shap_iter = 300 ! number of iteration of the shapiro filter ln_sssr_msk = .true. ! use a mask near the coast !___________!____________________!___________________!__________!_____________!________!___________!__________!__________!_______________! ! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask ! ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! filename ! sn_coast = 'eORCA025_distcoast_b0.2_v0.0' , 0. , 'Tcoast' , .false. , .true. , 'yearly' , '' , '' , '' sn_empc = 'eORCA025.L121-OPM021_y1999.1m.20y_flxT', -1. , 'sowafld' , .true. , .true. , 'yearly' , '' , '' , '' rn_dist = 150. ! distance to the coast / namelist_ice ------------ - in order to increase the sea ice production, we decrease the max sea ice fraction in the SH from 0.997 to 0.95 (same number as in IPCLCM6_LR) .. code-block:: console !------------------------------------------------------------------------------ &nampar ! Generic parameters !------------------------------------------------------------------------------ jpl = 5 ! number of ice categories nlay_i = 2 ! number of ice layers nlay_s = 1 ! number of snow layers (only 1 is working) ln_virtual_itd = .false. ! virtual ITD mono-category parameterization (jpl=1 only) ! i.e. enhanced thermal conductivity & virtual thin ice melting ln_icedyn = .true. ! ice dynamics (T) or not (F) ln_icethd = .true. ! ice thermo (T) or not (F) rn_amax_n = 0.997 ! maximum tolerated ice concentration NH rn_amax_s = 0.95 ! maximum tolerated ice concentration SH cn_icerst_in = "restart_ice" ! suffix of ice restart name (input) cn_icerst_out = "restart_ice" ! suffix of ice restart name (output) cn_icerst_indir = "" ! directory to read input ice restarts cn_icerst_outdir = "" ! directory to write output ice restarts / Input files =========== - Use eORCA025.L121_domain_cfg_b0.5_c3.0_d1.0.nc instead of eORCA025.L121_domain_cfg_b0.3_c3.0_d1.0.nc - Use eORCA025_shlat2d_v0.2.nc instead of eORCA025_shlat2d_v0.0.nc - Use eORCA025_bfr2d_v0.2.nc instead of eORCA025_bfr2d_v0.0.nc Code changes ============ - sbcssr.F90 has been changed to include possibilities to prescribed a correction on emp based on previous run sssr output. Full files are available here: https://github.com/pmathiot/NEMOCFG/tree/main/cfgs/eORCA025.L121-OPM026/src/ Monitoring ========== .. _eORCA025.L121-OPM026_monitoring: Global indicators ----------------- On these plot you can find a time series of: - ACC transport - AMOC at rapid array - AMHT at rapid array - Net global heat fluxes - mean sst in the southern ocean (see box in the map) - mean sst in the North West Corner (see box in the map) - sea ice extent (arctic/ant in summer/winter) .. image:: _static/VALGLO_PEN.png Regional indicators ------------------- On these plot, you can find time series of: - ACC transport - Maximum of the Weddell and Ross Gyre (box where the max compute show in the map) - Mean bottom salinity over the main dense water formation hot spot (West Ross and West FRIS) - Mean bottom temperature over East Ross and Amudsen sea to monitor CDW intrusion .. image:: _static/VALSO_PEN.png Local indicators ---------------- These plots monitor the evolution of ice shelf melting and the corresponding shelf properties (ROSS, FRIS, PINE, GETZ) .. image:: _static/VALSI_PEN.png Amundsen/Belingshausen seas --------------------------- These plot monitoring the evolution of temperature, salinity and ice shelf melt in Amundsen sea. .. image:: _static/VALAMU_PEN.png Ice shelves ----------- Amery ~~~~~ .. image:: _static/AMER_eORCA025.L121-OPM026.png Ross ~~~~ .. image:: _static/ROSS_eORCA025.L121-OPM026.png Getz ~~~~ .. image:: _static/GETZ_eORCA025.L121-OPM026.png Pine island ~~~~~~~~~~~ .. image:: _static/PINE_eORCA025.L121-OPM026.png George VI ~~~~~~~~~ .. image:: _static/GEVI_eORCA025.L121-OPM026.png Filschner Ronne ~~~~~~~~~~~~~~~ .. image:: _static/FRIS_eORCA025.L121-OPM026.png Riiser ~~~~~~ .. image:: _static/RIIS_eORCA025.L121-OPM026.png Fimbul ~~~~~~ .. image:: _static/FIMB_eORCA025.L121-OPM026.png Evaluation ========== Ice shelf melt: the pictures below are the climatological melt for all the ice shelves in NEMO for the OPM026 and OPM021 simulations over the climatology (1989-1993) by sectors. .. image:: _static/WAIS_OPM026_2009.png :scale: 27% .. image:: _static/EAIS_OPM026_2009.png :scale: 27% .. image:: _static/WEDD_OPM026_2009.png :scale: 27% Map of ice shelf melt for both simulations OPM021 and OPM026 over the last decade of the simulations (2009-2018) with bottom temperature. .. image:: _static/melt_sector_OPM026_2009-2018.png :scale: 20% .. image:: _static/melt_sector_OPM021_2009-2018.png :scale: 20% Map of barotropic stream function for both simulations OPM021 and OPM026 over climatology 2009-2018 .. image:: _static/BSF_SHLAT.png Map of bottomT and bottomS for both simulations OPM021 and OPM026 over the climatology 2009-2018 .. image:: _static/BOTS_REF_comb.png .. image:: _static/BOTT_REF_comb.png