From: "Dave"
Subject: [H390-MTS] CSMP on MTS
Date: January 7, 2012 10:24:34 AM EST
To: "H390 MTS Group" <H390-MTS@yahoogroups.com>

Folks,
 
Well I am very happy with this. I now have my sample CSMP program from “Applied Numerical Methods for Digital Computation with Fortran and CSMP” 2nd Edition by James, Smith & Wolford working after one minor typo. In case any one wants to try the source is below :-
 
INITIAL
CONSTANT W=10.0,C=0.00259,A=-0.6,B=1.0,K1=.255,K2=1.025
INCON X0=0.0,XDOT0=8.0
     M=W/386.
     COEF=C/M
DYNAMIC
NOSORT
      IF(X.GT.0)GOTO 2
      F0FX=K1*DEADSP(A,0.,X)
      GOTO 3
  2   F0FX=K2*DEADSP(0.,B,X)
  3   CONTINUE
SORT
      XDDOT=-COEF*XDOT-F0FX/M
      XDOT=INTGRL(XDOT0,XDDOT)
      X=INTGRL(X0,XDOT)
TERMINAL
      KE=0.5*M*XDOT**2
      WRITE(6,4)KE
  4   FORMAT(' ',E16.7)
TIMER DELT=0.05,OUTDEL=0.2,FINTIM=60.
METHOD RKSFX
LABEL SPRINT MASS VISCOUSLY DAMPED WITH DEAD SPACE
PRTPLOT XDOT,X,F0FX
END
STOP
ENDJOB
 
Giving the first few lines of print out :-
 
SPRINT MASS VISCOUSLY DAMPED WITH DEAD SPACE                       PAGE   1
 
                        MINIMUM             XDOT   VERSUS TIME           MAXIMUM
                      -7.6828E+00                                      8.0000E+00
TIME          XDOT        I                                                I
0.0           8.0000E+00   -------------------------------------------------+
2.0000E-01    7.0260E+00   ----------------------------------------------+
4.0000E-01   -1.2600E+00   --------------------+
6.0000E-01   -7.6546E+00   +
8.0000E-01   -7.5684E+00   +
1.0000E+00   -6.4429E+00   ---+
1.2000E+00   -2.9557E+00   ---------------+
1.4000E+00    1.5774E+00   -----------------------------+
1.6000E+00    5.4257E+00   -----------------------------------------+
1.8000E+00    7.1513E+00   -----------------------------------------------+
2.0000E+00    7.0579E+00   ----------------------------------------------+
2.2000E+00    4.4652E+00   --------------------------------------+
2.4000E+00   -3.6665E+00   ------------+
2.6000E+00   -6.8152E+00   --+
2.8000E+00   -6.6794E+00   ---+
3.0000E+00   -5.2724E+00   -------+
3.2000E+00   -1.9047E+00   ------------------+
3.4000E+00    2.1148E+00   -------------------------------+
3.6000E+00    5.2571E+00   -----------------------------------------+
3.8000E+00    6.3542E+00   --------------------------------------------+
4.0000E+00    6.2285E+00   --------------------------------------------+
4.2000E+00    3.7196E+00   ------------------------------------+
4.4000E+00   -3.4711E+00   -------------+
4.6000E+00   -6.0134E+00   -----+
4.8000E+00   -5.8943E+00   -----+
5.0000E+00   -4.9065E+00   --------+
 
Thanks to every one from the MTS team for making this possible. Now I need to start updating my CSMP web pages....
 
Dave
G4UGM