160 likes | 232 Views
This meeting, held on August 23, 2001 at NCAR, focused on comparing numerical elements in ocean models to enhance code development and user choices. The goal was to evaluate different schemes and parameterizations in POM and ROMS/TOMS based on criteria such as accuracy, stability, and computational cost.
E N D
TOMS Inaugural Meeting, NCAR, Aug. 23, 2001 Model Comparison and Evaluation Studies Tal Ezer (with H. Arango & A. Shchepetkin)
Goal: to evaluate and test new numerical elements in terrain-following ocean models. • Why: to help developers in improving codes and users in selecting options and parameterizations. • How: by comparing different schemes and parameterizations in POM and ROMS/TOMS. • Criteria: numerical errors, accuracy, stability, compatibility between schemes, computational cost.
Tested Elements • Advection schemes • Time stepping algorithms • Pressure gradient schemes • Others
Seamount configuration test Very steep case h=4050m, w=50km s=0.36, r=14.2 moderately steep h=2700m, w=100km s=0.07, r=2.7 S=max(DH/2H) r=max(sDH/HDs) grid=(64x64x20), Dx=8km
Zonal flow -50cm 0 +50cm Topography Sea surface height
Effect of advection scheme on model: Surface elevation anomaly ADV2- 2nd ord. cnt. POM ROMS ADV4- 4th ord. cnt. ADV3- 3rd ord. upst. ROMS ROMS
Time-stepping schemes(split mode: baroclinic/internal and barotropic/external)
Coupling of barotropic (external) and baroclinic (internal) modes in ROMS weights Un+1 = SamUm Un+½ = S bm’Um’ DTE 1<m<N, N=DTI/DTE DTI
Sensitivity to internal (DTI) & external (DTE) time steps ROMS UNSTABLE STABLE TDI/DTE POM CFL=13s
The adjustment process in POM and ROMS: forced case (zonal flow)
Roms- sensitivity of adjustment procesto time step choices DTE=12s DTE=24s DTI=360s DTI=720s
Structure ofV (cm/s) in ROMS for different PG schemes (medium seamount case) R-DJ (Vmax=3.7) R-PJQ (Vmax=0.03) R-FPJ (Vmax=30) R-DJC (Vmax=0.06) R-WDJ (Vmax=0.3)
More model elements and schemes need to be tested • Boundary conditions • Data assimilation • Nesting and air-sea coupling • Horizontal and vertical mixing schemes • Generalized coordinate system: comparisons between sigma, z-level & isopycnal grids.