READ_PHYPAR - could not find input file: ocean_frc.nc

[thomas.roc]

Hi everybody,

I am trying to run a very simple case (a simple rectangular channel with a constant inflow). My simulation is compiling without any trouble but when I am running the executable this message appears: "READ_PHYPAR - could not find input file: ocean_frc.nc"

I guess this is due to a consistency problem in my initial and analytical forcing because ROMS is asking me for ocean_frc.nc, yet I really cannot see which on is missing.

More surpringly, I did not have this problem before to update my ROMS trunk.

Does anybody have a clue how to figure out what is missing ?

Cheers.

FYU, my list of headers:
#define SOLVE3D
#define UV_ADV
#define UV_VIS2
#define MIX_S_UV
#define GLS_MIXING
#define CANUTO_B
#define UV_LOGDRAG

#define ANA_GRID
#define ANA_INITIAL
#define ANA_FSOBC
#define ANA_SMFLUx
#define ANA_STFLUX
#define ANA_BTFLUX
#define ANA_M3OBC
#define ANA_M2OBC

and the error message:
Operating system : Linux
CPU/hardware : x86_64
Compiler system : gfortran
Compiler command : /usr/bin/mpif90
Compiler flags : -frepack-arrays -O3 -ffast-math -ffree-form -ffree-line-length-none -ffree-form -ffree-line-length-none

Input Script : ./ocean_background_case.in

SVN Root URL : https://www.myroms.org/svn/src/trunk
SVN Revision : 634M

Local Root : /home/thomas/Desktop/ROMS_1
Header Dir : /home/thomas/Desktop/ROMS_1/ROMS/Include
Header file : background_case.h
Analytical Dir: /home/thomas/Desktop/ROMS_1/ROMS/Functionals

Resolution, Grid 01: 0200x0200x030, Parallel Nodes: 25, Tiling: 005x005


Physical Parameters, Grid: 01
=============================

1200 ntimes Number of timesteps for 3-D equations.
0.150 dt Timestep size (s) for 3-D equations.
30 ndtfast Number of timesteps for 2-D equations between
each 3D timestep.
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.
1200 nRST Number of timesteps between the writing of data
into restart fields.
120 ninfo Number of timesteps between print of information
to standard output.
T ldefout Switch to create a new output NetCDF file(s).
120 nHIS Number of timesteps between the writing fields
into history file.
1.9590E+01 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.
-1.000 gls_p GLS stability exponent.
0.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.
5.5500E-01 gls_c1 GLS shear production coefficient.
8.3300E-01 gls_c2 GLS dissipation coefficient.
-6.0000E-01 gls_c3m GLS stable buoyancy production coefficient.
1.0000E+00 gls_c3p GLS unstable buoyancy production coefficient.
2.0000E+00 gls_sigk GLS constant Schmidt number for TKE.
2.0000E+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.
0.0000E+00 rdrg Linear bottom drag coefficient (m/s).
0.0000E+00 rdrg2 Quadratic bottom drag coefficient.
0.0000E+00 Zob Bottom roughness (m).
0.0000E+00 Zos Surface roughness (m).
2 Vtransform S-coordinate transformation equation.
3 Vstretching S-coordinate stretching function.
6.0000E-01 theta_s S-coordinate surface control parameter.
6.0000E-01 theta_b S-coordinate bottom control parameter.
0.000 Tcline S-coordinate surface/bottom layer width (m) used
in vertical coordinate stretching.
1025.000 rho0 Mean density (kg/m3) for Boussinesq approximation.
0.000 dstart Time-stamp assigned to model initialization (days).
0.00 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
1.0000E-03 Znudg Nudging/relaxation time scale (days)
for free-surface.
1.0000E-03 M2nudg Nudging/relaxation time scale (days)
for 2D momentum.
1.0000E-03 M3nudg Nudging/relaxation time scale (days)
for 3D momentum.
5.0000E+00 obcfac Factor between passive and active
open boundary conditions.
F VolCons(1) NLM western edge boundary volume conservation.
F VolCons(2) NLM southern edge boundary volume conservation.
F VolCons(3) NLM eastern edge boundary volume conservation.
F VolCons(4) NLM northern edge boundary volume conservation.
20.000 T0 Background potential temperature (C) constant.
35.000 S0 Background salinity (PSU) constant.
1025.000 R0 Background density (kg/m3) used in linear Equation
of State.
1.7000E-04 Tcoef Thermal expansion coefficient (1/Celsius).
7.6000E-04 Scoef Saline contraction coefficient (1/PSU).
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(idOvel) Write out omega vertical velocity.
T Hout(idVvis) Write out vertical viscosity: AKv.
T Hout(idMtke) Write out turbulent kinetic energy.
T Hout(idMtls) Write out turbulent generic length-scale.

Output/Input Files:

Output Restart File: ocean_rst.nc
Output History File: ocean_his.nc

READ_PHYPAR - could not find input file: ocean_frc.nc

Elapsed CPU time (seconds):


ROMS/TOMS - Output NetCDF summary for Grid 01:

ROMS/TOMS - I/O error ............... exit_flag: 4

[kate]

#define ANA_SMFLUx

Did you cut and paste this? If so, make sure the X is capitalized in your file and try again.

[thomas.roc]

Yeah I spotted this one as well however the problem doesn't come from that...but thank you anyway.

I have actually cracked my problem though. In my input file (i.e. *.in) the open-boundary conditions for temperature and salinity were set to "Clamped" which imply to provide some values, analytic or data, for T and S. Switching to "Radiation" conditions solved my problem.

Cheers,

-T-