Dalhousie Geodynamics Group


Department of Oceanography
Dalhousie University
1355 Oxford Street
PO Box 15000
Halifax, NS, B3H 4R2
Canada




SOPALE

Sopale Nested

Processing SOPALE Output

Printing

MOZART

TMM (Thermo-mechanical Model)

old SOPALE documentation

Geodynamics home page

sopale_nested input

Introduction

The parameters here define what is written to the matlab output files. The files are written in the sub-directory "nameout1matlab_pnn", where nn is 00 for LS, and 01 for SS.

Some fields will only be written by the LS part of sopale_nested, and are flagged as such.

If any Lagrangian output is requested, the file "cell21<time-step>.ml" will also be written. cell21 contains the valid/invalid flags for Lagrangian particles.

The matlab output files DO NOT INCLUDE data for INJECTED Lagrangian particles. Only data for the grid Lagrangian particles is written. For output that does include injected lagrangian particles, use type 2 output files.

Section keywords:

matlab-free-format

[matlab-free-format]

Indicates start of section.

nout <value>

<value>
"matlab" output files will be written on the first, and every nout'th time step.

icolor1 <value>

<value>
value for icolor1. 0 or 1. See the fixed-format section for details.

And so on...

For any variable listed in the fixed-format section. The default for those not present is 0. The order of lines is not important.

end

terminates the free format matlab input section.

matlab

[matlab]

Indicates start of section.

In this section the input is fixed format; every line described below is required. In addition, order is crucial. The (perceived) advantage of free-format (described above) is that one only needs to include lines for the fields in which one is interested.

<ignored line>

ignored line
Can be used to comment the next line. There should be at least one item in this line..

nout

nout
"matlab" output files will be written on the first, and every nout'th time step.

<ignored line>

icolor1

icolor1
= 1 -   file "color1<time-step>.ml" will be written.
= anything else -   no file will be written for this parameter.

The mechanical material of each eulerian cell will be in the file. Same as record 16 of the type 1 output file. The format is binary integer*4.

<ignored line>

icolor2

icolor2
= 1 -   file "color2<time-step>.ml" will be written.
= anything else -   no file will be written for this parameter.
only effective in LS.

The mechanical material of each lagrangian grid particle will be in the file. Same as record 5 of the type 2 output file. The format is binary real*4.

<ignored line>

iecord

iecord
= 1 -   files "x1<time-step>.ml" and "y1<time-step>.ml" will be written.
= anything else -   no files will be written for this parameter.

The x and y coordinates of the eulerian grid will be in the files. Same as records 1 and 2 of the type 1 output file. Their formats are binary real*4.

<ignored line>

ilcord

ilcord
= 1 -   files "x2<time-step>.ml" and "y2<time-step>.ml" will be written.
= anything else -   no files will be written for this parameter.
only effective in LS.

The x and y coordinates of the lagrangian grid will be in the files. Same as records 1 and 2 of the type 2 output file. Their formats are binary real*4.

<ignored line>

istrain1

istrain1
= 1 -   file "strain1<time-step>.ml" will be written.
= anything else -   no file will be written for this parameter.

"strain1" for each eulerian node will be in the file. Same as record 18 of the type 1 output file. The format is binary real*4.

<ignored line>

istrain2

istrain2
= 1 -   file "strain2<time-step>.ml" will be written.
= anything else -   no file will be written for this parameter.
only effective in LS.

"strain2" for each lagrangian grid particle will be in the file. Same as record 7 of the type 2 output file. The format is binary real*4.

<ignored line>

istrainr

istrainr
= 1 -   file "e2d<time-step>.ml" will be written.
= anything else -   no file will be written for this parameter.

"f1_sr": strain rate 2nd invariant of the eulerian grid will be in the file. Same as record 13 of the type 1 output file. The format is binary real*4.

<ignored line>

itemper

itemper
= 1 -   file "temp1<time-step>.ml" will be written.
= anything else -   no file will be written for this parameter.

"t0" which is, I believe, the eulerian grid temperatures of the previous time step, will be in the file. The format is binary real*4.

<ignored line>

ievel

ievel
= 1 -   files "vx1<time-step>.ml" and "vy1<time-step>.ml" will be written.
= anything else -   no files will be written for this parameter.

