hallo
ich habe ein problem mit folgenden programm.
Das Programm ändert in jedem durchgang einen bestimmten wert der Paramaterdatei dirneu.parset.jedoch wenn man das programm startet wird der eintrag "theta=" in dirneu.parset nichts zugeordnet.
wo ist der fehler? :(
$level=1; #Anzahl der Gitter-level
$caunter=0; #Zählvariable
$theda=0.01;
while ($theda<=1.01) # theda läuft von 0.1 bis 1 in hunderstel-Schritten 0.01,0.02,0.03,0.04....
{
#print "Programm ist in theda for schliefe: $theda";
#print "\n"; #10
########################## Anfang Programmabschnitt 0
open(DATEI,">dirneu.parset"); #Parameterdatei wird geaendert und Programm wird ausgeführt
#überschreibe Datei dirneu.parset mit folgendem String.
print DATEI "# Number of grid levels
numLevels = $level
# Max number of domain decomposition iterations
numDDIt = 10000
# Number of multigrid iterations
numIt = 50
# Dirichlet-Neumann damping factor
theta = $theta
# Number of presmoothing steps
nu1 = 3
# Number of postsmoothing steps
nu2 = 3
# Number of coarse grid corrections
mu = 1
# Number of base solver iterations
baseIt = 100
# Tolerance of the DD iteration
ddTolerance = 1e-12
# Tolerance of the solver
tolerance = -1
# Tolerance of the base grid solver
baseTolerance = -1
####################################################
# Problem Geometry Specifications
####################################################
#path = /home/cocktail/sander/data/dd16/richards/twosquares/
#gridFile0 = twosquares0.grid
#gridFile1 = twosquares1.grid
#dirichletFile0 = twosquares0.dn
#dirichletFile1 = twosquares1.dn
#signoriniFile0 = twosquares0.nosig
#signoriniFile1 = twosquares1.nosig
#path = /home/haile/sander/data/richards/coffee/
#gridFile0 = coffee0.grid
#gridFile1 = coffee1.grid
#dirichletFile0 = coffee0.dn
#dirichletFile1 = coffee1.dn
#signoriniFile0 = coffee0.nosig
#signoriniFile1 = coffee1.nosig
#numGrids = 2
#path = /home/haile/sander/data/richards/unitsquare/
#gridFile0 = unitsquare_0.grid
#gridFile1 = unitsquare_1.grid
#gridFile2 = unitsquare_2.grid
#gridFile3 = unitsquare_3.grid
#dirichletFile0 = unitsquare_0.nn
##neumannNodesFile0 = unitsquare_0.nn
##neumannValuesFile0 = unitsquare_0.nv
#neumannNodesFile1 = unitsquare_1.nn
#neumannValuesFile1 = unitsquare_1.nv
numGrids = 2
path = /home/cocktail/jschreib/data/richards/yinyang/
gridFile0 = yin.grid
gridFile1 = yang.grid
dirichletFile0 = yin.dn
dirichletFile1 = yang.dn
materialFile0 = materials/very_coarse_sand.material
materialFile1 = materials/fine_clay.material";
close DATEI; #Ende der Überschreibung, schließe Datei. 90
$theda=$theda+0.01;
}
Datum: 04.03.2008-09:51
Problem mit Dateiveränderung
