## Copyright (C) 2013 fotios
## 
## This program is free software; you can redistribute it and/or modify
## it under the terms of the GNU General Public License as published by
## the Free Software Foundation; either version 3 of the License, or
## (at your option) any later version.
## 
## This program is distributed in the hope that it will be useful,
## but WITHOUT ANY WARRANTY; without even the implied warranty of
## MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
## GNU General Public License for more details.
## 
## You should have received a copy of the GNU General Public License
## along with Octave; see the file COPYING.  If not, see
## <http://www.gnu.org/licenses/>.

## fixmesh

## Author: fotios <address@hidden>
## Created: 2013-05-30

function [ s ] = fixmesh (noi, h)

  if ( !isscalar (noi) || !isscalar (h) )
    error ("both arguments must be scalars");
  endif

  fid = fopen ("ff2O_fix.idp", "w");
    for i = 1:noi
      fprintf (fid, ["Th = adaptmesh(Th, " num2str(h) ", IsMetric = 1, nbvx = 1000000);\n"]);
    endfor
	status = fclose (fid);

  if ( status == 0 )
    s = 1;
  else
    s = 0;
    warning ("mesh fix is not applied");
  endif
  
endfunction