[cig-commits] r6044 - geodyn/3D/MAG/trunk/doc

wei at geodynamics.org wei at geodynamics.org
Fri Feb 16 17:27:05 PST 2007 

Author: wei
Date: 2007-02-16 17:27:04 -0800 (Fri, 16 Feb 2007)
New Revision: 6044

Modified:
   geodyn/3D/MAG/trunk/doc/mag_book.lyx
Log:
More in manual

Modified: geodyn/3D/MAG/trunk/doc/mag_book.lyx
===================================================================
--- geodyn/3D/MAG/trunk/doc/mag_book.lyx	2007-02-17 01:25:10 UTC (rev 6043)
+++ geodyn/3D/MAG/trunk/doc/mag_book.lyx	2007-02-17 01:27:04 UTC (rev 6044)
@@ -1157,10 +1157,10 @@
 \end_layout
 
 \begin_layout Standard
-Procedure MAGTS.pro takes data from l-files generated by MAG and creates
- time series plots and statistics.
+Procedure MAGTS.pro and MAGXY.pro take data from l-files generated by MAG
+ and creates time series plots and statistics.
  This version reads an l-file consisting of 17 time series, the first record
- being dimensionless time.
+ being dimensionless time (see Appendix C for details on output file format).
  Energies and rms magnetic field and velocity are scaled as in MAG; tilt
  is dipole vector colatitude; pole longitude is dipole vector longitude.
  Figures 
@@ -1186,7 +1186,7 @@
 
 
 \begin_inset Float figure
-placement H
+placement h
 wide false
 sideways false
 status open
@@ -1222,7 +1222,7 @@
 
 
 \begin_inset Float figure
-placement H
+placement h
 wide false
 sideways false
 status open
@@ -1369,7 +1369,7 @@
 \begin_inset Float figure
 wide false
 sideways false
-status collapsed
+status open
 
 \begin_layout Standard
 \align center
@@ -1481,7 +1481,8 @@
 \begin_inset Formula $lmax=32$
 \end_inset
 
- and a four-fold symmetry is assumed in the longitudinal direction.
+ and a four-fold symmetry is assumed in the longitudinal direction (param.f
+ should be linked to param32s4.f when compile MAG).
  The input parameter files are par.bench0 for case 0 and par.bench1 in case
  1, both files resides in ~/src directory.
 \end_layout
@@ -1923,147 +1924,111 @@
 Reversal Dynamo Case
 \end_layout
 
-\begin_layout Subsection
-Creating MPEGs with IDL
-\end_layout
-
 \begin_layout Standard
-IDL provides an IDLgrMPEG class that allows the user to save an array of
- image frames as an MPEG movie.
- See the IDL Reference guide.
-\end_layout
+In this example, we produce a magnetic field reversal using MAG.
+ The input parameter for this case in the source directory is ~/src/par.Rev.
+ There is no longitudinal symmetry in this case, therefore when compile
+ MAG use param32s1.f linking to param.f.
+ The Ekman number is 
+\begin_inset Formula $E=0.02$
+\end_inset
 
-\begin_layout Standard
-The MPEG_PUT procedure stores the specified image array at a specified frame
- index in an MPEG sequence.
- The MPEG_SAVE procedure encodes and saves an open MPEG sequence.
- The images can be read using the READ_IMAGE function.
-\end_layout
+, the Prandtl number is 
+\begin_inset Formula $Pr=1$
+\end_inset
 
-\begin_layout Standard
-make_mpg, synopsis
+ and the magnetic Prandtl number is 
+\begin_inset Formula $Pm=10$
+\end_inset
+
+.
+ Rayleigh number is 
+\begin_inset Formula $Ea=12000$
+\end_inset
+
+.
 \end_layout
 
-\begin_layout Standard
-This module, mpg.pro reads a sequence of images (named prefix00001.suffix)
- and produces an MPEG1 or MPEG2 movie.
- Supported formats are: BMP, GIF, JPEG, PNG, PPM, SRF, TIFF, DICOM
+\begin_layout Subsection
+Results and Discussions
 \end_layout
 
 \begin_layout Standard
-Parameters
-\end_layout
+This case has been ran on 32 bit and 64 bit intel processors.
+ Figure 
+\begin_inset LatexCommand \ref{frevfield1}
 
-\begin_layout Itemize
-prefix string(required) A string specifying the prefix of the images to
- read included the path, for example "/home/data/images/Image"
-\end_layout
+\end_inset
 
