Changeset 1426 for trunk/UTIL/SPECTRA/readme

Timestamp:

May 11, 2015, 4:25:55 PM (10 years ago)

Author:

milmd

Message:

UTIL/SPECTRA. Correction of some bugs with spectra average. Now tested on Ciclad.

File:

• trunk/UTIL/SPECTRA/readme (modified) (2 diffs)

trunk/UTIL/SPECTRA/readme

@@ -51 +51 @@
 FC=pgf90 or ifort
 
+for some machine, may load netcdf module:
+module load netcdf
+
+for ciclad machine, if using netcdf4:
+in makefile change line 7
+LDFLAGS=-L${netcdfpath}/lib -lnetcdf -L${spherepackpath}/lib -lspherepack
+and replace by
+LDFLAGS=-L${netcdfpath}/lib -lnetcdf -lnetcdff -L${spherepackpath}/lib -lspherepack 
+
 make spectra_analysis
 
@@ -56 +65 @@
 
 > compile test program
-make test_analysis
+make test_harmonic
 
 > create harmonic winds
-./test_analysis
+./test_harmonic
 
-> compute kinetic energy spectrum for one harmonic
+>>>>> some examples
+
+> compute kinetic energy spectrum for one harmonic test file
+
 ./spectra_analysis harmonic_64x48_lmdz_22.nc -alt none -time none -o harmonic_64x48_lmdz_22_spectra
 
+output file harmonic_64x48_lmdz_22_spectra contains:
+first column -> harmonic numbers 
+second column -> amplitude of harmonic component, here only component 22 is non zero.
+
+
+> compute kinetic energy spectrum with temporal and vertical mean
+
+./spectra_analysis diagfi.nc -t 5 -mt 10 -z 15 -mz 2 -o spectra
+
+compute kinetic energy spectrum for t=5,6,...,15 and z=15,16,17 and average all spectra obtained.
+
+
+> compute divergence and rotational part of kinetic energy spectrum
+
+./spectra_analysis diagfi.nc -t 5 -z 15 -o spectra -divrot
+
+in spectra file, will find the rotationnal and divergent part of the decomposition of the velocity on the vectorial spherical harmonic basis:
+#Spherical  diagfi.nc                                         
+#wavenumber t=  150 z=    1                                   
+#           spec_tot        spec_div        spec_rot          
+    0       0.000000E+00    0.000000E+00    0.000000E+00     
+    1       0.119347E+01    0.101938E+01    0.174095E+00     
+    2       0.530388E+00    0.250037E+00    0.280350E+00     
+    3       0.106973E+01    0.589848E+00    0.479880E+00     
+...
+
+
+>>>>> output interpretation
+
+One can plot velocity projection against spherical wavenumber with gnuplot for instance.
+In geostrophic turbulence, a n^-3 slope must appear for wavenumber n=1..100(?) and a n^-(5/3) slope for higher n.
+
+
+>>>>> reference
+J. N. Koshyk, 2001, The Horizontal Kinetic Energy Spectrum and Spectral Budget Simulated by a High-Resolution Troposphere–Stratosphere–Mesosphere GCM
+
+
+
+
+
+