⁴We assume that Nt, NEQ are positive and small enough so that the requested memory is available. It is better to use the call allocate(T(Nt),STAT=IERR). The non zero value of IERR after the call indicates a successful allocation and the following test stops the program otherwise: IF(IERR .eq. 0) STOP ’Memory allocation for T failed’