Corrected a couple of small bugs that affect the Generalized Stability Theory (GST) algorithms for CPP options AFT_EIGENMODES, FORCING_SV, FT_EIGENMODES, OPT_PERTURBATION.
- Corrected the declaration of input/output arrays A and Aout in routine mp_gather_state (in distribute.F) to A(Nstr:Nend) and Aout(Asize), respectively. The declarion was reversed between these two arrays.
- Corrected a small bug in PARPACK routine pdnaupd.f in line 482. The
if ((ido .eq. 0) .or. (ido .eq. -2)) then
statement had an .eq. instead of .or.; I think that I introduced this bug when implementing checkpointing.
- Changed get_state.F so the initial time is set to dstart. This allows to use the zero-initial NetCDF with any starting time. Recall that in several adjoint applications zero-initial fields are used for the tangent linear and adjoint models.