I noticed two things from your ncdump output
1) NO3_TIME = UNLIMITED ; // (35 currently)
NO3 time as a RECORD variable, probable a harmless issue
2) It is NO3_TIME, and not NO3_time. And the same with O2_WEST, which should be O2_west (same for all the BGQ variables)
Not sure if this is the problem, but try renaming the variables, at least, to lower case.
Mines look like this:
double ZOO_north(zoo_time, s_rho, xi_rho) ;
ZOO_north:long_name = "northern boundary ZOO" ;
ZOO_north:units = "mMol N m-3" ;
double ZOO_south(zoo_time, s_rho, xi_rho) ;
ZOO_south:long_name = "southern boundary ZOO" ;
ZOO_south:units = "mMol N m-3" ;
double ZOO_west(zoo_time, s_rho, eta_rho) ;
ZOO_west:long_name = "western boundary ZOO" ;
ZOO_west:units = "mMol N m-3" ;
double bry_time(bry_time) ;
bry_time:long_name = "time for boundary climatology" ;
bry_time:units = "day" ;
bry_time:calendar = "360.0 days in every year" ;
bry_time:cycle_length = 360. ;
double chla_time(chla_time) ;
chla_time:long_name = "time for CHLA climatology" ;
chla_time:units = "day" ;
chla_time:cycle_length = 360. ;
double hc(one) ;
hc:long_name = "S-coordinate parameter, critical depth" ;
hc:units = "meter" ;
double no3_time(no3_time) ;
no3_time:long_name = "time for NO3 climatology" ;
no3_time:units = "day" ;
Hope this helps.
Andres
Also suggestion if any on adding the boundary for biology directly into croco_bry.nc (for running inter-annual runs boundary downloaded from CMEMS) are welcome. Tried and spent lot of time but unsuccessful.