/* rocot_check.f -- translated by f2c (version 20100827).
   You must link the resulting object file with libf2c:
	on Microsoft Windows system, link with libf2c.lib;
	on Linux or Unix systems, link with .../path/to/libf2c.a -lm
	or, if you install libf2c.a in a standard place, with -lf2c -lm
	-- in that order, at the end of the command line, as in
		cc *.o -lf2c -lm
	Source for libf2c is in /netlib/f2c/libf2c.zip, e.g.,

		http://www.netlib.org/f2c/libf2c.zip
*/

#include "f2c.h"

/* Table of constant values */

static integer c__1 = 1;
static integer c__4 = 4;
static integer c__9 = 9;
static real c_b595 = 0.f;
static real c_b597 = 1.f;

/* Main program */ int MAIN__(void)
{
    /* Format strings */
    static char fmt_100[] = "(a,3i6)";
    static char fmt_200[] = "(1x,a,f3.1,2a,2x,a)";
    static char fmt_105[] = "(a,1x,i4,1x,i2.2,1x,i2.2)";
    static char fmt_106[] = "(a,3(1x,i2.2))";
    static char fmt_108[] = "(f12.5,1x,f8.3,1x,f8.3)";
    static char fmt_109[] = "(f8.3,1x,f8.3)";
    static char fmt_110[] = "(a,f8.4)";
    static char fmt_112[] = "(a,3f10.3)";
    static char fmt_130[] = "(ss,32x,\002X=\002,f9.4,5x,\002r   =  \002,f9"
	    ".4,/,32x,\002Y=\002,f9.4,5x,\002teta=  \002,f7.2,/,32x,\002Z"
	    "=\002,f9.4,5x,\002phi =  \002,f7.2)";
    static char fmt_140[] = "(/,a44,\002 angle=\002,f8.3,\002    (deg.)\002,"
	    ":/,44x,\002 ratio=\002,f14.5)";
    static char fmt_150[] = "(/,a47,\002 angle= \002,f9.3,\002    (deg.)\002"
	    ",:/,47x,\002 ratio= \002,f11.5)";
    static char fmt_111[] = "(10x,a,f8.4,a)";

    /* System generated locals */
    real r__1, r__2;

    /* Builtin functions */
    double acos(doublereal);
    integer s_wsfe(cilist *), do_fio(integer *, char *, ftnlen), e_wsfe(void),
	     s_rsle(cilist *), do_lio(integer *, integer *, char *, ftnlen), 
	    e_rsle(void), s_wsle(cilist *), e_wsle(void);
    /* Subroutine */ int s_stop(char *, ftnlen);

    /* Local variables */
    extern /* Subroutine */ int g_gsm_dipole_tilt_angle__(real *), 
	    cv_date_to_jul2000__(integer *, integer *, integer *, integer *), 
	    cv_jul2000_to_date__(integer *, integer *, integer *, integer *), 
	    cv_date_to_jul1950__(integer *, integer *, integer *, integer *), 
	    cv_jul1950_to_date__(integer *, integer *, integer *, integer *);
    static real a, d__, r__, a2, b1, b2, b3, g1, g2, g3;
    extern /* Subroutine */ int g_gei_geo_sun_dir__(real *, real *, real *, 
	    real *, real *, real *);
    static real r2, u1, v1, v2, v3, w1, w2, w3, u2, u3;
    extern /* Subroutine */ int g_gei_geo_sun_rot__(real *, real *, real *, 
	    real *, real *, real *);
    static integer ih;
    static real bp, ra, br;
    static integer im;
    static real gp, bt, pi, gr, gt;
    static integer is, nv;
    static real up, nx, ny, px, py, pz, nz, tx, ty, tz, vr, vt, vp, wr, wt, 
	    wp, ur, ut, xo, yo, zo, xs, xm, ym, zm;
    static integer ih2;
    static real sa1, sa2, sa3;
    static integer im2;
    extern /* Subroutine */ int cp_angle_and_ratio__(real *, real *, real *, 
	    real *, real *, real *, real *, real *);
    static real bs1, bs2, bs3;
    static integer is2;
    static real sap, sar, bsr, sat, bst, bsp;
    static integer ily;
    extern /* Subroutine */ int t_sr2_to_sr__(real *, real *, real *, real *, 
	    real *, real *, real *), t_sr_to_sr2__(real *, real *, real *, 
	    real *, real *, real *, real *), t_geo_to_dm__(real *, real *, 
	    real *, real *, real *, real *, real *, real *), t_dm_to_geo__(
	    real *, real *, real *, real *, real *, real *, real *, real *), 
	    t_gei_to_sm__(real *, real *, real *, real *, real *, real *), 
	    t_sm_to_gei__(real *, real *, real *, real *, real *, real *), 
	    t_sm_to_mag__(real *, real *, real *, real *, real *, real *), 
	    t_mag_to_sm__(real *, real *, real *, real *, real *, real *), 
	    t_geo_to_sm__(real *, real *, real *, real *, real *, real *), 
	    t_sm_to_geo__(real *, real *, real *, real *, real *, real *), 
	    t_gsm_to_sm__(real *, real *, real *, real *, real *, real *), 
	    t_sm_to_gsm__(real *, real *, real *, real *, real *, real *), 
	    t_sr2_to_mfa__(real *, real *, real *, real *, real *, real *, 
	    real *, real *, real *, real *, real *, real *);
    static real sdec;
    static integer jud00, iday, jud50;
    static real obli;
    extern /* Subroutine */ int t_gse_to_sr2__(real *, real *, real *, real *,
	     real *, real *, real *, real *, real *), t_sr2_to_gse__(real *, 
	    real *, real *, real *, real *, real *, real *, real *, real *);
    static real gmst;
    extern /* Subroutine */ int t_gei_to_mag__(real *, real *, real *, real *,
	     real *, real *), t_mag_to_gei__(real *, real *, real *, real *, 
	    real *, real *), cp_leap_year__(integer *, integer *), 
	    t_geo_to_mag__(real *, real *, real *, real *, real *, real *), 
	    t_geo_to_gei__(real *, real *, real *, real *, real *, real *), 
	    t_mag_to_geo__(real *, real *, real *, real *, real *, real *), 
	    t_gei_to_gse__(real *, real *, real *, real *, real *, real *), 
	    t_gse_to_gei__(real *, real *, real *, real *, real *, real *), 
	    t_gei_to_geo__(real *, real *, real *, real *, real *, real *), 
	    t_gse_to_mfa__(real *, real *, real *, real *, real *, real *, 
	    real *, real *, real *);
    static real d1gei;
    extern /* Subroutine */ int g_gei_geo_dipole_dir__(real *, real *, real *,
	     real *, real *, real *);
    static real d2gei, d3gei, e1gei, e2gei, d1geo, d2geo, d3geo;
    extern /* Subroutine */ int t_sph_to_car__(real *, real *, real *, real *,
	     real *, real *);
    static real e3gei, e1geo, e2geo, e3geo, r1gei;
    static integer iday2, iday3, iyea3, iyea4, iday4;
    static real s1gei, s2gei, s3gei, s1geo, s2geo, s3geo, r2gei;
    static integer imon2, imon3, imon4;
    static real r3gei, r1geo, r2geo, r3geo;
    extern /* Subroutine */ int t_car_to_sph__(real *, real *, real *, real *,
	     real *, real *), t_gei_to_gsm__(real *, real *, real *, real *, 
	    real *, real *), t_geo_to_gsm__(real *, real *, real *, real *, 
	    real *, real *), t_geo_to_gse__(real *, real *, real *, real *, 
	    real *, real *), t_gsm_to_gei__(real *, real *, real *, real *, 
	    real *, real *), t_gse_to_gsm__(real *, real *, real *, real *, 
	    real *, real *), t_gsm_to_gse__(real *, real *, real *, real *, 
	    real *, real *), t_gsm_to_geo__(real *, real *, real *, real *, 
	    real *, real *), t_gse_to_geo__(real *, real *, real *, real *, 
	    real *, real *), t_geo_to_vdh__(real *, real *, real *, real *, 
	    real *, real *, real *, real *), t_vdh_to_geo__(real *, real *, 
	    real *, real *, real *, real *, real *, real *), cp_tpn_param__(
	    real *, real *, real *, real *, real *, real *, real *, real *, 
	    real *, real *, real *, real *, real *), t_gse_to_tpn__(real *, 
	    real *, real *, real *, real *, real *, real *, real *, real *, 
	    real *), t_gsm_to_tpn__(real *, real *, real *, real *, real *, 
	    real *, real *, real *, real *, real *);
    static real dpgei, epgei, drgei, ergei, dtgei, etgei, dpgeo, epgeo, drgeo,
	     ergeo, dtgeo, etgeo, rpgei, spgei, rrgei;
    static integer iyear;
    static real grlat, srgei, stgei, rtgei, spgeo, srgeo, rrgeo, grlon, slong,
	     stgeo, rtgeo, rpgeo, srasn;
    static integer idoty;
    extern /* Subroutine */ int cp_time_param__(integer *, integer *, integer 
	    *, integer *, integer *, integer *), t_gei_to_gseq__(real *, real 
	    *, real *, real *, real *, real *), t_gseq_to_gei__(real *, real *
	    , real *, real *, real *, real *), t_geo_to_gseq__(real *, real *,
	     real *, real *, real *, real *), t_gseq_to_geo__(real *, real *, 
	    real *, real *, real *, real *), t_gseq_to_gse__(real *, real *, 
	    real *, real *, real *, real *), t_gse_to_gseq__(real *, real *, 
	    real *, real *, real *, real *), r_date__(integer *, integer *, 
	    integer *);
    static real deltat;
    extern /* Subroutine */ int r_time__(integer *, integer *, integer *);
    static real diptan, spipha;
    static char verdat[14];
    static real spifre, houday;
    static integer imonth;
    extern /* Subroutine */ int cv_hms_to_dech__(integer *, integer *, 
	    integer *, real *), cv_dech_to_hms__(real *, integer *, integer *,
	     integer *), strlen_(char *, integer *, ftnlen), 
	    g_gei_geo_ecliptic_dir__(real *, real *, real *, real *, real *, 
	    real *);
    static real vernum;
    extern /* Subroutine */ int g_rocot_version_number__(real *, char *, 
	    ftnlen), g_gei_sun_param__(real *, real *, real *, real *, real *)
	    , cv_date_to_doty__(integer *, integer *, integer *, integer *), 
	    cv_doty_to_date__(integer *, integer *, integer *, integer *);

    /* Fortran I/O blocks */
    static cilist io___6 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___7 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___8 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___9 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___10 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___11 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___12 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___13 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___14 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___15 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___16 = { 0, 6, 0, fmt_200, 0 };
    static cilist io___17 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___18 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___19 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___20 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___21 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___22 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___23 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___24 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___25 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___26 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___27 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___28 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___29 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___30 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___31 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___32 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___33 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___34 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___35 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___36 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___37 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___38 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___39 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___40 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___41 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___42 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___43 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___44 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___45 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___46 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___47 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___48 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___49 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___50 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___51 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___52 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___53 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___54 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___55 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___56 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___57 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___58 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___59 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___60 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___61 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___62 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___63 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___64 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___65 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___66 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___67 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___68 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___69 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___70 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___71 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___72 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___73 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___74 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___81 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___82 = { 0, 6, 0, fmt_105, 0 };
    static cilist io___83 = { 0, 6, 0, fmt_106, 0 };
    static cilist io___84 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___85 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___86 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___87 = { 0, 5, 0, 0, 0 };
    static cilist io___91 = { 0, 6, 0, fmt_108, 0 };
    static cilist io___95 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___96 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___97 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___98 = { 0, 5, 0, 0, 0 };
    static cilist io___101 = { 0, 6, 0, fmt_109, 0 };
    static cilist io___102 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___103 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___104 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___105 = { 0, 5, 0, 0, 0 };
    static cilist io___109 = { 0, 6, 0, fmt_108, 0 };
    static cilist io___113 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___114 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___115 = { 0, 6, 0, 0, 0 };
    static cilist io___116 = { 0, 5, 0, 0, 0 };
    static cilist io___119 = { 0, 6, 0, fmt_108, 0 };
    static cilist io___120 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___121 = { 0, 5, 0, 0, 0 };
    static cilist io___123 = { 0, 6, 0, fmt_108, 0 };
    static cilist io___124 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___125 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___126 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___127 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___128 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___129 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___130 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___131 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___132 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___133 = { 0, 6, 0, fmt_200, 0 };
    static cilist io___134 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___135 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___136 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___137 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___138 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___139 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___140 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___141 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___142 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___148 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___149 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___150 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___151 = { 0, 6, 0, fmt_110, 0 };
    static cilist io___152 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___153 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___154 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___155 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___167 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___168 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___169 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___170 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___171 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___172 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___173 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___174 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___175 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___176 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___177 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___178 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___179 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___234 = { 0, 6, 0, fmt_112, 0 };
    static cilist io___235 = { 0, 6, 0, fmt_112, 0 };
    static cilist io___236 = { 0, 6, 0, fmt_112, 0 };
    static cilist io___237 = { 0, 6, 0, fmt_112, 0 };
    static cilist io___238 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___239 = { 0, 6, 0, fmt_112, 0 };
    static cilist io___240 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___241 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___242 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___243 = { 0, 6, 0, fmt_130, 0 };
    static cilist io___244 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___245 = { 0, 6, 0, fmt_130, 0 };
    static cilist io___246 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___247 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___248 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___249 = { 0, 6, 0, fmt_130, 0 };
    static cilist io___250 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___251 = { 0, 6, 0, fmt_130, 0 };
    static cilist io___252 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___253 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___254 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___255 = { 0, 6, 0, fmt_130, 0 };
    static cilist io___256 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___257 = { 0, 6, 0, fmt_130, 0 };
    static cilist io___258 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___259 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___260 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___261 = { 0, 6, 0, fmt_130, 0 };
    static cilist io___262 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___263 = { 0, 6, 0, fmt_130, 0 };
    static cilist io___264 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___265 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___266 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___267 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___268 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___273 = { 0, 6, 0, fmt_140, 0 };
    static cilist io___274 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___275 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___276 = { 0, 6, 0, fmt_140, 0 };
    static cilist io___277 = { 0, 6, 0, fmt_140, 0 };
    static cilist io___278 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___279 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___280 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___281 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___282 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___283 = { 0, 6, 0, fmt_130, 0 };
    static cilist io___284 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___285 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___286 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___287 = { 0, 6, 0, fmt_130, 0 };
    static cilist io___288 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___289 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___290 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___297 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___298 = { 0, 6, 0, fmt_130, 0 };
    static cilist io___299 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___300 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___301 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___302 = { 0, 6, 0, fmt_130, 0 };
    static cilist io___303 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___304 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___305 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___306 = { 0, 6, 0, fmt_130, 0 };
    static cilist io___307 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___308 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___309 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___310 = { 0, 6, 0, fmt_130, 0 };
    static cilist io___311 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___312 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___313 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___314 = { 0, 6, 0, fmt_130, 0 };
    static cilist io___315 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___316 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___317 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___318 = { 0, 6, 0, fmt_130, 0 };
    static cilist io___319 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___320 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___321 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___322 = { 0, 6, 0, fmt_130, 0 };
    static cilist io___323 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___324 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___325 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___326 = { 0, 6, 0, fmt_130, 0 };
    static cilist io___327 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___328 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___329 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___330 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___331 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___332 = { 0, 6, 0, fmt_130, 0 };
    static cilist io___333 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___334 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___335 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___336 = { 0, 6, 0, fmt_130, 0 };
    static cilist io___337 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___338 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___339 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___340 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___341 = { 0, 6, 0, fmt_130, 0 };
    static cilist io___342 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___343 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___344 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___345 = { 0, 6, 0, fmt_130, 0 };
    static cilist io___346 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___347 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___348 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___349 = { 0, 6, 0, fmt_130, 0 };
    static cilist io___350 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___351 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___352 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___353 = { 0, 6, 0, fmt_130, 0 };
    static cilist io___354 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___355 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___356 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___357 = { 0, 6, 0, fmt_130, 0 };
    static cilist io___358 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___359 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___360 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___361 = { 0, 6, 0, fmt_130, 0 };
    static cilist io___362 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___363 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___364 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___365 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___366 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___367 = { 0, 6, 0, fmt_130, 0 };
    static cilist io___368 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___369 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___370 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___371 = { 0, 6, 0, fmt_130, 0 };
    static cilist io___372 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___373 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___374 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___375 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___376 = { 0, 6, 0, fmt_130, 0 };
    static cilist io___377 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___378 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___379 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___380 = { 0, 6, 0, fmt_130, 0 };
    static cilist io___381 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___382 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___383 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___384 = { 0, 6, 0, fmt_130, 0 };
    static cilist io___385 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___386 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___387 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___388 = { 0, 6, 0, fmt_130, 0 };
    static cilist io___389 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___390 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___391 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___392 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___393 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___394 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___395 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___396 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___397 = { 0, 6, 0, fmt_130, 0 };
    static cilist io___398 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___399 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___400 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___401 = { 0, 6, 0, fmt_130, 0 };
    static cilist io___402 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___403 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___404 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___405 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___406 = { 0, 6, 0, fmt_130, 0 };
    static cilist io___407 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___408 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___409 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___410 = { 0, 6, 0, fmt_130, 0 };
    static cilist io___411 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___412 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___413 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___414 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___415 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___416 = { 0, 6, 0, fmt_130, 0 };
    static cilist io___417 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___418 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___419 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___426 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___427 = { 0, 6, 0, fmt_130, 0 };
    static cilist io___428 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___429 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___430 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___431 = { 0, 6, 0, fmt_130, 0 };
    static cilist io___432 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___433 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___434 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___435 = { 0, 6, 0, fmt_130, 0 };
    static cilist io___436 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___437 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___438 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___439 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___440 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___441 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___442 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___443 = { 0, 6, 0, fmt_130, 0 };
    static cilist io___444 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___445 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___452 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___453 = { 0, 6, 0, fmt_130, 0 };
    static cilist io___454 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___455 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___456 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___457 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___458 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___459 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___460 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___461 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___462 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___463 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___464 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___465 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___466 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___467 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___468 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___469 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___470 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___471 = { 0, 6, 0, fmt_130, 0 };
    static cilist io___472 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___473 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___474 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___475 = { 0, 6, 0, fmt_130, 0 };
    static cilist io___477 = { 0, 6, 0, fmt_150, 0 };
    static cilist io___478 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___479 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___480 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___481 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___482 = { 0, 6, 0, fmt_130, 0 };
    static cilist io___483 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___484 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___485 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___486 = { 0, 6, 0, fmt_130, 0 };
    static cilist io___487 = { 0, 6, 0, fmt_150, 0 };
    static cilist io___488 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___489 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___490 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___491 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___492 = { 0, 6, 0, fmt_130, 0 };
    static cilist io___493 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___494 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___495 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___496 = { 0, 6, 0, fmt_130, 0 };
    static cilist io___497 = { 0, 6, 0, fmt_150, 0 };
    static cilist io___498 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___499 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___500 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___501 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___502 = { 0, 6, 0, fmt_130, 0 };
    static cilist io___503 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___504 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___505 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___506 = { 0, 6, 0, fmt_130, 0 };
    static cilist io___507 = { 0, 6, 0, fmt_150, 0 };
    static cilist io___508 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___509 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___510 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___511 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___512 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___513 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___514 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___515 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___516 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___517 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___518 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___519 = { 0, 6, 0, fmt_130, 0 };
    static cilist io___520 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___521 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___522 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___523 = { 0, 6, 0, fmt_130, 0 };
    static cilist io___524 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___525 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___526 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___527 = { 0, 6, 0, fmt_130, 0 };
    static cilist io___528 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___529 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___530 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___531 = { 0, 6, 0, fmt_130, 0 };
    static cilist io___532 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___533 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___534 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___535 = { 0, 6, 0, fmt_130, 0 };
    static cilist io___536 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___537 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___538 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___539 = { 0, 6, 0, fmt_130, 0 };
    static cilist io___540 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___541 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___542 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___543 = { 0, 6, 0, fmt_130, 0 };
    static cilist io___544 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___545 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___546 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___547 = { 0, 6, 0, fmt_130, 0 };
    static cilist io___548 = { 0, 6, 0, fmt_150, 0 };
    static cilist io___549 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___550 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___551 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___552 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___553 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___554 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___555 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___556 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___557 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___558 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___559 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___560 = { 0, 6, 0, fmt_130, 0 };
    static cilist io___561 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___562 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___563 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___564 = { 0, 6, 0, fmt_130, 0 };
    static cilist io___565 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___566 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___567 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___568 = { 0, 6, 0, fmt_130, 0 };
    static cilist io___569 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___570 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___571 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___572 = { 0, 6, 0, fmt_130, 0 };
    static cilist io___573 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___574 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___575 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___576 = { 0, 6, 0, fmt_130, 0 };
    static cilist io___577 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___578 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___579 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___580 = { 0, 6, 0, fmt_130, 0 };
    static cilist io___581 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___582 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___583 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___584 = { 0, 6, 0, fmt_130, 0 };
    static cilist io___585 = { 0, 6, 0, fmt_150, 0 };
    static cilist io___586 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___587 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___588 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___589 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___590 = { 0, 6, 0, fmt_130, 0 };
    static cilist io___591 = { 0, 6, 0, fmt_150, 0 };
    static cilist io___592 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___593 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___594 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___595 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___596 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___597 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___598 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___599 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___600 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___601 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___602 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___603 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___604 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___605 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___606 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___607 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___608 = { 0, 6, 0, fmt_130, 0 };
    static cilist io___609 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___610 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___611 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___612 = { 0, 6, 0, fmt_130, 0 };
    static cilist io___613 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___614 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___615 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___616 = { 0, 6, 0, fmt_130, 0 };
    static cilist io___617 = { 0, 6, 0, fmt_150, 0 };
    static cilist io___618 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___619 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___620 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___621 = { 0, 6, 0, fmt_130, 0 };
    static cilist io___622 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___623 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___624 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___625 = { 0, 6, 0, fmt_130, 0 };
    static cilist io___626 = { 0, 6, 0, fmt_150, 0 };
    static cilist io___627 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___628 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___629 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___630 = { 0, 6, 0, fmt_130, 0 };
    static cilist io___631 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___632 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___633 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___634 = { 0, 6, 0, fmt_130, 0 };
    static cilist io___635 = { 0, 6, 0, fmt_150, 0 };
    static cilist io___636 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___637 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___638 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___639 = { 0, 6, 0, fmt_130, 0 };
    static cilist io___640 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___641 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___642 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___643 = { 0, 6, 0, fmt_130, 0 };
    static cilist io___644 = { 0, 6, 0, fmt_150, 0 };
    static cilist io___645 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___646 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___647 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___648 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___649 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___650 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___651 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___652 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___653 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___654 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___655 = { 0, 6, 0, fmt_112, 0 };
    static cilist io___656 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___657 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___658 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___659 = { 0, 6, 0, fmt_130, 0 };
    static cilist io___660 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___661 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___662 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___663 = { 0, 6, 0, fmt_130, 0 };
    static cilist io___664 = { 0, 6, 0, fmt_150, 0 };
    static cilist io___665 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___666 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___667 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___668 = { 0, 6, 0, fmt_130, 0 };
    static cilist io___669 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___670 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___671 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___672 = { 0, 6, 0, fmt_130, 0 };
    static cilist io___673 = { 0, 6, 0, fmt_150, 0 };
    static cilist io___674 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___675 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___676 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___677 = { 0, 6, 0, fmt_112, 0 };
    static cilist io___678 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___679 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___680 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___681 = { 0, 6, 0, fmt_130, 0 };
    static cilist io___682 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___683 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___684 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___685 = { 0, 6, 0, fmt_130, 0 };
    static cilist io___686 = { 0, 6, 0, fmt_150, 0 };
    static cilist io___687 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___688 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___689 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___690 = { 0, 6, 0, fmt_130, 0 };
    static cilist io___691 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___692 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___693 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___694 = { 0, 6, 0, fmt_130, 0 };
    static cilist io___695 = { 0, 6, 0, fmt_150, 0 };
    static cilist io___696 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___697 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___698 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___699 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___700 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___701 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___702 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___703 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___704 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___705 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___706 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___707 = { 0, 6, 0, fmt_130, 0 };
    static cilist io___708 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___709 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___710 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___711 = { 0, 6, 0, fmt_130, 0 };
    static cilist io___712 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___713 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___714 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___715 = { 0, 6, 0, fmt_111, 0 };
    static cilist io___716 = { 0, 6, 0, fmt_111, 0 };
    static cilist io___717 = { 0, 6, 0, fmt_111, 0 };
    static cilist io___718 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___719 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___720 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___721 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___722 = { 0, 6, 0, fmt_130, 0 };
    static cilist io___723 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___724 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___725 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___726 = { 0, 6, 0, fmt_130, 0 };
    static cilist io___727 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___728 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___729 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___730 = { 0, 6, 0, fmt_130, 0 };
    static cilist io___731 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___732 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___733 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___734 = { 0, 6, 0, fmt_130, 0 };
    static cilist io___735 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___736 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___737 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___738 = { 0, 6, 0, fmt_130, 0 };
    static cilist io___739 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___740 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___741 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___742 = { 0, 6, 0, fmt_130, 0 };
    static cilist io___743 = { 0, 6, 0, fmt_150, 0 };
    static cilist io___744 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___745 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___746 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___747 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___748 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___749 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___750 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___751 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___752 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___753 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___754 = { 0, 6, 0, fmt_130, 0 };
    static cilist io___755 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___756 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___763 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___764 = { 0, 6, 0, fmt_130, 0 };
    static cilist io___765 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___766 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___767 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___768 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___769 = { 0, 6, 0, fmt_130, 0 };
    static cilist io___770 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___771 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___772 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___773 = { 0, 6, 0, fmt_130, 0 };
    static cilist io___774 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___775 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___782 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___783 = { 0, 6, 0, fmt_130, 0 };
    static cilist io___784 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___785 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___786 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___787 = { 0, 6, 0, fmt_130, 0 };
    static cilist io___788 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___789 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___790 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___791 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___792 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___793 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___794 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___795 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___796 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___797 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___798 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___799 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___803 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___804 = { 0, 6, 0, fmt_112, 0 };
    static cilist io___806 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___807 = { 0, 6, 0, fmt_112, 0 };
    static cilist io___817 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___818 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___819 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___820 = { 0, 6, 0, fmt_112, 0 };
    static cilist io___821 = { 0, 6, 0, fmt_112, 0 };
    static cilist io___822 = { 0, 6, 0, fmt_112, 0 };
    static cilist io___823 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___824 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___825 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___826 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___827 = { 0, 6, 0, fmt_130, 0 };
    static cilist io___828 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___829 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___830 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___831 = { 0, 6, 0, fmt_130, 0 };
    static cilist io___832 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___833 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___834 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___835 = { 0, 6, 0, fmt_130, 0 };
    static cilist io___836 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___837 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___841 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___842 = { 0, 6, 0, fmt_112, 0 };
    static cilist io___843 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___844 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___845 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___846 = { 0, 6, 0, fmt_130, 0 };
    static cilist io___847 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___848 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___849 = { 0, 6, 0, fmt_100, 0 };
    static cilist io___850 = { 0, 6, 0, fmt_100, 0 };



/*     ***************************************************************0** */

/* *   Rocot_Check program:  check subroutines of the rocotlib library; */

/*     Ask date and time and check calendar conversions and coordinate */
/*     transformations by walking across all path */
/*                       ___________________ */

/*     For details, see schematic diagram of section III the orginal */
/*     document untitled: */
/*    "CLUSTER Software Tools, part I: Coordinate Transformation Library" */
/*                     Document de travail DT/CRPE/1231 */
/*                  Patrick Robert, CRPE/TID, Juillet 1993 */
/*                       ___________________ */

/*     Adaptation for successives version of rocotlib: */

/*               version 1.0, November 1992 */
/*               version 1.1, July     1993 */
/*               version 1.2, January  1995 */
/*               version 1.3, July     2000 */
/*               version 1.4, June     2001 */
/*               version 1.5, December 2001 */
/*               version 1.6, Juin     2002 (upgrade IGRF -> 2005) */
/*               version 1.7, December 2002 (Version for CDPP) */
/*               version 1.8, November 2003 (adding t_sr2_to_sr tests) */
/*               version 1.9, March    2004 (compatibility  V 1.9 IDL) */
/*               version 2.0, May      2016 (test GSE->MFA et TPN) */

/*    Copyright 1992, Patrick ROBERT, CNRS-ESA, All Rights reserved */

/*     ***************************************************************0** */


    pi = acos(-1.f);
    d__ = 180.f / pi;


    g_rocot_version_number__(&vernum, verdat, (ftnlen)14);
    strlen_(verdat, &nv, (ftnlen)14);

    s_wsfe(&io___6);
    do_fio(&c__1, " ", (ftnlen)1);
    e_wsfe();
    s_wsfe(&io___7);
    do_fio(&c__1, "XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX", (
	    ftnlen)53);
    e_wsfe();
    s_wsfe(&io___8);
    do_fio(&c__1, " ", (ftnlen)1);
    e_wsfe();
    s_wsfe(&io___9);
    do_fio(&c__1, "     Coordinates Transformation Library  ROCOTLIB", (
	    ftnlen)49);
    e_wsfe();
    s_wsfe(&io___10);
    do_fio(&c__1, "                  ___________________", (ftnlen)37);
    e_wsfe();
    s_wsfe(&io___11);
    do_fio(&c__1, " ", (ftnlen)1);
    e_wsfe();
    s_wsfe(&io___12);
    do_fio(&c__1, " Rocotche program: Check transformations of Rocotlib", (
	    ftnlen)52);
    e_wsfe();
    s_wsfe(&io___13);
    do_fio(&c__1, " ", (ftnlen)1);
    e_wsfe();
    s_wsfe(&io___14);
    do_fio(&c__1, "          Patrick ROBERT, CRPE, November 1992", (ftnlen)45)
	    ;
    e_wsfe();
    s_wsfe(&io___15);
    do_fio(&c__1, " ", (ftnlen)1);
    e_wsfe();
    s_wsfe(&io___16);
    do_fio(&c__1, "         Revised Version ", (ftnlen)25);
    do_fio(&c__1, (char *)&vernum, (ftnlen)sizeof(real));
    do_fio(&c__1, " - ", (ftnlen)3);
    do_fio(&c__1, verdat, nv);
    e_wsfe();
    s_wsfe(&io___17);
    do_fio(&c__1, " ", (ftnlen)1);
    e_wsfe();
    s_wsfe(&io___18);
    do_fio(&c__1, "XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX", (
	    ftnlen)53);
    e_wsfe();
    s_wsfe(&io___19);
    e_wsfe();
    s_wsfe(&io___20);
    e_wsfe();
    s_wsfe(&io___21);
    e_wsfe();
    s_wsfe(&io___22);
    do_fio(&c__1, "SUMMARY OF CHECK OPERATIONS :", (ftnlen)29);
    e_wsfe();
    s_wsfe(&io___23);
    do_fio(&c__1, "***************************", (ftnlen)27);
    e_wsfe();
    s_wsfe(&io___24);
    e_wsfe();
    s_wsfe(&io___25);
    do_fio(&c__1, "  I) TEST OF ROCOTLIB CALENDAR CONVERSIONS", (ftnlen)42);
    e_wsfe();
    s_wsfe(&io___26);
    e_wsfe();
    s_wsfe(&io___27);
    do_fio(&c__1, " II) CHECK BASIC TIME PARAMETERS", (ftnlen)32);
    e_wsfe();
    s_wsfe(&io___28);
    e_wsfe();
    s_wsfe(&io___29);
    do_fio(&c__1, "III) CHECK BASIC TRANSFORMATIONS", (ftnlen)32);
    e_wsfe();
    s_wsfe(&io___30);
    e_wsfe();
    s_wsfe(&io___31);
    do_fio(&c__1, " IV) TEST OF ROCOTLIB TRANSFORMATION SUBROUTINES", (ftnlen)
	    48);
    e_wsfe();
    s_wsfe(&io___32);
    e_wsfe();
    s_wsfe(&io___33);
    do_fio(&c__1, "    A) STAR TRANSFORMATIONS AROUND  GEI  SYSTEM", (ftnlen)
	    47);
    e_wsfe();
    s_wsfe(&io___34);
    e_wsfe();
    s_wsfe(&io___35);
    do_fio(&c__1, "       GEI -----------------------------> MAG", (ftnlen)45)
	    ;
    e_wsfe();
    s_wsfe(&io___36);
    do_fio(&c__1, "       GEI ------------------------> SM", (ftnlen)39);
    e_wsfe();
    s_wsfe(&io___37);
    do_fio(&c__1, "       GEI ------------------> GSM", (ftnlen)34);
    e_wsfe();
    s_wsfe(&io___38);
    do_fio(&c__1, "       GEI ------------> GSE", (ftnlen)28);
    e_wsfe();
    s_wsfe(&io___39);
    e_wsfe();
    s_wsfe(&io___40);
    do_fio(&c__1, "    B) COUNTERCLOCKWISE RING CUMULATIVE TRANSFORM.", (
	    ftnlen)50);
    e_wsfe();
    s_wsfe(&io___41);
    e_wsfe();
    s_wsfe(&io___42);
    do_fio(&c__1, "       GEO ->GEI ->GSQ ->GSE ->GSM ->SM ->MAG ->GEO", (
	    ftnlen)51);
    e_wsfe();
    s_wsfe(&io___43);
    e_wsfe();
    s_wsfe(&io___44);
    do_fio(&c__1, "    C) CLOCKWISE RING CUMULATIVE TRANSFORMATIONS", (ftnlen)
	    48);
    e_wsfe();
    s_wsfe(&io___45);
    e_wsfe();
    s_wsfe(&io___46);
    do_fio(&c__1, "       GEO ->MAG ->SM ->GSM ->GSE ->GSQ ->GEI ->GEO", (
	    ftnlen)51);
    e_wsfe();
    s_wsfe(&io___47);
    e_wsfe();
    s_wsfe(&io___48);
    do_fio(&c__1, "    D) STAR TRANSFORMATIONS AROUND  GEO  SYSTEM", (ftnlen)
	    47);
    e_wsfe();
    s_wsfe(&io___49);
    e_wsfe();
    s_wsfe(&io___50);
    do_fio(&c__1, "       GEO ------> SM", (ftnlen)21);
    e_wsfe();
    s_wsfe(&io___51);
    do_fio(&c__1, "       GEO -----------> GSM", (ftnlen)27);
    e_wsfe();
    s_wsfe(&io___52);
    do_fio(&c__1, "       GEO -----------------> GSE", (ftnlen)33);
    e_wsfe();
    s_wsfe(&io___53);
    do_fio(&c__1, "       GEO -----------------------> GSQ", (ftnlen)39);
    e_wsfe();
    s_wsfe(&io___54);
    e_wsfe();
    s_wsfe(&io___55);
    do_fio(&c__1, "    E) LOCAL SYSTEMS", (ftnlen)20);
    e_wsfe();
    s_wsfe(&io___56);
    e_wsfe();
    s_wsfe(&io___57);
    do_fio(&c__1, "       GEO ----> DM", (ftnlen)19);
    e_wsfe();
    s_wsfe(&io___58);
    do_fio(&c__1, "       GEO ----> VDH", (ftnlen)20);
    e_wsfe();
    s_wsfe(&io___59);
    e_wsfe();
    s_wsfe(&io___60);
    do_fio(&c__1, "    F) SPACECRAFT SYSTEMS", (ftnlen)25);
    e_wsfe();
    s_wsfe(&io___61);
    e_wsfe();
    s_wsfe(&io___62);
    do_fio(&c__1, "       GEO -> GSE -> SR2 -> SR -> SR2 -> GSE -> GEO", (
	    ftnlen)51);
    e_wsfe();
    s_wsfe(&io___63);
    e_wsfe();
    s_wsfe(&io___64);
    do_fio(&c__1, "    G) PSEUDO-MAGNETOPAUSE SYSTEMS", (ftnlen)34);
    e_wsfe();
    s_wsfe(&io___65);
    e_wsfe();
    s_wsfe(&io___66);
    do_fio(&c__1, "       GSE ----> TPN", (ftnlen)20);
    e_wsfe();
    s_wsfe(&io___67);
    do_fio(&c__1, "       GSM ----> TPN", (ftnlen)20);
    e_wsfe();

    s_wsfe(&io___68);
    do_fio(&c__1, "\f", (ftnlen)1);
    e_wsfe();
    s_wsfe(&io___69);
    e_wsfe();
    s_wsfe(&io___70);
    e_wsfe();
    s_wsfe(&io___71);
    do_fio(&c__1, "INPUT DATA TEST :", (ftnlen)17);
    e_wsfe();
    s_wsfe(&io___72);
    do_fio(&c__1, "***************", (ftnlen)15);
    e_wsfe();
    s_wsfe(&io___73);
    e_wsfe();
    s_wsfe(&io___74);
    do_fio(&c__1, "INPUT DATE AND TIME:", (ftnlen)20);
    e_wsfe();

    r_date__(&iyear, &imonth, &iday);
    r_time__(&ih, &im, &is);

    s_wsfe(&io___81);
    e_wsfe();
    s_wsfe(&io___82);
    do_fio(&c__1, "iyear,imonth,iday: ", (ftnlen)19);
    do_fio(&c__1, (char *)&iyear, (ftnlen)sizeof(integer));
    do_fio(&c__1, (char *)&imonth, (ftnlen)sizeof(integer));
    do_fio(&c__1, (char *)&iday, (ftnlen)sizeof(integer));
    e_wsfe();
    s_wsfe(&io___83);
    do_fio(&c__1, "ih, im, is       : ", (ftnlen)19);
    do_fio(&c__1, (char *)&ih, (ftnlen)sizeof(integer));
    do_fio(&c__1, (char *)&im, (ftnlen)sizeof(integer));
    do_fio(&c__1, (char *)&is, (ftnlen)sizeof(integer));
    e_wsfe();

    s_wsfe(&io___84);
    e_wsfe();
    s_wsfe(&io___85);
    do_fio(&c__1, "INPUT VECTOR TEST:", (ftnlen)18);
    e_wsfe();

    s_wsfe(&io___86);
    do_fio(&c__1, "r, teta,phi of any vector in GEO system ? (deg.)", (ftnlen)
	    48);
    e_wsfe();
    s_rsle(&io___87);
    do_lio(&c__4, &c__1, (char *)&gr, (ftnlen)sizeof(real));
    do_lio(&c__4, &c__1, (char *)&gt, (ftnlen)sizeof(real));
    do_lio(&c__4, &c__1, (char *)&gp, (ftnlen)sizeof(real));
    e_rsle();
    s_wsfe(&io___91);
    do_fio(&c__1, (char *)&gr, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&gt, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&gp, (ftnlen)sizeof(real));
    e_wsfe();
    gt /= d__;
    gp /= d__;
    t_sph_to_car__(&gr, &gt, &gp, &g1, &g2, &g3);

    s_wsfe(&io___95);
    e_wsfe();
    s_wsfe(&io___96);
    do_fio(&c__1, "INPUT GEOG. LAT AND LONG. FOR LOCAL SYSTEMS TEST:", (
	    ftnlen)49);
    e_wsfe();
    s_wsfe(&io___97);
    do_fio(&c__1, "lat. and long. of local direct. of observation ? (d.)", (
	    ftnlen)53);
    e_wsfe();
    s_rsle(&io___98);
    do_lio(&c__4, &c__1, (char *)&grlat, (ftnlen)sizeof(real));
    do_lio(&c__4, &c__1, (char *)&grlon, (ftnlen)sizeof(real));
    e_rsle();
    s_wsfe(&io___101);
    do_fio(&c__1, (char *)&grlat, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&grlon, (ftnlen)sizeof(real));
    e_wsfe();
    grlat /= d__;
    grlon /= d__;

    s_wsfe(&io___102);
    e_wsfe();
    s_wsfe(&io___103);
    do_fio(&c__1, "INPUT SPIN DIRECTION for SR and SR2 SYSTEMS:", (ftnlen)44);
    e_wsfe();
    s_wsfe(&io___104);
    do_fio(&c__1, "r, teta,phi of spin axis in GSE system ? (deg.)", (ftnlen)
	    47);
    e_wsfe();
    s_rsle(&io___105);
    do_lio(&c__4, &c__1, (char *)&sar, (ftnlen)sizeof(real));
    do_lio(&c__4, &c__1, (char *)&sat, (ftnlen)sizeof(real));
    do_lio(&c__4, &c__1, (char *)&sap, (ftnlen)sizeof(real));
    e_rsle();
    s_wsfe(&io___109);
    do_fio(&c__1, (char *)&sar, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&sat, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&sap, (ftnlen)sizeof(real));
    e_wsfe();
    sat /= d__;
    sap /= d__;
    t_sph_to_car__(&sar, &sat, &sap, &sa1, &sa2, &sa3);

    s_wsfe(&io___113);
    e_wsfe();
    s_wsfe(&io___114);
    do_fio(&c__1, "SPIN FREQUENCY, SPIN PHASE & Deltat T :", (ftnlen)39);
    e_wsfe();
    s_wsle(&io___115);
    do_lio(&c__9, &c__1, "Spin frequency, Spin phase ? (Hz, deg.)", (ftnlen)
	    39);
    e_wsle();
    s_rsle(&io___116);
    do_lio(&c__4, &c__1, (char *)&spifre, (ftnlen)sizeof(real));
    do_lio(&c__4, &c__1, (char *)&spipha, (ftnlen)sizeof(real));
    e_rsle();
    s_wsfe(&io___119);
    do_fio(&c__1, (char *)&spifre, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&spipha, (ftnlen)sizeof(real));
    e_wsfe();
    s_wsfe(&io___120);
    do_fio(&c__1, "delta T ? (sec. from time of spin phase)", (ftnlen)40);
    e_wsfe();
    s_rsle(&io___121);
    do_lio(&c__4, &c__1, (char *)&deltat, (ftnlen)sizeof(real));
    e_rsle();
    s_wsfe(&io___123);
    do_fio(&c__1, (char *)&deltat, (ftnlen)sizeof(real));
    e_wsfe();
    spipha /= d__;

    s_wsfe(&io___124);
    e_wsfe();
    s_wsfe(&io___125);
    e_wsfe();
    s_wsfe(&io___126);
    e_wsfe();
    s_wsfe(&io___127);
    do_fio(&c__1, "RUN OF CHECK PROGRAM", (ftnlen)20);
    e_wsfe();
    s_wsfe(&io___128);
    do_fio(&c__1, "********************", (ftnlen)20);
    e_wsfe();
    s_wsfe(&io___129);
    e_wsfe();
    s_wsfe(&io___130);
    do_fio(&c__1, "Please wait...", (ftnlen)14);
    e_wsfe();
    s_wsfe(&io___131);
    do_fio(&c__1, "\f", (ftnlen)1);
    e_wsfe();

    s_wsfe(&io___132);
    e_wsfe();
    s_wsfe(&io___133);
    do_fio(&c__1, "            T E S T   O F   R O C O T L I B  Version ", (
	    ftnlen)53);
    do_fio(&c__1, (char *)&vernum, (ftnlen)sizeof(real));
    e_wsfe();
    s_wsfe(&io___134);
    e_wsfe();
    s_wsfe(&io___135);
    e_wsfe();
    s_wsfe(&io___136);
    do_fio(&c__1, "I) TEST OF ROCOTLIB CALENDAR CONVERSIONS", (ftnlen)40);
    e_wsfe();
    s_wsfe(&io___137);
    do_fio(&c__1, "   =====================================", (ftnlen)40);
    e_wsfe();
    s_wsfe(&io___138);
    e_wsfe();
    s_wsfe(&io___139);
    do_fio(&c__1, "INPUT DATE AND TIME:", (ftnlen)20);
    e_wsfe();
    s_wsfe(&io___140);
    e_wsfe();
    s_wsfe(&io___141);
    do_fio(&c__1, "year, month,  day   :     ", (ftnlen)26);
    do_fio(&c__1, (char *)&iyear, (ftnlen)sizeof(integer));
    do_fio(&c__1, (char *)&imonth, (ftnlen)sizeof(integer));
    do_fio(&c__1, (char *)&iday, (ftnlen)sizeof(integer));
    e_wsfe();
    s_wsfe(&io___142);
    do_fio(&c__1, "hour, minute, second:     ", (ftnlen)26);
    do_fio(&c__1, (char *)&ih, (ftnlen)sizeof(integer));
    do_fio(&c__1, (char *)&im, (ftnlen)sizeof(integer));
    do_fio(&c__1, (char *)&is, (ftnlen)sizeof(integer));
    e_wsfe();

    cv_hms_to_dech__(&ih, &im, &is, &houday);
    cv_date_to_doty__(&iyear, &imonth, &iday, &idoty);
    cv_date_to_jul1950__(&iyear, &imonth, &iday, &jud50);
    cv_date_to_jul2000__(&iyear, &imonth, &iday, &jud00);
    cp_leap_year__(&iyear, &ily);

    s_wsfe(&io___148);
    e_wsfe();
    s_wsfe(&io___149);
    do_fio(&c__1, "computed values:", (ftnlen)16);
    e_wsfe();
    s_wsfe(&io___150);
    e_wsfe();
    s_wsfe(&io___151);
    do_fio(&c__1, "decimal hour of the day:  ", (ftnlen)26);
    do_fio(&c__1, (char *)&houday, (ftnlen)sizeof(real));
    e_wsfe();
    s_wsfe(&io___152);
    do_fio(&c__1, "day of the year:          ", (ftnlen)26);
    do_fio(&c__1, (char *)&idoty, (ftnlen)sizeof(integer));
    e_wsfe();
    s_wsfe(&io___153);
    do_fio(&c__1, "Julian day from 1-1-1950: ", (ftnlen)26);
    do_fio(&c__1, (char *)&jud50, (ftnlen)sizeof(integer));
    e_wsfe();
    s_wsfe(&io___154);
    do_fio(&c__1, "Julian day from 1-1-2000: ", (ftnlen)26);
    do_fio(&c__1, (char *)&jud00, (ftnlen)sizeof(integer));
    e_wsfe();
    s_wsfe(&io___155);
    do_fio(&c__1, "leap year (1=yes,0=no):   ", (ftnlen)26);
    do_fio(&c__1, (char *)&ily, (ftnlen)sizeof(integer));
    e_wsfe();

    cv_doty_to_date__(&idoty, &iyear, &imon2, &iday2);
    cv_jul1950_to_date__(&jud50, &iyea3, &imon3, &iday3);
    cv_jul2000_to_date__(&jud00, &iyea4, &imon4, &iday4);
    cv_dech_to_hms__(&houday, &ih2, &im2, &is2);

    s_wsfe(&io___167);
    e_wsfe();
    s_wsfe(&io___168);
    do_fio(&c__1, "recompute date from year and day of the year:", (ftnlen)45)
	    ;
    do_fio(&c__1, (char *)&iyear, (ftnlen)sizeof(integer));
    do_fio(&c__1, (char *)&imon2, (ftnlen)sizeof(integer));
    do_fio(&c__1, (char *)&iday2, (ftnlen)sizeof(integer));
    e_wsfe();
    s_wsfe(&io___169);
    do_fio(&c__1, "recompute date from Julian day 1950         :", (ftnlen)45)
	    ;
    do_fio(&c__1, (char *)&iyea3, (ftnlen)sizeof(integer));
    do_fio(&c__1, (char *)&imon3, (ftnlen)sizeof(integer));
    do_fio(&c__1, (char *)&iday3, (ftnlen)sizeof(integer));
    e_wsfe();
    s_wsfe(&io___170);
    do_fio(&c__1, "recompute date from Julian day 2000         :", (ftnlen)45)
	    ;
    do_fio(&c__1, (char *)&iyea4, (ftnlen)sizeof(integer));
    do_fio(&c__1, (char *)&imon4, (ftnlen)sizeof(integer));
    do_fio(&c__1, (char *)&iday4, (ftnlen)sizeof(integer));
    e_wsfe();
    s_wsfe(&io___171);
    do_fio(&c__1, "recompute time from decimal hour:", (ftnlen)33);
    do_fio(&c__1, (char *)&ih2, (ftnlen)sizeof(integer));
    do_fio(&c__1, (char *)&im2, (ftnlen)sizeof(integer));
    do_fio(&c__1, (char *)&is2, (ftnlen)sizeof(integer));
    e_wsfe();

    s_wsfe(&io___172);
    e_wsfe();
    s_wsfe(&io___173);
    e_wsfe();
    s_wsfe(&io___174);
    do_fio(&c__1, "II) CHECK BASIC TIME PARAMETERS", (ftnlen)31);
    e_wsfe();
    s_wsfe(&io___175);
    do_fio(&c__1, "    ===========================", (ftnlen)31);
    e_wsfe();
    s_wsfe(&io___176);
    e_wsfe();
    s_wsfe(&io___177);
    e_wsfe();
    s_wsfe(&io___178);
    do_fio(&c__1, "SUN PARAMETERS IN GEI SYSTEM", (ftnlen)28);
    e_wsfe();
    s_wsfe(&io___179);
    e_wsfe();

    cp_time_param__(&iyear, &imonth, &iday, &ih, &im, &is);

    g_gei_geo_sun_dir__(&s1gei, &s2gei, &s3gei, &s1geo, &s2geo, &s3geo);
    g_gei_geo_dipole_dir__(&d1gei, &d2gei, &d3gei, &d1geo, &d2geo, &d3geo);
    g_gei_geo_ecliptic_dir__(&e1gei, &e2gei, &e3gei, &e1geo, &e2geo, &e3geo);
    g_gei_geo_sun_rot__(&r1gei, &r2gei, &r3gei, &r1geo, &r2geo, &r3geo);
    g_gei_sun_param__(&gmst, &slong, &srasn, &sdec, &obli);
    g_gsm_dipole_tilt_angle__(&diptan);

    t_car_to_sph__(&s1gei, &s2gei, &s3gei, &srgei, &stgei, &spgei);
    t_car_to_sph__(&d1gei, &d2gei, &d3gei, &drgei, &dtgei, &dpgei);
    t_car_to_sph__(&e1gei, &e2gei, &e3gei, &ergei, &etgei, &epgei);
    t_car_to_sph__(&r1gei, &r2gei, &r3gei, &rrgei, &rtgei, &rpgei);

    t_car_to_sph__(&s1geo, &s2geo, &s3geo, &srgeo, &stgeo, &spgeo);
    t_car_to_sph__(&d1geo, &d2geo, &d3geo, &drgeo, &dtgeo, &dpgeo);
    t_car_to_sph__(&e1geo, &e2geo, &e3geo, &ergeo, &etgeo, &epgeo);
    t_car_to_sph__(&r1geo, &r2geo, &r3geo, &rrgeo, &rtgeo, &rpgeo);


    s_wsfe(&io___234);
    do_fio(&c__1, "Greenwich Sideral Time (deg.):", (ftnlen)30);
    r__1 = gmst * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    e_wsfe();
    s_wsfe(&io___235);
    do_fio(&c__1, "ecliptic longitude     (deg.):", (ftnlen)30);
    r__1 = slong * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    e_wsfe();
    s_wsfe(&io___236);
    do_fio(&c__1, "right ascention        (deg.):", (ftnlen)30);
    r__1 = srasn * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    e_wsfe();
    s_wsfe(&io___237);
    do_fio(&c__1, "declination            (deg.):", (ftnlen)30);
    r__1 = sdec * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    e_wsfe();
    s_wsfe(&io___238);
    e_wsfe();
    s_wsfe(&io___239);
    do_fio(&c__1, "DIPOLE TILT ANGLE      (deg.):", (ftnlen)30);
    r__1 = diptan * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    e_wsfe();
    s_wsfe(&io___240);
    e_wsfe();
    s_wsfe(&io___241);
    do_fio(&c__1, "SUN DIRECTION", (ftnlen)13);
    e_wsfe();
    s_wsfe(&io___242);
    do_fio(&c__1, "                     in GEI:", (ftnlen)28);
    e_wsfe();
    s_wsfe(&io___243);
    do_fio(&c__1, (char *)&s1gei, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&srgei, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&s2gei, (ftnlen)sizeof(real));
    r__1 = stgei * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&s3gei, (ftnlen)sizeof(real));
    r__2 = spgei * d__;
    do_fio(&c__1, (char *)&r__2, (ftnlen)sizeof(real));
    e_wsfe();
    s_wsfe(&io___244);
    do_fio(&c__1, "                     in GEO:", (ftnlen)28);
    e_wsfe();
    s_wsfe(&io___245);
    do_fio(&c__1, (char *)&s1geo, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&srgeo, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&s2geo, (ftnlen)sizeof(real));
    r__1 = stgeo * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&s3geo, (ftnlen)sizeof(real));
    r__2 = spgeo * d__;
    do_fio(&c__1, (char *)&r__2, (ftnlen)sizeof(real));
    e_wsfe();
    s_wsfe(&io___246);
    e_wsfe();
    s_wsfe(&io___247);
    do_fio(&c__1, "DIPOLE DIRECTION", (ftnlen)16);
    e_wsfe();
    s_wsfe(&io___248);
    do_fio(&c__1, "                     in GEI:", (ftnlen)28);
    e_wsfe();
    s_wsfe(&io___249);
    do_fio(&c__1, (char *)&d1gei, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&drgei, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&d2gei, (ftnlen)sizeof(real));
    r__1 = dtgei * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&d3gei, (ftnlen)sizeof(real));
    r__2 = dpgei * d__;
    do_fio(&c__1, (char *)&r__2, (ftnlen)sizeof(real));
    e_wsfe();
    s_wsfe(&io___250);
    do_fio(&c__1, "                     in GEO:", (ftnlen)28);
    e_wsfe();
    s_wsfe(&io___251);
    do_fio(&c__1, (char *)&d1geo, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&drgeo, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&d2geo, (ftnlen)sizeof(real));
    r__1 = dtgeo * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&d3geo, (ftnlen)sizeof(real));
    r__2 = dpgeo * d__;
    do_fio(&c__1, (char *)&r__2, (ftnlen)sizeof(real));
    e_wsfe();
    s_wsfe(&io___252);
    e_wsfe();
    s_wsfe(&io___253);
    do_fio(&c__1, "ECLIPTIC DIRECTION", (ftnlen)18);
    e_wsfe();
    s_wsfe(&io___254);
    do_fio(&c__1, "                     in GEI:", (ftnlen)28);
    e_wsfe();
    s_wsfe(&io___255);
    do_fio(&c__1, (char *)&e1gei, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&ergei, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&e2gei, (ftnlen)sizeof(real));
    r__1 = etgei * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&e3gei, (ftnlen)sizeof(real));
    r__2 = epgei * d__;
    do_fio(&c__1, (char *)&r__2, (ftnlen)sizeof(real));
    e_wsfe();
    s_wsfe(&io___256);
    do_fio(&c__1, "                     in GEO:", (ftnlen)28);
    e_wsfe();
    s_wsfe(&io___257);
    do_fio(&c__1, (char *)&e1geo, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&ergeo, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&e2geo, (ftnlen)sizeof(real));
    r__1 = etgeo * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&e3geo, (ftnlen)sizeof(real));
    r__2 = epgeo * d__;
    do_fio(&c__1, (char *)&r__2, (ftnlen)sizeof(real));
    e_wsfe();
    s_wsfe(&io___258);
    e_wsfe();
    s_wsfe(&io___259);
    do_fio(&c__1, "SUN EQUATOR DIRECTION", (ftnlen)21);
    e_wsfe();
    s_wsfe(&io___260);
    do_fio(&c__1, "                     in GEI:", (ftnlen)28);
    e_wsfe();
    s_wsfe(&io___261);
    do_fio(&c__1, (char *)&r1gei, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&rrgei, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&r2gei, (ftnlen)sizeof(real));
    r__1 = rtgei * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&r3gei, (ftnlen)sizeof(real));
    r__2 = rpgei * d__;
    do_fio(&c__1, (char *)&r__2, (ftnlen)sizeof(real));
    e_wsfe();
    s_wsfe(&io___262);
    do_fio(&c__1, "                     in GEO:", (ftnlen)28);
    e_wsfe();
    s_wsfe(&io___263);
    do_fio(&c__1, (char *)&r1geo, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&rrgeo, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&r2geo, (ftnlen)sizeof(real));
    r__1 = rtgeo * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&r3geo, (ftnlen)sizeof(real));
    r__2 = rpgeo * d__;
    do_fio(&c__1, (char *)&r__2, (ftnlen)sizeof(real));
    e_wsfe();


    s_wsfe(&io___264);
    do_fio(&c__1, "\f", (ftnlen)1);
    e_wsfe();
    s_wsfe(&io___265);
    e_wsfe();
    s_wsfe(&io___266);
    do_fio(&c__1, "III) CHECK BASIC TRANSFORMATIONS", (ftnlen)32);
    e_wsfe();
    s_wsfe(&io___267);
    do_fio(&c__1, "     ===========================", (ftnlen)32);
    e_wsfe();
    s_wsfe(&io___268);
    e_wsfe();

    cp_angle_and_ratio__(&s1gei, &s2gei, &s3gei, &e1gei, &e2gei, &e3gei, &a, &
	    r__);
    cp_angle_and_ratio__(&s1geo, &s2geo, &s3geo, &e1geo, &e2geo, &e3geo, &a2, 
	    &r2);

    s_wsfe(&io___273);
    do_fio(&c__1, "Check (S,E) angle is equal to 90 deg. :", (ftnlen)39);
    r__1 = a * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    e_wsfe();
    if ((r__1 = a - a2, dabs(r__1)) > 1e-4f / d__) {
	s_wsfe(&io___274);
	do_fio(&c__1, " *** angle computed in GEI or GEO system", (ftnlen)40);
	e_wsfe();
	s_wsfe(&io___275);
	do_fio(&c__1, "     give not the same result:", (ftnlen)30);
	e_wsfe();
	s_wsfe(&io___276);
	do_fio(&c__1, "in GEI:", (ftnlen)7);
	r__1 = a * d__;
	do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
	e_wsfe();
	s_wsfe(&io___277);
	do_fio(&c__1, "in GEO:", (ftnlen)7);
	r__1 = a2 * d__;
	do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
	e_wsfe();
    }

    s_wsfe(&io___278);
    e_wsfe();
    s_wsfe(&io___279);
    e_wsfe();
    s_wsfe(&io___280);
    do_fio(&c__1, "1) input vector: Sun direction in GEI system", (ftnlen)44);
    e_wsfe();
    s_wsfe(&io___281);
    do_fio(&c__1, "   -----------------------------------------", (ftnlen)44);
    e_wsfe();
    s_wsfe(&io___282);
    e_wsfe();
    s_wsfe(&io___283);
    do_fio(&c__1, (char *)&s1gei, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&srgei, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&s2gei, (ftnlen)sizeof(real));
    r__1 = stgei * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&s3gei, (ftnlen)sizeof(real));
    r__2 = spgei * d__;
    do_fio(&c__1, (char *)&r__2, (ftnlen)sizeof(real));
    e_wsfe();
    s_wsfe(&io___284);
    e_wsfe();
    s_wsfe(&io___285);
    do_fio(&c__1, "                               in GEO system", (ftnlen)44);
    e_wsfe();
    s_wsfe(&io___286);
    e_wsfe();
    s_wsfe(&io___287);
    do_fio(&c__1, (char *)&s1geo, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&srgeo, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&s2geo, (ftnlen)sizeof(real));
    r__1 = stgeo * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&s3geo, (ftnlen)sizeof(real));
    r__2 = spgeo * d__;
    do_fio(&c__1, (char *)&r__2, (ftnlen)sizeof(real));
    e_wsfe();

    s_wsfe(&io___288);
    e_wsfe();
    s_wsfe(&io___289);
    e_wsfe();
    s_wsfe(&io___290);
    do_fio(&c__1, "   check Y=0 in SM  system with t_gei_to_sm:", (ftnlen)44);
    e_wsfe();
    t_gei_to_sm__(&s1gei, &s2gei, &s3gei, &v1, &v2, &v3);
    t_car_to_sph__(&v1, &v2, &v3, &vr, &vt, &vp);
    s_wsfe(&io___297);
    e_wsfe();
    s_wsfe(&io___298);
    do_fio(&c__1, (char *)&v1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&vr, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&v2, (ftnlen)sizeof(real));
    r__1 = vt * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&v3, (ftnlen)sizeof(real));
    r__2 = vp * d__;
    do_fio(&c__1, (char *)&r__2, (ftnlen)sizeof(real));
    e_wsfe();
    s_wsfe(&io___299);
    e_wsfe();
    s_wsfe(&io___300);
    do_fio(&c__1, "                           with t_geo_to_sm:", (ftnlen)44);
    e_wsfe();
    t_geo_to_sm__(&s1geo, &s2geo, &s3geo, &v1, &v2, &v3);
    t_car_to_sph__(&v1, &v2, &v3, &vr, &vt, &vp);
    s_wsfe(&io___301);
    e_wsfe();
    s_wsfe(&io___302);
    do_fio(&c__1, (char *)&v1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&vr, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&v2, (ftnlen)sizeof(real));
    r__1 = vt * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&v3, (ftnlen)sizeof(real));
    r__2 = vp * d__;
    do_fio(&c__1, (char *)&r__2, (ftnlen)sizeof(real));
    e_wsfe();

    s_wsfe(&io___303);
    e_wsfe();
    s_wsfe(&io___304);
    do_fio(&c__1, "   check Y=0 and Z=0 in GSM system with t_gei_to_gsm:", (
	    ftnlen)53);
    e_wsfe();
    t_gei_to_gsm__(&s1gei, &s2gei, &s3gei, &v1, &v2, &v3);
    t_car_to_sph__(&v1, &v2, &v3, &vr, &vt, &vp);
    s_wsfe(&io___305);
    e_wsfe();
    s_wsfe(&io___306);
    do_fio(&c__1, (char *)&v1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&vr, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&v2, (ftnlen)sizeof(real));
    r__1 = vt * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&v3, (ftnlen)sizeof(real));
    r__2 = vp * d__;
    do_fio(&c__1, (char *)&r__2, (ftnlen)sizeof(real));
    e_wsfe();
    s_wsfe(&io___307);
    e_wsfe();
    s_wsfe(&io___308);
    do_fio(&c__1, "                           with t_geo_to_gsm:", (ftnlen)45)
	    ;
    e_wsfe();
    t_geo_to_gsm__(&s1geo, &s2geo, &s3geo, &v1, &v2, &v3);
    t_car_to_sph__(&v1, &v2, &v3, &vr, &vt, &vp);
    s_wsfe(&io___309);
    e_wsfe();
    s_wsfe(&io___310);
    do_fio(&c__1, (char *)&v1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&vr, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&v2, (ftnlen)sizeof(real));
    r__1 = vt * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&v3, (ftnlen)sizeof(real));
    r__2 = vp * d__;
    do_fio(&c__1, (char *)&r__2, (ftnlen)sizeof(real));
    e_wsfe();

    s_wsfe(&io___311);
    e_wsfe();
    s_wsfe(&io___312);
    do_fio(&c__1, "   check Y=0 and Z=0 in GSE system with t_gei_to_gse:", (
	    ftnlen)53);
    e_wsfe();
    t_gei_to_gse__(&s1gei, &s2gei, &s3gei, &v1, &v2, &v3);
    t_car_to_sph__(&v1, &v2, &v3, &vr, &vt, &vp);
    s_wsfe(&io___313);
    e_wsfe();
    s_wsfe(&io___314);
    do_fio(&c__1, (char *)&v1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&vr, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&v2, (ftnlen)sizeof(real));
    r__1 = vt * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&v3, (ftnlen)sizeof(real));
    r__2 = vp * d__;
    do_fio(&c__1, (char *)&r__2, (ftnlen)sizeof(real));
    e_wsfe();
    s_wsfe(&io___315);
    e_wsfe();
    s_wsfe(&io___316);
    do_fio(&c__1, "                           with t_geo_to_gse:", (ftnlen)45)
	    ;
    e_wsfe();
    t_geo_to_gse__(&s1geo, &s2geo, &s3geo, &v1, &v2, &v3);
    t_car_to_sph__(&v1, &v2, &v3, &vr, &vt, &vp);
    s_wsfe(&io___317);
    e_wsfe();
    s_wsfe(&io___318);
    do_fio(&c__1, (char *)&v1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&vr, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&v2, (ftnlen)sizeof(real));
    r__1 = vt * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&v3, (ftnlen)sizeof(real));
    r__2 = vp * d__;
    do_fio(&c__1, (char *)&r__2, (ftnlen)sizeof(real));
    e_wsfe();

    s_wsfe(&io___319);
    e_wsfe();
    s_wsfe(&io___320);
    do_fio(&c__1, "   check Y=0 and Z=0 in GSQ system with t_gei_to_gseq:", (
	    ftnlen)54);
    e_wsfe();
    t_gei_to_gseq__(&s1gei, &s2gei, &s3gei, &v1, &v2, &v3);
    t_car_to_sph__(&v1, &v2, &v3, &vr, &vt, &vp);
    s_wsfe(&io___321);
    e_wsfe();
    s_wsfe(&io___322);
    do_fio(&c__1, (char *)&v1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&vr, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&v2, (ftnlen)sizeof(real));
    r__1 = vt * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&v3, (ftnlen)sizeof(real));
    r__2 = vp * d__;
    do_fio(&c__1, (char *)&r__2, (ftnlen)sizeof(real));
    e_wsfe();
    s_wsfe(&io___323);
    e_wsfe();
    s_wsfe(&io___324);
    do_fio(&c__1, "                           with t_geo_to_gseq:", (ftnlen)
	    46);
    e_wsfe();
    t_geo_to_gseq__(&s1geo, &s2geo, &s3geo, &v1, &v2, &v3);
    t_car_to_sph__(&v1, &v2, &v3, &vr, &vt, &vp);
    s_wsfe(&io___325);
    e_wsfe();
    s_wsfe(&io___326);
    do_fio(&c__1, (char *)&v1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&vr, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&v2, (ftnlen)sizeof(real));
    r__1 = vt * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&v3, (ftnlen)sizeof(real));
    r__2 = vp * d__;
    do_fio(&c__1, (char *)&r__2, (ftnlen)sizeof(real));
    e_wsfe();

    s_wsfe(&io___327);
    do_fio(&c__1, "\f", (ftnlen)1);
    e_wsfe();
    s_wsfe(&io___328);
    e_wsfe();
    s_wsfe(&io___329);
    do_fio(&c__1, "2) input vector: Dipole direction in GEI system", (ftnlen)
	    47);
    e_wsfe();
    s_wsfe(&io___330);
    do_fio(&c__1, "   --------------------------------------------", (ftnlen)
	    47);
    e_wsfe();
    s_wsfe(&io___331);
    e_wsfe();
    s_wsfe(&io___332);
    do_fio(&c__1, (char *)&d1gei, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&drgei, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&d2gei, (ftnlen)sizeof(real));
    r__1 = dtgei * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&d3gei, (ftnlen)sizeof(real));
    r__2 = dpgei * d__;
    do_fio(&c__1, (char *)&r__2, (ftnlen)sizeof(real));
    e_wsfe();
    s_wsfe(&io___333);
    e_wsfe();
    s_wsfe(&io___334);
    do_fio(&c__1, "                                  in GEO system", (ftnlen)
	    47);
    e_wsfe();
    s_wsfe(&io___335);
    e_wsfe();
    s_wsfe(&io___336);
    do_fio(&c__1, (char *)&d1geo, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&drgeo, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&d2geo, (ftnlen)sizeof(real));
    r__1 = dtgeo * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&d3geo, (ftnlen)sizeof(real));
    r__2 = dpgeo * d__;
    do_fio(&c__1, (char *)&r__2, (ftnlen)sizeof(real));
    e_wsfe();


    s_wsfe(&io___337);
    e_wsfe();
    s_wsfe(&io___338);
    e_wsfe();
    s_wsfe(&io___339);
    do_fio(&c__1, "   check X=0 and Y=0 in MAG system with t_gei_to_mag:", (
	    ftnlen)53);
    e_wsfe();
    t_gei_to_mag__(&d1gei, &d2gei, &d3gei, &v1, &v2, &v3);
    t_car_to_sph__(&v1, &v2, &v3, &vr, &vt, &vp);
    s_wsfe(&io___340);
    e_wsfe();
    s_wsfe(&io___341);
    do_fio(&c__1, (char *)&v1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&vr, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&v2, (ftnlen)sizeof(real));
    r__1 = vt * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&v3, (ftnlen)sizeof(real));
    r__2 = vp * d__;
    do_fio(&c__1, (char *)&r__2, (ftnlen)sizeof(real));
    e_wsfe();
    s_wsfe(&io___342);
    e_wsfe();
    s_wsfe(&io___343);
    do_fio(&c__1, "                           with t_geo_to_mag:", (ftnlen)45)
	    ;
    e_wsfe();
    t_geo_to_mag__(&d1geo, &d2geo, &d3geo, &v1, &v2, &v3);
    t_car_to_sph__(&v1, &v2, &v3, &vr, &vt, &vp);
    s_wsfe(&io___344);
    e_wsfe();
    s_wsfe(&io___345);
    do_fio(&c__1, (char *)&v1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&vr, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&v2, (ftnlen)sizeof(real));
    r__1 = vt * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&v3, (ftnlen)sizeof(real));
    r__2 = vp * d__;
    do_fio(&c__1, (char *)&r__2, (ftnlen)sizeof(real));
    e_wsfe();

    s_wsfe(&io___346);
    e_wsfe();
    s_wsfe(&io___347);
    do_fio(&c__1, "   check X=0 and Y=0 in SM  system with t_gei_to_sm:", (
	    ftnlen)52);
    e_wsfe();
    t_gei_to_sm__(&d1gei, &d2gei, &d3gei, &v1, &v2, &v3);
    t_car_to_sph__(&v1, &v2, &v3, &vr, &vt, &vp);
    s_wsfe(&io___348);
    e_wsfe();
    s_wsfe(&io___349);
    do_fio(&c__1, (char *)&v1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&vr, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&v2, (ftnlen)sizeof(real));
    r__1 = vt * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&v3, (ftnlen)sizeof(real));
    r__2 = vp * d__;
    do_fio(&c__1, (char *)&r__2, (ftnlen)sizeof(real));
    e_wsfe();
    s_wsfe(&io___350);
    e_wsfe();
    s_wsfe(&io___351);
    do_fio(&c__1, "                           with t_geo_to_sm:", (ftnlen)44);
    e_wsfe();
    t_geo_to_sm__(&d1geo, &d2geo, &d3geo, &v1, &v2, &v3);
    t_car_to_sph__(&v1, &v2, &v3, &vr, &vt, &vp);
    s_wsfe(&io___352);
    e_wsfe();
    s_wsfe(&io___353);
    do_fio(&c__1, (char *)&v1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&vr, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&v2, (ftnlen)sizeof(real));
    r__1 = vt * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&v3, (ftnlen)sizeof(real));
    r__2 = vp * d__;
    do_fio(&c__1, (char *)&r__2, (ftnlen)sizeof(real));
    e_wsfe();

    s_wsfe(&io___354);
    e_wsfe();
    s_wsfe(&io___355);
    do_fio(&c__1, "   check Y=0 in GSM system with t_gei_to_gsm:", (ftnlen)45)
	    ;
    e_wsfe();
    t_gei_to_gsm__(&d1gei, &d2gei, &d3gei, &v1, &v2, &v3);
    t_car_to_sph__(&v1, &v2, &v3, &vr, &vt, &vp);
    s_wsfe(&io___356);
    e_wsfe();
    s_wsfe(&io___357);
    do_fio(&c__1, (char *)&v1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&vr, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&v2, (ftnlen)sizeof(real));
    r__1 = vt * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&v3, (ftnlen)sizeof(real));
    r__2 = vp * d__;
    do_fio(&c__1, (char *)&r__2, (ftnlen)sizeof(real));
    e_wsfe();
    s_wsfe(&io___358);
    e_wsfe();
    s_wsfe(&io___359);
    do_fio(&c__1, "                           with t_geo_to_gsm:", (ftnlen)45)
	    ;
    e_wsfe();
    t_geo_to_gsm__(&d1geo, &d2geo, &d3geo, &v1, &v2, &v3);
    t_car_to_sph__(&v1, &v2, &v3, &vr, &vt, &vp);
    s_wsfe(&io___360);
    e_wsfe();
    s_wsfe(&io___361);
    do_fio(&c__1, (char *)&v1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&vr, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&v2, (ftnlen)sizeof(real));
    r__1 = vt * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&v3, (ftnlen)sizeof(real));
    r__2 = vp * d__;
    do_fio(&c__1, (char *)&r__2, (ftnlen)sizeof(real));
    e_wsfe();


    s_wsfe(&io___362);
    do_fio(&c__1, "\f", (ftnlen)1);
    e_wsfe();
    s_wsfe(&io___363);
    e_wsfe();
    s_wsfe(&io___364);
    do_fio(&c__1, "3) input vector: Ecliptic direction in GEI system", (
	    ftnlen)49);
    e_wsfe();
    s_wsfe(&io___365);
    do_fio(&c__1, "   ----------------------------------------------", (
	    ftnlen)49);
    e_wsfe();
    s_wsfe(&io___366);
    e_wsfe();
    s_wsfe(&io___367);
    do_fio(&c__1, (char *)&e1gei, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&ergei, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&e2gei, (ftnlen)sizeof(real));
    r__1 = etgei * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&e3gei, (ftnlen)sizeof(real));
    r__2 = epgei * d__;
    do_fio(&c__1, (char *)&r__2, (ftnlen)sizeof(real));
    e_wsfe();
    s_wsfe(&io___368);
    e_wsfe();
    s_wsfe(&io___369);
    do_fio(&c__1, "                                    in GEO system", (
	    ftnlen)49);
    e_wsfe();
    s_wsfe(&io___370);
    e_wsfe();
    s_wsfe(&io___371);
    do_fio(&c__1, (char *)&e1geo, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&ergeo, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&e2geo, (ftnlen)sizeof(real));
    r__1 = etgeo * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&e3geo, (ftnlen)sizeof(real));
    r__2 = epgeo * d__;
    do_fio(&c__1, (char *)&r__2, (ftnlen)sizeof(real));
    e_wsfe();


    s_wsfe(&io___372);
    e_wsfe();
    s_wsfe(&io___373);
    e_wsfe();
    s_wsfe(&io___374);
    do_fio(&c__1, "   check X=0 and Y=0 in GSE system with t_gei_to_gse:", (
	    ftnlen)53);
    e_wsfe();
    t_gei_to_gse__(&e1gei, &e2gei, &e3gei, &v1, &v2, &v3);
    t_car_to_sph__(&v1, &v2, &v3, &vr, &vt, &vp);
    s_wsfe(&io___375);
    e_wsfe();
    s_wsfe(&io___376);
    do_fio(&c__1, (char *)&v1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&vr, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&v2, (ftnlen)sizeof(real));
    r__1 = vt * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&v3, (ftnlen)sizeof(real));
    r__2 = vp * d__;
    do_fio(&c__1, (char *)&r__2, (ftnlen)sizeof(real));
    e_wsfe();
    s_wsfe(&io___377);
    e_wsfe();
    s_wsfe(&io___378);
    do_fio(&c__1, "                           with t_geo_to_gse:", (ftnlen)45)
	    ;
    e_wsfe();
    t_geo_to_gse__(&e1geo, &e2geo, &e3geo, &v1, &v2, &v3);
    t_car_to_sph__(&v1, &v2, &v3, &vr, &vt, &vp);
    s_wsfe(&io___379);
    e_wsfe();
    s_wsfe(&io___380);
    do_fio(&c__1, (char *)&v1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&vr, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&v2, (ftnlen)sizeof(real));
    r__1 = vt * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&v3, (ftnlen)sizeof(real));
    r__2 = vp * d__;
    do_fio(&c__1, (char *)&r__2, (ftnlen)sizeof(real));
    e_wsfe();

    s_wsfe(&io___381);
    e_wsfe();
    s_wsfe(&io___382);
    do_fio(&c__1, "   check X=0 in GSQ system with t_gei_to_gseq:", (ftnlen)
	    46);
    e_wsfe();
    t_gei_to_gseq__(&e1gei, &e2gei, &e3gei, &v1, &v2, &v3);
    t_car_to_sph__(&v1, &v2, &v3, &vr, &vt, &vp);
    s_wsfe(&io___383);
    e_wsfe();
    s_wsfe(&io___384);
    do_fio(&c__1, (char *)&v1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&vr, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&v2, (ftnlen)sizeof(real));
    r__1 = vt * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&v3, (ftnlen)sizeof(real));
    r__2 = vp * d__;
    do_fio(&c__1, (char *)&r__2, (ftnlen)sizeof(real));
    e_wsfe();
    s_wsfe(&io___385);
    e_wsfe();
    s_wsfe(&io___386);
    do_fio(&c__1, "                           with t_geo_to_gseq:", (ftnlen)
	    46);
    e_wsfe();
    t_geo_to_gseq__(&e1geo, &e2geo, &e3geo, &v1, &v2, &v3);
    t_car_to_sph__(&v1, &v2, &v3, &vr, &vt, &vp);
    s_wsfe(&io___387);
    e_wsfe();
    s_wsfe(&io___388);
    do_fio(&c__1, (char *)&v1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&vr, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&v2, (ftnlen)sizeof(real));
    r__1 = vt * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&v3, (ftnlen)sizeof(real));
    r__2 = vp * d__;
    do_fio(&c__1, (char *)&r__2, (ftnlen)sizeof(real));
    e_wsfe();


    s_wsfe(&io___389);
    e_wsfe();
    s_wsfe(&io___390);
    e_wsfe();
    s_wsfe(&io___391);
    e_wsfe();
    s_wsfe(&io___392);
    e_wsfe();
    s_wsfe(&io___393);
    e_wsfe();
    s_wsfe(&io___394);
    do_fio(&c__1, "4) input vector: Sun equator in GEI system", (ftnlen)42);
    e_wsfe();
    s_wsfe(&io___395);
    do_fio(&c__1, "   ---------------------------------------", (ftnlen)42);
    e_wsfe();
    s_wsfe(&io___396);
    e_wsfe();
    s_wsfe(&io___397);
    do_fio(&c__1, (char *)&r1gei, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&rrgei, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&r2gei, (ftnlen)sizeof(real));
    r__1 = rtgei * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&r3gei, (ftnlen)sizeof(real));
    r__2 = rpgei * d__;
    do_fio(&c__1, (char *)&r__2, (ftnlen)sizeof(real));
    e_wsfe();
    s_wsfe(&io___398);
    e_wsfe();
    s_wsfe(&io___399);
    do_fio(&c__1, "                             in GEO system", (ftnlen)42);
    e_wsfe();
    s_wsfe(&io___400);
    e_wsfe();
    s_wsfe(&io___401);
    do_fio(&c__1, (char *)&r1geo, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&rrgeo, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&r2geo, (ftnlen)sizeof(real));
    r__1 = rtgeo * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&r3geo, (ftnlen)sizeof(real));
    r__2 = rpgeo * d__;
    do_fio(&c__1, (char *)&r__2, (ftnlen)sizeof(real));
    e_wsfe();


    s_wsfe(&io___402);
    e_wsfe();
    s_wsfe(&io___403);
    e_wsfe();
    s_wsfe(&io___404);
    do_fio(&c__1, "   check and Y=0 in GSQ system with t_gei_to_gseq:", (
	    ftnlen)50);
    e_wsfe();
    t_gei_to_gseq__(&r1gei, &r2gei, &r3gei, &v1, &v2, &v3);
    t_car_to_sph__(&v1, &v2, &v3, &vr, &vt, &vp);
    s_wsfe(&io___405);
    e_wsfe();
    s_wsfe(&io___406);
    do_fio(&c__1, (char *)&v1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&vr, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&v2, (ftnlen)sizeof(real));
    r__1 = vt * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&v3, (ftnlen)sizeof(real));
    r__2 = vp * d__;
    do_fio(&c__1, (char *)&r__2, (ftnlen)sizeof(real));
    e_wsfe();
    s_wsfe(&io___407);
    e_wsfe();
    s_wsfe(&io___408);
    do_fio(&c__1, "                           with t_geo_to_gseq:", (ftnlen)
	    46);
    e_wsfe();
    t_geo_to_gseq__(&r1geo, &r2geo, &r3geo, &v1, &v2, &v3);
    t_car_to_sph__(&v1, &v2, &v3, &vr, &vt, &vp);
    s_wsfe(&io___409);
    e_wsfe();
    s_wsfe(&io___410);
    do_fio(&c__1, (char *)&v1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&vr, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&v2, (ftnlen)sizeof(real));
    r__1 = vt * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&v3, (ftnlen)sizeof(real));
    r__2 = vp * d__;
    do_fio(&c__1, (char *)&r__2, (ftnlen)sizeof(real));
    e_wsfe();


    s_wsfe(&io___411);
    do_fio(&c__1, "\f", (ftnlen)1);
    e_wsfe();
    s_wsfe(&io___412);
    e_wsfe();
    s_wsfe(&io___413);
    do_fio(&c__1, "5) input vector: North geographic in GEO system", (ftnlen)
	    47);
    e_wsfe();
    s_wsfe(&io___414);
    do_fio(&c__1, "   --------------------------------------------", (ftnlen)
	    47);
    e_wsfe();
    t_car_to_sph__(&c_b595, &c_b595, &c_b597, &vr, &vt, &vp);
    s_wsfe(&io___415);
    e_wsfe();
    s_wsfe(&io___416);
    do_fio(&c__1, (char *)&c_b595, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&vr, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&c_b595, (ftnlen)sizeof(real));
    r__1 = vt * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&c_b597, (ftnlen)sizeof(real));
    r__2 = vp * d__;
    do_fio(&c__1, (char *)&r__2, (ftnlen)sizeof(real));
    e_wsfe();
    s_wsfe(&io___417);
    e_wsfe();
    s_wsfe(&io___418);
    e_wsfe();
    s_wsfe(&io___419);
    do_fio(&c__1, "   check X=0 and Y=0 in GEI system with t_geo_to_gei:", (
	    ftnlen)53);
    e_wsfe();
    t_geo_to_gei__(&c_b595, &c_b595, &c_b597, &w1, &w2, &w3);
    t_car_to_sph__(&w1, &w2, &w3, &wr, &wt, &wp);
    s_wsfe(&io___426);
    e_wsfe();
    s_wsfe(&io___427);
    do_fio(&c__1, (char *)&w1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&wr, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&w2, (ftnlen)sizeof(real));
    r__1 = wt * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&w3, (ftnlen)sizeof(real));
    r__2 = wp * d__;
    do_fio(&c__1, (char *)&r__2, (ftnlen)sizeof(real));
    e_wsfe();

    s_wsfe(&io___428);
    e_wsfe();
    s_wsfe(&io___429);
    do_fio(&c__1, "   check Y=0 in MAG system with t_gei_to_mag:", (ftnlen)45)
	    ;
    e_wsfe();
    t_gei_to_mag__(&w1, &w2, &w3, &v1, &v2, &v3);
    t_car_to_sph__(&v1, &v2, &v3, &vr, &vt, &vp);
    s_wsfe(&io___430);
    e_wsfe();
    s_wsfe(&io___431);
    do_fio(&c__1, (char *)&v1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&vr, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&v2, (ftnlen)sizeof(real));
    r__1 = vt * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&v3, (ftnlen)sizeof(real));
    r__2 = vp * d__;
    do_fio(&c__1, (char *)&r__2, (ftnlen)sizeof(real));
    e_wsfe();
    s_wsfe(&io___432);
    e_wsfe();
    s_wsfe(&io___433);
    do_fio(&c__1, "                           with t_geo_to_mag:", (ftnlen)45)
	    ;
    e_wsfe();
    t_geo_to_mag__(&w1, &w2, &w3, &v1, &v2, &v3);
    t_car_to_sph__(&v1, &v2, &v3, &vr, &vt, &vp);
    s_wsfe(&io___434);
    e_wsfe();
    s_wsfe(&io___435);
    do_fio(&c__1, (char *)&v1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&vr, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&v2, (ftnlen)sizeof(real));
    r__1 = vt * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&v3, (ftnlen)sizeof(real));
    r__2 = vp * d__;
    do_fio(&c__1, (char *)&r__2, (ftnlen)sizeof(real));
    e_wsfe();


    s_wsfe(&io___436);
    e_wsfe();
    s_wsfe(&io___437);
    e_wsfe();
    s_wsfe(&io___438);
    do_fio(&c__1, "IV) TEST OF ROCOTLIB TRANSFORMATION SUBROUTINES", (ftnlen)
	    47);
    e_wsfe();
    s_wsfe(&io___439);
    do_fio(&c__1, "    ===========================================", (ftnlen)
	    47);
    e_wsfe();
    s_wsfe(&io___440);
    e_wsfe();
    s_wsfe(&io___441);
    do_fio(&c__1, "1) input choosed test vector in GEO system :", (ftnlen)44);
    e_wsfe();
    s_wsfe(&io___442);
    e_wsfe();
    s_wsfe(&io___443);
    do_fio(&c__1, (char *)&g1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&gr, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&g2, (ftnlen)sizeof(real));
    r__1 = gt * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&g3, (ftnlen)sizeof(real));
    r__2 = gp * d__;
    do_fio(&c__1, (char *)&r__2, (ftnlen)sizeof(real));
    e_wsfe();

    s_wsfe(&io___444);
    e_wsfe();
    s_wsfe(&io___445);
    do_fio(&c__1, "2) converts GEO system to GEI system by t_geo_to_gei", (
	    ftnlen)52);
    e_wsfe();
    t_geo_to_gei__(&g1, &g2, &g3, &u1, &u2, &u3);
    t_car_to_sph__(&u1, &u2, &u3, &ur, &ut, &up);
    s_wsfe(&io___452);
    e_wsfe();
    s_wsfe(&io___453);
    do_fio(&c__1, (char *)&u1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&ur, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&u2, (ftnlen)sizeof(real));
    r__1 = ut * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&u3, (ftnlen)sizeof(real));
    r__2 = up * d__;
    do_fio(&c__1, (char *)&r__2, (ftnlen)sizeof(real));
    e_wsfe();


    s_wsfe(&io___454);
    do_fio(&c__1, "\f", (ftnlen)1);
    e_wsfe();
    s_wsfe(&io___455);
    e_wsfe();
    s_wsfe(&io___456);
    do_fio(&c__1, "A) STAR TRANSFORMATIONS AROUND  GEI  SYSTEM", (ftnlen)43);
    e_wsfe();
    s_wsfe(&io___457);
    do_fio(&c__1, "   ========================================", (ftnlen)43);
    e_wsfe();
    s_wsfe(&io___458);
    e_wsfe();
    s_wsfe(&io___459);
    do_fio(&c__1, "            # 3", (ftnlen)15);
    e_wsfe();
    s_wsfe(&io___460);
    do_fio(&c__1, "       GEI --------> MAG", (ftnlen)24);
    e_wsfe();
    s_wsfe(&io___461);
    do_fio(&c__1, "            # 4", (ftnlen)15);
    e_wsfe();
    s_wsfe(&io___462);
    do_fio(&c__1, "       GEI --------> SM", (ftnlen)23);
    e_wsfe();
    s_wsfe(&io___463);
    do_fio(&c__1, "            # 5", (ftnlen)15);
    e_wsfe();
    s_wsfe(&io___464);
    do_fio(&c__1, "       GEI --------> GSM", (ftnlen)24);
    e_wsfe();
    s_wsfe(&io___465);
    do_fio(&c__1, "            # 6", (ftnlen)15);
    e_wsfe();
    s_wsfe(&io___466);
    do_fio(&c__1, "       GEI --------> GSE", (ftnlen)24);
    e_wsfe();
    s_wsfe(&io___467);
    e_wsfe();

    s_wsfe(&io___468);
    e_wsfe();
    s_wsfe(&io___469);
    do_fio(&c__1, "1) t# 3 converts GEI system to MAG system byt_gei_to_mag", 
	    (ftnlen)56);
    e_wsfe();
    t_gei_to_mag__(&u1, &u2, &u3, &v1, &v2, &v3);
    t_car_to_sph__(&v1, &v2, &v3, &vr, &vt, &vp);
    s_wsfe(&io___470);
    e_wsfe();
    s_wsfe(&io___471);
    do_fio(&c__1, (char *)&v1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&vr, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&v2, (ftnlen)sizeof(real));
    r__1 = vt * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&v3, (ftnlen)sizeof(real));
    r__2 = vp * d__;
    do_fio(&c__1, (char *)&r__2, (ftnlen)sizeof(real));
    e_wsfe();

    s_wsfe(&io___472);
    e_wsfe();
    s_wsfe(&io___473);
    do_fio(&c__1, "              come back to GEI system by t_mag_to_gei", (
	    ftnlen)53);
    e_wsfe();
    t_mag_to_gei__(&v1, &v2, &v3, &w1, &w2, &w3);
    t_car_to_sph__(&w1, &w2, &w3, &wr, &wt, &wp);
    s_wsfe(&io___474);
    e_wsfe();
    s_wsfe(&io___475);
    do_fio(&c__1, (char *)&w1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&wr, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&w2, (ftnlen)sizeof(real));
    r__1 = wt * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&w3, (ftnlen)sizeof(real));
    r__2 = wp * d__;
    do_fio(&c__1, (char *)&r__2, (ftnlen)sizeof(real));
    e_wsfe();

    cp_angle_and_ratio__(&u1, &u2, &u3, &w1, &w2, &w3, &a, &ra);
    s_wsfe(&io___477);
    do_fio(&c__1, "difference beetween first GEI vector:", (ftnlen)37);
    r__1 = a * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&ra, (ftnlen)sizeof(real));
    e_wsfe();


    s_wsfe(&io___478);
    e_wsfe();
    s_wsfe(&io___479);
    e_wsfe();
    s_wsfe(&io___480);
    do_fio(&c__1, "2) t# 4 converts GEI system to SM  system byt_gei_to_sm", (
	    ftnlen)55);
    e_wsfe();
    t_gei_to_sm__(&u1, &u2, &u3, &v1, &v2, &v3);
    t_car_to_sph__(&v1, &v2, &v3, &vr, &vt, &vp);
    s_wsfe(&io___481);
    e_wsfe();
    s_wsfe(&io___482);
    do_fio(&c__1, (char *)&v1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&vr, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&v2, (ftnlen)sizeof(real));
    r__1 = vt * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&v3, (ftnlen)sizeof(real));
    r__2 = vp * d__;
    do_fio(&c__1, (char *)&r__2, (ftnlen)sizeof(real));
    e_wsfe();

    s_wsfe(&io___483);
    e_wsfe();
    s_wsfe(&io___484);
    do_fio(&c__1, "              come back to GEI system by t_sm_to_gei", (
	    ftnlen)52);
    e_wsfe();
    t_sm_to_gei__(&v1, &v2, &v3, &w1, &w2, &w3);
    t_car_to_sph__(&w1, &w2, &w3, &wr, &wt, &wp);
    s_wsfe(&io___485);
    e_wsfe();
    s_wsfe(&io___486);
    do_fio(&c__1, (char *)&w1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&wr, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&w2, (ftnlen)sizeof(real));
    r__1 = wt * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&w3, (ftnlen)sizeof(real));
    r__2 = wp * d__;
    do_fio(&c__1, (char *)&r__2, (ftnlen)sizeof(real));
    e_wsfe();

    cp_angle_and_ratio__(&u1, &u2, &u3, &w1, &w2, &w3, &a, &ra);
    s_wsfe(&io___487);
    do_fio(&c__1, "difference beetween first GEI vector:", (ftnlen)37);
    r__1 = a * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&ra, (ftnlen)sizeof(real));
    e_wsfe();


    s_wsfe(&io___488);
    e_wsfe();
    s_wsfe(&io___489);
    e_wsfe();
    s_wsfe(&io___490);
    do_fio(&c__1, "3) t# 5 converts GEI system to GSM system byt_gei_to_gsm", 
	    (ftnlen)56);
    e_wsfe();
    t_gei_to_gsm__(&u1, &u2, &u3, &v1, &v2, &v3);
    t_car_to_sph__(&v1, &v2, &v3, &vr, &vt, &vp);
    s_wsfe(&io___491);
    e_wsfe();
    s_wsfe(&io___492);
    do_fio(&c__1, (char *)&v1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&vr, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&v2, (ftnlen)sizeof(real));
    r__1 = vt * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&v3, (ftnlen)sizeof(real));
    r__2 = vp * d__;
    do_fio(&c__1, (char *)&r__2, (ftnlen)sizeof(real));
    e_wsfe();

    s_wsfe(&io___493);
    e_wsfe();
    s_wsfe(&io___494);
    do_fio(&c__1, "              come back to GEI system by t_gsm_to_gei", (
	    ftnlen)53);
    e_wsfe();
    t_gsm_to_gei__(&v1, &v2, &v3, &w1, &w2, &w3);
    t_car_to_sph__(&w1, &w2, &w3, &wr, &wt, &wp);
    s_wsfe(&io___495);
    e_wsfe();
    s_wsfe(&io___496);
    do_fio(&c__1, (char *)&w1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&wr, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&w2, (ftnlen)sizeof(real));
    r__1 = wt * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&w3, (ftnlen)sizeof(real));
    r__2 = wp * d__;
    do_fio(&c__1, (char *)&r__2, (ftnlen)sizeof(real));
    e_wsfe();

    cp_angle_and_ratio__(&u1, &u2, &u3, &w1, &w2, &w3, &a, &ra);
    s_wsfe(&io___497);
    do_fio(&c__1, "difference beetween first GEI vector:", (ftnlen)37);
    r__1 = a * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&ra, (ftnlen)sizeof(real));
    e_wsfe();


    s_wsfe(&io___498);
    do_fio(&c__1, "\f", (ftnlen)1);
    e_wsfe();
    s_wsfe(&io___499);
    e_wsfe();
    s_wsfe(&io___500);
    do_fio(&c__1, "4) t# 6 converts GEI system to GSE system by t_gei_to_gse",
	     (ftnlen)57);
    e_wsfe();
    t_gei_to_gse__(&u1, &u2, &u3, &v1, &v2, &v3);
    t_car_to_sph__(&v1, &v2, &v3, &vr, &vt, &vp);
    s_wsfe(&io___501);
    e_wsfe();
    s_wsfe(&io___502);
    do_fio(&c__1, (char *)&v1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&vr, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&v2, (ftnlen)sizeof(real));
    r__1 = vt * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&v3, (ftnlen)sizeof(real));
    r__2 = vp * d__;
    do_fio(&c__1, (char *)&r__2, (ftnlen)sizeof(real));
    e_wsfe();

    s_wsfe(&io___503);
    e_wsfe();
    s_wsfe(&io___504);
    do_fio(&c__1, "              come back to GEI system by t_gse_to_gei", (
	    ftnlen)53);
    e_wsfe();
    t_gse_to_gei__(&v1, &v2, &v3, &w1, &w2, &w3);
    t_car_to_sph__(&w1, &w2, &w3, &wr, &wt, &wp);
    s_wsfe(&io___505);
    e_wsfe();
    s_wsfe(&io___506);
    do_fio(&c__1, (char *)&w1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&wr, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&w2, (ftnlen)sizeof(real));
    r__1 = wt * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&w3, (ftnlen)sizeof(real));
    r__2 = wp * d__;
    do_fio(&c__1, (char *)&r__2, (ftnlen)sizeof(real));
    e_wsfe();

    cp_angle_and_ratio__(&u1, &u2, &u3, &w1, &w2, &w3, &a, &ra);
    s_wsfe(&io___507);
    do_fio(&c__1, "difference beetween first GEI vector:", (ftnlen)37);
    r__1 = a * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&ra, (ftnlen)sizeof(real));
    e_wsfe();

    s_wsfe(&io___508);
    do_fio(&c__1, "\f", (ftnlen)1);
    e_wsfe();
    s_wsfe(&io___509);
    e_wsfe();
    s_wsfe(&io___510);
    do_fio(&c__1, "B) COUNTERCLOCKWISE RING CUMULATIVE TRANSFORMATIONS", (
	    ftnlen)51);
    e_wsfe();
    s_wsfe(&io___511);
    do_fio(&c__1, "   ================================================", (
	    ftnlen)51);
    e_wsfe();
    s_wsfe(&io___512);
    e_wsfe();
    s_wsfe(&io___513);
    do_fio(&c__1, "    # 2      7      8      9     10    11     12", (ftnlen)
	    48);
    e_wsfe();
    s_wsfe(&io___514);
    do_fio(&c__1, "  GEO -> GEI -> GSQ -> GSE -> GSM -> SM -> MAG -> GEO", (
	    ftnlen)53);
    e_wsfe();
    s_wsfe(&io___515);
    e_wsfe();
    s_wsfe(&io___516);
    e_wsfe();
    s_wsfe(&io___517);
    do_fio(&c__1, "   input choosed test vector in GEO system :", (ftnlen)44);
    e_wsfe();
    s_wsfe(&io___518);
    e_wsfe();
    s_wsfe(&io___519);
    do_fio(&c__1, (char *)&g1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&gr, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&g2, (ftnlen)sizeof(real));
    r__1 = gt * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&g3, (ftnlen)sizeof(real));
    r__2 = gp * d__;
    do_fio(&c__1, (char *)&r__2, (ftnlen)sizeof(real));
    e_wsfe();

    s_wsfe(&io___520);
    e_wsfe();
    s_wsfe(&io___521);
    do_fio(&c__1, "1) t# 2 converts GEO system to GEI system by t_geo_to_gei",
	     (ftnlen)57);
    e_wsfe();
    t_geo_to_gei__(&g1, &g2, &g3, &u1, &u2, &u3);
    t_car_to_sph__(&u1, &u2, &u3, &ur, &ut, &up);
    s_wsfe(&io___522);
    e_wsfe();
    s_wsfe(&io___523);
    do_fio(&c__1, (char *)&u1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&ur, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&u2, (ftnlen)sizeof(real));
    r__1 = ut * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&u3, (ftnlen)sizeof(real));
    r__2 = up * d__;
    do_fio(&c__1, (char *)&r__2, (ftnlen)sizeof(real));
    e_wsfe();

    s_wsfe(&io___524);
    e_wsfe();
    s_wsfe(&io___525);
    do_fio(&c__1, "2) t# 7 converts GEI system to GSQ system by t_gei_to_gseq"
	    , (ftnlen)58);
    e_wsfe();
    t_gei_to_gseq__(&u1, &u2, &u3, &v1, &v2, &v3);
    t_car_to_sph__(&v1, &v2, &v3, &vr, &vt, &vp);
    s_wsfe(&io___526);
    e_wsfe();
    s_wsfe(&io___527);
    do_fio(&c__1, (char *)&v1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&vr, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&v2, (ftnlen)sizeof(real));
    r__1 = vt * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&v3, (ftnlen)sizeof(real));
    r__2 = vp * d__;
    do_fio(&c__1, (char *)&r__2, (ftnlen)sizeof(real));
    e_wsfe();


    s_wsfe(&io___528);
    e_wsfe();
    s_wsfe(&io___529);
    do_fio(&c__1, "3) t# 8 converts GSQ system to GSE system by t_gseq_to_gse"
	    , (ftnlen)58);
    e_wsfe();
    t_gseq_to_gse__(&v1, &v2, &v3, &w1, &w2, &w3);
    t_car_to_sph__(&w1, &w2, &w3, &wr, &wt, &wp);
    s_wsfe(&io___530);
    e_wsfe();
    s_wsfe(&io___531);
    do_fio(&c__1, (char *)&w1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&wr, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&w2, (ftnlen)sizeof(real));
    r__1 = wt * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&w3, (ftnlen)sizeof(real));
    r__2 = wp * d__;
    do_fio(&c__1, (char *)&r__2, (ftnlen)sizeof(real));
    e_wsfe();


    s_wsfe(&io___532);
    e_wsfe();
    s_wsfe(&io___533);
    do_fio(&c__1, "4) \303\276# 9 converts GSE system to GSM system by t_gse"
	    "_to_gsm", (ftnlen)58);
    e_wsfe();
    t_gse_to_gsm__(&w1, &w2, &w3, &v1, &v2, &v3);
    t_car_to_sph__(&v1, &v2, &v3, &vr, &vt, &vp);
    s_wsfe(&io___534);
    e_wsfe();
    s_wsfe(&io___535);
    do_fio(&c__1, (char *)&v1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&vr, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&v2, (ftnlen)sizeof(real));
    r__1 = vt * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&v3, (ftnlen)sizeof(real));
    r__2 = vp * d__;
    do_fio(&c__1, (char *)&r__2, (ftnlen)sizeof(real));
    e_wsfe();


    s_wsfe(&io___536);
    e_wsfe();
    s_wsfe(&io___537);
    do_fio(&c__1, "5) t# 10 converts GSM system to SM  system by t_gsm_to_sm",
	     (ftnlen)57);
    e_wsfe();
    t_gsm_to_sm__(&v1, &v2, &v3, &w1, &w2, &w3);
    t_car_to_sph__(&w1, &w2, &w3, &wr, &wt, &wp);
    s_wsfe(&io___538);
    e_wsfe();
    s_wsfe(&io___539);
    do_fio(&c__1, (char *)&w1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&wr, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&w2, (ftnlen)sizeof(real));
    r__1 = wt * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&w3, (ftnlen)sizeof(real));
    r__2 = wp * d__;
    do_fio(&c__1, (char *)&r__2, (ftnlen)sizeof(real));
    e_wsfe();


    s_wsfe(&io___540);
    e_wsfe();
    s_wsfe(&io___541);
    do_fio(&c__1, "6) t# 11 converts SM  system to MAG system by t_sm_to_mag",
	     (ftnlen)57);
    e_wsfe();
    t_sm_to_mag__(&w1, &w2, &w3, &v1, &v2, &v3);
    t_car_to_sph__(&v1, &v2, &v3, &vr, &vt, &vp);
    s_wsfe(&io___542);
    e_wsfe();
    s_wsfe(&io___543);
    do_fio(&c__1, (char *)&v1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&vr, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&v2, (ftnlen)sizeof(real));
    r__1 = vt * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&v3, (ftnlen)sizeof(real));
    r__2 = vp * d__;
    do_fio(&c__1, (char *)&r__2, (ftnlen)sizeof(real));
    e_wsfe();


    s_wsfe(&io___544);
    e_wsfe();
    s_wsfe(&io___545);
    do_fio(&c__1, "7) t# 12 converts MAG system to GEO system by t_mag_to_geo"
	    , (ftnlen)58);
    e_wsfe();
    t_mag_to_geo__(&v1, &v2, &v3, &w1, &w2, &w3);
    t_car_to_sph__(&w1, &w2, &w3, &wr, &wt, &wp);
    s_wsfe(&io___546);
    e_wsfe();
    s_wsfe(&io___547);
    do_fio(&c__1, (char *)&w1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&wr, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&w2, (ftnlen)sizeof(real));
    r__1 = wt * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&w3, (ftnlen)sizeof(real));
    r__2 = wp * d__;
    do_fio(&c__1, (char *)&r__2, (ftnlen)sizeof(real));
    e_wsfe();

    cp_angle_and_ratio__(&g1, &g2, &g3, &w1, &w2, &w3, &a, &ra);
    s_wsfe(&io___548);
    do_fio(&c__1, "difference beetween first GEO vector:", (ftnlen)37);
    r__1 = a * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&ra, (ftnlen)sizeof(real));
    e_wsfe();


    s_wsfe(&io___549);
    do_fio(&c__1, "\f", (ftnlen)1);
    e_wsfe();
    s_wsfe(&io___550);
    e_wsfe();
    s_wsfe(&io___551);
    do_fio(&c__1, "C) CLOCKWISE RING CUMULATIVE TRANSFORMATIONS", (ftnlen)44);
    e_wsfe();
    s_wsfe(&io___552);
    do_fio(&c__1, "   =========================================", (ftnlen)44);
    e_wsfe();
    s_wsfe(&io___553);
    e_wsfe();
    s_wsfe(&io___554);
    do_fio(&c__1, "    # 12     11    10     9      8      7      2", (ftnlen)
	    48);
    e_wsfe();
    s_wsfe(&io___555);
    do_fio(&c__1, "  GEO -> MAG -> SM -> GSM -> GSE -> GSQ -> GEI -> GEO", (
	    ftnlen)53);
    e_wsfe();
    s_wsfe(&io___556);
    e_wsfe();
    s_wsfe(&io___557);
    e_wsfe();
    s_wsfe(&io___558);
    do_fio(&c__1, "   input choosed test vector in GEO system :", (ftnlen)44);
    e_wsfe();
    s_wsfe(&io___559);
    e_wsfe();
    s_wsfe(&io___560);
    do_fio(&c__1, (char *)&g1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&gr, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&g2, (ftnlen)sizeof(real));
    r__1 = gt * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&g3, (ftnlen)sizeof(real));
    r__2 = gp * d__;
    do_fio(&c__1, (char *)&r__2, (ftnlen)sizeof(real));
    e_wsfe();

    s_wsfe(&io___561);
    e_wsfe();
    s_wsfe(&io___562);
    do_fio(&c__1, "1) t# 12 converts GEO system to MAG system by t_geo_to_mag"
	    , (ftnlen)58);
    e_wsfe();
    t_geo_to_mag__(&g1, &g2, &g3, &w1, &w2, &w3);
    t_car_to_sph__(&w1, &w2, &w3, &wr, &wt, &wp);
    s_wsfe(&io___563);
    e_wsfe();
    s_wsfe(&io___564);
    do_fio(&c__1, (char *)&w1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&wr, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&w2, (ftnlen)sizeof(real));
    r__1 = wt * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&w3, (ftnlen)sizeof(real));
    r__2 = wp * d__;
    do_fio(&c__1, (char *)&r__2, (ftnlen)sizeof(real));
    e_wsfe();


    s_wsfe(&io___565);
    e_wsfe();
    s_wsfe(&io___566);
    do_fio(&c__1, "2) t# 11 converts MAG system to SM  system by t_mag_to_sm",
	     (ftnlen)57);
    e_wsfe();
    t_mag_to_sm__(&w1, &w2, &w3, &v1, &v2, &v3);
    t_car_to_sph__(&v1, &v2, &v3, &vr, &vt, &vp);
    s_wsfe(&io___567);
    e_wsfe();
    s_wsfe(&io___568);
    do_fio(&c__1, (char *)&v1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&vr, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&v2, (ftnlen)sizeof(real));
    r__1 = vt * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&v3, (ftnlen)sizeof(real));
    r__2 = vp * d__;
    do_fio(&c__1, (char *)&r__2, (ftnlen)sizeof(real));
    e_wsfe();


    s_wsfe(&io___569);
    e_wsfe();
    s_wsfe(&io___570);
    do_fio(&c__1, "3) t# 10 converts SM  system to GSM system by t_sm_to_gsm",
	     (ftnlen)57);
    e_wsfe();
    t_sm_to_gsm__(&v1, &v2, &v3, &w1, &w2, &w3);
    t_car_to_sph__(&w1, &w2, &w3, &wr, &wt, &wp);
    s_wsfe(&io___571);
    e_wsfe();
    s_wsfe(&io___572);
    do_fio(&c__1, (char *)&w1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&wr, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&w2, (ftnlen)sizeof(real));
    r__1 = wt * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&w3, (ftnlen)sizeof(real));
    r__2 = wp * d__;
    do_fio(&c__1, (char *)&r__2, (ftnlen)sizeof(real));
    e_wsfe();


    s_wsfe(&io___573);
    e_wsfe();
    s_wsfe(&io___574);
    do_fio(&c__1, "4) t# 9 converts GSM system to GSE system by t_gsm_to_gse",
	     (ftnlen)57);
    e_wsfe();
    t_gsm_to_gse__(&w1, &w2, &w3, &v1, &v2, &v3);
    t_car_to_sph__(&v1, &v2, &v3, &vr, &vt, &vp);
    s_wsfe(&io___575);
    e_wsfe();
    s_wsfe(&io___576);
    do_fio(&c__1, (char *)&v1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&vr, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&v2, (ftnlen)sizeof(real));
    r__1 = vt * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&v3, (ftnlen)sizeof(real));
    r__2 = vp * d__;
    do_fio(&c__1, (char *)&r__2, (ftnlen)sizeof(real));
    e_wsfe();


    s_wsfe(&io___577);
    e_wsfe();
    s_wsfe(&io___578);
    do_fio(&c__1, "5) t# 8 converts GSE system to GSQ system by t_gse_to_gseq"
	    , (ftnlen)58);
    e_wsfe();
    t_gse_to_gseq__(&v1, &v2, &v3, &w1, &w2, &w3);
    t_car_to_sph__(&w1, &w2, &w3, &wr, &wt, &wp);
    s_wsfe(&io___579);
    e_wsfe();
    s_wsfe(&io___580);
    do_fio(&c__1, (char *)&w1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&wr, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&w2, (ftnlen)sizeof(real));
    r__1 = wt * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&w3, (ftnlen)sizeof(real));
    r__2 = wp * d__;
    do_fio(&c__1, (char *)&r__2, (ftnlen)sizeof(real));
    e_wsfe();


    s_wsfe(&io___581);
    e_wsfe();
    s_wsfe(&io___582);
    do_fio(&c__1, "6) t# 7 converts GSQ system to GEI system by t_gseq_to_gei"
	    , (ftnlen)58);
    e_wsfe();
    t_gseq_to_gei__(&w1, &w2, &w3, &v1, &v2, &v3);
    t_car_to_sph__(&v1, &v2, &v3, &vr, &vt, &vp);
    s_wsfe(&io___583);
    e_wsfe();
    s_wsfe(&io___584);
    do_fio(&c__1, (char *)&v1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&vr, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&v2, (ftnlen)sizeof(real));
    r__1 = vt * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&v3, (ftnlen)sizeof(real));
    r__2 = vp * d__;
    do_fio(&c__1, (char *)&r__2, (ftnlen)sizeof(real));
    e_wsfe();

    cp_angle_and_ratio__(&u1, &u2, &u3, &v1, &v2, &v3, &a, &ra);
    s_wsfe(&io___585);
    do_fio(&c__1, "difference beetween first GEI vector:", (ftnlen)37);
    r__1 = a * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&ra, (ftnlen)sizeof(real));
    e_wsfe();


    s_wsfe(&io___586);
    e_wsfe();
    s_wsfe(&io___587);
    e_wsfe();
    s_wsfe(&io___588);
    do_fio(&c__1, "7) t# 2 converts GEI system to GEO system by t_gei_to_geo",
	     (ftnlen)57);
    e_wsfe();
    t_gei_to_geo__(&v1, &v2, &v3, &w1, &w2, &w3);
    t_car_to_sph__(&w1, &w2, &w3, &wr, &wt, &wp);
    s_wsfe(&io___589);
    e_wsfe();
    s_wsfe(&io___590);
    do_fio(&c__1, (char *)&w1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&wr, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&w2, (ftnlen)sizeof(real));
    r__1 = wt * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&w3, (ftnlen)sizeof(real));
    r__2 = wp * d__;
    do_fio(&c__1, (char *)&r__2, (ftnlen)sizeof(real));
    e_wsfe();

    cp_angle_and_ratio__(&g1, &g2, &g3, &w1, &w2, &w3, &a, &ra);
    s_wsfe(&io___591);
    do_fio(&c__1, "difference beetween first GEO vector:", (ftnlen)37);
    r__1 = a * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&ra, (ftnlen)sizeof(real));
    e_wsfe();

    s_wsfe(&io___592);
    do_fio(&c__1, "\f", (ftnlen)1);
    e_wsfe();
    s_wsfe(&io___593);
    do_fio(&c__1, "D) STAR TRANSFORMATIONS AROUND  GEO  SYSTEM", (ftnlen)43);
    e_wsfe();
    s_wsfe(&io___594);
    do_fio(&c__1, "   ========================================", (ftnlen)43);
    e_wsfe();
    s_wsfe(&io___595);
    e_wsfe();
    s_wsfe(&io___596);
    do_fio(&c__1, "          # 13", (ftnlen)14);
    e_wsfe();
    s_wsfe(&io___597);
    do_fio(&c__1, "     GEO --------> SM", (ftnlen)21);
    e_wsfe();
    s_wsfe(&io___598);
    do_fio(&c__1, "          # 14", (ftnlen)14);
    e_wsfe();
    s_wsfe(&io___599);
    do_fio(&c__1, "     GEO --------> GSM", (ftnlen)22);
    e_wsfe();
    s_wsfe(&io___600);
    do_fio(&c__1, "          # 15", (ftnlen)14);
    e_wsfe();
    s_wsfe(&io___601);
    do_fio(&c__1, "     GEO --------> GSE", (ftnlen)22);
    e_wsfe();
    s_wsfe(&io___602);
    do_fio(&c__1, "          # 16", (ftnlen)14);
    e_wsfe();
    s_wsfe(&io___603);
    do_fio(&c__1, "     GEO --------> GSQ", (ftnlen)22);
    e_wsfe();

    s_wsfe(&io___604);
    e_wsfe();
    s_wsfe(&io___605);
    e_wsfe();
    s_wsfe(&io___606);
    do_fio(&c__1, "   input choosed test vector in GEO system :", (ftnlen)44);
    e_wsfe();
    s_wsfe(&io___607);
    e_wsfe();
    s_wsfe(&io___608);
    do_fio(&c__1, (char *)&g1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&gr, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&g2, (ftnlen)sizeof(real));
    r__1 = gt * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&g3, (ftnlen)sizeof(real));
    r__2 = gp * d__;
    do_fio(&c__1, (char *)&r__2, (ftnlen)sizeof(real));
    e_wsfe();

    s_wsfe(&io___609);
    e_wsfe();
    s_wsfe(&io___610);
    do_fio(&c__1, "1) t# 13 converts GEO system to SM  system by t_geo_to_sm",
	     (ftnlen)57);
    e_wsfe();
    t_geo_to_sm__(&g1, &g2, &g3, &v1, &v2, &v3);
    t_car_to_sph__(&v1, &v2, &v3, &vr, &vt, &vp);
    s_wsfe(&io___611);
    e_wsfe();
    s_wsfe(&io___612);
    do_fio(&c__1, (char *)&v1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&vr, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&v2, (ftnlen)sizeof(real));
    r__1 = vt * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&v3, (ftnlen)sizeof(real));
    r__2 = vp * d__;
    do_fio(&c__1, (char *)&r__2, (ftnlen)sizeof(real));
    e_wsfe();

    s_wsfe(&io___613);
    e_wsfe();
    s_wsfe(&io___614);
    do_fio(&c__1, "              come back to GEO system by t_sm_to_geo", (
	    ftnlen)52);
    e_wsfe();
    t_sm_to_geo__(&v1, &v2, &v3, &w1, &w2, &w3);
    t_car_to_sph__(&w1, &w2, &w3, &wr, &wt, &wp);
    s_wsfe(&io___615);
    e_wsfe();
    s_wsfe(&io___616);
    do_fio(&c__1, (char *)&w1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&wr, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&w2, (ftnlen)sizeof(real));
    r__1 = wt * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&w3, (ftnlen)sizeof(real));
    r__2 = wp * d__;
    do_fio(&c__1, (char *)&r__2, (ftnlen)sizeof(real));
    e_wsfe();

    cp_angle_and_ratio__(&g1, &g2, &g3, &w1, &w2, &w3, &a, &ra);
    s_wsfe(&io___617);
    do_fio(&c__1, "difference beetween first GEO vector:", (ftnlen)37);
    r__1 = a * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&ra, (ftnlen)sizeof(real));
    e_wsfe();


    s_wsfe(&io___618);
    e_wsfe();
    s_wsfe(&io___619);
    do_fio(&c__1, "2) t# 14 converts GEO system to GSM system by t_geo_to_gsm"
	    , (ftnlen)58);
    e_wsfe();
    t_geo_to_gsm__(&g1, &g2, &g3, &v1, &v2, &v3);
    t_car_to_sph__(&v1, &v2, &v3, &vr, &vt, &vp);
    s_wsfe(&io___620);
    e_wsfe();
    s_wsfe(&io___621);
    do_fio(&c__1, (char *)&v1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&vr, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&v2, (ftnlen)sizeof(real));
    r__1 = vt * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&v3, (ftnlen)sizeof(real));
    r__2 = vp * d__;
    do_fio(&c__1, (char *)&r__2, (ftnlen)sizeof(real));
    e_wsfe();

    s_wsfe(&io___622);
    e_wsfe();
    s_wsfe(&io___623);
    do_fio(&c__1, "              come back to GEO system by t_gsm_to_geo", (
	    ftnlen)53);
    e_wsfe();
    t_gsm_to_geo__(&v1, &v2, &v3, &w1, &w2, &w3);
    t_car_to_sph__(&w1, &w2, &w3, &wr, &wt, &wp);
    s_wsfe(&io___624);
    e_wsfe();
    s_wsfe(&io___625);
    do_fio(&c__1, (char *)&w1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&wr, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&w2, (ftnlen)sizeof(real));
    r__1 = wt * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&w3, (ftnlen)sizeof(real));
    r__2 = wp * d__;
    do_fio(&c__1, (char *)&r__2, (ftnlen)sizeof(real));
    e_wsfe();

    cp_angle_and_ratio__(&g1, &g2, &g3, &w1, &w2, &w3, &a, &ra);
    s_wsfe(&io___626);
    do_fio(&c__1, "difference beetween first GEO vector:", (ftnlen)37);
    r__1 = a * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&ra, (ftnlen)sizeof(real));
    e_wsfe();


    s_wsfe(&io___627);
    e_wsfe();
    s_wsfe(&io___628);
    do_fio(&c__1, "3) t# 15 converts GEO system to GSE system by t_geo_to_gse"
	    , (ftnlen)58);
    e_wsfe();
    t_geo_to_gse__(&g1, &g2, &g3, &v1, &v2, &v3);
    t_car_to_sph__(&v1, &v2, &v3, &vr, &vt, &vp);
    s_wsfe(&io___629);
    e_wsfe();
    s_wsfe(&io___630);
    do_fio(&c__1, (char *)&v1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&vr, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&v2, (ftnlen)sizeof(real));
    r__1 = vt * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&v3, (ftnlen)sizeof(real));
    r__2 = vp * d__;
    do_fio(&c__1, (char *)&r__2, (ftnlen)sizeof(real));
    e_wsfe();

    s_wsfe(&io___631);
    e_wsfe();
    s_wsfe(&io___632);
    do_fio(&c__1, "              come back to GEO system by t_gse_to_geo", (
	    ftnlen)53);
    e_wsfe();
    t_gse_to_geo__(&v1, &v2, &v3, &w1, &w2, &w3);
    t_car_to_sph__(&w1, &w2, &w3, &wr, &wt, &wp);
    s_wsfe(&io___633);
    e_wsfe();
    s_wsfe(&io___634);
    do_fio(&c__1, (char *)&w1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&wr, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&w2, (ftnlen)sizeof(real));
    r__1 = wt * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&w3, (ftnlen)sizeof(real));
    r__2 = wp * d__;
    do_fio(&c__1, (char *)&r__2, (ftnlen)sizeof(real));
    e_wsfe();

    cp_angle_and_ratio__(&g1, &g2, &g3, &w1, &w2, &w3, &a, &ra);
    s_wsfe(&io___635);
    do_fio(&c__1, "difference beetween first GEO vector:", (ftnlen)37);
    r__1 = a * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&ra, (ftnlen)sizeof(real));
    e_wsfe();


    s_wsfe(&io___636);
    e_wsfe();
    s_wsfe(&io___637);
    do_fio(&c__1, "4) t# 16 converts GEO system to GSQ system by t_geo_to_gs"
	    "eq", (ftnlen)59);
    e_wsfe();
    t_geo_to_gseq__(&g1, &g2, &g3, &v1, &v2, &v3);
    t_car_to_sph__(&v1, &v2, &v3, &vr, &vt, &vp);
    s_wsfe(&io___638);
    e_wsfe();
    s_wsfe(&io___639);
    do_fio(&c__1, (char *)&v1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&vr, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&v2, (ftnlen)sizeof(real));
    r__1 = vt * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&v3, (ftnlen)sizeof(real));
    r__2 = vp * d__;
    do_fio(&c__1, (char *)&r__2, (ftnlen)sizeof(real));
    e_wsfe();

    s_wsfe(&io___640);
    e_wsfe();
    s_wsfe(&io___641);
    do_fio(&c__1, "              come back to GEO system by t_gseq_to_geo", (
	    ftnlen)54);
    e_wsfe();
    t_gseq_to_geo__(&v1, &v2, &v3, &w1, &w2, &w3);
    t_car_to_sph__(&w1, &w2, &w3, &wr, &wt, &wp);
    s_wsfe(&io___642);
    e_wsfe();
    s_wsfe(&io___643);
    do_fio(&c__1, (char *)&w1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&wr, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&w2, (ftnlen)sizeof(real));
    r__1 = wt * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&w3, (ftnlen)sizeof(real));
    r__2 = wp * d__;
    do_fio(&c__1, (char *)&r__2, (ftnlen)sizeof(real));
    e_wsfe();

    cp_angle_and_ratio__(&g1, &g2, &g3, &w1, &w2, &w3, &a, &ra);
    s_wsfe(&io___644);
    do_fio(&c__1, "difference beetween first GEO vector:", (ftnlen)37);
    r__1 = a * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&ra, (ftnlen)sizeof(real));
    e_wsfe();

    s_wsfe(&io___645);
    do_fio(&c__1, "\f", (ftnlen)1);
    e_wsfe();
    s_wsfe(&io___646);
    do_fio(&c__1, "E) LOCAL SYSTEMS TRANSFORMATIONS", (ftnlen)32);
    e_wsfe();
    s_wsfe(&io___647);
    do_fio(&c__1, "   =============================", (ftnlen)32);
    e_wsfe();
    s_wsfe(&io___648);
    e_wsfe();
    s_wsfe(&io___649);
    do_fio(&c__1, "          # 17", (ftnlen)14);
    e_wsfe();
    s_wsfe(&io___650);
    do_fio(&c__1, "     GEO --------> DM", (ftnlen)21);
    e_wsfe();
    s_wsfe(&io___651);
    do_fio(&c__1, "          # 18", (ftnlen)14);
    e_wsfe();
    s_wsfe(&io___652);
    do_fio(&c__1, "     GEO --------> VDH", (ftnlen)22);
    e_wsfe();
    s_wsfe(&io___653);
    e_wsfe();
    s_wsfe(&io___654);
    do_fio(&c__1, "   first direction of observation in GEO system:", (ftnlen)
	    48);
    e_wsfe();
    s_wsfe(&io___655);
    do_fio(&c__1, "        geog. lat., long. =", (ftnlen)27);
    r__1 = grlat * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    r__2 = grlon * d__;
    do_fio(&c__1, (char *)&r__2, (ftnlen)sizeof(real));
    e_wsfe();

    s_wsfe(&io___656);
    e_wsfe();
    s_wsfe(&io___657);
    do_fio(&c__1, "1) t# 17 converts GEO system to DM  system by t_geo_to_dm",
	     (ftnlen)57);
    e_wsfe();
    t_geo_to_dm__(&g1, &g2, &g3, &grlat, &grlon, &v1, &v2, &v3);
    t_car_to_sph__(&v1, &v2, &v3, &vr, &vt, &vp);
    s_wsfe(&io___658);
    e_wsfe();
    s_wsfe(&io___659);
    do_fio(&c__1, (char *)&v1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&vr, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&v2, (ftnlen)sizeof(real));
    r__1 = vt * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&v3, (ftnlen)sizeof(real));
    r__2 = vp * d__;
    do_fio(&c__1, (char *)&r__2, (ftnlen)sizeof(real));
    e_wsfe();

    s_wsfe(&io___660);
    e_wsfe();
    s_wsfe(&io___661);
    do_fio(&c__1, "              come back to GEO system by t_dm_to_geo", (
	    ftnlen)52);
    e_wsfe();
    t_dm_to_geo__(&v1, &v2, &v3, &grlat, &grlon, &w1, &w2, &w3);
    t_car_to_sph__(&w1, &w2, &w3, &wr, &wt, &wp);
    s_wsfe(&io___662);
    e_wsfe();
    s_wsfe(&io___663);
    do_fio(&c__1, (char *)&w1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&wr, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&w2, (ftnlen)sizeof(real));
    r__1 = wt * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&w3, (ftnlen)sizeof(real));
    r__2 = wp * d__;
    do_fio(&c__1, (char *)&r__2, (ftnlen)sizeof(real));
    e_wsfe();

    cp_angle_and_ratio__(&g1, &g2, &g3, &w1, &w2, &w3, &a, &ra);
    s_wsfe(&io___664);
    do_fio(&c__1, "difference beetween first GEO vector:", (ftnlen)37);
    r__1 = a * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&ra, (ftnlen)sizeof(real));
    e_wsfe();


    s_wsfe(&io___665);
    e_wsfe();
    s_wsfe(&io___666);
    do_fio(&c__1, "2) t# 18 converts GEO system to VDH system by t_geo_to_vdh"
	    , (ftnlen)58);
    e_wsfe();
    t_geo_to_vdh__(&g1, &g2, &g3, &grlat, &grlon, &v1, &v2, &v3);
    t_car_to_sph__(&v1, &v2, &v3, &vr, &vt, &vp);
    s_wsfe(&io___667);
    e_wsfe();
    s_wsfe(&io___668);
    do_fio(&c__1, (char *)&v1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&vr, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&v2, (ftnlen)sizeof(real));
    r__1 = vt * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&v3, (ftnlen)sizeof(real));
    r__2 = vp * d__;
    do_fio(&c__1, (char *)&r__2, (ftnlen)sizeof(real));
    e_wsfe();

    s_wsfe(&io___669);
    e_wsfe();
    s_wsfe(&io___670);
    do_fio(&c__1, "              come back to GEO system by t_vdh_to_geo", (
	    ftnlen)53);
    e_wsfe();
    t_vdh_to_geo__(&v1, &v2, &v3, &grlat, &grlon, &w1, &w2, &w3);
    t_car_to_sph__(&w1, &w2, &w3, &wr, &wt, &wp);
    s_wsfe(&io___671);
    e_wsfe();
    s_wsfe(&io___672);
    do_fio(&c__1, (char *)&w1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&wr, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&w2, (ftnlen)sizeof(real));
    r__1 = wt * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&w3, (ftnlen)sizeof(real));
    r__2 = wp * d__;
    do_fio(&c__1, (char *)&r__2, (ftnlen)sizeof(real));
    e_wsfe();

    cp_angle_and_ratio__(&g1, &g2, &g3, &w1, &w2, &w3, &a, &ra);
    s_wsfe(&io___673);
    do_fio(&c__1, "difference beetween first GEO vector:", (ftnlen)37);
    r__1 = a * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&ra, (ftnlen)sizeof(real));
    e_wsfe();

    s_wsfe(&io___674);
    e_wsfe();
    s_wsfe(&io___675);
    e_wsfe();
    s_wsfe(&io___676);
    do_fio(&c__1, "   second direction of observation in GEO system:", (
	    ftnlen)49);
    e_wsfe();
    s_wsfe(&io___677);
    do_fio(&c__1, "        geog. lat., long. =", (ftnlen)27);
    r__1 = grlat * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    r__2 = grlon * d__;
    do_fio(&c__1, (char *)&r__2, (ftnlen)sizeof(real));
    e_wsfe();

    s_wsfe(&io___678);
    e_wsfe();
    s_wsfe(&io___679);
    do_fio(&c__1, "3) t# 17 converts GEO system to DM  system by t_geo_to_dm",
	     (ftnlen)57);
    e_wsfe();
    t_geo_to_dm__(&g1, &g2, &g3, &grlat, &grlon, &v1, &v2, &v3);
    t_car_to_sph__(&v1, &v2, &v3, &vr, &vt, &vp);
    s_wsfe(&io___680);
    e_wsfe();
    s_wsfe(&io___681);
    do_fio(&c__1, (char *)&v1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&vr, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&v2, (ftnlen)sizeof(real));
    r__1 = vt * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&v3, (ftnlen)sizeof(real));
    r__2 = vp * d__;
    do_fio(&c__1, (char *)&r__2, (ftnlen)sizeof(real));
    e_wsfe();

    s_wsfe(&io___682);
    e_wsfe();
    s_wsfe(&io___683);
    do_fio(&c__1, "              come back to GEO system by t_dm_to_geo", (
	    ftnlen)52);
    e_wsfe();
    t_dm_to_geo__(&v1, &v2, &v3, &grlat, &grlon, &w1, &w2, &w3);
    t_car_to_sph__(&w1, &w2, &w3, &wr, &wt, &wp);
    s_wsfe(&io___684);
    e_wsfe();
    s_wsfe(&io___685);
    do_fio(&c__1, (char *)&w1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&wr, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&w2, (ftnlen)sizeof(real));
    r__1 = wt * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&w3, (ftnlen)sizeof(real));
    r__2 = wp * d__;
    do_fio(&c__1, (char *)&r__2, (ftnlen)sizeof(real));
    e_wsfe();

    cp_angle_and_ratio__(&g1, &g2, &g3, &w1, &w2, &w3, &a, &ra);
    s_wsfe(&io___686);
    do_fio(&c__1, "difference beetween first GEO vector:", (ftnlen)37);
    r__1 = a * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&ra, (ftnlen)sizeof(real));
    e_wsfe();


    s_wsfe(&io___687);
    e_wsfe();
    s_wsfe(&io___688);
    do_fio(&c__1, "4) t# 18 converts GEO system to VDH system by t_geo_to_vdh"
	    , (ftnlen)58);
    e_wsfe();
    t_geo_to_vdh__(&g1, &g2, &g3, &grlat, &grlon, &v1, &v2, &v3);
    t_car_to_sph__(&v1, &v2, &v3, &vr, &vt, &vp);
    s_wsfe(&io___689);
    e_wsfe();
    s_wsfe(&io___690);
    do_fio(&c__1, (char *)&v1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&vr, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&v2, (ftnlen)sizeof(real));
    r__1 = vt * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&v3, (ftnlen)sizeof(real));
    r__2 = vp * d__;
    do_fio(&c__1, (char *)&r__2, (ftnlen)sizeof(real));
    e_wsfe();

    s_wsfe(&io___691);
    e_wsfe();
    s_wsfe(&io___692);
    do_fio(&c__1, "              come back to GEO system by t_vdh_to_geo", (
	    ftnlen)53);
    e_wsfe();
    t_vdh_to_geo__(&v1, &v2, &v3, &grlat, &grlon, &w1, &w2, &w3);
    t_car_to_sph__(&w1, &w2, &w3, &wr, &wt, &wp);
    s_wsfe(&io___693);
    e_wsfe();
    s_wsfe(&io___694);
    do_fio(&c__1, (char *)&w1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&wr, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&w2, (ftnlen)sizeof(real));
    r__1 = wt * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&w3, (ftnlen)sizeof(real));
    r__2 = wp * d__;
    do_fio(&c__1, (char *)&r__2, (ftnlen)sizeof(real));
    e_wsfe();

    cp_angle_and_ratio__(&g1, &g2, &g3, &w1, &w2, &w3, &a, &ra);
    s_wsfe(&io___695);
    do_fio(&c__1, "difference beetween first GEO vector:", (ftnlen)37);
    r__1 = a * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&ra, (ftnlen)sizeof(real));
    e_wsfe();

    s_wsfe(&io___696);
    do_fio(&c__1, "\f", (ftnlen)1);
    e_wsfe();
    s_wsfe(&io___697);
    e_wsfe();
    s_wsfe(&io___698);
    do_fio(&c__1, "F) LOCAL SPACECRAFT SYSTEMS TRANSFORMATIONS", (ftnlen)43);
    e_wsfe();
    s_wsfe(&io___699);
    do_fio(&c__1, "   ========================================", (ftnlen)43);
    e_wsfe();
    s_wsfe(&io___700);
    e_wsfe();
    s_wsfe(&io___701);
    do_fio(&c__1, "    # 15    19     21    -21    -19    -15", (ftnlen)42);
    e_wsfe();
    s_wsfe(&io___702);
    do_fio(&c__1, "  GEO -> GSE -> SR2 -> SR -> SR2 -> GSE -> GEO", (ftnlen)
	    46);
    e_wsfe();
    s_wsfe(&io___703);
    e_wsfe();
    s_wsfe(&io___704);
    e_wsfe();
    s_wsfe(&io___705);
    do_fio(&c__1, "   input choosed test vector in GEO system :", (ftnlen)44);
    e_wsfe();
    s_wsfe(&io___706);
    e_wsfe();
    s_wsfe(&io___707);
    do_fio(&c__1, (char *)&g1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&gr, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&g2, (ftnlen)sizeof(real));
    r__1 = gt * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&g3, (ftnlen)sizeof(real));
    r__2 = gp * d__;
    do_fio(&c__1, (char *)&r__2, (ftnlen)sizeof(real));
    e_wsfe();

    s_wsfe(&io___708);
    e_wsfe();
    s_wsfe(&io___709);
    do_fio(&c__1, "   input choosed spin axis in GSE system :", (ftnlen)42);
    e_wsfe();
    s_wsfe(&io___710);
    e_wsfe();
    s_wsfe(&io___711);
    do_fio(&c__1, (char *)&sa1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&sar, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&sa2, (ftnlen)sizeof(real));
    r__1 = sat * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&sa3, (ftnlen)sizeof(real));
    r__2 = sap * d__;
    do_fio(&c__1, (char *)&r__2, (ftnlen)sizeof(real));
    e_wsfe();

    s_wsfe(&io___712);
    e_wsfe();
    s_wsfe(&io___713);
    do_fio(&c__1, "   input choosed Spin frequency, Spin phase, dt:", (ftnlen)
	    48);
    e_wsfe();
    s_wsfe(&io___714);
    e_wsfe();
    s_wsfe(&io___715);
    do_fio(&c__1, "         Spin frequency= ", (ftnlen)25);
    do_fio(&c__1, (char *)&spifre, (ftnlen)sizeof(real));
    do_fio(&c__1, " Hz", (ftnlen)3);
    e_wsfe();
    s_wsfe(&io___716);
    do_fio(&c__1, "         Spin phase    = ", (ftnlen)25);
    r__1 = spipha * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, " deg.", (ftnlen)5);
    e_wsfe();
    s_wsfe(&io___717);
    do_fio(&c__1, "         delta T       = ", (ftnlen)25);
    do_fio(&c__1, (char *)&deltat, (ftnlen)sizeof(real));
    do_fio(&c__1, " sec.", (ftnlen)5);
    e_wsfe();
    s_wsfe(&io___718);
    e_wsfe();

    s_wsfe(&io___719);
    e_wsfe();
    s_wsfe(&io___720);
    do_fio(&c__1, "1) t# 15 converts GEO system to GSE system by t_geo_to_gse"
	    , (ftnlen)58);
    e_wsfe();
    t_geo_to_gse__(&g1, &g2, &g3, &w1, &w2, &w3);
    t_car_to_sph__(&w1, &w2, &w3, &wr, &wt, &wp);
    s_wsfe(&io___721);
    e_wsfe();
    s_wsfe(&io___722);
    do_fio(&c__1, (char *)&w1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&wr, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&w2, (ftnlen)sizeof(real));
    r__1 = wt * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&w3, (ftnlen)sizeof(real));
    r__2 = wp * d__;
    do_fio(&c__1, (char *)&r__2, (ftnlen)sizeof(real));
    e_wsfe();

    s_wsfe(&io___723);
    e_wsfe();
    s_wsfe(&io___724);
    do_fio(&c__1, "2) t# 19 converts GSE system to SR2 system by t_gse_to_sr2"
	    , (ftnlen)58);
    e_wsfe();
    t_gse_to_sr2__(&w1, &w2, &w3, &sa1, &sa2, &sa3, &v1, &v2, &v3);
    t_car_to_sph__(&v1, &v2, &v3, &vr, &vt, &vp);
    s_wsfe(&io___725);
    e_wsfe();
    s_wsfe(&io___726);
    do_fio(&c__1, (char *)&v1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&vr, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&v2, (ftnlen)sizeof(real));
    r__1 = vt * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&v3, (ftnlen)sizeof(real));
    r__2 = vp * d__;
    do_fio(&c__1, (char *)&r__2, (ftnlen)sizeof(real));
    e_wsfe();


    s_wsfe(&io___727);
    e_wsfe();
    s_wsfe(&io___728);
    do_fio(&c__1, "3) t#-21 converts SR2 system to SR  system by t_sr2_to_sr",
	     (ftnlen)57);
    e_wsfe();
    t_sr2_to_sr__(&v1, &v2, &spifre, &spipha, &deltat, &w1, &w2);
    w3 = v3;
    t_car_to_sph__(&w1, &w2, &w3, &wr, &wt, &wp);
    s_wsfe(&io___729);
    e_wsfe();
    s_wsfe(&io___730);
    do_fio(&c__1, (char *)&w1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&wr, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&w2, (ftnlen)sizeof(real));
    r__1 = wt * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&w3, (ftnlen)sizeof(real));
    r__2 = wp * d__;
    do_fio(&c__1, (char *)&r__2, (ftnlen)sizeof(real));
    e_wsfe();


    s_wsfe(&io___731);
    e_wsfe();
    s_wsfe(&io___732);
    do_fio(&c__1, "4) t# 21 converts SR  system to SR2 system by t_sr_to_sr2",
	     (ftnlen)57);
    e_wsfe();
    t_sr_to_sr2__(&w1, &w2, &spifre, &spipha, &deltat, &v1, &v2);
    v3 = w3;
    t_car_to_sph__(&v1, &v2, &v3, &vr, &vt, &vp);
    s_wsfe(&io___733);
    e_wsfe();
    s_wsfe(&io___734);
    do_fio(&c__1, (char *)&v1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&vr, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&v2, (ftnlen)sizeof(real));
    r__1 = vt * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&v3, (ftnlen)sizeof(real));
    r__2 = vp * d__;
    do_fio(&c__1, (char *)&r__2, (ftnlen)sizeof(real));
    e_wsfe();


    s_wsfe(&io___735);
    e_wsfe();
    s_wsfe(&io___736);
    do_fio(&c__1, "5) t#-19 converts SR2 system to GSE system by t_sr2_to_gse"
	    , (ftnlen)58);
    e_wsfe();
    t_sr2_to_gse__(&v1, &v2, &v3, &sa1, &sa2, &sa3, &w1, &w2, &w3);
    t_car_to_sph__(&w1, &w2, &w3, &wr, &wt, &wp);
    s_wsfe(&io___737);
    e_wsfe();
    s_wsfe(&io___738);
    do_fio(&c__1, (char *)&w1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&wr, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&w2, (ftnlen)sizeof(real));
    r__1 = wt * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&w3, (ftnlen)sizeof(real));
    r__2 = wp * d__;
    do_fio(&c__1, (char *)&r__2, (ftnlen)sizeof(real));
    e_wsfe();


    s_wsfe(&io___739);
    e_wsfe();
    s_wsfe(&io___740);
    do_fio(&c__1, "6) t#-15 converts GSE system to GEO system by t_gse_to_geo"
	    , (ftnlen)58);
    e_wsfe();
    t_gse_to_geo__(&w1, &w2, &w3, &v1, &v2, &v3);
    t_car_to_sph__(&v1, &v2, &v3, &vr, &vt, &vp);
    s_wsfe(&io___741);
    e_wsfe();
    s_wsfe(&io___742);
    do_fio(&c__1, (char *)&v1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&vr, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&v2, (ftnlen)sizeof(real));
    r__1 = vt * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&v3, (ftnlen)sizeof(real));
    r__2 = vp * d__;
    do_fio(&c__1, (char *)&r__2, (ftnlen)sizeof(real));
    e_wsfe();

    cp_angle_and_ratio__(&g1, &g2, &g3, &v1, &v2, &v3, &a, &ra);
    s_wsfe(&io___743);
    do_fio(&c__1, "difference beetveen first GEO wector:", (ftnlen)37);
    r__1 = a * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&ra, (ftnlen)sizeof(real));
    e_wsfe();

    s_wsfe(&io___744);
    e_wsfe();
    s_wsfe(&io___745);
    e_wsfe();
    s_wsfe(&io___746);
    do_fio(&c__1, "G) MAGNETIC FIELD ALIGNED TRANSFORMATION", (ftnlen)40);
    e_wsfe();
    s_wsfe(&io___747);
    do_fio(&c__1, "   ======================================", (ftnlen)41);
    e_wsfe();
    s_wsfe(&io___748);
    e_wsfe();
    s_wsfe(&io___749);
    do_fio(&c__1, "  GSE -> MFA", (ftnlen)12);
    e_wsfe();
    s_wsfe(&io___750);
    do_fio(&c__1, "  SR2 -> MFA", (ftnlen)12);
    e_wsfe();

    s_wsfe(&io___751);
    e_wsfe();
    s_wsfe(&io___752);
    do_fio(&c__1, "   input choosed test vector in GSE system :", (ftnlen)44);
    e_wsfe();
    s_wsfe(&io___753);
    e_wsfe();
    s_wsfe(&io___754);
    do_fio(&c__1, (char *)&w1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&wr, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&w2, (ftnlen)sizeof(real));
    r__1 = wt * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&w3, (ftnlen)sizeof(real));
    r__2 = wp * d__;
    do_fio(&c__1, (char *)&r__2, (ftnlen)sizeof(real));
    e_wsfe();

    s_wsfe(&io___755);
    e_wsfe();
    s_wsfe(&io___756);
    do_fio(&c__1, "   input choosed DC mag. field in GSE system :", (ftnlen)
	    46);
    e_wsfe();
    b1 = 30.f;
    b2 = 50.f;
    b3 = 200.f;
    t_car_to_sph__(&b1, &b2, &b3, &br, &bt, &bp);
    s_wsfe(&io___763);
    e_wsfe();
    s_wsfe(&io___764);
    do_fio(&c__1, (char *)&b1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&br, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&b2, (ftnlen)sizeof(real));
    r__1 = bt * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&b3, (ftnlen)sizeof(real));
    r__2 = bp * d__;
    do_fio(&c__1, (char *)&r__2, (ftnlen)sizeof(real));
    e_wsfe();

    s_wsfe(&io___765);
    e_wsfe();
    s_wsfe(&io___766);
    e_wsfe();
    s_wsfe(&io___767);
    do_fio(&c__1, "1) t# 22 converts GSE system to MFA system by t_gse_to_mfa"
	    , (ftnlen)58);
    e_wsfe();
    t_gse_to_mfa__(&w1, &w2, &w3, &b1, &b2, &b3, &v1, &v2, &v3);
    t_car_to_sph__(&v1, &v2, &v3, &vr, &vt, &vp);
    s_wsfe(&io___768);
    e_wsfe();
    s_wsfe(&io___769);
    do_fio(&c__1, (char *)&v1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&vr, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&v2, (ftnlen)sizeof(real));
    r__1 = vt * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&v3, (ftnlen)sizeof(real));
    r__2 = vp * d__;
    do_fio(&c__1, (char *)&r__2, (ftnlen)sizeof(real));
    e_wsfe();

    s_wsfe(&io___770);
    e_wsfe();
    s_wsfe(&io___771);
    do_fio(&c__1, "   input choosed test vector in SR2 system :", (ftnlen)44);
    e_wsfe();
    t_gse_to_sr2__(&w1, &w2, &w3, &sa1, &sa2, &sa3, &v1, &v2, &v3);
    t_car_to_sph__(&v1, &v2, &v3, &vr, &vt, &vp);
    s_wsfe(&io___772);
    e_wsfe();
    s_wsfe(&io___773);
    do_fio(&c__1, (char *)&v1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&vr, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&v2, (ftnlen)sizeof(real));
    r__1 = vt * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&v3, (ftnlen)sizeof(real));
    r__2 = vp * d__;
    do_fio(&c__1, (char *)&r__2, (ftnlen)sizeof(real));
    e_wsfe();

    s_wsfe(&io___774);
    e_wsfe();
    s_wsfe(&io___775);
    do_fio(&c__1, "   input choosed DC mag. field in SR2 system :", (ftnlen)
	    46);
    e_wsfe();
    t_gse_to_sr2__(&b1, &b2, &b3, &sa1, &sa2, &sa3, &bs1, &bs2, &bs3);
    t_car_to_sph__(&bs1, &bs2, &bs3, &bsr, &bst, &bsp);
    s_wsfe(&io___782);
    e_wsfe();
    s_wsfe(&io___783);
    do_fio(&c__1, (char *)&bs1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&bsr, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&bs2, (ftnlen)sizeof(real));
    r__1 = bst * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&bs3, (ftnlen)sizeof(real));
    r__2 = bsp * d__;
    do_fio(&c__1, (char *)&r__2, (ftnlen)sizeof(real));
    e_wsfe();

    s_wsfe(&io___784);
    e_wsfe();
    s_wsfe(&io___785);
    do_fio(&c__1, "2) t# 20 converts SR2 system to MFA system by t_sr2_to_mfa"
	    , (ftnlen)58);
    e_wsfe();
    t_sr2_to_mfa__(&v1, &v2, &v3, &bs1, &bs2, &bs3, &sa1, &sa2, &sa3, &w1, &
	    w2, &w3);
    t_car_to_sph__(&w1, &w2, &w3, &wr, &wt, &wp);
    s_wsfe(&io___786);
    e_wsfe();
    s_wsfe(&io___787);
    do_fio(&c__1, (char *)&w1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&wr, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&w2, (ftnlen)sizeof(real));
    r__1 = wt * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&w3, (ftnlen)sizeof(real));
    r__2 = wp * d__;
    do_fio(&c__1, (char *)&r__2, (ftnlen)sizeof(real));
    e_wsfe();

    s_wsfe(&io___788);
    do_fio(&c__1, "\f", (ftnlen)1);
    e_wsfe();
    s_wsfe(&io___789);
    e_wsfe();
    s_wsfe(&io___790);
    do_fio(&c__1, "H) PSEUDO-MAGNETOPAUSE SYSTEMS", (ftnlen)30);
    e_wsfe();
    s_wsfe(&io___791);
    do_fio(&c__1, "   ===========================", (ftnlen)30);
    e_wsfe();
    s_wsfe(&io___792);
    e_wsfe();
    s_wsfe(&io___793);
    do_fio(&c__1, "       GSE ----> TPN", (ftnlen)20);
    e_wsfe();
    s_wsfe(&io___794);
    do_fio(&c__1, "       GSM ----> TPN", (ftnlen)20);
    e_wsfe();

    s_wsfe(&io___795);
    e_wsfe();
    s_wsfe(&io___796);
    do_fio(&c__1, "1) Computation of the TPN system", (ftnlen)32);
    e_wsfe();

    s_wsfe(&io___797);
    e_wsfe();
    s_wsfe(&io___798);
    do_fio(&c__1, "   Input position in GSE system, in RE:", (ftnlen)39);
    e_wsfe();
    s_wsfe(&io___799);
    do_fio(&c__1, "   (supposed to be on the magnetopause)", (ftnlen)39);
    e_wsfe();

    xo = 5.72f;
    yo = 8.23f;
    zo = 9.22f;

    s_wsfe(&io___803);
    e_wsfe();
    s_wsfe(&io___804);
    do_fio(&c__1, "   Xo, Yo, Zo in GSE =", (ftnlen)22);
    do_fio(&c__1, (char *)&xo, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&yo, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&zo, (ftnlen)sizeof(real));
    e_wsfe();

    xs = 12.f;

    s_wsfe(&io___806);
    e_wsfe();
    s_wsfe(&io___807);
    do_fio(&c__1, "   abcisse of subsolar point for paraboloid model:", (
	    ftnlen)50);
    do_fio(&c__1, (char *)&xs, (ftnlen)sizeof(real));
    e_wsfe();

    cp_tpn_param__(&xo, &yo, &zo, &xs, &tx, &ty, &tz, &px, &py, &pz, &nx, &ny,
	     &nz);

    s_wsfe(&io___817);
    e_wsfe();
    s_wsfe(&io___818);
    do_fio(&c__1, "   TPN vectors in GSE:", (ftnlen)22);
    e_wsfe();
    s_wsfe(&io___819);
    e_wsfe();
    s_wsfe(&io___820);
    do_fio(&c__1, "   Tx,Ty,Tz=", (ftnlen)12);
    do_fio(&c__1, (char *)&tx, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&ty, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&tz, (ftnlen)sizeof(real));
    e_wsfe();
    s_wsfe(&io___821);
    do_fio(&c__1, "   Px,Py,Pz=", (ftnlen)12);
    do_fio(&c__1, (char *)&px, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&py, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&pz, (ftnlen)sizeof(real));
    e_wsfe();
    s_wsfe(&io___822);
    do_fio(&c__1, "   Nx,Ny,Nz=", (ftnlen)12);
    do_fio(&c__1, (char *)&nx, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&ny, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&nz, (ftnlen)sizeof(real));
    e_wsfe();
    s_wsfe(&io___823);
    do_fio(&c__1, " ", (ftnlen)1);
    e_wsfe();

    s_wsfe(&io___824);
    e_wsfe();
    s_wsfe(&io___825);
    do_fio(&c__1, "   input choosed test vector in GSE system :", (ftnlen)44);
    e_wsfe();

    v1 = nx;
    v2 = ny;
    v3 = nz;

    t_car_to_sph__(&v1, &v2, &v3, &vr, &vt, &vp);
    s_wsfe(&io___826);
    e_wsfe();
    s_wsfe(&io___827);
    do_fio(&c__1, (char *)&v1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&vr, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&v2, (ftnlen)sizeof(real));
    r__1 = vt * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&v3, (ftnlen)sizeof(real));
    r__2 = vp * d__;
    do_fio(&c__1, (char *)&r__2, (ftnlen)sizeof(real));
    e_wsfe();

    s_wsfe(&io___828);
    do_fio(&c__1, " ", (ftnlen)1);
    e_wsfe();
    s_wsfe(&io___829);
    do_fio(&c__1, "2) t# 23 convert GSE system to TPN system by t_gse_to_tpn",
	     (ftnlen)57);
    e_wsfe();

    t_gse_to_tpn__(&v1, &v2, &v3, &xo, &yo, &zo, &xs, &w1, &w2, &w3);

    t_car_to_sph__(&w1, &w2, &w3, &wr, &wt, &wp);
    s_wsfe(&io___830);
    e_wsfe();
    s_wsfe(&io___831);
    do_fio(&c__1, (char *)&w1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&wr, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&w2, (ftnlen)sizeof(real));
    r__1 = wt * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&w3, (ftnlen)sizeof(real));
    r__2 = wp * d__;
    do_fio(&c__1, (char *)&r__2, (ftnlen)sizeof(real));
    e_wsfe();

    s_wsfe(&io___832);
    do_fio(&c__1, " ", (ftnlen)1);
    e_wsfe();
    s_wsfe(&io___833);
    do_fio(&c__1, "3) t# 9 convert GSE system to GSM system by t_gse_to_gsm", 
	    (ftnlen)56);
    e_wsfe();

    t_gse_to_gsm__(&v1, &v2, &v3, &w1, &w2, &w3);
    t_car_to_sph__(&w1, &w2, &w3, &wr, &wt, &wp);
    s_wsfe(&io___834);
    e_wsfe();
    s_wsfe(&io___835);
    do_fio(&c__1, (char *)&w1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&wr, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&w2, (ftnlen)sizeof(real));
    r__1 = wt * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&w3, (ftnlen)sizeof(real));
    r__2 = wp * d__;
    do_fio(&c__1, (char *)&r__2, (ftnlen)sizeof(real));
    e_wsfe();

    s_wsfe(&io___836);
    do_fio(&c__1, " ", (ftnlen)1);
    e_wsfe();
    s_wsfe(&io___837);
    do_fio(&c__1, "4) t# 9 convert GSE position of S/C to GSM position", (
	    ftnlen)51);
    e_wsfe();

    t_gse_to_gsm__(&xo, &yo, &zo, &xm, &ym, &zm);

    s_wsfe(&io___841);
    e_wsfe();
    s_wsfe(&io___842);
    do_fio(&c__1, "   Xo, Yo, Zo in GSM =", (ftnlen)22);
    do_fio(&c__1, (char *)&xm, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&ym, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&zm, (ftnlen)sizeof(real));
    e_wsfe();

    s_wsfe(&io___843);
    do_fio(&c__1, " ", (ftnlen)1);
    e_wsfe();
    s_wsfe(&io___844);
    do_fio(&c__1, "5) t# 24 convert GSM system to TPN system by t_gsm_to_tpn",
	     (ftnlen)57);
    e_wsfe();

    t_gsm_to_tpn__(&w1, &w2, &w3, &xm, &ym, &zm, &xs, &v1, &v2, &v3);

    t_car_to_sph__(&v1, &v2, &v3, &vr, &vt, &vp);
    s_wsfe(&io___845);
    e_wsfe();
    s_wsfe(&io___846);
    do_fio(&c__1, (char *)&v1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&vr, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&v2, (ftnlen)sizeof(real));
    r__1 = vt * d__;
    do_fio(&c__1, (char *)&r__1, (ftnlen)sizeof(real));
    do_fio(&c__1, (char *)&v3, (ftnlen)sizeof(real));
    r__2 = vp * d__;
    do_fio(&c__1, (char *)&r__2, (ftnlen)sizeof(real));
    e_wsfe();


    s_wsfe(&io___847);
    e_wsfe();
    s_wsfe(&io___848);
    do_fio(&c__1, "===============================", (ftnlen)31);
    e_wsfe();
    s_wsfe(&io___849);
    do_fio(&c__1, "END OF rocot_check test program", (ftnlen)31);
    e_wsfe();
    s_wsfe(&io___850);
    do_fio(&c__1, "===============================", (ftnlen)31);
    e_wsfe();



    s_stop("rocot_check.exe: normal termination", (ftnlen)35);
    return 0;
} /* MAIN__ */


/*     XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX0XX */

/* Subroutine */ int strlen_(char *str, integer *nc, ftnlen str_len)
{
    /* System generated locals */
    integer i__1;

    /* Builtin functions */
    integer i_len(char *, ftnlen);

    /* Local variables */
    static integer i__, ii, nbc;


/* ---------------------------------------------------------------------- */

/* *   Class  : basic compute modules of Rocotlib Software */
/* *   Object : compute usefull lengh of a string (ignore last blanks) */
/* *   Author : P. Robert, CRPE, 1992 */

/* ---------------------------------------------------------------------- */


    nbc = i_len(str, str_len);

/* *** nc= longueur au de la de laquelle il n'y a que des blancs */
/*         (mais la chaine 1:nc peu contenir des blancs) */

    i__1 = nbc;
    for (i__ = 1; i__ <= i__1; ++i__) {
	ii = nbc - i__ + 1;
	if (*(unsigned char *)&str[ii - 1] != ' ') {
	    goto L12;
	}
/* L10: */
    }
    ii = 0;
L12:
    *nc = ii;

    return 0;
} /* strlen_ */

/* Main program alias */ int rocot_check__ () { MAIN__ (); return 0; }
