Changes between Version 1 and Version 2 of Ticket #681


Ignore:
Timestamp:
08/24/15 17:51:35 (9 years ago)
Author:
arango
Comment:

Legend:

Unmodified
Added
Removed
Modified

  • Ticket #681 – Description

    v1 v2  
    11The C-preprocessing option '''SPLINES''' have been in ROMS since the beginning.  It is used to activate conservative, parabolic splines reconstruction of vertical derivatives in the following routines:
    22
    3  * '''gls_corstep.F''': vertical velocity shear for Richardson Number
     3 * '''gls_corstep.F''': vertical velocity shear for Richardson Number (if option '''GLS_MIXING''')
    44
    5  * '''my25_corstep.F''': vertical velocity shear for Richardson Number
     5 * '''my25_corstep.F''': vertical velocity shear for Richardson Number (if option '''MY25_MIXING''')
    66
    7  * '''lmd_bkpp.F''': bulk Richardson Number
     7 * '''lmd_bkpp.F''': bulk Richardson Number (if options '''LMD_MIXING''' and '''LMD_BKPP''')
    88
    9  * '''lmd_skpp.F''': bulk Richardson Number
     9 * '''lmd_skpp.F''': bulk Richardson Number (if options '''LMD_MIXING''' and '''LMD_SKPP''')
    1010
    11  * '''lmd_vmix.F''': bulk Richardson Number
     11
     12 * '''lmd_vmix.F''': bulk Richardson Number (if option '''LMD_MIXING''')
     13
    1214
    1315 * '''step3d_t.F''': vertical diffusion
    … …  
    1517 * '''step3d_uv.F''': vertical viscosity
    1618
    17 There has been several messages in the ROMS forum discussing the problems when '''SPLINES''' is used to compute vertical viscosity and vertical diffusion. This was first brought to our attention by [https://www.myroms.org/forum/viewtopic.php?f=14&t=764&p=2419&hilit=logvinov#p2419 Evgeny Logvinov] when checking the stress condition ('''sustr = Akv * du/dz''' and '''svstr = Akv * dv/dz'''). He noticed that the stress condition was not holding and got a very sharp change in '''du/dz''' and '''dv/dz''' 1m below the surface.  Everything worked fine when '''SPLINES''' was not activated. Similar behavior was discussed [https://www.myroms.org/forum/viewtopic.php?f=33&t=3126 Chris Sherwood].
     19There has been several messages in the ROMS forum discussing the problems when '''SPLINES''' is used to compute vertical viscosity and vertical diffusion. This was first brought to our attention by [https://www.myroms.org/forum/viewtopic.php?f=14&t=764&p=2419&hilit=logvinov#p2419 Evgeny Logvinov] when checking the stress condition ('''sustr = Akv * du/dz''' and '''svstr = Akv * dv/dz'''). He noticed that the stress condition was not holding and he got a very sharp change in '''du/dz''' and '''dv/dz''' 1m below the surface.  Everything worked fine when '''SPLINES''' was '''not''' activated. Similar behavior was discussed [https://www.myroms.org/forum/viewtopic.php?f=33&t=3126 Chris Sherwood].
    1820
    19 I know several users that avoid using the '''SPLINES''' option because the of the surface and bottom stress condition needs to be preserved. However, it is nice is nice to have the parabolic splines reconstruction when computing either Richardson Number or Bulk Richardson Number. The vertical velocity shear (used as a squared quantity in the denominator) is a very noisy quantity so using parabolic splines yields a smoother values.
     21I know of several users that avoid using the '''SPLINES''' option because the of the surface and bottom stress condition needs to be preserved. However, it is nice is nice to have the parabolic splines reconstruction when computing either Richardson Number or Bulk Richardson Number. The vertical velocity shear (used as a squared quantity in the denominator) is a very noisy quantity so using parabolic splines yields a smoother values.
    2022
    2123Therefore, the C-preprocessing option '''SPLINES''' is removed and replaced with following three options for more flexibility:
    … …  
    2931----
    3032
    31 '''WARNING''': ROMS will issue an error message if the option SPLINES is activated after this update.  The user will need to activate the '''SPLINES_VVISC''' and '''SPLINES_VDIFF'''.  Also you will need to activate '''RI_SPLINES''' if using '''GLS_MIXING''', '''LMD_MIXING''', '''MY25_MIXING'''.  Therefore, users need to update the header file for their applications.
     33'''WARNING''': ROMS will issue an error message if the option '''SPLINES''' is activated after this update.  The user will need to activate the '''SPLINES_VVISC''' and '''SPLINES_VDIFF'''.  Also you will need to activate '''RI_SPLINES''' if using '''GLS_MIXING''', '''LMD_MIXING''', '''MY25_MIXING'''.  Therefore, users need to update the header file for their applications.
    3234
    33 Several header files were update in the code repository and the test cases repository.  I highly recommend users to update their repositories.
     35Several header files were updated in the '''trunk''' repository and '''test''' repositories.  I highly recommend users to update their repositories.
    3436
    3537----