ESTUARY_TEST sediment-related bugs

[m.hadfield]

In sed_bed.F and sed_fluxes.F, when ASSUMED_SHAPE is not defined the following appear:

   real(r8), intent(inout) :: ero_flux(LBi:UBi,LBj:UBj,:) 
   real(r8), intent(inout) :: settling_flux(LBi:UBi,LBj:UBj,:) 

That should be

   real(r8), intent(inout) :: ero_flux(LBi:UBi,LBj:UBj,NST) 
   real(r8), intent(inout) :: settling_flux(LBi:UBi,LBj:UBj,NST) 

Corrected files attached.

Also in subroutine readPhyPar, file in_par.F, I'm getting bounds-checking areas in this area of the code

          ELSE IF (TRIM(KeyWord).eq.'Hout(idBott)') THEN
            IF (MAXVAL(idBott).eq.0) THEN
              WRITE (out,280) 'idBott'
              STOP
            END IF
            Npts=load_l(Nval, Cval, MBOTP*Ngrids, Lbottom)
            DO ng=1,Ngrids
              DO itrc=1,MBOTP
                i=idBott(itrc)
>>error here>>  Hout(i,ng)=Lbottom(itrc,ng)
              END DO
            END DO

because for itrc=17, idBott(itrc) is equal to zero.

[arango]

Yes, thank you. The problem with idBott(itrc) for itrc=17 is a typo in mod_ncparam.F at line 1494. We need idefx instead idrfx.