Version 17 (modified by flavoni, 5 years ago) (diff)

Robustness-TestCases

    TEST CASE: EQUATORIAL BOX

    Motivations:

    This equatorial box is done to :

    show how to create in NEMO a test case without modifying code,

    illustrate equatorial circulation, and to illustrate what appens when we remove Coriolis force, stratification, …

    The structure of this test case is done like in CONFIGURATION directory:

    For this test case we created a directory TEST_EQ_BOX01 in which we give:

    namelist_cfg

    iodef.xml

    cpp keys

    We give files:

    coordinates.nc

    bathy_r025.nc TScst_init_r025.nc flx_r025.nc

    (We give also python scripts to create these files in CREATE_FILES directory)

    NEMO code read input files, all is setted in namelist and iodef.xml files, and we give ferret scripts to produce plots.

    This test case wants to show:

    Outils:

    in https://forge.ipsl.jussieu.fr/nemo/browser/branches/2015/dev_r5102_CNRS11_TestCase/NEMOGCM/TOOLS/TEST_CASES

    There are these directories:

    CREATE_FILES:

    In this directory there are PYTHON scripts to create these files :

    bathymetry coordinates with e1,e2 constant coordinates with gphi0=constant (to do not have Coriolis) constant fluxes input file for salinity and temperature constant input file for salinity constant and a temperature profile (given by Levitus profile, averaged for pacific ocean)

    CREATE_PLOTS:

    In this directory there are FERRET scripts to create plots of :

    vertical slides of velocities (between à and 1000 meters) over all bassin (0° - 180°), average between 1°S and 1°N. vertical slides of temperature (between à and 1000 meters) over all bassin (0° - 180°), average between 1°S and 1°N.

    TEST_EQ_BOX01:

    TEST_EQ_BOX01_cst:

    TEST_EQ_BOX01_noCoriolis:

    TEST_EQ_BOX01_tprof:

    TEST_EQ_BOX01_kelvin:

    there are some examples to create a test case of this type :

    Equatorial BOX, grid of 0.25 ° resolution, between latitude: -30°,30°, and longitude 0°,180°,

    with zps coordinates,

    with ubs schema on dynamics and on tracers, with constant TS initial state, with flux only wind (u composant), with flat bathymetry.

    If you want to create a test case what you have to do?

    1. go into (branches/2015/dev_r5102_CNRS11_TestCase/)NEMOGCM/TOOLS/TEST_CASES/CREATE_FILES
    2. read README file
    3. run python make_coordinates.py to generate coordinates_r025.nc
    4. link coordinates.nc to coordinates_r025.nc
    5. run python make_bathy.py to generate bathy_meter.nc → bathy_r025.nc
    6. link bathy_meter.nc to bathy_r025.nc
    7. run python make_TScst_init.py to generate TScst_init_r025.nc (you do not need link, filename is written in namelist_cfg file)
    8. run python make_FLX.py to generate flx_r025.nc (you do not need link, filename is written in namelist_cfg file)

    NOTA BENE: if you change resolution (in python scripts) you need to change also in namelist_cfg:

    jpidta, jpjdta and jpkdta jpiglo and jpjglo '

    Example 1: Equatorial box without Coriolis

    1. TEST_EQ_BOX01_noCoriolis : Equatorial box (180°,60°) : with NO stratification / constant wind trades / NO Coriolis
    options values
    resolution 0.25° (¼)
    longitude lon_T=[ 0,180 ]
    latitude lat_T=[-31,31] (to have one point at 0° Equator)
    vertical levels 31
    batymetry 5000m
    forcing files utau=-0.1 / vtau=0 / qtot=0 / qsr=0 / emp=0
    NO stratification so=35 / thetao=10
    total time steps 3600 (=50days)
    output frequency 5 days

    output plot:https://forge.ipsl.jussieu.fr/nemo/attachment/wiki/WorkingGroups/Robustness-TestCases/ImplementationStrategyAndTemplate/current-temp-1S1N-all.pdf gphi0-current-temp-1S1N-all.pdf

    Example 2: Equatorial box e1, e2 constant values

    2. TEST_EQ_BOX01_cst : Equatorial box (180°,60°) : with NO stratification / constant wind trades / e1=e2=constant
    options values
    resolution 0.25° (¼)
    longitude lon_T=[ 0,180 ]
    latitude lat_T=[-31,31] (to have one point at 0° Equator)
    vertical levels 31
    batymetry 5000m
    forcing files utau=-0.1 / vtau=0 / qtot=0 / qsr=0 / emp=0
    NO stratification so=35 / thetao=10
    total time steps 3600(=50days)
    output frequency 5 days

    output plot: e1e2cst-current-temp-1S1N-all.pdf

    Example 3: Equatorial box with e1, e2 variables with latitudes

    3. TEST_EQ_BOX01_cst : Equatorial box (180°,60°) : with NO stratification / constant wind trades / e1, e2 variable with longitude
    options values
    resolution 0.25° (¼)
    longitude lon_T=[ 0,180 ]
    latitude lat_T=[-31,31] (to have one point at 0° Equator)
    vertical levels 31
    batymetry 5000m
    forcing files utau=-0.1 / vtau=0 / qtot=0 / qsr=0 / emp=0
    NO stratification so=35 / thetao=10
    total time steps 3600(=50days)
    output frequency 5 days

    output plot:https://forge.ipsl.jussieu.fr/nemo/attachment/wiki/WorkingGroups/Robustness-TestCases/ImplementationStrategyAndTemplate/current-temp-1S1N-all.pdf TScst-current-temp-all.pdf

    Example 4: Equatorial box with stratification and constant winds

    4. TEST_EQ_BOX01_cst : Equatorial box (180°,60°) : with stratification / constant wind trades
    options values
    resolution 0.25° (¼)
    longitude lon_T=[ 0,180 ]
    latitude lat_T=[-31,31] (to have one point at 0° Equator)
    vertical levels 31
    batymetry 5000m
    forcing files utau=-0.1 / vtau=0 / qtot=0 / qsr=0 / emp=0
    NO stratification so=35 / thetao=10
    total time steps 3600(=50days)
    output frequency 5 days

    output plot:https://forge.ipsl.jussieu.fr/nemo/attachment/wiki/WorkingGroups/Robustness-TestCases/ImplementationStrategyAndTemplate/current-temp-1S1N-all.pdf tprof-current-temp-1S1N-all.pdf

    Example 5: Equatorial box stratification and variable wind

    5. TEST_EQ_BOX01_cst : Equatorial box (180°,60°) : with stratification / variable wind (to see Kelvin waves)
    options values
    resolution 0.25° (¼)
    longitude lon_T=[ 0,180 ]
    latitude lat_T=[-31,31] (to have one point at 0° Equator)
    vertical levels 31
    batymetry 5000m
    forcing files utau=-0.1 / vtau=0 / qtot=0 / qsr=0 / emp=0
    NO stratification so=35 / thetao=10
    total time steps 3600 (restart from a 1801 time step, restart file from test 6) with winds = 0
    output frequency 5 days

    output plot:https://forge.ipsl.jussieu.fr/nemo/attachment/wiki/WorkingGroups/Robustness-TestCases/ImplementationStrategyAndTemplate/current-temp-1S1N-all.pdf Kelv-current-temp-1S1N-all.pdf

    Attachments (5)

    Download all attachments as: .zip