Yannos (Federica, 2/22/2005)
------
Yannos computes eigenfrequencies and eigenfunctions for a specified 1D model.


-> Run Yannos:
--------------
/data/08/federica/SCRIPTS/kernel/yannos_ucb_new.s $model_name

This script create an input file (yannos.dat) for both the spheroidal and toroidal mode case according to the parameters specified by the user in the first part of the script and then runs the new version of yannos (correct also when the switch for gravity is on) for both spheroidal and toroidal modes.

Federica, 9/27/05
In the directory /a/perry/data/26/yann/programmes/yannos you can find more information on yannos, in particular the source codes, makefiles and examples.
However the new version of Yannos currently used is /data/08/barbara/,codes04/,yannos/src/yannos1. This new version of Yannos allows the computation of the eigenfunction for a limited portion of the mantle. Therefore it has an additional line with the information on the portion of the mantle for which eigenfucntions are stored.


-> Run the reformatting files:
------------------------------
/data/08/federica/SCRIPTS/kernel/yannos_ST_format.s $model_name $max_freq
/data/08/federica/SCRIPTS/kernel/yannos_ST_format700km.s $model_name $max_freq

These scripts correct some mislabeling of the modes and set to zero the 
eigenfunctions of core and Stoneley modes which are not computed accurately by 
Yannos.
For planets other than the Earth, these scripts should be adapted (Mark Panning has a version for Europa).

The new version of yannos_ST_format800km.s (2/25/2005) is calling the new code yannos_ST_format800km_3_ellrot.f (instead of yannos_ST_format800km_3.f) which computes the ellipticity and rotational splitting parameters and stores them in the eigenfunction headers. For reading this reformatted eigenfunction files into other codes use the modified routine geteigys_ellrot.c (instead of geteigys.c).

Federica, 9/27/05
For the new version of Yannos (see above) use:
/data/08/federica/SCRIPTS/kernel/yannos_ST_format_rad.s
/data/08/federica/SCRIPTS/kernel/yannos_ST_format700km_rad.s

-> Ellipticity and rotational splitting parameters:
---------------------------------------------------
The ellipticity and rotational splitting parameters are computed in the subroutine ellcoup_rot (in /a/perry/data/08/federica/SRC/KERNEL/ANISOZ/ellgravy_rot.c). This routine takes as input two sets of eigenfunctions for two modes and the corresponding structure model and it computes and returns the ellipticity (ellterm) and rotational (rotterm) splitting parameters. 

To obtain the same ellipticity term as the one stored in the original XD eigenfuctions (e.g. prem.ell) and stored in the structure sprem_st (/data/07/gung/,eigfun/premy_st.h), the value return by the ellcoup_rot subroutine (ellterm) needs to be multiplied by some normalization factors (as e.g. done in the help routine ell). 
The value for the rotational slitting parameter returned by the ellcoup_rot subroutine (rotterm) has already the same normalization as the one stored in the original XD eigenfunctions.
The rotational splitting parameter is computed according to the formulas in Woodhouse & Dahlen (1978) and Dahlen (1968) and follows the normalization in Dahlen (1968). The detailed comments in the source code (ellgravy_rot.c) give more information on the used normalizations.