I got some questions on how to plot the spectra. I am sure there are many ways to do this, but I use the method from Jean-Daniel Gril at Meteo-France using the "Ectoplasm" program which allows you to compute the spectral energy on ellipses referenced by integer ( between 0 to NMSMAX ). Then I usually plot the output using Matlab.

These are the steps:

1) Copy the binary "ALTO" to "ECTO" (ALTO is compiled by default in HARMONIE and should be in your bin/ directory).

2) Link to your history file:

ln -s $DIR1/$FILE infile${EXP}.fa ($FILE is ICMSHHARM+00${LL})

3) ./ECTO namlist infile${EXP}.fa (See below how to construct the namlist)

4) Plot wavenumbers 2:end, I usually plot using loglog plot in Matlab, e.g:

loglog((2:lenl),LINEAR(2:end,5),'b','linewidth',2);

* namlist:

&NAMNAM

C(1)='S***TEMPE???????'

L(:,1)=24,41,50

CKE(1)='S0**WIND.?.PHYS'

LKE(:,1)=24,41,50

CVD(1)='S0**WIND.?.PHYS'

LVD(:,1)=24,41,50

&END

With this namelist KE for lev 41 ends up in column 5, that is why I select "5" in the plot command

Here is the documentation from Jean-Daniel on how to construct the namelist:

www.cnrm.meteo.fr/gmapdoc/meshtml/ecto3.0