Changes between Version 16 and Version 17 of Lorenz_96_model
 Timestamp:
 Sep 1, 2018, 12:20:43 PM (3 months ago)
Legend:
 Unmodified
 Added
 Removed
 Modified

Lorenz_96_model
v16 v17 104 104 === 4. Plot output from the assimilation experiments === 105 105 106 To display the output of the assimilation experiments we provide several plotting scripts for Matlab and Octave. To use them do106 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 107 107 {{{ 108 108 cd ../src/lorenz96/plotting/ … … 110 110 and see the file README there for a description of the available plotting scripts. 111 111 112 `plot_example .m` plots the true and estimated mean RMS errors as a function of the112 `plot_example` plots the true and estimated mean RMS errors as a function of the 113 113 forgetting factor. 114 114 115 115 The other plotting scripts require the specification 116 116 of the directory and name of the file to be read. Sometimes, 117 there are additional arguments like the time step index. Use118 'help' in Matlabto display the information about required input.117 there are additional arguments like the time step index. In Matlab use 118 'help' to display the information about required input. 119 119 120 Plotting examples:[[BR]]120 '''!Matlab/Octave''' plotting examples:[[BR]] 121 121 `plotobs('../../../bin/obs.nc',100)` plots the observation at time step 100[[BR]] 122 122 `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]] 125 125 126 '''Python plotting''' 127 128 For Python the scripts are provided by `plot_l96.py`. 129 130 One can either import the module plot_l96 and use its functions interactively, e.g[[BR]] 131 `>>> import plot_l96`[[BR]] 132 `>>> plot_l96.plot_obs('../../../bin/obs.nc', 4)`[[BR]] 133 134 Alternatively the script can be run from the command line, providing the function 135 name and its argument as command line parameters:[[BR]] 136 `./plot_l96.py plot_obs ../../../bin/obs.nc 4`[[BR]] 137 (If this fails you can also try to run the script as 138 `python plot_l96.py plot_obs ../../../bin/obs.nc 4`) 126 139 127 140