GMSH translator (umesh1.f) usage unsuccessful

[huxiang_shi]

Dear Professors, dear FEAP community,
I am using Version 8.6.1o  19/11/2022, and I am using the GMSH translator which is made by Prof. Taylor in this linkage http://feap.berkeley.edu/forum/index.php?topic=1842.msg7326#msg7326. I have attached the files in the linkage that I used.
I replaced umesh1.f by the new file made by Prof. Taylor, and implemented the input file Isquare for calculation, but the results came out with errors, here is the content of Lsquare file:

   --------------------------------- START OF FEAP LOG ----------------------------------

 *ERROR* ENDCLR: End of file encountered in TINPUT
         Unit Number =  15

       0    0   0  0.000E+00  0.00E+00  0.00E+00            0.00E+00                0.05   

 Total        0        0
  ---------------------------------- END OF FEAP LOG -----------------------------------

Can you point out where the places I should correct. Thanks very much.

[Prof. S. Govindjee]

Can you say more about what you are doing?  I compiled the umesh1.f into 8.6.1o (November 19, 2022) and it ran your input file just fine.  Attached are the outputs.  On the screen, I get:

Code: [Select]

MESH:square             
 MESH:square.msh                    10
     Number Nodes    =     564
     Number Elements =    1130
        Type =  15 Part   =   0 Region =   1 NEL =   1 Node Point       
        Type =  15 Part   =   0 Region =   2 NEL =   1 Node Point       
        Type =  15 Part   =   0 Region =   3 NEL =   1 Node Point       
        Type =  15 Part   =   0 Region =   4 NEL =   1 Node Point       
        Type =   1 Part   =   0 Region =   1 NEL =   2 Node Line         
        Type =   1 Part   =   0 Region =   2 NEL =   2 Node Line         
        Type =   1 Part   =   0 Region =   3 NEL =   2 Node Line         
        Type =   1 Part   =   0 Region =   4 NEL =   2 Node Line         
        Type =   2 Part   =   0 Region =   6 NEL =   3 Node Triangle     

     Material    Element Tag   Element Type  NH1/2 History    NH3 History
          1              1      Solid                  0              0

     P a r t i t i o n   1

     E q u a t i o n / P r o b l e m   S u m m a r y:

       Mesh dimension  (ndm) =         2 :  Number nodes     =             564
       Number dof/node (ndf) =         2 :  Number elements  =            1046
       Element eqs.   (nadd) =         0 :  Number materials =               1
       Number rigid bodies   =         0 :  Number equations =            1128
                                         :  Number global eq =               0
       Number joints         =         0 :  Number tang terms=          427936
       Transient integrator  = Static    :  Average col. ht. =             380



[huxiang_shi]

Dear Professor,
I debugged the process, and found that at the line 108 and 176 of umesh1.f, do while (.not.pcomp(tr,'   ',4)), in the '   ', only three blank spaces are used, I think there may be four blank spaces same with the number 4, after I changed it, I get the same result as yours, but a window showing the mesh as shown in your picture did't come out.
Thanks very much for your reply.

[Prof. S. Govindjee]

I made the plot separately from the interactive prompt; plot,mesh.  I agree with your correction to the file.  The blank string should have 4 spaces in it. It is interesting that your computer failed on the string check whereas mine did not care; I guess I got lucky.