Version 7 (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 show what appens when in equatorial channel we remove Coriolis force

    show what appens when we remove 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/TEST_EQ_BOX01

    and 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)

    You choose in namelist_cfg cn_exp name and nitend.


    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

    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

    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

    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

    Example 5: Equatorial box stratification and constant wind profile

    5. TEST_EQ_BOX01_cst : Equatorial box (180°,60°) : with stratification / constant winds (profile)
    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

    Example 6: Equatorial box stratification and variable wind

    6. 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:current-temp-1S1N-all.pdf

    Attachments (5)

    Download all attachments as: .zip