Error when running with EBUS after compilation, 2nd nesting

Im using the script below to run the model (second nesting), but Im getting this error:

ERROR: keyword not found: diagbioFlux_history_fields
ERROR: keyword not found: diagbioVSink_history_fields
ERROR: keyword not found: diagbioGasExc_history_fields
ERROR: keyword not found: diagbioFlux_average_fields
ERROR: keyword not found: diagbioVSink_average_fields
ERROR: keyword not found: diagbioGasExc_average_fields

script:

title:
BENGUELA TEST MODEL ZOOM LEVEL #1
time_stepping: NTIMES dt[sec] NDTFAST NINFO
1 25 60 1
S-coord: THETA_S, THETA_B, Hc (m)
5.0d0 0.88d0 70.0d0
grid: filename
ROMS_FILES/roms_grd.nc.2
forcing: filename
ROMS_FILES/roms_frc.nc.2
bulk_forcing: filename
ROMS_FILES/roms_blk.nc.2
climatology: filename
ROMS_FILES/roms_clm.nc
boundary: filename
ROMS_FILES/roms_bry.nc
initial: NRREC filename
1
ROMS_FILES/roms_ini.nc.1
restart: NRST, NRPFRST / filename
103680 -1
ROMS_FILES/roms_rst.nc.2
history: LDEFHIS, NWRT, NRPFHIS / filename
T 10368 0
ROMS_FILES/roms_his.nc.2
averages: NTSAVG, NAVG, NRPFAVG / filename
1 576 0
ROMS_FILES/roms_avg.nc.2

primary_history_fields: zeta UBAR VBAR U V wrtT(1:NT)
T T T T T 30T
auxiliary_history_fields: rho Omega W Akv Akt Aks HBL HBBL Bostr Wstr Ustr Vstr rsw rlw lat sen HEL
F F T F T F T T T T T T 10
F

primary_averages: zeta UBAR VBAR U V wrtT(1:NT)
T T T T T 30T
auxiliary_averages: rho Omega W Akv Akt Aks HBL HBBL Bostr Wstr Ustr Vstr rsw rlw lat sen HEL
F T T F T F T T T T T T 10
F

rho0:
1025.d0

lateral_visc: VISC2, VISC4 [m^2/sec for all]
0. 0.

tracer_diff2: TNU2(1:NT) [m^2/sec for all]
30*0.0d0

tracer_diff4: TNU4(1:NT) [m^4/sec for all]
30*0.0d8

vertical_mixing: Akv_bak, Akt_bak [m^2/sec]
0.d0 30*0.d0

bottom_drag: RDRG [m/s], RDRG2, Zob [m], Cdb_min, Cdb_max
0.0d-04 1.25d-3 0.d-3 1.d-4 1.d-1

gamma2:
1.0d0

sponge: X_SPONGE [m], V_SPONGE [m^2/sec]
0.5e5 10.0d0

nudg_cof: TauT_in, TauT_out, TauM_in, TauM_out [days for all]
0.1 5. 0.1 5.

diagnostics: ldefdia nwrtdia nrpfdia /filename
T 5760 0
ROMS_FILES/roms_dia.nc.1

diag_avg: ldefdia_avg ntsdia_avg nwrtdia_avg nprfdia_avg /filename
T 1 432 0
ROMS_FILES/roms_dia_avg.nc.1

diag3D_history_fields: diag_tracers3D(1:NT)
30*T

diag2D_history_fields: diag_tracers2D(1:NT)
30*T

diag3D_average_fields: diag_tracers3D_avg(1:NT)
30*T

diag2D_average_fields: diag_tracers2D_avg(1:NT)
30*T

diagnosticsM: ldefdiaM nwrtdiaM nrpfdiaM /filename
F 5760 0
ROMS_FILES/roms_diaM.nc.2

diagM_avg: ldefdiaM_avg ntsdiaM_avg nwrtdiaM_avg nprfdiaM_avg /filename
F 1 5760 0
ROMS_FILES/roms_diaM_avg.nc.2

diagM_history_fields: diag_momentum(1:2)
T T

diagM_average_fields: diag_momentum_avg(1:2)
T T

diagnostics_bio: ldefdiabio nwrtdiabio nrpfdiabio /filename
T 5760 0
ROMS_FILES/roms_diabio.nc.2

