Changes between Version 1 and Version 2 of Ticket #681
- Timestamp:
- 08/24/15 17:51:35 (9 years ago)
Legend:
- Unmodified
- Added
- Removed
- Modified
-
Ticket #681 – Description
v1 v2 1 1 The 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: 2 2 3 * '''gls_corstep.F''': vertical velocity shear for Richardson Number 3 * '''gls_corstep.F''': vertical velocity shear for Richardson Number (if option '''GLS_MIXING''') 4 4 5 * '''my25_corstep.F''': vertical velocity shear for Richardson Number 5 * '''my25_corstep.F''': vertical velocity shear for Richardson Number (if option '''MY25_MIXING''') 6 6 7 * '''lmd_bkpp.F''': bulk Richardson Number 7 * '''lmd_bkpp.F''': bulk Richardson Number (if options '''LMD_MIXING''' and '''LMD_BKPP''') 8 8 9 * '''lmd_skpp.F''': bulk Richardson Number 9 * '''lmd_skpp.F''': bulk Richardson Number (if options '''LMD_MIXING''' and '''LMD_SKPP''') 10 10 11 * '''lmd_vmix.F''': bulk Richardson Number 11 12 * '''lmd_vmix.F''': bulk Richardson Number (if option '''LMD_MIXING''') 13 12 14 13 15 * '''step3d_t.F''': vertical diffusion … … 15 17 * '''step3d_uv.F''': vertical viscosity 16 18 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 notactivated. Similar behavior was discussed [https://www.myroms.org/forum/viewtopic.php?f=33&t=3126 Chris Sherwood].19 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 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]. 18 20 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.21 I 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. 20 22 21 23 Therefore, the C-preprocessing option '''SPLINES''' is removed and replaced with following three options for more flexibility: … … 29 31 ---- 30 32 31 '''WARNING''': ROMS will issue an error message if the option SPLINESis 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. 32 34 33 Several header files were update in the code repository and the test cases repository. I highly recommend users to update their repositories.35 Several header files were updated in the '''trunk''' repository and '''test''' repositories. I highly recommend users to update their repositories. 34 36 35 37 ----