The Lorenz-96 model is a popular testbed for numerical weather
prediction algorithms, reproducing coarse features of the mid-latitude
atmosphere [1]}, [2]}. The
state at time \(t\)  is a 40-dimensional vector,
\(Z(t)=(Z_1(t),\ldots ,Z_{40}(t))\) .
The state dynamics are defined by a set of nonlinear ODEs that
represent the spatial discretization of a time-dependent PDE,
i.e.,
\( \frac{{\rm d}Z_j}{{\rm d}t} =(Z_{j+1} - Z_{j-2}) Z_{j-1} - Z_j + F,\qquad j=1,\ldots ,40,\) 