diagbio_avg: ldefdiabio_avg ntsdiabio_avg nwrtdiabio_avg nprfdiabio_avg /filename
T 1 1152 0
ROMS_FILES/roms_diabio_avg.nc.2

biology: forcing file
roms_frcbio.nc.1

sediments: input file
sediment.in
sediment_history_fields: bed_thick bed_poros bed_fra(sand,silt)
T F T T

bbl_history_fields: Abed Hripple Lripple Zbnot Zbapp Bostrw
T F F T F T

floats: LDEFFLT, NFLT, NRPFFLT / inpname, hisname
T 6 0
floats.in
ROMS_FILES/floats.nc
float_fields: Grdvar Temp Salt Rho Vel
T T T T T

stations: LDEFSTA, NSTA, NRPFSTA / inpname, hisname
T 400 0
stations.in
ROMS_FILES/stations.nc
station_fields: Grdvar Temp Salt Rho Vel
T T T T T

psource: Nsrc Isrc Jsrc Dsrc Qbar [m3/s] Lsrc Tsrc
0
129 174 1 -28797.3034 F T 17.2152 2 6.797
40 145 0 -6960.7174 F T 25.1539 2 0.2943
44 148 1 -1334.8814 F T 25.2814 2 0.2848
45 147 0 -741.8801 F T 25.2814 2 0.2848
psource_ncfile: Nsrc Isrc Jsrc Dsrc rivdir Lsrc Tsrc runoff file name
ROMS_FILES/roms_runoff.nc
4
129 174 1 -1 30T 17.2152 2
40 145 0 0 30
T 25.1539 2
44 148 1 -1 30T 25.2814 2
45 147 0 0 30
T 25.2814 2

online: byear bmonth recordsperday byearend bmonthend / data path
2000 1 4 2000 3
DATA/CFSR_Benguela_LR/

Hi, I don’t know BioEBUS, but did you try to copy this in your croco_inter.in.2 ?

diagbioFlux_history_fields: wrtdiabioFlux
50*T

diagbioVSink_history_fields: wrtdiabioVSink
50*T

diagbioGasExc_history_fields: wrtdiabioGasExc
50*T

diagbioFlux_average_fields: wrtdiabioFlux_avg
50*T

diagbioVSink_average_fields: wrtdiabioVSink_avg
50*T

diagbioGasExc_average_fields: wrtdiabioGasExc_avg
50*T

Guillaume

Hi!, I’ve tried running again the model with that, but now Im getting this error:
Do you have any idea why this might be happening?
Thank you in advance

Computing YEAR 8 MONTH 4
forrtl: severe (59): list-directed I/O syntax error, unit 7, file /Volumes/Elements/ROMS_bio_gandy/RUN_AgrifEBUS/AGRIF_FixedGrids.in
Image PC Routine Line Source
romsEBUS1 000000010AA4D071 for__io_return Unknown Unknown
romsEBUS1 000000010AA7F85D for_read_seq_lis_ Unknown Unknown
romsEBUS1 000000010AA7E2AD for_read_seq_lis Unknown Unknown
romsEBUS1 000000010A8E0878 agrif_read_fixed_ Unknown Unknown
romsEBUS1 000000010A8E0064 computenbmaxtimes Unknown Unknown
romsEBUS1 000000010A7A4257 MAIN__ Unknown Unknown
romsEBUS1 000000010A7A296E main Unknown Unknown

Hi, Can you post here the file AGRIF_FixedGrids.in you try to read , please ?

Im using this :

2
75    239    15    125    3    3    3    3
120    440    30    215    3    3    4    4
0

number of children per parent

imin imax jmin jmax spacerefx spacerefy timerefx timerefy

[all coordinates are relative to each parent grid!]

number of children : 2

imin imax jmin jmax spacerefx spacerefy timerefx timerefy : 75 239 15 125 3 3 3 3

imin imax jmin jmax spacerefx spacerefy timerefx timerefy :120 440 30 215 3 3 4 4

Can send your file directly with a link in the forum instead ? Thanks

This is the AGRIF im using for the second nesting:

AGRIF_FixedGrids.in

I don t know but you have “~” in last line, can you remove it to see what happens ? thanks

I removed that, but nope. It doesn’t run. Actually I compile the model again just to double check all was ok. and now Im back to the first sets of errors Ive posted here. I just wanted to get Net Primary Production (NPP) from this second nesting, since I couldn’t get it from the first one.