! ---------------------------------
! -- Diagnostics
!     INTEGER           :: unitnum
!     CHARACTER(len=30) :: ctemp
!     CHARACTER(len=3)  :: cagain
!     DATA iagain /0/
!     SAVE iagain
!     iagain=iagain+1
!     write(cagain,'(i3.3)') iagain
! ---------------------------------
! ---------------------------------------------------------------
! -- Diagnostics
!        envstring = 'WORKDIR'
!        unitnum=89
!        ctemp = 'diag'//trim(adjustl(get_cmytid()))//"."//cagain
!        call openf(unitnum,ctemp,envstring)
! ---------------------------------------------------------------
!Diag close(89)