I am trying to simulate tide by 2D-ROMS.
I only retain sustr and sms_time in the frc.nc. At the same time, tide is added to the frc.nc. ANA_INITIAL and ANA_BMFLUX is defined.Compiling is ok, but something is wrong when running

Can anyone help me?
Output info:
INITIAL: Configuring and initializing forward nonlinear model ...
Minimum X-grid spacing, DXmin = 2.75248895E+00 km
Maximum X-grid spacing, DXmax = 3.20812877E+00 km
Minimum Y-grid spacing, DYmin = 2.75409673E+00 km
Maximum Y-grid spacing, DYmax = 3.20766198E+00 km
Minimum barotropic Courant Number = 1.85253133E-01
Maximum barotropic Courant Number = 9.71708342E-01
Maximum Coriolis Courant Number = 5.84041068E-03
GET_NGFLD - tidal period
(Min = 4.30819210E+04 Max = 2.38071599E+06)
GET_2DFLD - tidal elevation amplitude
(Min = 3.26393686E-05 Max = 6.68537894E+00)
GET_2DFLD - tidal elevation phase angle
(Min = 2.41898787E-06 Max = 6.28318020E+00)
GET_2DFLD - tidal current inclination angle
(Min = 4.97589188E-07 Max = 6.28318253E+00)
GET_2DFLD - tidal current phase angle
(Min = 2.59098604E-06 Max = 6.28318291E+00)
GET_2DFLD - maximum tidal current, ellipse major axis
(Min = 1.27754016E-06 Max = 1.42023097E+00)
GET_2DFLD - minimum tidal current, ellipse minor axis
(Min = -5.23102760E-01 Max = 4.38586697E-01)
GET_2DFLD - surface u-momentum stress, t = 345 00:00:00
(Rec=0000012, Index=1, File: barotropic_tide_frc.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_2DFLD - surface v-momentum stress, t = 345 00:00:00
(Rec=0000012, Index=1, File: barotropic_tide_frc.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
INQUIRE - unable to assign file counter, Fcount = 0
while processing structure:
and variable; zeta_east
NETCDF_OPEN - unable to open existing NetCDF file:
call from: inquire.FQUIRE - unable to assign file counter, Fcount = 0
while processing structure:
and variable; zeta_east
NETCDF_OPEN - unable to open existing NetCDF file:
call from: inquire.F
INITIAL: Configuring and initializing forward nonlinear model ...
Minimum X-grid spacing, DXmin = 2.75248895E+00 km
Maximum X-grid spacing, DXmax = 3.20812877E+00 km
Minimum Y-grid spacing, DYmin = 2.75409673E+00 km
Maximum Y-grid spacing, DYmax = 3.20766198E+00 km
Minimum barotropic Courant Number = 1.85253133E-01
Maximum barotropic Courant Number = 9.71708342E-01
Maximum Coriolis Courant Number = 5.84041068E-03
GET_NGFLD - tidal period
(Min = 4.30819210E+04 Max = 2.38071599E+06)
GET_2DFLD - tidal elevation amplitude
(Min = 3.26393686E-05 Max = 6.68537894E+00)
GET_2DFLD - tidal elevation phase angle
(Min = 2.41898787E-06 Max = 6.28318020E+00)
GET_2DFLD - tidal current inclination angle
(Min = 4.97589188E-07 Max = 6.28318253E+00)
GET_2DFLD - tidal current phase angle
(Min = 2.59098604E-06 Max = 6.28318291E+00)
GET_2DFLD - maximum tidal current, ellipse major axis
(Min = 1.27754016E-06 Max = 1.42023097E+00)
GET_2DFLD - minimum tidal current, ellipse minor axis
(Min = -5.23102760E-01 Max = 4.38586697E-01)
GET_2DFLD - surface u-momentum stress, t = 345 00:00:00
(Rec=0000012, Index=1, File: barotropic_tide_frc.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_2DFLD - surface v-momentum stress, t = 345 00:00:00
(Rec=0000012, Index=1, File: barotropic_tide_frc.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
INQUIRE - unable to assign file counter, Fcount = 0
while processing structure:
and variable; zeta_east
NETCDF_OPEN - unable to open existing NetCDF file:
call from: inquire.FQUIRE