Changes between Version 16 and Version 17 of Lorenz_96_model

Timestamp:

Sep 1, 2018, 12:20:43 PM (6 years ago)

Author:

lnerger

Comment:

--

Lorenz_96_model

@@ -104 +104 @@
 === 4. Plot output from the assimilation experiments ===
 
-To display the output of the assimilation experiments we provide several plotting scripts for Matlab and Octave. To use them do
+To display the output of the assimilation experiments we provide plotting scripts for Matlab, Octave and (from PDAF 1.13.2) Python. To use them do
 {{{
 cd ../src/lorenz96/plotting/
@@ -110 +110 @@
 and see the file README there for a description of the available plotting scripts. 
 
-`plot_example.m` plots the true and estimated mean RMS errors as a function of the
+`plot_example` plots the true and estimated mean RMS errors as a function of the
 forgetting factor. 
 
 The other plotting scripts require the specification 
 of the directory and name of the file to be read. Sometimes, 
-there are additional arguments like the time step index. Use
-'help' in Matlab to display the information about required input. 
+there are additional arguments like the time step index. In Matlab use
+'help' to display the information about required input. 
 
-Plotting examples:[[BR]]
+'''!Matlab/Octave''' plotting examples:[[BR]]
 `plotobs('../../../bin/obs.nc',100)` plots the observation at time step 100[[BR]]
 `plotstate('../../../bin/t1_N30_f0.97.nc',100,'f')` plots the forecast state estimate at the 100th analysis step[[BR]]
@@ -124 +124 @@
 `plotstate('../../../bin/state.nc',1101)` plots the true state at model time step 1101 (= analysis step 100)[[BR]] 
 
+'''Python plotting'''
+
+For Python the scripts are provided by `plot_l96.py`.
+
+One can either import the module plot_l96 and use its functions interactively, e.g[[BR]]
+`>>> import plot_l96`[[BR]]
+`>>> plot_l96.plot_obs('../../../bin/obs.nc', 4)`[[BR]]
+
+Alternatively the script can be run from the command line, providing the function
+name and its argument as command line parameters:[[BR]]
+`./plot_l96.py plot_obs ../../../bin/obs.nc 4`[[BR]]
+(If this fails you can also try to run the script as
+`python plot_l96.py plot_obs ../../../bin/obs.nc 4`)