Custom Query (964 matches)

When running FLT_TEST with SOLVE3D defined, with MPI and bounds-checking enabled, an out of bounds error is triggered in step_floats_tile in the code that keeps the floats at their initial depth. The specific statement where the out-of-bounds error occurs is at line 292 in step_floats.F and the failure occurs for float l=231 (i.e. near the top of the grid). The reason it occurs is the horizontal locations, track(ixgrd,nfp1,l) and track(ixgrd,nfp1,l), appear not to have been initialised.

I will look into this further when I can.

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.

In the Kelvin test case, a spike in sea surface height develops in the SW corner of the domain, reaching a value of 46 m at t = 2 hours. The model blows up at t = 3 hours.