| RUNNING CODE: (How to make
movies of poloidal field evolution in Dikpati's kinematic solar
dynamo model)
Go to working directory on zev,
and set up directories for your runs.
Write aliases into /home/zita/.cshrc
zev05 (/net/zev/i/zita/dyn05)
goz (/net/zev/i/zita/dyn05/dt001/pb68/3z1_1z2)
goz2 (/net/zev/i/zita/dyn05/dt001/pb68/3z1_02z2)
[Results can later be moved to http://download.hao.ucar.edu/pub/green/...]
green05 (/toshi/ftp/pub/green/dynamo/2005/etath)
got (/toshi/ftp/pub/green/dynamo/2005/etath/dt001/pb68/3z1_1z2)
got2 (/toshi/ftp/pub/green/dynamo/2005/etath/dt001/pb68/3z1_02z2)
Compile and run code to generate
data:
$ pgf77 -o evolveta evolveta.f
$ ./evolveta > log & (& puts run in background
so you can still work on terminal)
Or run using input parameter file in.dat:
$ runevo (./evolveta <in.dat >log &)
------------------------------------------------------------------------------
Run using scripts: (examples; modify
for new runs)
1. Start ss run:
check in.dat, directory structure, delete excess files
$ getrun ("cp /net/zev/i/zita/dyn05/dt001/pb68/runss
.")
$ runss cp ../../evolveta .
cp ../../mov1 .
mv params.dat in.dat
./evolveta <in.dat >log &
2. Start movie run:
$ cat log (note it=number of iterations to convergence,
and cfl limit)
(decrease dt for movie run if necessary)
$ mov1 mkdir mov
cp final.dat mov/init.dat
mv evolveta mov
cp params.dat mov/in.dat
cp in.dat mov/in.dat
cp ../../getplt mov/.
cp ../../mkall mov/.
cd mov
$ cd mov
(mov) $ vi in.dat (itmax=100->1, irelax=0->1, imovie=0->1,
check other params)
(mov) $ runevo ("./evolveta <in.dat >log &)
3. Generate plots: (still in mov
directory)
$ getplt cp ../../../p* (get plot files)
cp ../../../mk* . (get idl scripts)
cp ../../../cleanp .
./pwrit68 (separate final.dat into a.dat, r.dat, th.dat)
mv params.dat in.dat (just in case old name was used)
$ mkall idl (type in idl script names with @ below. ORDER
IS CRUCIAL)
@mkplt plte68, pfly
pltsnap,infile='movbtor.dat'
mv snap~.eps -> bt~.eps
@mketa pltsnap,infile='moveta.dat'
mv snap~.eps -> e~.eps
setenv DISPLAY pan:0.0 (or replace "pan" with
your local machine)
@mkmov pltsnap,infile='movie.dat'
mv snap~.eps -> bp~.eps
pltmov,infile='movie.dat',outfile='movie.gif',fb=0.68
4. Check results, move plots to toshi, and cleanup
$ ls
$ cps or cpslog, edited appropriately for all runs
Ex: cp ./q0s1c3e2/mov/*.eps /toshi/ftp/pub/green/dynamo/2005/etath/dt001/pb68/3z1_1z2/q0s1c3e2/.
$ cleanp rm files no longer needed (but KEEP files.DAT)
------------------------------------------------------------------------------
View movies:
$ xanim movie.gif
OR
- open a web browser,
- navigate to toshi directory on
http://download.hao.ucar.edu/pub/green
- and click on movie.gif or file.eps
Browsers generally know how to display these files.
|