"vx1" and "vy1", the velocities of the eulerian grid, will be in the files. Same as records 3 and 4 of the type 1 output file. The format is binary real*4.

<ignored line>

ilvel

ilvel
= 1 -   files "vx2<time-step>.ml" and "vy2<time-step>.ml" will be written.
= anything else -   no files will be written for this parameter.
only effective in LS.

"vx2" and "vy2", the velocities of the lagrangian grid, will be in the file. Same as records 3 and 4 of the type 2 output file. The format is binary real*4.

<ignored line>

isxx

isxx
= 1 -   file "sxx<time-step>.ml" will be written.
= anything else -   no file will be written for this parameter.

"sigxx" will be in the file. This is calculated in compute_viscosities. The format is binary real*4.

<ignored line>

isyy

isyy
= 1 -   file "syy<time-step>.ml" will be written.
= anything else -   no file will be written for this parameter.

"sigyy" will be in the file. This is calculated in compute_viscosities. The format is binary real*4.

<ignored line>

isxy

isxy
= 1 -   file "sxy<time-step>.ml" will be written.
= anything else -   no file will be written for this parameter.

"sigxy" will be in the file. This is calculated in compute_viscosities. The format is binary real*4.

<ignored line>

inodpres

inodpres
= 1 -   file "p<time-step>.ml" will be written.
= anything else -   no file will be written for this parameter.

"nodpress" will be in the file. Same as record 6 of the type 1 output file. The format is binary real*4.

<ignored line>

iefx1

iefx1
= 1 -   file "efx1<time-step>.ml" will be written.
= anything else -   no file will be written for this parameter.

"e_fx1" the horizontal strain rate for the eulerian grid, will be in the file. Same as record 14 of the type 1 output file. The format is binary real*4.

<ignored line>

iefy1

iefy1
= 1 -   file "efy1<time-step>.ml" will be written.
= anything else -   no file will be written for this parameter.

"e_fy1" the shear strain rate (not vertical strain rate, as one might conclude from the name) for the eulerian grid, will be in the file. Same as record 15 of the type 1 output file. The format is binary real*4.

<ignored line>

irpress

irpress
= 1 -   file "rpress<time-step>.ml" will be written.
= anything else -   no file will be written for this parameter.

"rpress" = epress1(i) - plitho(i). The data is elemental. The format is binary real*4.

<ignored line>

ilag_dens

Since version 2009-07-26_19:

ilag_dens
= 1 -   file "lagdens<time-step>.ml" will be written.
= anything else -   no file will be written for this parameter.
only effective in LS.

<ignored line>

iheat_flux

Since version 2009-07-26_22.7:

iheat_flux
= 1 -   files "heatflux_x<time-step>.ml" and "heatflux_y<time-step>.ml" will be written.
= anything else -   no file will be written for this parameter.

"heatflux" is the heat flux passing through each Eulerian cell. The format is binary real*4.

<ignored line>

ipresdiff

Since version 2009-07-26_29.0:

ipresdiff
= 1 -   file "presdiff<time-step>.ml" will be written.
= anything else -   no file will be written for this parameter.

"presdiff" = nodpres(i) - plithon(i). The data is nodal. The format is binary real*4. "plithon" is plitho calculated for each Eulerian node.

<ignored line>

iplitho

Since version 2012-05-02_1.0:

iplitho
= 1 -   file "plitho<time-step>.ml" will be written.
= anything else -   no file will be written for this parameter.

plitho is the lithostatic pressure, calculated for each cell. The data is elemental. The format is binary real*4.

<ignored line>

iplithon

Since version 2012-05-02_1.0:

iplithon
= 1 -   file "plithon<time-step>.ml" will be written.
= anything else -   no file will be written for this parameter.

plithon is the lithostatic pressure, calculated for each node. The data is nodal. The format is binary real*4.

<ignored line>

ivy1r

Since version 2013-01-22_1.0:

ivy1r
= 1 -   file "vy1r<time-step>.ml" will be written.
= anything else -   no file will be written for this parameter.

vy1r is viscosity, calculated for each node. The data is nodal. The format is binary real*4.



This page was last modified on Monday, 17-Mar-2014 12:06:11 ADT
Comments to geodynam at dal dot ca