Declaration problem in cgradient.h

[m.hadfield]

In file cgradient.h, subroutine read_state.F, the explicit-shape version of the declaration for argument s_t is

real(r8), intent(inout) :: s_t(LBi:UBi,LBj:UBj,N(ng),NT(ng))

This has only 4 dimensions, whereas the corresponding assumed-shape declaration has 5:

real(r8), intent(inout) :: s_t(LBi:,LBj:,:,:,:)

I think the correct explicit-shape declaration is

real(r8), intent(inout) :: s_t(LBi:UBi,LBj:UBj,N(ng),2,NT(ng))

Corrected file attached.

[arango]

This array needs to be declared as:

real(r8), intent(inout) :: s_t(LBi:UBi,LBj:UBj,N(ng),3,NT(ng))