#!./mcsh

### ##########################################################################
### MC = < Manifold Code >
### Copyright (C) 1994-- Michael Holst
###
### This library is free software; you can redistribute it and/or
### modify it under the terms of the GNU Lesser General Public
### License as published by the Free Software Foundation; either
### version 2.1 of the License, or (at your option) any later version.
###
### This library 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
### Lesser General Public License for more details.
###
### You should have received a copy of the GNU Lesser General Public
### License along with this library; if not, write to the Free Software
### Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
###
### rcsid="$Id: runme,v 1.3 2010/08/12 05:16:38 fetk Exp $"
### ##########################################################################

##############################################################################
# File:     runexample (interactive)
# Purpose:  Backward Euler time-stepping
# Authors:  Kaihsu Tai, Stephen Bond, Michael Holst
# rcsid="$Id: runme,v 1.3 2010/08/12 05:16:38 fetk Exp $"
##############################################################################

set ISKEY file
set IFNAM mcin-cube.m
set ISNAM 0
set IHVAL unix

set LMAX 1000
set LTOL 1e-10

# write/plot
set FKEY 0
set GVAL 1.0

qread
qinit

mark
refine
qinit

mark
refine
qinit

mark
refine
qinit

mark
refine
qinit

mark
refine
qinit

mark
refine
qinit

set TMETH 0
set TSTEP 0.025
set NSTEP 4000
set PTYPE 2 # (0=Matlab,1=GV,2=GMV,3=DX,4=UCD,5=TEC)
set PFREQ 200
set EFREQ 20

set OSKEY file
set OSFMT asc
set OFNAM output/stepNUM.gmv
# GV     # set OFNAM output/stepNUM.gv
# OPENDX # set OFNAM output/stepNUM.dx
# UCD    # set OFNAM output/stepNUM.ucd
# TECPLOT# set OFNAM output/stepNUM.plt

qstep
