Synchronization with mat_ini_cubep3m.F90. The common parts should really
be in a common module!!