Changes between Version 16 and Version 17 of Lorenz_96_model


Ignore:
Timestamp:
Sep 1, 2018, 12:20:43 PM (3 months ago)
Author:
lnerger
Comment:

--

Legend:

Unmodified
Added
Removed
Modified
  • Lorenz_96_model

    v16 v17  
    104104=== 4. Plot output from the assimilation experiments ===
    105105
    106 To display the output of the assimilation experiments we provide several plotting scripts for Matlab and Octave. To use them do
     106To 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
    107107{{{
    108108cd ../src/lorenz96/plotting/
     
    110110and see the file README there for a description of the available plotting scripts.
    111111
    112 `plot_example.m` plots the true and estimated mean RMS errors as a function of the
     112`plot_example` plots the true and estimated mean RMS errors as a function of the
    113113forgetting factor.
    114114
    115115The other plotting scripts require the specification
    116116of the directory and name of the file to be read. Sometimes,
    117 there are additional arguments like the time step index. Use
    118 'help' in Matlab to display the information about required input.
     117there are additional arguments like the time step index. In Matlab use
     118'help' to display the information about required input.
    119119
    120 Plotting examples:[[BR]]
     120'''!Matlab/Octave''' plotting examples:[[BR]]
    121121`plotobs('../../../bin/obs.nc',100)` plots the observation at time step 100[[BR]]
    122122`plotstate('../../../bin/t1_N30_f0.97.nc',100,'f')` plots the forecast state estimate at the 100th analysis step[[BR]]
     
    124124`plotstate('../../../bin/state.nc',1101)` plots the true state at model time step 1101 (= analysis step 100)[[BR]]
    125125
     126'''Python plotting'''
     127
     128For Python the scripts are provided by `plot_l96.py`.
     129
     130One 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
     134Alternatively the script can be run from the command line, providing the function
     135name 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`)
    126139
    127140