diff --git a/fem/src/MeshRemeshing.F90 b/fem/src/MeshRemeshing.F90
index 809fe80330..6760e51e19 100644
--- a/fem/src/MeshRemeshing.F90
+++ b/fem/src/MeshRemeshing.F90
@@ -1634,12 +1634,12 @@ SUBROUTINE RemeshMMG3D(Model, InMesh,OutMesh,EdgePairs,PairCount,&
     END IF
     
     IF(SaveMMGMeshes) THEN
-      WRITE(MeshName, '(A,i0,A)') TRIM(premmg_meshfile), time, '.mesh'
+      MeshName = TRIM(premmg_meshfile) // I2S(time) // '.mesh'
       CALL MMG3D_SaveMesh(mmgMesh,MeshName,LEN(TRIM(MeshName)),ierr)
     END IF
     
     IF(SaveMMGSols) THEN
-      WRITE(SolName, '(A,i0,A)') TRIM(premmg_solfile), time, '.sol'
+      SolName =  TRIM(premmg_solfile) // I2S(time) // '.sol'
       CALL MMG3D_SaveSol(mmgMesh, mmgSol,SolName,LEN(TRIM(SolName)),ierr)
     END IF
     
@@ -3072,7 +3072,7 @@ FUNCTION GET_MMG2D_MESH(MeshNumber,OutputFilename) RESULT(NewMesh)
     ! Initialize the new mesh stucture
     NewMesh => AllocateMesh(nt + nt0,na + na0,np + np0,ParEnv%PEs > 1)
     IF (MeshNumber > 0 ) THEN
-      WRITE(NewMesh % Name,'(A,A,I0)') TRIM(OutPutFileName),'_N',MeshNumber
+      NewMesh % Name = TRIM(OutPutFileName) // '_N' // I2S(MeshNumber)
     ELSE
       NewMesh % Name = TRIM(OutPutFileName)
     END IF
diff --git a/fem/src/ModelDescription.F90 b/fem/src/ModelDescription.F90
index bbb3296da2..a419c0f498 100644
--- a/fem/src/ModelDescription.F90
+++ b/fem/src/ModelDescription.F90
@@ -3257,6 +3257,7 @@ FUNCTION LoadModel( ModelName,BoundariesOnly,numprocs,mype,MeshIndex) RESULT( Mo
           i = i + 1
         END DO
 
+        IF(ALLOCATED(Solver % Mesh % Name)) DEALLOCATE(Solver % Mesh % Name)
         ALLOCATE(CHARACTER(j)::Solver % Mesh % Name)
         i = k
         k = 1