I am attempting to run a coupled WRF/ROMS model simulation with open boundary conditions, but I am experiencing problems initializing ROMS, not related to WRF. Currently, I have ROMS fully compiled and I have already generated ROMS model grid, initialization, climatology, and boundary files. When I attempt to run the coupled system WRF initializes correctly, but ROMS fails because it is missing a file ocean_frc.nc. If I assume closed-boundary conditions, all my input files feed into ROMS correctly, but the run fails after 3 hours due to model instability associated with the closed boundaries. Is there any way to modify my configuration to allow for open boundary conditions yet work around this forcing condition file? Otherwise, if I must generate it, are there tools available when I can generate these files and from where may I find the need data?
Thanks
roms.error output file
Coupled Input File name = Projects/NOREASTER/coupling_coawst.in
Model Coupling:
Ocean Model MPI nodes: 000 - 007
Atmos Model MPI nodes: 008 - 059
Process Information:
Node # 0 (pid= 20114) is active.
Node # 3 (pid= 20118) is active.
Node # 2 (pid= 20113) is active.
Model Input Parameters: ROMS/TOMS version 3.4
Wednesday - June 6, 2012 - 6:35:23 PM
-----------------------------------------------------------------------------
Node # 1 (pid= 20108) is active.
Node # 5 (pid= 20117) is active.
Node # 4 (pid= 20119) is active.
Node # 7 (pid= 20116) is active.
Node # 6 (pid= 20115) is active.
Idealized Coastal shelf
Operating system : Linux
CPU/hardware : x86_64
Compiler system : ifort
Compiler command : /usr/local/intel/mpi/3.2.2.006/bin64/mpif90
Compiler flags : -heap-arrays -fp-model precise -m64 -ip -O3 -xW -I/usr/local/other/SLES11/MCT/include -I/usr/local/other/SLES11/MCT/include -free
Input Script : Projects/NOREASTER/roms.namelist.in
SVN Root URL : https://www.myroms.org/svn/omlab/branches/jcwarner
SVN Revision : 583M
Local Root : /gpfsm/dnb32/snicholl/WRF_COAWST/WRF_ROMS_open/COAWST
Header Dir : /discover/nobackup/snicholl/WRF_COAWST/WRF_ROMS_open/COAWST/Projects/NOREASTER
Header file : noreaster.h
Analytical Dir: /discover/nobackup/snicholl/WRF_COAWST/WRF_ROMS_open/COAWST/Projects/NOREASTER
Resolution, Grid 01: 0158x0338x016, Parallel Nodes: 8, Tiling: 002x004
Physical Parameters, Grid: 01
=============================
64800 ntimes Number of timesteps for 3-D equations.
10.000 dt Timestep size (s) for 3-D equations.
30 ndtfast Number of timesteps for 2-D equations between
each 3D timestep.
600.000 TI_ATM_OCN Time interval (s) between coupling ATM-OCN models.
60 nOCN_ATM Number of OCN timesteps between coupling to ATM.
1 WRF_CPLG_GRID WRF grid number that roms/swan will couple to.
1 ERstr Starting ensemble/perturbation run number.
1 ERend Ending ensemble/perturbation run number.
0 nrrec Number of restart records to read from disk.
T LcycleRST Switch to recycle time-records in restart file.
2160 nRST Number of timesteps between the writing of data
into restart fields.
1 ninfo Number of timesteps between print of information
to standard output.
T ldefout Switch to create a new output NetCDF file(s).
2160 nHIS Number of timesteps between the writing fields
into history file.
1 ntsAVG Starting timestep for the accumulation of output
time-averaged data.
2160 nAVG Number of timesteps between the writing of
time-averaged data into averages file.
1 ntsDIA Starting timestep for the accumulation of output
time-averaged diagnostics data.
2160 nDIA Number of timesteps between the writing of
time-averaged data into diagnostics file.
0.0000E+00 nl_tnu2(01) NLM Horizontal, harmonic mixing coefficient
(m2/s) for tracer 01: temp
0.0000E+00 nl_tnu2(02) NLM Horizontal, harmonic mixing coefficient
(m2/s) for tracer 02: salt
5.0000E-02 nl_visc2 NLM Horizontal, harmonic mixing coefficient
(m2/s) for momentum.
1.0000E-06 Akt_bak(01) Background vertical mixing coefficient (m2/s)
for tracer 01: temp
1.0000E-06 Akt_bak(02) Background vertical mixing coefficient (m2/s)
for tracer 02: salt
1.0000E-05 Akv_bak Background vertical mixing coefficient (m2/s)
for momentum.
5.0000E-06 Akk_bak Background vertical mixing coefficient (m2/s)
for turbulent energy.
5.0000E-06 Akp_bak Background vertical mixing coefficient (m2/s)
for turbulent generic statistical field.
3.000 gls_p GLS stability exponent.
1.500 gls_m GLS turbulent kinetic energy exponent.
-1.000 gls_n GLS turbulent length scale exponent.
7.6000E-06 gls_Kmin GLS minimum value of turbulent kinetic energy.
1.0000E-12 gls_Pmin GLS minimum value of dissipation.
5.4770E-01 gls_cmu0 GLS stability coefficient.
1.4400E+00 gls_c1 GLS shear production coefficient.
1.9200E+00 gls_c2 GLS dissipation coefficient.
-4.0000E-01 gls_c3m GLS stable buoyancy production coefficient.
1.0000E+00 gls_c3p GLS unstable buoyancy production coefficient.
1.0000E+00 gls_sigk GLS constant Schmidt number for TKE.
1.3000E+00 gls_sigp GLS constant Schmidt number for PSI.
1400.000 charnok_alpha Charnok factor for Zos calculation.
0.500 zos_hsig_alpha Factor for Zos calculation using Hsig(Awave).
0.250 sz_alpha Factor for Wave dissipation surface tke flux .
100.000 crgban_cw Factor for Craig/Banner surface tke flux.
3.0000E-04 rdrg Linear bottom drag coefficient (m/s).
3.0000E-03 rdrg2 Quadratic bottom drag coefficient.
2.0000E-02 Zob Bottom roughness (m).
2.0000E-02 Zos Surface roughness (m).
1 Vtransform S-coordinate transformation equation.
1 Vstretching S-coordinate stretching function.
5.0000E+00 theta_s S-coordinate surface control parameter.
4.0000E-01 theta_b S-coordinate bottom control parameter.
50.000 Tcline S-coordinate surface/bottom layer width (m) used
in vertical coordinate stretching.
1025.000 rho0 Mean density (kg/m3) for Boussinesq approximation.
54058.000 dstart Time-stamp assigned to model initialization (days).
0.000 tide_start Reference time origin for tidal forcing (days).
20061119.50 time_ref Reference time for units attribute (yyyymmdd.dd)
0.0000E+00 Tnudg(01) Nudging/relaxation time scale (days)
for tracer 01: temp
0.0000E+00 Tnudg(02) Nudging/relaxation time scale (days)
for tracer 02: salt
0.0000E+00 Znudg Nudging/relaxation time scale (days)
for free-surface.
0.0000E+00 M2nudg Nudging/relaxation time scale (days)
for 2D momentum.
0.0000E+00 M3nudg Nudging/relaxation time scale (days)
for 3D momentum.
0.0000E+00 obcfac Factor between passive and active
open boundary conditions.
18.000 T0 Background potential temperature (C) constant.
35.000 S0 Background salinity (PSU) constant.
1.000 gamma2 Slipperiness variable: free-slip (1.0) or
no-slip (-1.0).
T Hout(idFsur) Write out free-surface.
T Hout(idUbar) Write out 2D U-momentum component.
T Hout(idVbar) Write out 2D V-momentum component.
T Hout(idUvel) Write out 3D U-momentum component.
T Hout(idVvel) Write out 3D V-momentum component.
T Hout(idWvel) Write out W-momentum component.
T Hout(idTvar) Write out tracer 01: temp
T Hout(idUsms) Write out surface U-momentum stress.
T Hout(idVsms) Write out surface V-momentum stress.
T Hout(idTsur) Write out surface net heat flux.
T Hout(idTsur) Write out surface net salt flux.
T Hout(idSrad) Write out shortwave radiation flux.
Output/Input Files:
Output Restart File: ocean_rst.nc
Output History File: ocean_noreaster_his.nc
Output Averages File: ocean_avg.nc
Output Diagnostics File: ocean_dia.nc
Physical parameters File: Projects/NOREASTER/roms.namelist.in
Input Grid File: Projects/NOREASTER/ROMS_coawst_grd_final.nc
Input Nonlinear Initial File: Projects/NOREASTER/ROMS_coawst_ini_final.nc
Tile partition information for Grid 01: 0158x0338x0016 tiling: 002x004
tile Istr Iend Jstr Jend Npts
0 1 79 1 84 106176
1 80 158 1 84 106176
2 1 79 85 169 107440
3 80 158 85 169 107440
4 1 79 170 254 107440
5 80 158 170 254 107440
6 1 79 255 338 106176
7 80 158 255 338 106176
Tile minimum and maximum fractional grid coordinates:
(interior points only)
tile Xmin Xmax Ymin Ymax grid
0 0.50 79.50 0.50 84.50 RHO-points
1 79.50 158.50 0.50 84.50 RHO-points
2 0.50 79.50 84.50 169.50 RHO-points
3 79.50 158.50 84.50 169.50 RHO-points
4 0.50 79.50 169.50 254.50 RHO-points
5 79.50 158.50 169.50 254.50 RHO-points
6 0.50 79.50 254.50 338.50 RHO-points
7 79.50 158.50 254.50 338.50 RHO-points
0 1.00 79.50 0.50 84.50 U-points
1 79.50 158.00 0.50 84.50 U-points
2 1.00 79.50 84.50 169.50 U-points
3 79.50 158.00 84.50 169.50 U-points
4 1.00 79.50 169.50 254.50 U-points
5 79.50 158.00 169.50 254.50 U-points
6 1.00 79.50 254.50 338.50 U-points
7 79.50 158.00 254.50 338.50 U-points
0 0.50 79.50 1.00 84.50 V-points
1 79.50 158.50 1.00 84.50 V-points
2 0.50 79.50 84.50 169.50 V-points
3 79.50 158.50 84.50 169.50 V-points
4 0.50 79.50 169.50 254.50 V-points
5 79.50 158.50 169.50 254.50 V-points
6 0.50 79.50 254.50 338.00 V-points
7 79.50 158.50 254.50 338.00 V-points
Maximum halo size in XI and ETA directions:
HaloSizeI(1) = 182
HaloSizeJ(1) = 194
TileSide(1) = 91
TileSize(1) = 7735
Activated C-preprocessing Options:
NOREASTER Idealized Coastal shelf
ADD_FSOBC Adding tidal elevation to proccesed OBC data.
ADD_M2OBC Adding tidal currents to proccesed OBC data.
ATM_PRESS Impose atmospheric pressure onto sea surface.
ANA_BSFLUX Analytical kinematic bottom salinity flux.
ANA_BTFLUX Analytical kinematic bottom temperature flux.
ANA_FSOBC Analytical free-surface boundary conditions.
ANA_M2OBC Analytical 2D momentum boundary conditions.
ANA_SSFLUX Analytical kinematic surface salinity flux.
ASSUMED_SHAPE Using assumed-shape arrays.
ATM2OCN_FLUXES Surface fluxes from ATM Model.
AVERAGES Writing out time-averaged fields.
DIAGNOSTICS_TS Computing and writing tracer diagnostic terms.
DIAGNOSTICS_UV Computing and writing momentum diagnostic terms.
DJ_GRADPS Parabolic Splines density Jacobian (Shchepetkin, 2002).
DOUBLE_PRECISION Double precision arithmetic.
EAST_FSCHAPMAN Eastern edge, free-surface, Chapman condition.
EAST_M2FLATHER Eastern edge, 2D momentum, Flather condition.
EAST_M3RADIATION Eastern edge, 3D momentum, radiation condition.
EAST_TRADIATION Eastern edge, tracers, radiation condition.
GLS_MIXING Generic Length-Scale turbulence closure.
KANTHA_CLAYSON Kantha and Clayson stability function formulation.
MASKING Land/Sea masking.
MCT_LIB Using Model Coupling Toolkit library.
MIX_GEO_TS Mixing of tracers along geopotential surfaces.
MIX_S_UV Mixing of momentum along constant S-surfaces.
MPI MPI distributed-memory configuration.
NONLINEAR Nonlinear Model.
NONLIN_EOS Nonlinear Equation of State for seawater.
NORTHERN_WALL Wall boundary at Northern edge.
N2S2_HORAVG Horizontal smoothing of buoyancy and shear.
POWER_LAW Power-law shape time-averaging barotropic filter.
PROFILE Time profiling activated .
K_GSCHEME Third-order upstream advection of TKE fields.
RADIATION_2D Use tangential phase speed in radiation conditions.
RAMP_TIDES Ramping tidal forcing for one day.
!RST_SINGLE Double precision fields in restart NetCDF file.
SALINITY Using salinity.
SOLVE3D Solving 3D Primitive Equations.
SOUTH_FSCHAPMAN Southern edge, free-surface, Chapman condition.
SOUTH_M2FLATHER Southern edge, 2D momentum, Flather condition.
SOUTH_M3RADIATION Southern edge, 3D momentum, radiation condition.
SOUTH_TRADIATION Southern edge, tracers, radiation condition.
SPLINES Conservative parabolic spline reconstruction.
SSH_TIDES Add tidal elevation to SSH climatology.
TS_U3HADVECTION Third-order upstream horizontal advection of tracers.
TS_C4VADVECTION Fourth-order centered vertical advection of tracers.
TS_DIF2 Harmonic mixing of tracers.
UV_ADV Advection of momentum.
UV_COR Coriolis term.
UV_U3HADVECTION Third-order upstream horizontal advection of 3D momentum.
UV_C4VADVECTION Fourth-order centered vertical advection of momentum.
UV_LOGDRAG Logarithmic bottom stress.
UV_TIDES Add tidal currents to 2D momentum climatologies.
UV_VIS2 Harmonic mixing of momentum.
VAR_RHO_2D Variable density barotropic mode.
WESTERN_WALL Wall boundary at Western edge.
WRF_COUPLING WRF model coupling.
NETCDF_OPEN - unable to open existing NetCDF file:
ocean_frc.nc
call from:
rank 6 in job 1 borg01v021_48088 caused collective abort of all ranks
exit status of rank 6: killed by signal 9