Custom Query (986 matches)
Results (106 - 108 of 986)
| Ticket | Owner | Reporter | Resolution | Summary |
|---|---|---|---|---|
| #144 | Done | Updated model coupling using the MCT library | ||
| Description |
The code was updated to consolidate the model coupling with either the ESMF or MCT libraries. Both libraries interfaces have the same design. They use input script coupling.in and metadata file coupling.dat for model coupling. The design is generic so it is easier to add new fields to interchange them between models. I cleaned all the files that use the MCT library and I removed the Lib/MCT_WRF directory. It is now obsolete. I also built the drivers for ROMS/WRF/SWAN coupling. However, the WRF model is not released a this time since it requires further development. It was loaded to my branch in the omlab repository which have limited access. At this time all the coupled models are in the same grid. However, the roms_import.F and roms_export.F routines can be expanded in the future to include the spatial interpolation between coupled models with different grids. I also updated the way that time is reported in the standard output. Now, it is reported as string having the following form: day hh:mm:ss. This allows more detailed information that just a float number. It is also more accurate. Many thanks to John Warner for his work and great help. |
|||
| #146 | Done | Updated SWAN model to version 40.51 | ||
| Description |
Updated SWAN model to the latest released version (40.51) taken from their website. The patches 40.51.A and 40.51.B were applied. Our modifications to SWAN were redesigned so it is easy to update this model in the future. They are now done with c-preprocessing. This allows to run SWAN as originally distributed by just not defining MODEL_COUPLING. There are few modifications to SWAN kernel so it is not possible to get identical solutions when compared with previous version. However, these differences are hard to see by eye during plotting. Added two new directories SWAN/Programs and SWAN/Docs. Also, few bugs were corrected in the coupling with the MCT library. |
|||
| #149 | Done | Time dependent horizontal diffusion/viscosity coefficients | ||
| Description |
The horizontal diffusion (t3dmix2_s.h, t3dmix4_s.h, t3dmix2_geo.h, t3dmix4_geo.h, t3dmix2_iso.h, t3dmix4_iso.h) and viscosity (uv3dmix2_s.h, uv3dmix4_s.h, uv3dmix2_geo.h, uv3dmix4_geo.h) routines were updated to include time-dependent mixing coefficients. Few new CPP options are introduced: TS_U3ADV_SPLIT use if 3rd-order upstream split tracer advection TS_SMAGORINSKY use to turn ON or OFF Smagorinsky-like diffusion UV_U3ADV_SPLIT use if 3rd-order upstream split momentum advection UV_SMAGORINSKY use to turn ON or OFF Smagorinsky-like diffusion The primary goal of this update is to correct for the spurious numerical diapycnal diffusion of the tracer advection in coarse basin scale applications. Since this requires a velocity-dependent 3D horizontal diffusion coefficient, a Smagorinsky-like mixing coefficient is also coded. However, its not recommended to use this option on tracers since it is not clear, physically, how to compute the Smagorinsky tracer diffusion from momentum viscosity, as discussed by Griffies and Hallberg (2000). It is added for pedagogical reasons. The Smagorisky mixing on momentum makes more sense. The 3rd-order upstream split advection (TS_U3ADV_SPLIT) can be used to correct for the spurious diapycnal diffusion of the 3rd-order upstream advection operator in terrain-following coordinates. The advection operator is split between advective and diffusive components, as shown by Holland et al. (1998), Webb et al. (1998), and Marchesiello et al. (2008). They show that the split operator has a 4th-order centered differences term plus a 3rd-order hyper-diffusion to correct for spurious diapycnal mixing. The total (time dependent plus time independent) horizontal mixing coefficient are computed in new routine hmixing.F. Many thanks to Patrick Marchesiello for providing his codes from which hmixing.F is based. Notice that several internal cpp options are activated in globaldefs.h to facilitate the split of the advection operator: #ifdef TS_U3ADV_SPLIT # define DIFF_3DCOEF # ifdef TS_U3HADVECTION # undef TS_U3HADVECTION # endif # ifndef TS_C4HADVECTION # define TS_C4HADVECTION # endif # ifndef TS_C4VADVECTION # define TS_C4VADVECTION # endif # ifndef TS_DIF4 # define TS_DIF4 # endif # ifdef TS_DIF2 # undef TS_DIF2 # endif # ifdef TS_SMAGORINSKY # undef TS_SMAGORINSKY # endif #endif It is also recommended to use the rotated mixing tensor along geopotentials (MIX_GEO_TS) for the biharmonic operator. The modifications to the tangent linear, representer, and adjoint models will be updated soon. |
|||