-\begin_layout Itemize
-suffix string(required) A string specifying the suffix of the images to
- read, for example "jpg".
-\end_layout
+ shows a plot of mean velocity Vrms, mean magnetic field Brms, axial dipole
+ and dipole tilt on outer boundary.
+ It indicated a magnetic field reversal between time 25-30.
+ Figure 
+\begin_inset LatexCommand \ref{frevfield2}
 
-\begin_layout Itemize
-n_start integer(required) An integer specifying the first image in the sequence.
-\end_layout
+\end_inset
 
-\begin_layout Itemize
-n_end integer(required) An integer specifying the last image in the sequence.
-\end_layout
+ shows a longer run of MAG, we see the magnetic field reversed again.
+ At this time, magnetic field has weakened substantially.
+ On the pole plot
+\begin_inset LatexCommand \ref{fpole1}
 
-\begin_layout Itemize
-digits integer(required) An integer specifying the number of digits of the
- number field in the sequence, for example image00001.jpg would require digits=5.
-\end_layout
+\end_inset
 
-\begin_layout Itemize
-dims integer(optional) An integer array specifying the size of the output
- image.
- If not specified the size of the first image is used.
+ the top 
 \end_layout
 
-\begin_layout Itemize
-format integer(optional) An integer with values 0 for MPEG1 and 1 for MPEG2.
- Default is MPEG1.
-\end_layout
+\begin_layout Standard
+\begin_inset LatexCommand \label{frevfield1}
 
-\begin_layout Itemize
-frame_rate integer(optional) An integer with values.
-\end_layout
+\end_inset
 
-\begin_layout LyX-Code
-1 = 23.976 frames/sec: NTSC encapsulated film rate.
-\end_layout
 
-\begin_layout LyX-Code
-2 = 24 frames/sec: standard film rate.
-\end_layout
+\begin_inset Float figure
+wide false
+sideways false
+status open
 
-\begin_layout LyX-Code
-3 = 25 frames/sec: PAL video frame rate
-\end_layout
+\begin_layout Standard
+\InsetSpace ~
+\InsetSpace ~
 
-\begin_layout LyX-Code
-4 = 29.97 frames/sec: NTSC video frame rate
-\end_layout
+\begin_inset Graphics
+	filename images/field-64.ps
+	lyxscale 50
+	scale 75
 
-\begin_layout LyX-Code
-5 = 30 frames/sec: NTSC drop frame video rate (the default).
-\end_layout
+\end_inset
 
-\begin_layout LyX-Code
-6 = 50 frames/sec: double frame rate/progressive PAL
-\end_layout
 
-\begin_layout LyX-Code
-7 = 59.94 frames/sec: double frame rate NTSC
 \end_layout
 
-\begin_layout LyX-Code
-8 = 60 frames/sec: double frame rate NTSC
+\begin_layout Caption
+Field Plot for Reversal Dynamo Case
 \end_layout
 
-\begin_layout Itemize
-mpeg_file string(optional) A string specifying the output MPEG file (default
- outfile.mpg)
-\end_layout
+\end_inset
 
-\begin_layout Itemize
-tmp_dir string(optional) A string specifying temporary directory to use
- for the temporary image files.
-\end_layout
 
-\begin_layout Standard
-To run the IDL movie generator type in a shell:
 \end_layout
 
 \begin_layout Standard
-%idl
-\end_layout
+\begin_inset LatexCommand \label{frevfield2}
 
-\begin_layout Standard
-> .compile mpg.pro
-\end_layout
+\end_inset
 
-\begin_layout Standard
-> make_mpg, prefix='image', suffix='.bmp', n_start=0, n_end=100, digits=5,
- dims=[512,512], mpeg_file='mympeg.mpg'
-\end_layout
 
-\begin_layout Standard
-\noindent
-\align center
 \begin_inset Float figure
 wide false
 sideways false
 status open
 
 \begin_layout Standard
-\noindent
-\align center
+\InsetSpace ~
+\InsetSpace ~
+
 \begin_inset Graphics
-	filename images/g7revR.ps
+	filename images/field-64-revR.ps
 	lyxscale 50
 	scale 75
 
@@ -2073,7 +2038,7 @@
 \end_layout
 
 \begin_layout Caption
-Pole Plot 3
+Field Plot for Reversal Dynamo Case
 \end_layout
 
 \end_inset
