CALL MDSFINDUNIT( mUnit, myThid )

        OPEN( mUnit, FILE=mFileName, STATUS='old',

         READ( mUnit, FMT='(A)', END=1001 ) lineBuf

           READ( mUnit, FMT='(A)', END=1001 ) lineBuf

          READ(  mUnit, FMT='(A)', END=1001 ) lineBuf

           READ( mUnit, FMT='(20(2X,A8,1X))', ERR=1004, END=1004 )

          READ(  mUnit, FMT='(A)', END=1001 ) lineBuf

        CLOSE(mUnit)

      INTEGER mUnit, errIO

       WRITE(mUnit,'(3A)') " simulation = { '",simulName(1:iL),"' };"

      WRITE(mUnit,'(1X,A,I3,A)') 'nDims = [ ',nDims,' ];'

      WRITE(mUnit,'(1X,A)') 'dimList = ['

         WRITE(mUnit,'(1X,3(I5,","))')    (dimList(i,j),i=1,3)

         WRITE(mUnit,'(1X,2(I5,","),I5)') (dimList(i,j),i=1,3)

         WRITE(mUnit,'(1X,3(I10,","))')     (dimList(i,j),i=1,3)

         WRITE(mUnit,'(1X,2(I10,","),I10)') (dimList(i,j),i=1,3)

      WRITE(mUnit,'(1X,A)') '];'

        WRITE(mUnit,'(1X,2(A,I5),A)') 'map2glob = [ ',

       WRITE(mUnit,'(1X,A)') "dataprec = [ 'float32' ];"

       WRITE(mUnit,'(1X,A)') "dataprec = [ 'float64' ];"

      WRITE(mUnit,'(1X,A,I5,A)') 'nrecords = [ ',nrecords,' ];'

     & WRITE(mUnit,'(1X,A,I10,A)') 'timeStepNumber = [ ',myIter,' ];'

       WRITE(mUnit,'(1X,3A)') 'timeInterval = [', msgBuf(1:20*ii),' ];'

       WRITE(mUnit,'(1X,A,1PE21.14,A)')

        WRITE(mUnit,'(1X,A,I4,A)') 'nFlds = [ ', nFlds, ' ];'

        WRITE(mUnit,'(1X,A)') 'fldList = {'

        WRITE(mUnit,'(20(A2,A8,A1))')

        WRITE(mUnit,'(1X,A)') '};'

       WRITE(mUnit,'(3A)')' /* ', titleLine(1:iL), ' */'

      CLOSE(mUnit)

      INTEGER mUnit

      CALL MDSFINDUNIT( mUnit, myThid )

      OPEN( mUnit, file=mFileName, status='unknown',