Custom Query (986 matches)

In mod_param.F, NFV is dimensioned by 10+NT or 12+NT, depending. Later, in step_floats.F, the two optional values are overwritten by tracers, plus it is assumed in wrt_floats.F that the tracers are after position 10. It seems to me you need to get that 10/12 thing to propagate into step_floats.F so that the vertical mixing can be output if desired.

ROMS/Nonlinear/sed_bed.F needs the ifdef SUSPLOAD around the ero_flux and settling_flux pointers for argument list, such as:

# ifdef SUSPLOAD
     &                   OCEAN(ng) % ero_flux,                          &
     &                   OCEAN(ng) % settling_flux,                     &
# endif

(this is needed at 4 places at the top of the file, in the argument list).

The code is okay, but I needed to make some modifications to the supplied TEST_HEAD input files to get bedload to be non-zero. These modifications are required because of the change in R224, where bedload was split from the other routines, and the looping was changed to DO ised=NCS+1,NST, yet the input files still used NCS=1, NNT=0 (NST=1), so the loop never gets executed.

I needed to modify the ocean_test_head.in file to set NCS=0, NNS=1, modify the sediment_test_head.in file to specify parameters for SAND instead of MUD, and also change the variable names from MUDFRAC_01 to SANDFRAC_01, MUD_01 to SAND_01, etc in the initialization NetCDF file test_head_ini.nc. Once I made these changes, the TEST_HEAD test case run with the R271 code (current trunk code) produced non-zero bedload values equivalent to versions prior to R224.