@@ -2082,46 +2047,245 @@
 \end_layout
 
 \begin_layout Standard
-\noindent
-\align center
+\begin_inset LatexCommand \label{fpole1}
+
+\end_inset
+
+
 \begin_inset Float figure
 wide false
 sideways false
 status open
 
 \begin_layout Standard
-\noindent
-\align center
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+
 \begin_inset Graphics
 	filename images/g1revR.ps
 	lyxscale 50
-	scale 75
+	scale 40
+	clip
 
 \end_inset
 
 
 \end_layout
 
+\begin_layout Standard
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+
+\begin_inset Graphics
+	filename images/g7revR.ps
+	lyxscale 50
+	scale 40
+
 \end_inset
 
 
 \end_layout
 
+\begin_layout Caption
+The top is pole plot at the begining of the second field reversal
+\end_layout
+
 \begin_layout Standard
-\noindent
-\align center
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+bottom is the pole plot at the end of the second field reversal
+\end_layout
+
+\end_inset
+
+
+\end_layout
+
+\begin_layout Subsection
+Creating Reversal Dynamo Movie
+\end_layout
+
+\begin_layout Subsubsection
+Generate Movie Files
+\end_layout
+
+\begin_layout Standard
+MAG has the function to record movie.
+ The input paramter imovopt gives option to write 
+\begin_inset Formula $B_{r}$
+\end_inset
+
+ at the outer surface, the output file is produced with prefix 
+\begin_inset Quotes sld
+\end_inset
+
+mm.
+\begin_inset Quotes srd
+\end_inset
+
+ when imovopt=100, or the digit at the hundredth is larger than zero.
+ For our reversal dynamo case we exam the field plot and decide to record
+ the first field reversal.
+ As shown the figure 
+\begin_inset LatexCommand \ref{fmovie1}
+
+\end_inset
+
+, we pick the time interval to generate the movie file.
+ The input parameter for movie recording is in ~/src/par.revRmv.
+ We choose restart file d5 as starting point, and set to record movie at
+ time 
+\begin_inset Formula $t=22$
+\end_inset
+
+.
+ We recorded 400 frames over 8 time unit.
+ The sample output files are in ~/rev-data.
+ 
+\end_layout
+
+\begin_layout Standard
+This version of MAG provides an IDL routine magmovieCIG.pro (~/idl/magmovieCIG.pro
+), it read in the movie file 
+\begin_inset Quotes sld
+\end_inset
+
+mm.
+\begin_inset Quotes srd
+\end_inset
+
+ and displays the magnetic field at the outer surface.
+ This procedure can also create JPG files of the movie-frame images.
+ Figure 
+\begin_inset LatexCommand \ref{fmframe}
+
+\end_inset
+
+shows the first and the last frames of the reversal dynamo movie.
+ 
+\end_layout
+
+\begin_layout Standard
+\begin_inset LatexCommand \label{fmovie1}
+
+\end_inset
+
+
 \begin_inset Float figure
 wide false
 sideways false
 status open
 
 \begin_layout Standard
-\noindent
-\align center
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+
 \begin_inset Graphics
-	filename images/poler.ps
+	filename images/ex3fp.eps
 	lyxscale 50
-	scale 75
+	scale 50
 
 \end_inset
 
@@ -2129,7 +2293,7 @@
 \end_layout
 
 \begin_layout Caption
-Pole Plot
+Interval to record movie
 \end_layout
 
 \end_inset
@@ -2138,28 +2302,135 @@
 \end_layout
 
 \begin_layout Standard
-\noindent
-\align center
+\begin_inset LatexCommand \label{fmframe}
+
+\end_inset
+
+
 \begin_inset Float figure
 wide false
 sideways false
 status open
 
 \begin_layout Standard
-\noindent
-\align center
 \begin_inset Graphics
-	filename images/field-64-revR.ps
+	filename images/revmv006.JPG
 	lyxscale 50
-	scale 75
+	scale 50
 
 \end_inset
 
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
 
+\begin_inset Graphics
+	filename images/revmv398.JPG
+	lyxscale 50
+	scale 50
+
+\end_inset
+
+
 \end_layout
 
+\begin_layout Standard
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+a\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+\InsetSpace ~
+b 
+\end_layout
+
 \begin_layout Caption
-Field Plot
+a: first frame of movie, b: last frame of movie
 \end_layout
 
 \end_inset
