missing line for PISCES ln_river in p4zsed.F90 (master branch)
! Add the external input of nutrients from the river **(phosphate and silicate)**
! ----------------------------------------------------------
IF( ln_river ) THEN
DO jj = JRANGE
DO ji = IRANGE
tra(ji,jj,1,jppo4) = tra(ji,jj,1,jppo4) + rivdip(ji,jj) * rfact2 ! **for PO4**
tra(ji,jj,1,jpsil) = tra(ji,jj,1,jpsil) + rivdsi(ji,jj) * rfact2 ! **for Si**
ENDDO
ENDDO
ENDIF