Sediment bottom properties initial conditions logic

[arango]

Corrected a problem in checkvars.F in sediment or BBL applications when the bottom properties are specified with analytical expressions (ANA_SEDIMENT).

# if defined SEDIMENT || defined BBL_MODEL
#  ifndef ANA_SEDIMENT
         DO itrc=1,MBOTP
           get_var(idBott(itrc))=.TRUE.
         END DO
#  endif
# endif

The conditional # ifndef ANA_SEDIMENT was missing. This forced to look for this variables in the initial conditions NetCDF file.

Recall that it is possible to combine analytical expression with physical initial conditions (zeta, u, v, T, S) in the NetCDF file in sediment or biology models applications. If either ANA_SEDIMENT or ANA_BIOLOGY are activated, the user needs to provide analytical expressions for all the state variables associated with these models.

Many thanks to Laura Bianucci for reporting this in the ​forum. Also many thanks to John Wilkin and Chris Sherwood for bringing this to my attention.