Using rwo files as input to orbfit

warthogastronomer Mar 13, 2014

Hi All

If you have any experience with orbfit (4.2) and think you might be able to help, please read on!

As a simple test of the code, I am trying to reproduce the orbital elements of object 88611 contained in the astdys database. I have downloaded 88611.rwo (copied below) from the astdys page, as well as 88611.cat (copied below) to be used as input orbital elements for the fit, and have created 88611.oop for options to orbfit.

When I run > orbfit.x and provide 88611 as the name, it succeeds in finding a fit. But on inspection, it overwrites the rwo file, and ignored the BIAS and RMS values contained in the rwo file, rather than using them in the fit. It seems it is setting BIAS=0 and RMS=1 for all entries in the rwo file.

Does anyone know how to get orbfit to use the BIAS and RMS values in the input .rwo file?

Note: -changing the value of .print_res of reject.auto has no affect on the screen output or fit results.
Thanks

Wes

88611.oop:
object1.
.name = '88611'
.obs_fname='88611'
input_files.
.incond=88611.cat[OEF]
operations.print_res = 1
reject.auto = F

88611.cat:
format  = 'OEF1.1'       ! file format
rectype = 'ML'           ! record type (1L/ML)
refsys  = ECLM J2000     ! default reference system
END_OF_HEADER
88611
! Equinoctial elements: a, e*sin(LP), e*cos(LP), tan(i/2)*sin(LN), tan(i/2)*cos(LN), mean long.
 EQU   4.38401434576208E+01  -0.000023322945265  -0.031656954489436   -0.018542187353306    0.012881306692064 325.5869998303674
 MJD     53245.984115890 TDT
 MAG   5.609  0.150
! RMS    1.02718E-02   2.41615E-04   9.69533E-05   2.20197E-06   3.44174E-06   2.15932E-02
! EIG   2.48003E-07   2.98351E-07   3.93499E-06   4.91329E-06   2.67184E-04   1.02785E-02
! WEA   0.99934  -0.02138   0.00729  -0.00000   0.00001  -0.02830
 COV   1.055098072121325E-04 -2.255990021951618E-06  7.701782146877127E-07
 COV  -4.092556660651663E-10  5.811367808792076E-10 -1.710475317183444E-04
 COV   5.837761634225183E-08 -1.030414813050828E-08  1.296712438419433E-11
 COV  -1.873658231300820E-11  5.041125374767595E-06  9.399935164594918E-09
 COV   1.411563135977395E-12 -2.488459635842624E-12 -4.054502896438521E-07
 COV   4.848693583168296E-12 -7.502004305630544E-12  1.347632097175391E-09
 COV   1.184558288861564E-11 -1.971820096594674E-09  4.662660329835423E-04
 NOR   1.343350999389266E+07  7.041845944463621E+08 -4.619712491975347E+08
 NOR   7.811075870347890E+07 -1.067815715215493E+08 -3.087800084880433E+06
 NOR   7.190884053201854E+12  4.657358661712934E+12  1.404857923289951E+12
 NOR   9.832430322081732E+11 -7.343731613182695E+10  3.079631275556257E+12
 NOR   9.140654106469681E+11  6.508671883871273E+11 -4.784534456932166E+10
 NOR   1.052818618979664E+13  6.686351797772004E+12 -1.436754889260029E+10
 NOR   4.334162490267313E+12 -1.010471716514754E+10  7.512451285687269E+08

