INTEGER diagLoc_ioUnit
& diagLoc_ioUnit
IF ( diagLoc_ioUnit.EQ.0 ) THEN
CALL MDSFINDUNIT( diagLoc_ioUnit, myThid )
OPEN( diagLoc_ioUnit, FILE=dataFName, STATUS='unknown' )
& 'open unit=',diagLoc_ioUnit, ', file: ',dataFName(1:iL)
IF ( diagLoc_ioUnit.GT.0 ) THEN
WRITE(diagLoc_ioUnit,'(1P2E18.10,A,2I4,I8,A,2I4,I6,2A)')
diagLoc_ioUnit = 0
CLOSE(diagLoc_ioUnit)