This shows you the differences between two versions of the page.
Both sides previous revision Previous revision Next revision | Previous revision | ||
tools:ocean_nudging2 [2018/06/21 08:56] vsicardi |
tools:ocean_nudging2 [2019/08/12 12:53] (current) vlapin [How to run a nudging experiment with EC-Earth 3.2.3] |
||
---|---|---|---|
Line 1: | Line 1: | ||
======PAGE UNDER CONSTRUCTION. PLEASE BE PATIENT================= | ======PAGE UNDER CONSTRUCTION. PLEASE BE PATIENT================= | ||
- | ==== How to run a nudging experiment with EC-Earth 3.2==== | + | ==== How to run a nudging experiment with EC-Earth 3.2.3 ==== |
- | You can run experiments activating surface restoring or Newtonian Dumping (3D nudging). | ||
- | Both feature are already part of NEMO code, nevertheless with auto-ecearth there are some parameters to activate in the configuration files and some variables to take care of in the namelist. | ||
- | First of all, in both cases, please | + | You can run experiments activating surface restoring or Newtonian Damping (3D nudging). |
- | For surface restoring the location is: | + | Both features are already part of NEMO code, nevertheless with auto-ecearth there are some parameters to activate |
- | If you would like to prepare new reference files, their naming should follow the convention : | + | |
+ | First of all, check that the nudging | ||
+ | |||
+ | In MN4, data for surface restoring | ||
+ | |||
+ | Any new surface restoring files should be located in the above path following this structure: ocean/ | ||
+ | |||
+ | For 3d nudging, | ||
+ | |||
+ | |||
+ | If you would like to prepare new reference files, their naming should follow the convention : REFERENCE_fcMEMB_YYYYMM_grid_T.nc.gz = one file containing 3-dimensional monthly mean temperature and salinity. For a given year the previous year and the following year are linked. So be sure you have all the files before you start. | ||
+ | |||
+ | In theory the user should only specify the name of the file to be used, then auto-ecearth via templates (ecearth3.ini and nemo3.ini) will take care of creating the right link and put them in the right places. If not, then there is a bug to be fixed. No panic. We do it. | ||
+ | |||
- | For 3d nudging, the location is : HSM_DIR/ | ||
- | If you would like to prepare new reference files, their naming should follow the convention : REFERENCE_fcMEMB_YYYYMM_grid_T.nc.gz = one file containing 3-dimensional monthly mean temperature and salinity. | ||
========Surface Restoring======== | ========Surface Restoring======== | ||
Line 24: | Line 33: | ||
- | **Proj.conf settings: | + | STEPS to activate surface restoring: |
- | In proj.conf you activate | + | |
- | You should specify | + | **1. Proj.conf settings: |
+ | |||
+ | - Ordered List Item Activate | ||
+ | - Define | ||
+ | |||
So far you can use oras4, | So far you can use oras4, | ||
+ | |||
Please note that is responsibility of the user to be sure that the surface restoring data to be used are in the right format and in the right place. In case of doubts speak with Valentina or Yohan. | Please note that is responsibility of the user to be sure that the surface restoring data to be used are in the right format and in the right place. In case of doubts speak with Valentina or Yohan. | ||
- | **Namelist settings** | + | **2.Namelist settings** |
- | You can use nemo-only or coupled version | + | Note that, depending on the branch the user wants to use this part is already set-up. See section about available branches |
- | Depending on this the namelist | + | |
+ | You can use nemo-only or coupled version for your nudged experiments. Depending on this option, | ||
+ | |||
+ | |||
+ | Here is an example to activate SST restoring in **namelist_nemo-ORCA1L75-standalone.cfg.sh**. | ||
+ | |||
+ | in section nambsrc you will find: | ||
+ | |||
+ | | ||
+ | |||
+ | This is switched on from configuration file. So do not worry about it. If SURF_RESTO = TRUE in proj.conf, | ||
+ | |||
+ | |||
+ | Here there are some parameters that the user should check: | ||
+ | |||
+ | | ||
+ | & | ||
+ | | ||
+ | ! filename | ||
+ | | ||
+ | | ||
+ | | ||
+ | | ||
+ | | ||
+ | ! | ||
+ | | ||
+ | | ||
+ | | ||
+ | | ||
+ | | ||
+ | | ||
+ | | ||
+ | | ||
+ | | ||
+ | | ||
+ | / | ||
+ | | ||
- | Here is an example to activate SST restoring : | ||
- | ln_ssr = .true. ! Sea Surface Restoring on T and/or S (T => fill namsbc_ssr) | ||
- | !----------------------------------------------------------------------- | ||
- | & | ||
- | !----------------------------------------------------------------------- | ||
- | ! ! file name ! frequency (hours) ! variable ! time interpol. ! clim ! ' | ||
- | ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! ' | ||
- | sn_sst = ' | ||
- | sn_sss = ' | ||
- | cn_dir = ' | ||
- | nn_sstr = 1 ! add a retroaction term in the surface heat flux (=1) or not (=0) | ||
- | nn_sssr = 0 !lolo only for ocean standalone! add a damping term in the surface freshwater flux (=2) | ||
- | ! or to SSS only (=1) or no damping term (=0) | ||
- | rn_dqdt = -40. ! magnitude of the retroaction on temperature [W/m2/K] | ||
- | rn_deds = -864. !lolo 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 = 5.e0 !lolo ABS(Max/ | ||
- | / | ||
Reference data | Reference data | ||
- | Nudging files, as defined in fldread.F90. They can be one file per year (even if there are one record per month in it), or one separate file for each and every month. In the following example it's one file with 12 months, for each variable. | ||
Note that the name of the file contains " | Note that the name of the file contains " | ||
- | sst_data_yXXXX.nc | ||
- | sss_data_yXXXX.nc | ||
- | SST parameters : | + | SST parameters |
nn_sstr | nn_sstr | ||
- | sn_sst | + | sn_sst file name = name of the file, in the format: sst_restore_data. if you use s4 or one of your exp, do not modify it. auto-ecearth will do the job. The name of the nudging file, is defined in fldread.F90, |
- | + | ||
- | file name : Name of the nudging file, as defined in fldread.F90, | + | |
frequency : the frequency of the data, usually 1 month, i.e. -1 | frequency : the frequency of the data, usually 1 month, i.e. -1 | ||
Line 94: | Line 121: | ||
rn_deds (-864 mm/day/psu is the standard value consistent with rn_dqdt = -40) | rn_deds (-864 mm/day/psu is the standard value consistent with rn_dqdt = -40) | ||
- | The user must create the nudging files, according to experiment needs. These files can be created either all the once at the beginning of the simulation or after each leg during the simulation. | + | The user must create the nudging files, according to experiment needs. |
- | The script to create such files will be uploaded here soon. | + | |
The user can activate SST restoring and SSS restoring independently with nn_sstr and nn_sssr | The user can activate SST restoring and SSS restoring independently with nn_sstr and nn_sssr | ||
- | To switch it on, just add nemo: | + | =========== 3D Nudging ========== |
+ | |||
+ | 3D nudging in the ocean is done by using the tracer damping functionality of NEMO (tradmp.F90). The nudging can be don only under the mixed layer or in the entire water column. The user has to nudge both the temperature and salinity, while activating | ||
+ | |||
+ | In order to activate | ||
+ | |||
+ | |||
+ | | ||
+ | & | ||
+ | | ||
+ | | ||
+ | | ||
+ | | ||
+ | | ||
+ | | ||
+ | |||
+ | / | ||
+ | Th | ||
+ | The namelist parameter nn_zdmp sets whether the damping should be applied in the whole water column or only below the mixed layer (defined either on a density | ||
+ | or So criterion). The user should modify it the namelist according to the experiment necessity. | ||
+ | |||
+ | **resto.nc** files contains the restoration coefficient. A default one is located in: / | ||
+ | |||
+ | In proj.conf the use should specify | ||
+ | |||
+ | Example of proj.conf session: | ||
+ | OCE_NUDG = TRUE | ||
+ | |||
+ | # Toward which target ocean is nuged to; essential if " | ||
+ | # A directory | ||
+ | DATA_OCE_NUDG = s4 | ||
+ | # A different DATA_OCE_NUDG member per EXPID member [" | ||
+ | DATA_OCE_NUDG_member = fc0 | ||
+ | |||
+ | |||
+ | At this point the user must manually edit the corresponding namelist coupled or standalone, in the following session: | ||
+ | |||
+ | |||
+ | | ||
+ | & | ||
+ | | ||
+ | | ||
+ | ! (if <0 months) | ||
+ | | ||
+ | | ||
+ | | ||
+ | | ||
+ | | ||
+ | / | ||
+ | |||
+ | |||
+ | The 3D nudging files are genrically renamed to salinity_yXXXX and temperature_yXXXX. The user must be sure that the correct links are created in the inidata folder via eceart3.ini and nemo3.ini templates. | ||
- | Please note that this is an ongoing work. | + | For any information |
- | At BSC we work on this development with the branch bsc_nemo_nudging. No changes have been so far in classic. | + | |
- | For further informations | + |