| 46 | | Consequently, the header file '''ana_hmixcoef.h''' is deprecated and replaced with '''ana_sponge.h'''. So if you were used this capability before in older ROMS versions, you need to code your sponge areas in the new header file '''ana_sponge.h'''. This is a much simpler file. Notice that the routine '''ini_hmixcoef''' is called before than '''ana_sponge''' is called. Therefore, users have two options: ('''1''') completely overwrite the horizontal mixing coefficient assigned earlier in '''ini_hmixcoef''', or ('''2''') only increase horizontal mixing in the sponge areas. The sponge areas in '''ana_sponge''' can be coded as a non-dimensional factor as it is done in '''ini_hmixcoef''' with '''visc_factor''' and '''diff_factor'''. Or a complete overwrite, which is not recommended if either '''VISC_GRID''' and '''DIFF_GRID''' are activated. |
| | 48 | Consequently, the header file '''ana_hmixcoef.h''' is deprecated and replaced with '''ana_sponge.h'''. So if you were using this capability before with older ROMS versions, you need to code your sponge areas in the new header file '''ana_sponge.h'''. This is a much simpler file. Notice that the routine '''ini_hmixcoef''' is called before than '''ana_sponge''' is called. Therefore, users have two options: ('''1''') completely overwrite the horizontal mixing coefficient assigned earlier in '''ini_hmixcoef''', or ('''2''') only increase horizontal mixing in the sponge areas. The sponge areas in '''ana_sponge''' can be coded as a non-dimensional factor as it is done in '''ini_hmixcoef''' with '''visc_factor''' and '''diff_factor'''. Or a complete overwrite, which is not recommended if either '''VISC_GRID''' and '''DIFF_GRID''' are activated. |
| 146 | | A new routine '''get_nudgcoef.F''' is added to read these nudging inverse time scales. This routine will check the '''units''' attribute to convert the scales to 1/second. If the nudging scales for a specific tracer are available (say '''salt_NudgCoef''') it will read that NetCDF variable. If '''not and''' the generic scales are available ('''tracer_NudgeCoef'''), it will process those values instead. This strategy will give a lot of flexibility when setting nudging for a particular application. The generic '''tracer_NudgeCoef''' variable is useful when nudging passive (biology and sediment) tracers. |
| | 148 | A template script '''d_nudgcoef.m''' is provided in the matlab/initial repository to create this new NetCDF file. |
| | 149 | |
| | 150 | A new routine '''get_nudgcoef.F''' is added to read these nudging inverse time scales. This routine will check the '''units''' attribute to convert the scales to 1/second. If the nudging scales for a specific tracer are available (say '''salt_NudgCoef''') it will read that NetCDF variable. If '''not and''' the generic scales are available ('''tracer_NudgeCoef'''), it will process those values instead. This strategy give us a lot of flexibility when setting nudging for a particular application. The generic '''tracer_NudgeCoef''' variable is useful when nudging passive (biology and sediment) tracers. |
| | 151 | |
| | 152 | ''' WARNING:''' The header file '''ana_nudgcoef.h''' was modified to include depth dependency in the nudging inverse time scales. The routine is simpler now. However, please use the NetCDF file instead to avoid parallel coding errors. Recall that you can plot and fine tune the variables outside of ROMS. |
| | 153 | |
| | 154 | If passive/active radiation open boundary conditions and activated climatology nudging, the inverse time scales are used directly in the open boundary conditions routines instead of the uniform variables '''*obc_out''' and '''*obc_in'''. Notice that this logic was removed from '''ana_nudgcoef.h'''. |
| | 155 | |
| | 156 | ---- |
| | 157 | |
| | 158 | Again, please use the NetCDF options discussed above instead of coding with analytical functions. We keep these analytical functions as part of the legacy code. Coding such routines require extensive parallel expertise. So if you are in doubt, set a NetCDF file instead. It is much easier. I will ignore forum messages asking for help or guidance how to code such analytical functions. |
| | 159 | |
| | 160 | This constitutes the final piece in the nesting algorithms. I have seen some nesting inquires in the forum. Nesting requires expertise to set-up and a lot of patience. There are a lot things that you need to think when setting a nesting application: dynamical regimes, atmospheric forcing, volume/mass conservation, grid topology, bathymetry, land/sea masking, and so on. You just cannot put a nesting grid anywhere! Many users think about nesting as a grid generation problem, but it is much more than that. We cannot provide you with a cooking recipe for nesting so you just need to get your hands dirty and try various strategies until you find one that works for your particular application. |