@@ -2167,6 +2438,137 @@
 
 \end_layout
 
+\begin_layout Subsubsection
+Creating MPEGs with IDL
+\end_layout
+
+\begin_layout Standard
+IDL provides an IDLgrMPEG class that allows the user to save an array of
+ image frames as an MPEG movie.
+ See the IDL Reference guide.
+\end_layout
+
+\begin_layout Standard
+The MPEG_PUT procedure stores the specified image array at a specified frame
+ index in an MPEG sequence.
+ The MPEG_SAVE procedure encodes and saves an open MPEG sequence.
+ The images can be read using the READ_IMAGE function.
+\end_layout
+
+\begin_layout Standard
+This module, mpg.pro reads a sequence of images (named prefix00001.suffix)
+ and produces an MPEG1 or MPEG2 movie.
+ Supported formats are: BMP, GIF, JPEG, PNG, PPM, SRF, TIFF, DICOM.
+ mpg.pro runs with the following parameters:
+\end_layout
+
+\begin_layout List
+\labelwidthstring 00.00.0000
+prefix string(required) A string specifying the prefix of the images to
+ read included the path, for example "/home/data/images/Image"
+\end_layout
+
+\begin_layout List
+\labelwidthstring 00.00.0000
+suffix string(required) A string specifying the suffix of the images to
+ read, for example "jpg".
+\end_layout
+
+\begin_layout List
+\labelwidthstring 00.00.0000
+n_start integer(required) An integer specifying the first image in the sequence.
+\end_layout
+
+\begin_layout List
+\labelwidthstring 00.00.0000
+n_end integer(required) An integer specifying the last image in the sequence.
+\end_layout
+
+\begin_layout List
+\labelwidthstring 00.00.0000
+digits integer(required) An integer specifying the number of digits of the
+ number field in the sequence, for example image00001.jpg would require digits=5.
+\end_layout
+
+\begin_layout List
+\labelwidthstring 00.00.0000
+dims integer(optional) An integer array specifying the size of the output
+ image.
+ If not specified the size of the first image is used.
+\end_layout
+
+\begin_layout List
+\labelwidthstring 00.00.0000
+format integer(optional) An integer with values 0 for MPEG1 and 1 for MPEG2.
+ Default is MPEG1.
+\end_layout
+
+\begin_layout List
+\labelwidthstring 00.00.0000
+frame_rate integer(optional) An integer with values.
+\end_layout
+
+\begin_layout LyX-Code
+1 = 23.976 frames/sec: NTSC encapsulated film rate.
+\end_layout
+
+\begin_layout LyX-Code
+2 = 24 frames/sec: standard film rate.
+\end_layout
+
+\begin_layout LyX-Code
+3 = 25 frames/sec: PAL video frame rate
+\end_layout
+
+\begin_layout LyX-Code
+4 = 29.97 frames/sec: NTSC video frame rate
+\end_layout
+
+\begin_layout LyX-Code
+5 = 30 frames/sec: NTSC drop frame video rate (the default).
+\end_layout
+
+\begin_layout LyX-Code
+6 = 50 frames/sec: double frame rate/progressive PAL
+\end_layout
+
+\begin_layout LyX-Code
+7 = 59.94 frames/sec: double frame rate NTSC
+\end_layout
+
+\begin_layout LyX-Code
+8 = 60 frames/sec: double frame rate NTSC
+\end_layout
+
+\begin_layout List
+\labelwidthstring 00.00.0000
+mpeg_file string(optional) A string specifying the output MPEG file (default
+ outfile.mpg).
+\end_layout
+
+\begin_layout List
+\labelwidthstring 00.00.0000
+tmp_dir string(optional) A string specifying temporary directory to use
+ for the temporary image files.
+\end_layout
+
+\begin_layout Standard
+To run the IDL movie generator type in a shell:
+\end_layout
+
+\begin_layout LyX-Code
+%idl
+\end_layout
+
+\begin_layout LyX-Code
+> .compile mpg.pro
+\end_layout
+
+\begin_layout LyX-Code
+> make_mpg, prefix='image', suffix='JPG', n_start=0, n_end=100, digits=3,
+   dims=[512,512], frame_rate=2,mpeg_file='mympeg.mpg'
+\end_layout
+
 \begin_layout Part
 Appendices
 \end_layout

More information about the cig-commits mailing list