88611.rwo:
version =   2
errmod  = 'cbm10'
RMSast  =   2.92790E-01
RMSmag  =   7.01925E-01
END_OF_HEADER
! Object   Obser ============= Date ============= ================== Right Ascension =================  ================= Declination ===================== ==== Magnitude ==== Ast Obs  Residual SEL
! Design   K T N YYYY MM DD.dddddddddd   Accuracy HH MM SS.sss  Accuracy      RMS  F     Bias    Resid sDD MM SS.ss  Accuracy      RMS  F     Bias    Resid Val  B   RMS  Resid Cat Cod       Chi A M
 88611     O C 1 2000 08 01.29094       1.000E-05 21 40 29.660  1.500E-01    1.500 F    0.000    0.255 -13 03 36.80  1.000E-01    1.500 F    0.000    0.097 21.8 R  0.70   0.04     807      0.19 1 1
 88611     O C 1 2000 08 01.36718       1.000E-05 21 40 29.360  1.500E-01    1.500 F    0.000    0.683 -13 03 38.30  1.000E-01    1.500 F    0.000    0.162                         807      0.49 1 0
 88611     O C 1 2001 08 20.12602       1.000E-05 21 43 45.440  1.500E-01    1.020 F   -0.115    0.170 -12 43 33.90  1.000E-01    0.790 F    0.256    0.037 21.0 R  0.70  -0.72   c 807      0.18 1 1
 88611     O C 1 2001 08 20.20431       1.000E-05 21 43 45.090  1.500E-01    1.020 F   -0.115    0.151 -12 43 35.80  1.000E-01    0.790 F    0.256   -0.144                       c 807      0.24 1 0
 88611     O C 1 2001 08 21.13965       1.000E-05 21 43 41.000  1.500E-01    1.020 F   -0.115    0.160 -12 43 56.30  1.000E-01    0.790 F    0.256   -0.143 21.0 R  0.70  -0.72   c 807      0.24 1 1
 88611     O C 1 2001 08 21.21386       1.000E-05 21 43 40.670  1.500E-01    1.020 F   -0.115    0.161 -12 43 58.00  1.000E-01    0.790 F    0.256   -0.216                       c 807      0.32 1 0
 88611     O C 1 2001 09 12.23025       1.000E-05 21 42 08.900  1.500E-01    1.020 F   -0.014    0.245 -12 51 41.20  1.000E-01    0.790 F    0.232    0.053 21.8 R  0.70   0.00   c 695      0.26 1 1
 88611     O C 1 2001 09 12.27730       1.000E-05 21 42 08.710  1.500E-01    1.020 F   -0.014    0.162 -12 51 42.10  1.000E-01    0.790 F    0.232    0.072                       c 695      0.19 1 0
 88611     O C 1 2001 09 18.18911       1.000E-05 21 41 46.770  1.500E-01    1.020 F   -0.014    0.569 -12 53 34.00  1.000E-01    0.790 F    0.232    0.274 22.6 R  0.70   0.78   c 688      0.68 1 1
 88611     O C 1 2001 09 18.27980       1.000E-05 21 41 46.430  1.500E-01    1.020 F   -0.014    0.434 -12 53 36.40  1.000E-01    0.790 F    0.232   -0.468                       c 688      0.74 1 0
 88611     O C 1 2001 09 19.16569       1.000E-05 21 41 43.270  1.500E-01    1.020 F   -0.014    0.157 -12 53 52.10  1.000E-01    0.790 F    0.232   -0.051 21.8 R  0.70  -0.02   c 688      0.17 1 1
 88611     O C 1 2001 09 19.29448       1.000E-05 21 41 42.800  1.500E-01    1.020 F   -0.014    0.059 -12 53 54.30  1.000E-01    0.790 F    0.232    0.075                       c 688      0.11 1 0
 88611     O C   2001 11 03.01134       1.000E-05 21 40 11.740  1.500E-01    1.020 F   -0.004   -0.162 -13 01 57.40  1.000E-01    0.790 F    0.293   -0.390 22.0    0.70  -1.10   c 304      0.53 1 1
 88611     O C   2001 11 03.02280       1.000E-05 21 40 11.730  1.500E-01    1.020 F   -0.004   -0.252 -13 01 57.20  1.000E-01    0.790 F    0.293   -0.167                       c 304      0.34 1 0
 88611     O C 5 2002 08 03.98333       1.000E-05 21 49 31.800  1.500E-01    1.020 F   -0.009   -0.325 -12 10 41.00  1.000E-01    0.790 F    0.334    0.160 22.5 R  0.70   0.74   c 950      0.38 1 1
 88611     O C 5 2002 08 03.99122       1.000E-05 21 49 31.750  1.500E-01    1.020 F   -0.009   -0.562 -12 10 41.80  1.000E-01    0.790 F    0.334   -0.471 21.8 R  0.70   0.04   c 950      0.82 1 1
 88611     O C 5 2002 08 04.00021       1.000E-05 21 49 31.720  1.500E-01    1.020 F   -0.009   -0.437 -12 10 42.10  1.000E-01    0.790 F    0.334   -0.580 21.9 R  0.70   0.14   c 950      0.86 1 1
 88611     O C 5 2002 08 04.00723       1.000E-05 21 49 31.670  1.500E-01    1.020 F   -0.009   -0.729 -12 10 42.10  1.000E-01    0.790 F    0.334   -0.430 21.9 R  0.70   0.14   c 950      0.91 1 1
 88611     O C 3 2002 09 05.35971       1.000E-05 21 47 12.430  1.500E-01    1.170 F   -0.305    0.047 -12 22 30.40  1.000E-01    1.020 F    0.052   -0.549                       b 568      0.54 1 0
 88611     O C 3 2002 09 07.35630       1.000E-05 21 47 04.270  1.500E-01    1.170 F   -0.305    0.149 -12 23 12.10  1.000E-01    1.020 F    0.052   -0.171                       b 568      0.21 1 0
 88611     O C 3 2002 09 07.42880       1.000E-05 21 47 03.970  1.500E-01    1.170 F   -0.305    0.150 -12 23 13.60  1.000E-01    1.020 F    0.052   -0.155                       b 568      0.20 1 0
 88611     O C   2002 09 09.20428       1.000E-05 21 46 56.790  1.500E-01    1.020 F   -0.009   -0.662 -12 23 50.20  1.000E-01    0.790 F    0.334   -0.074 19.5 R  0.70  -2.29   c 662      0.67 1 0
 88611     O C   2002 09 09.29464       1.000E-05 21 46 56.480  1.500E-01    1.020 F   -0.009    0.176 -12 23 50.70  1.000E-01    0.790 F    0.334    1.295                       c 662      1.67 1 0
 88611     O C   2003 09 03.07624       1.000E-05 21 51 56.960  1.500E-01    0.990 F   -0.108   -0.240 -11 54 31.00  1.000E-01    0.810 F    0.070    0.323 21.8 R  0.70   0.03   o 304      0.47 1 1
 88611     O C   2003 09 03.12404       1.000E-05 21 51 56.760  1.500E-01    0.990 F   -0.108   -0.182 -11 54 32.70  1.000E-01    0.810 F    0.070   -0.326                       o 304      0.45 1 0
 88611     O C   2003 10 23.02215       1.000E-05 21 49 23.950  1.500E-01    0.990 F   -0.093   -0.159 -12 08 10.60  1.000E-01    0.810 F    0.051    0.150 22.0 R  0.70   0.12   o 304      0.25 1 1
 88611     O C   2003 10 23.15818       1.000E-05 21 49 23.760  1.500E-01    0.990 F   -0.093    0.038 -12 08 11.70  1.000E-01    0.810 F    0.051    0.157                       o 304      0.20 1 0
 88611     O C   2003 10 24.01169       1.000E-05 21 49 22.600  1.500E-01    0.990 F   -0.093    0.288 -12 08 18.70  1.000E-01    0.810 F    0.051    0.021 21.9 R  0.70   0.01   o 304      0.31 1 1
 88611     O C   2003 10 24.13110       1.000E-05 21 49 22.430  1.500E-01    0.990 F   -0.093    0.279 -12 08 19.60  1.000E-01    0.810 F    0.051    0.047                       o 304      0.31 1 0
 88611     O C   2004 05 25.30093       1.000E-05 22 01 19.280  1.500E-01    0.990 F   -0.176   -0.019 -11 03 21.80  1.000E-01    0.810 F    0.094    0.146 22.3 R  0.70   0.39   o 304      0.18 1 1
 88611     O C   2004 05 25.30653       1.000E-05 22 01 19.280  1.500E-01    0.990 F   -0.176   -0.035 -11 03 21.90  1.000E-01    0.810 F    0.094    0.029                       o 304      0.05 1 0
 88611     O C   2004 05 25.31029       1.000E-05 22 01 19.280  1.500E-01    0.990 F   -0.176   -0.045 -11 03 21.80  1.000E-01    0.810 F    0.094    0.117                       o 304      0.15 1 0
 88611     O C   2004 05 25.31405       1.000E-05 22 01 19.260  1.500E-01    0.990 F   -0.176   -0.350 -11 03 22.10  1.000E-01    0.810 F    0.094   -0.194                       o 304      0.44 1 0
 88611     O C   2004 05 25.31781       1.000E-05 22 01 19.280  1.500E-01    0.990 F   -0.176   -0.066 -11 03 21.70  1.000E-01    0.810 F    0.094    0.194                       o 304      0.25 1 0
 88611     O C   2004 05 25.32157       1.000E-05 22 01 19.280  1.500E-01    0.990 F   -0.176   -0.077 -11 03 21.80  1.000E-01    0.810 F    0.094    0.082                       o 304      0.13 1 0
 88611     O C   2004 05 25.32534       1.000E-05 22 01 19.300  1.500E-01    0.990 F   -0.176    0.208 -11 03 21.60  1.000E-01    0.810 F    0.094    0.271                       o 304      0.40 1 0
 88611     O C   2004 05 25.32910       1.000E-05 22 01 19.290  1.500E-01    0.990 F   -0.176    0.050 -11 03 21.70  1.000E-01    0.810 F    0.094    0.159                       o 304      0.21 1 0
 88611     O C   2004 05 25.33286       1.000E-05 22 01 19.280  1.500E-01    0.990 F   -0.176   -0.107 -11 03 21.80  1.000E-01    0.810 F    0.094    0.047                       o 304      0.12 1 0
 88611     O C   2004 05 25.33699       1.000E-05 22 01 19.270  1.500E-01    0.990 F   -0.176   -0.265 -11 03 21.60  1.000E-01    0.810 F    0.094    0.234                       o 304      0.40 1 0
 88611     O C   2004 05 25.34873       1.000E-05 22 01 19.280  1.500E-01    0.990 F   -0.176   -0.149 -11 03 21.80  1.000E-01    0.810 F    0.094   -0.002                       o 304      0.15 1 0
 88611     O C   2004 05 25.35302       1.000E-05 22 01 19.290  1.500E-01    0.990 F   -0.176   -0.013 -11 03 21.60  1.000E-01    0.810 F    0.094    0.184                       o 304      0.23 1 0
 88611     O C   2004 05 26.29405       1.000E-05 22 01 19.470  1.500E-01    0.990 F   -0.176   -0.153 -11 03 18.90  1.000E-01    0.810 F    0.094    0.037                       o 304      0.16 1 0
 88611     O C   2004 05 26.30231       1.000E-05 22 01 19.470  1.500E-01    0.990 F   -0.176   -0.166 -11 03 18.60  1.000E-01    0.810 F    0.094    0.315                       o 304      0.43 1 0
 88611     O C   2004 05 27.26691       1.000E-05 22 01 19.600  1.500E-01    0.990 F   -0.176    0.056 -11 03 16.30  1.000E-01    0.810 F    0.094    0.122 21.6 R  0.70  -0.31   o 304      0.16 1 1
 88611     O C   2004 05 27.37124       1.000E-05 22 01 19.590  1.500E-01    0.990 F   -0.176   -0.119 -11 03 16.10  1.000E-01    0.810 F    0.094    0.092                       o 304      0.17 1 0
 88611     O C   2004 05 29.31226       1.000E-05 22 01 19.570  1.500E-01    0.990 F   -0.176   -0.112 -11 03 12.50  1.000E-01    0.810 F    0.094    0.055                       o 304      0.13 1 0
 88611     O C   2004 05 29.32355       1.000E-05 22 01 19.570  1.500E-01    0.990 F   -0.176   -0.087 -11 03 12.40  1.000E-01    0.810 F    0.094    0.141                       o 304      0.20 1 0
 88611     O P   2004 09 13.02868       1.000E-05 21 55 46.340  1.500E-01    0.990 F   -0.100   -0.184 -11 31 04.00  1.000E-01    0.810 F    0.101   -0.290 22.2 R  0.70   0.40   o 304      0.42 1 1
 88611     O P   2004 09 13.26652       1.000E-05 21 55 45.400  1.500E-01    0.990 F   -0.100   -0.119 -11 31 08.90  1.000E-01    0.810 F    0.101   -0.224                       o 304      0.31 1 0
 88611     O C   2005 06 07.31697       1.000E-05 22 05 42.610  1.500E-01    0.990 F   -0.091    0.351 -10 35 52.50  1.000E-01    0.810 F    0.122   -0.009 22.5 R  0.70   0.60   o 304      0.36 1 1
 88611     O C   2005 06 07.34192       1.000E-05 22 05 42.570  1.500E-01    0.990 F   -0.091    0.047 -10 35 52.40  1.000E-01    0.810 F    0.122    0.142                       o 304      0.18 1 0
 88611     O C   2005 07 10.35703       1.000E-05 22 04 38.160  1.500E-01    0.990 F   -0.091    0.236 -10 40 46.90  1.000E-01    0.810 F    0.122    0.030 22.9 R  0.70   1.06   o 304      0.25 1 1
 88611     O C   2005 07 10.43502       1.000E-05 22 04 37.920  1.500E-01    0.990 F   -0.091    0.292 -10 40 48.10  1.000E-01    0.810 F    0.122    0.011                       o 304      0.30 1 0
 88611     O C   2005 07 11.21186       1.000E-05 22 04 35.540  1.500E-01    0.990 F   -0.091    0.327 -10 41 00.00  1.000E-01    0.810 F    0.122    0.062                       o 304      0.35 1 0
 88611     O C   2005 07 11.27812       1.000E-05 22 04 35.330  1.500E-01    0.990 F   -0.091    0.335 -10 41 01.10  1.000E-01    0.810 F    0.122   -0.004                       o 304      0.34 1 0
 88611     O C 5 2009 08 20.19006       1.000E-05 22 20 00.630  1.500E-01    0.990 F   -0.174    0.220 -09 01 25.80  1.000E-01    0.810 F    0.176   -0.017 22.2 R  0.70   0.45   o 688      0.23 1 1
 88611     O C 5 2009 08 20.23589       1.000E-05 22 20 00.400  1.500E-01    0.990 F   -0.174   -0.263 -09 01 26.70  1.000E-01    0.810 F    0.176    0.173                       o 688      0.35 1 0
 88611     O C   2012 10 19.00542       1.000E-05 22 29 39.300  1.500E-01    0.990 F   -0.069   -0.032 -07 56 43.60  1.000E-01    0.810 F    0.247   -0.263 22.1 R  0.70   0.23   o 304      0.35 1 1
 88611     O C   2012 10 19.06667       1.000E-05 22 29 39.150  1.500E-01    0.990 F   -0.069    0.008 -07 56 44.30  1.000E-01    0.810 F    0.247   -0.022                       o 304      0.03 1 0
 88611     O C   2012 10 19.20360       1.000E-05 22 29 38.820  1.500E-01    0.990 F   -0.069    0.144 -07 56 46.60  1.000E-01    0.810 F    0.247   -0.235                       o 304      0.35 1 0
 88611     O C   2012 10 20.00112       1.000E-05 22 29 36.920  1.500E-01    0.990 F   -0.069   -0.059 -07 56 58.50  1.000E-01    0.810 F    0.247   -0.061                       o 304      0.10 1 0
 88611     O C   2012 10 20.11332       1.000E-05 22 29 36.650  1.500E-01    0.990 F   -0.069   -0.032 -07 57 00.20  1.000E-01    0.810 F    0.247   -0.084                       o 304      0.12 1 0
 88611     O C   2012 10 20.20082       1.000E-05 22 29 36.450  1.500E-01    0.990 F   -0.069    0.118 -07 57 01.30  1.000E-01    0.810 F    0.247    0.113                       o 304      0.20 1 0