Out-of-bounds error in mod_ncparam.F, variable idTTav

[m.hadfield]

When I run the UPWELLING case with ECOSIM defined and with bounds-checking enabled, I get the following error

Code: Select all

0: Subscript out of range for array idttav (mod_ncparam.f90: 1938)
    subscript=3, lower bound=1, upper bound=2, dimension=1

The error arises in this section of code:

Code: Select all

            SELECT CASE (Vinfo(1))
              CASE ('tracer2')
                varid=varid-1
                DO i=1,MT
                  varid=varid+1
                  idTTav(i)=varid
                  ...

Array idTTav is described as "<t*t> for active tracers" and has been allocated with an upper bound of NAT

Code: Select all

      allocate ( idTTav(NAT) )

The solution seems to be to allocate idTTav (and similar variables iHUTav, iHVTav, idUTav & idVTav) with an upper bound of MT.

[arango]

Oops, I missed that one. It is amazing what compilers allow nowadays. I did ran the upwelling test case with biology inside of TotalView with ifort. The TotalView and ifort have a lot of problems due to multiple directives in the assemble code. I have been struggling with this for two years. Each one is blaming the other. I think that it is TotalView fault. I haven't being able to run TotalView and PGI on the iMac.

Thank you for reporting this bug The code is now corrected.