enod's home: April 2008

Monday, April 14, 2008

Is there a way to count the nubmer of files in many many directories conviently?

Without writing scripts to do so.
Answer awaiting...

The keyboard does not feel good now.

I just got my fan of notebook computer changed. It ran down for a long time. The computer got so hot after several minutes running.

However, the keyboard doesn't feel good now. It makes much louder noise to me. Stupid HP nx6330 keyboard quality!

Sunday, April 13, 2008

Very unstable internet in office.

The internet in office is very unstable. Sometimes, it is very slow, with only tens of bits per second. And the speed can reach as high as several hundred bytes per second; but don't last for long time. The most weird thing is even the file wasn't downloaded completely, FlashFXP stop downloading and think it has been finished, and also removed it from the queue. That's the worst thing for a correct downloading.

Is ftp.cdc.noaa.gov down?

I cannot login it for several days.

File which stores Solaris Calendar information.

/var/spool/calendar/callog.user
It's a text file.

Cannot get access to Google Blog from yesterday.

It presents great difficulties.

Saturday, April 12, 2008

Comparison between rate solution between CMONOC center and GAG@ICD.

SITE N Rate E Rate U Rate
GAG CMONOC GAG CMONOC GAG CMONOC
----|-------------|---------|------------|--------|--------------|------
bjsh -0.01481 -0.01258 0.02700 0.02820 -0.00240 -0.00267
dlha -0.00686 -0.00137 0.02811 0.03352 0.00091 -0.00069
dxin -0.00961 -0.00713 0.02714 0.02917 -0.00461 -0.00114
guan -0.01639 -0.01384 0.03028 0.03138 -0.00634 -0.00559
hrbn -0.00535 -0.01272 0.02512 0.02496 -0.02066 -0.00281
jixn -0.01425 -0.01139 0.02662 0.02756 -0.00327 -0.00076
kmin -0.01952 -0.02142 0.03409 0.03319 -0.00559 -0.00497
luzh -0.01453 -0.01234 0.03135 0.03323 -0.00058 -0.00100
qdao -0.00785 -0.01038 0.03269 0.03079 0.02293 0.00729
qion -0.01425 -0.01200 0.02882 0.03021 0.00060 -0.00076
suiy -0.01657 -0.01404 0.02224 0.02408 -0.00126 0.00049
tain -0.01252 -0.01426 0.03184 0.03078 -0.00387 -0.00315
tash *********** 0.02110 ********* 0.02471 2602.83188 -0.00236
wush 0.01127 0.01315 0.02781 0.02947 -0.00436 -0.00202
xiaa -0.00969 -0.00855 0.02793 0.02810 -0.00571 -0.00709
xiag -0.02227 -0.02021 0.02482 0.02690 -0.00098 -0.00106
xiam -0.01800 -0.01501 0.02967 0.03114 -0.00220 -0.00222
xnin -0.00719 -0.00588 0.02896 0.03621 0.00172 -0.00198
zhnz -0.01292 -0.01271 0.03158 0.02959 0.00438 -0.00501

Plot of XIAA positions time series.

Fig. Coordinates time series of XIAA. (top) North; (middle) East; (bottom) Up.
Note, there are a data gap between 140~365/2005. GAMIT failed to get solutions for those days in the 2nd batch run. I don't know why. If it is run day by day, it seems there is no problems.

Offsets in JPLM time series.

There are obviously several jumps in the coordinates of JPLM time series. This series is obtained by combining SIO igs1, igs2 and igs3 H-files.

ASH701945.02B is missing from the latest igs05.atx file.

The latest IGS absolute phase center corrections file [IGSMAIL #5742] came out. ftp://igscb.jpl.nasa.gov/igscb/station/general/igs05.atx

The gamit cannot get information for antanna (ASH701945.02B ) of XIAM site, and stopped proceeding. However, there are records in the old igs_05.atx files. Did they delete them?

Can I install two versions of IDL on the same Linux machine?

RHEL-5:
1. IDL V7.0 (Already installed);
2. IDL V6.4 (Not installed).
I wander whether there is a way to make those two versions of IDL co-exist in one Linux machine. For different user, I think there is no problem. But how to install those two for the same user?

[IGSMAIL-5757]: IGS Analysis Center Workshop 2008 - Call for abstracts for the Reanalysis Session

******************************************************************************
IGS Electronic Mail 10 Apr 23:08:52 PDT 2008 Message Number 5757
******************************************************************************

Author: Gerd Gendt

Dear colleagues,

we would like to remind you to submit contributions (posters) for the
session "Reanalysis Campaign" at the IGS Analysis Center Workshop, 2-6
June 2008, Miami Beach, Florida, USA (see
http://www.ngs.noaa.gov/IGSWorkshop2008)

The deadline for the submissions is 1 May 2008.

You can find a session description below.
Please send authors and title of your contribution to the session chairs:

G. Gendt: gendt@gfz-potsdam.de
R. Ferland: rferland@nrcan.gc.ca
I. Romero: nacho@nachoromero.com
P. Steigenberger: steigenberger@bv.tum.de

and a copy to Jim Ray: jimr@ngs.noaa.gov

Best regards
Gerd Gendt

Session Description:

During the IGS history several TRFs were applied and the ACs have
implemented many changes for software and technology improvements; all
led to various inconsistencies in the IGS products. A major change was
the introduction of the absolute antenna models, which entered only with
some approximation into the recent IRTF2005. The IGS has to overcome
these inconsistencies in their time series with a reprocessing of its
historical data since 1994 using the latest models and the latest AC's
software implementations. The improved homogeneous products will allow
better geophysical interpretations and will result in a significantly
better contribution to the next realization of the ITRF. The reanalysis
will also include more stations to better address the requirements of
the user community.

The session will concentrate on all activities related to the IGS
reprocessing. Evaluations of the reprocessing test campaign, which shall
identify problems in an early stage, and first results from the IGS
reprocessing, available up to this IGS Workshop, are welcome. Also
welcome are experiences from any reprocessing activities outside the
IGS, which can give valuable feedback to the IGS reprocessing project.

--
**********************************************************************
* Gerd Gendt *
* _/_/_/ _/_/_/_/_/_/ GeoForschungsZentrum *
* _/ _/ _/ Department Geodesy and Remote Sensing *
* _/ _/_/ _/_/ _/ Telegrafenberg A17 *
* _/ _/ _/ _/ D-14473 Potsdam , GERMANY *
* _/_/_/ _/ _/_/_/ Phone : ++49 331 288 1114 *
* Fax :++49 331 288 1759 *
* E-Mail: gendt@gfz-potsdam.de *
* IGS GFZ AC : AFTP://ftp.gfz-potsdam.de/pub/igsprod *
* http://www.gfz-potsdam.de/igs *
* IGS Analysis Coord.: AFTP://ftp.gfz-potsdam.de/pub/igsacc *
* http://www.gfz-potsdam.de/igsacc *
**********************************************************************

Problem with zcat, zgrep, gzcat, gzgrep.

The zcat and zgrep works well in RHEL-5. However, zcat (or gzgrep and gzcat) in Sol10X86 behave differently with SIO h-files and local CMON h-files. Moreover, in RHEL-5, I failed to get information from the compressed h-files of CMON experiment. This cross-platform problem seems so complicated. I can't handle this. :(

There are problems with grep of RHEL-5. It said "Binary file (standard input) matches" when using zcat hbjgpa.05001.Z | grep "Number of stations" or zgrep "Number of stations" hbjgpa.05001.Z.

IERS Bulletin A

**********************************************************************
* *
* I E R S B U L L E T I N - A *
* *
* Rapid Service/Prediction of Earth Orientation *
**********************************************************************
10 April 2008 Vol. XXI No. 015
______________________________________________________________________
GENERAL INFORMATION:
To receive this information electronically, contact:
ser7@maia.usno.navy.mil or use <http://maia.usno.navy.mil/>
MJD = Julian Date - 2 400 000.5 days
UT2-UT1 = 0.022 sin(2*pi*T) - 0.012 cos(2*pi*T)
- 0.006 sin(4*pi*T) + 0.007 cos(4*pi*T)
where pi = 3.14159265... and T is the date in Besselian years.
TT = TAI + 32.184 seconds
DUT1= (UT1-UTC) transmitted with time signals
= -0.4 seconds beginning 13 March 2008 at 0000 UTC
Beginning 1 January 2006:
TAI-UTC(BIPM) = 33.000 000 seconds
________________________________________________________________________
The contributed observations used in the preparation of this Bulletin
are available at <http://maia.usno.navy.mil/bulla-data.html>. The
contributed analysis results are based on data from Very Long Baseline
Interferometry (VLBI), Satellite Laser Ranging (SLR), the Global
Positioning System (GPS) satellites, Lunar Laser Ranging (LLR), and
meteorological predictions of variations in Atmospheric Angular
Momentum (AAM).
________________________________________________________________________

COMBINED EARTH ORIENTATION PARAMETERS:

IERS Rapid Service
MJD x error y error UT1-UTC error
" " " " s s
8 4 4 54560 -.05614 .00009 .49535 .00009 -.365281 .000010
8 4 5 54561 -.05359 .00009 .49704 .00009 -.367124 .000009
8 4 6 54562 -.05088 .00009 .49872 .00009 -.369052 .000009
8 4 7 54563 -.04795 .00009 .50046 .00009 -.370810 .000010
8 4 8 54564 -.04517 .00009 .50243 .00010 -.372215 .000010
8 4 9 54565 -.04292 .00009 .50468 .00010 -.373256 .000048
8 4 10 54566 -.04110 .00009 .50674 .00010 -.374044 .000057

_______________________________________________________________________

PREDICTIONS:
The following formulas will not reproduce the predictions given below,
but may be used to extend the predictions beyond the end of this table.

x = .0718 - .0799 cos A + .0631 sin A - .0266 cos C + .0970 sin C
y = .3543 + .0561 cos A + .0778 sin A + .0970 cos C + .0266 sin C
UT1-UTC = -.3427 - .00080 (MJD - 54574) - (UT2-UT1)

where A = 2*pi*(MJD-54566)/365.25 and C = 2*pi*(MJD-54566)/435.

TAI-UTC(MJD 54567) = 33.0
The accuracy may be estimated from the expressions:
S x,y = 0.00068 (MJD-54566)**0.80 S t = 0.00025 (MJD-54566)**0.75
Estimated accuracies are: Predictions 10 d 20 d 30 d 40 d
Polar coord's 0.004 0.007 0.010 0.013
UT1-UTC 0.0014 0.0024 0.0032 0.0040

MJD x(arcsec) y(arcsec) UT1-UTC(sec)
2008 4 11 54567 -0.0390 0.5086 -0.37469
2008 4 12 54568 -0.0368 0.5102 -0.37532
2008 4 13 54569 -0.0345 0.5117 -0.37606
2008 4 14 54570 -0.0320 0.5132 -0.37698
2008 4 15 54571 -0.0294 0.5147 -0.37808
2008 4 16 54572 -0.0267 0.5161 -0.37934
2008 4 17 54573 -0.0239 0.5175 -0.38069
2008 4 18 54574 -0.0211 0.5189 -0.38205
2008 4 19 54575 -0.0183 0.5202 -0.38331
2008 4 20 54576 -0.0154 0.5215 -0.38441
2008 4 21 54577 -0.0125 0.5228 -0.38532
2008 4 22 54578 -0.0096 0.5240 -0.38604
2008 4 23 54579 -0.0067 0.5252 -0.38658
2008 4 24 54580 -0.0037 0.5263 -0.38702
2008 4 25 54581 -0.0007 0.5274 -0.38743
2008 4 26 54582 0.0023 0.5284 -0.38783
2008 4 27 54583 0.0054 0.5294 -0.38831
2008 4 28 54584 0.0084 0.5303 -0.38894
2008 4 29 54585 0.0115 0.5312 -0.38979
2008 4 30 54586 0.0145 0.5321 -0.39089
2008 5 1 54587 0.0176 0.5329 -0.39224
2008 5 2 54588 0.0208 0.5336 -0.39376
2008 5 3 54589 0.0239 0.5343 -0.39533
2008 5 4 54590 0.0270 0.5350 -0.39679
2008 5 5 54591 0.0302 0.5356 -0.39801
2008 5 6 54592 0.0333 0.5362 -0.39890
2008 5 7 54593 0.0365 0.5367 -0.39949
2008 5 8 54594 0.0397 0.5372 -0.39987
2008 5 9 54595 0.0429 0.5376 -0.40020
2008 5 10 54596 0.0460 0.5380 -0.40060
2008 5 11 54597 0.0492 0.5383 -0.40116
2008 5 12 54598 0.0524 0.5386 -0.40189
2008 5 13 54599 0.0557 0.5388 -0.40276
2008 5 14 54600 0.0589 0.5390 -0.40371
2008 5 15 54601 0.0621 0.5391 -0.40464
2008 5 16 54602 0.0653 0.5392 -0.40550
2008 5 17 54603 0.0685 0.5392 -0.40621
2008 5 18 54604 0.0717 0.5392 -0.40675
2008 5 19 54605 0.0750 0.5391 -0.40710
2008 5 20 54606 0.0782 0.5390 -0.40727
2008 5 21 54607 0.0814 0.5389 -0.40731
2008 5 22 54608 0.0846 0.5387 -0.40727
2008 5 23 54609 0.0878 0.5384 -0.40721
2008 5 24 54610 0.0910 0.5381 -0.40722
2008 5 25 54611 0.0942 0.5377 -0.40736
2008 5 26 54612 0.0974 0.5373 -0.40768
2008 5 27 54613 0.1006 0.5369 -0.40822
2008 5 28 54614 0.1038 0.5364 -0.40899
2008 5 29 54615 0.1070 0.5358 -0.40994
2008 5 30 54616 0.1101 0.5352 -0.41100
2008 5 31 54617 0.1133 0.5346 -0.41206
2008 6 1 54618 0.1164 0.5339 -0.41296
2008 6 2 54619 0.1196 0.5331 -0.41361
2008 6 3 54620 0.1227 0.5324 -0.41397
2008 6 4 54621 0.1258 0.5315 -0.41411
2008 6 5 54622 0.1289 0.5306 -0.41414
2008 6 6 54623 0.1320 0.5297 -0.41421
2008 6 7 54624 0.1350 0.5287 -0.41445
2008 6 8 54625 0.1381 0.5277 -0.41490
2008 6 9 54626 0.1411 0.5266 -0.41553
2008 6 10 54627 0.1441 0.5255 -0.41626
2008 6 11 54628 0.1471 0.5243 -0.41699
2008 6 12 54629 0.1500 0.5231 -0.41764
2008 6 13 54630 0.1530 0.5219 -0.41814
2008 6 14 54631 0.1559 0.5206 -0.41845
2008 6 15 54632 0.1588 0.5192 -0.41856
2008 6 16 54633 0.1617 0.5179 -0.41847
2008 6 17 54634 0.1646 0.5164 -0.41821
2008 6 18 54635 0.1674 0.5149 -0.41785
2008 6 19 54636 0.1702 0.5134 -0.41743
2008 6 20 54637 0.1730 0.5119 -0.41705
2008 6 21 54638 0.1757 0.5103 -0.41677
2008 6 22 54639 0.1785 0.5086 -0.41664
2008 6 23 54640 0.1812 0.5069 -0.41670
2008 6 24 54641 0.1838 0.5052 -0.41696
2008 6 25 54642 0.1865 0.5034 -0.41738
2008 6 26 54643 0.1891 0.5016 -0.41792
2008 6 27 54644 0.1917 0.4998 -0.41847
2008 6 28 54645 0.1942 0.4979 -0.41892
2008 6 29 54646 0.1967 0.4960 -0.41918
2008 6 30 54647 0.1992 0.4940 -0.41919
2008 7 1 54648 0.2017 0.4920 -0.41898
2008 7 2 54649 0.2041 0.4900 -0.41864
2008 7 3 54650 0.2065 0.4879 -0.41830
2008 7 4 54651 0.2088 0.4858 -0.41810
2008 7 5 54652 0.2111 0.4836 -0.41811
2008 7 6 54653 0.2134 0.4814 -0.41835
2008 7 7 54654 0.2156 0.4792 -0.41872
2008 7 8 54655 0.2178 0.4769 -0.41915
2008 7 9 54656 0.2200 0.4747 -0.41950
2008 7 10 54657 0.2221 0.4723 -0.41971
2008 7 11 54658 0.2242 0.4700 -0.41973
2008 7 12 54659 0.2262 0.4676 -0.41954
2008 7 13 54660 0.2282 0.4652 -0.41915
2008 7 14 54661 0.2302 0.4627 -0.41858
2008 7 15 54662 0.2321 0.4603 -0.41790
2008 7 16 54663 0.2340 0.4578 -0.41717
2008 7 17 54664 0.2358 0.4552 -0.41646
2008 7 18 54665 0.2376 0.4527 -0.41586
2008 7 19 54666 0.2394 0.4501 -0.41543
2008 7 20 54667 0.2411 0.4474 -0.41520
2008 7 21 54668 0.2427 0.4448 -0.41519
2008 7 22 54669 0.2443 0.4421 -0.41536
2008 7 23 54670 0.2459 0.4394 -0.41565
2008 7 24 54671 0.2474 0.4367 -0.41595
2008 7 25 54672 0.2489 0.4340 -0.41618
2008 7 26 54673 0.2503 0.4312 -0.41624
2008 7 27 54674 0.2517 0.4284 -0.41609
2008 7 28 54675 0.2531 0.4256 -0.41569
2008 7 29 54676 0.2544 0.4228 -0.41516
2008 7 30 54677 0.2556 0.4199 -0.41460
2008 7 31 54678 0.2568 0.4171 -0.41415
2008 8 1 54679 0.2579 0.4142 -0.41391
2008 8 2 54680 0.2590 0.4113 -0.41393
2008 8 3 54681 0.2601 0.4084 -0.41418
2008 8 4 54682 0.2611 0.4054 -0.41455
2008 8 5 54683 0.2620 0.4025 -0.41492
2008 8 6 54684 0.2629 0.3995 -0.41521
2008 8 7 54685 0.2638 0.3965 -0.41532
2008 8 8 54686 0.2646 0.3936 -0.41525
2008 8 9 54687 0.2653 0.3906 -0.41497
2008 8 10 54688 0.2660 0.3875 -0.41454
2008 8 11 54689 0.2666 0.3845 -0.41398
2008 8 12 54690 0.2672 0.3815 -0.41337
2008 8 13 54691 0.2678 0.3784 -0.41279
2008 8 14 54692 0.2683 0.3754 -0.41231
2008 8 15 54693 0.2687 0.3723 -0.41200
2008 8 16 54694 0.2691 0.3692 -0.41190
2008 8 17 54695 0.2694 0.3662 -0.41206
2008 8 18 54696 0.2697 0.3631 -0.41243
2008 8 19 54697 0.2699 0.3600 -0.41296
2008 8 20 54698 0.2701 0.3569 -0.41356
2008 8 21 54699 0.2702 0.3538 -0.41410
2008 8 22 54700 0.2703 0.3507 -0.41449
2008 8 23 54701 0.2703 0.3476 -0.41467
2008 8 24 54702 0.2703 0.3445 -0.41464
2008 8 25 54703 0.2702 0.3414 -0.41445
2008 8 26 54704 0.2701 0.3383 -0.41421
2008 8 27 54705 0.2699 0.3353 -0.41404
2008 8 28 54706 0.2696 0.3322 -0.41407
2008 8 29 54707 0.2693 0.3291 -0.41436
2008 8 30 54708 0.2690 0.3260 -0.41489
2008 8 31 54709 0.2686 0.3229 -0.41561
2008 9 1 54710 0.2682 0.3198 -0.41641
2008 9 2 54711 0.2676 0.3168 -0.41717
2008 9 3 54712 0.2671 0.3137 -0.41779
2008 9 4 54713 0.2665 0.3107 -0.41824
2008 9 5 54714 0.2658 0.3076 -0.41849
2008 9 6 54715 0.2651 0.3046 -0.41856
2008 9 7 54716 0.2644 0.3016 -0.41850
2008 9 8 54717 0.2636 0.2986 -0.41838
2008 9 9 54718 0.2627 0.2956 -0.41825
2008 9 10 54719 0.2618 0.2926 -0.41820
2008 9 11 54720 0.2608 0.2896 -0.41831
2008 9 12 54721 0.2598 0.2867 -0.41863
2008 9 13 54722 0.2588 0.2837 -0.41919
2008 9 14 54723 0.2576 0.2808 -0.42000
2008 9 15 54724 0.2565 0.2779 -0.42102
2008 9 16 54725 0.2553 0.2750 -0.42215
2008 9 17 54726 0.2540 0.2721 -0.42327
2008 9 18 54727 0.2527 0.2693 -0.42427
2008 9 19 54728 0.2514 0.2665 -0.42505
2008 9 20 54729 0.2500 0.2637 -0.42559
2008 9 21 54730 0.2485 0.2609 -0.42595
2008 9 22 54731 0.2470 0.2581 -0.42623
2008 9 23 54732 0.2455 0.2553 -0.42654
2008 9 24 54733 0.2439 0.2526 -0.42701
2008 9 25 54734 0.2423 0.2499 -0.42771
2008 9 26 54735 0.2406 0.2472 -0.42865
2008 9 27 54736 0.2389 0.2446 -0.42978
2008 9 28 54737 0.2371 0.2420 -0.43104
2008 9 29 54738 0.2353 0.2394 -0.43229
2008 9 30 54739 0.2334 0.2368 -0.43346
2008 10 1 54740 0.2315 0.2343 -0.43446
2008 10 2 54741 0.2296 0.2318 -0.43527
2008 10 3 54742 0.2276 0.2293 -0.43589
2008 10 4 54743 0.2256 0.2268 -0.43635
2008 10 5 54744 0.2235 0.2244 -0.43671
2008 10 6 54745 0.2214 0.2220 -0.43705
2008 10 7 54746 0.2193 0.2197 -0.43743
2008 10 8 54747 0.2171 0.2173 -0.43793
2008 10 9 54748 0.2149 0.2150 -0.43860
2008 10 10 54749 0.2126 0.2128 -0.43950
2008 10 11 54750 0.2103 0.2106 -0.44063
2008 10 12 54751 0.2080 0.2084 -0.44199
2008 10 13 54752 0.2056 0.2062 -0.44349
2008 10 14 54753 0.2032 0.2041 -0.44504
2008 10 15 54754 0.2008 0.2020 -0.44649
2008 10 16 54755 0.1983 0.2000 -0.44774
2008 10 17 54756 0.1958 0.1980 -0.44872
2008 10 18 54757 0.1933 0.1960 -0.44944
2008 10 19 54758 0.1907 0.1941 -0.44998
2008 10 20 54759 0.1881 0.1922 -0.45048
2008 10 21 54760 0.1854 0.1903 -0.45104
2008 10 22 54761 0.1828 0.1885 -0.45176
2008 10 23 54762 0.1801 0.1867 -0.45266
2008 10 24 54763 0.1774 0.1850 -0.45372
2008 10 25 54764 0.1746 0.1833 -0.45487
2008 10 26 54765 0.1718 0.1817 -0.45605
2008 10 27 54766 0.1690 0.1801 -0.45717
2008 10 28 54767 0.1662 0.1785 -0.45817
2008 10 29 54768 0.1633 0.1770 -0.45901
2008 10 30 54769 0.1604 0.1755 -0.45967
2008 10 31 54770 0.1575 0.1741 -0.46018
2008 11 1 54771 0.1546 0.1727 -0.46057
2008 11 2 54772 0.1516 0.1713 -0.46092
2008 11 3 54773 0.1487 0.1700 -0.46131
2008 11 4 54774 0.1457 0.1688 -0.46181
2008 11 5 54775 0.1427 0.1676 -0.46250
2008 11 6 54776 0.1396 0.1664 -0.46340
2008 11 7 54777 0.1366 0.1653 -0.46455
2008 11 8 54778 0.1335 0.1642 -0.46590
2008 11 9 54779 0.1304 0.1632 -0.46743
2008 11 10 54780 0.1273 0.1623 -0.46904
2008 11 11 54781 0.1242 0.1613 -0.47061
2008 11 12 54782 0.1211 0.1604 -0.47204
2008 11 13 54783 0.1179 0.1596 -0.47323
2008 11 14 54784 0.1147 0.1588 -0.47420
2008 11 15 54785 0.1116 0.1581 -0.47498
2008 11 16 54786 0.1084 0.1574 -0.47568
2008 11 17 54787 0.1052 0.1568 -0.47646
2008 11 18 54788 0.1020 0.1562 -0.47742
2008 11 19 54789 0.0988 0.1556 -0.47859
2008 11 20 54790 0.0956 0.1552 -0.47993
2008 11 21 54791 0.0923 0.1547 -0.48137
2008 11 22 54792 0.0891 0.1543 -0.48280
2008 11 23 54793 0.0859 0.1540 -0.48414
2008 11 24 54794 0.0826 0.1537 -0.48528
2008 11 25 54795 0.0794 0.1534 -0.48624
2008 11 26 54796 0.0761 0.1532 -0.48704
2008 11 27 54797 0.0729 0.1531 -0.48765
2008 11 28 54798 0.0696 0.1530 -0.48811
2008 11 29 54799 0.0664 0.1529 -0.48850
2008 11 30 54800 0.0631 0.1529 -0.48888
2008 12 1 54801 0.0599 0.1530 -0.48932
2008 12 2 54802 0.0566 0.1531 -0.48989
2008 12 3 54803 0.0534 0.1533 -0.49063
2008 12 4 54804 0.0501 0.1535 -0.49151
2008 12 5 54805 0.0469 0.1537 -0.49251
2008 12 6 54806 0.0437 0.1540 -0.49356
2008 12 7 54807 0.0405 0.1544 -0.49463
2008 12 8 54808 0.0372 0.1548 -0.49567
2008 12 9 54809 0.0340 0.1552 -0.49659
2008 12 10 54810 0.0308 0.1557 -0.49734
2008 12 11 54811 0.0277 0.1562 -0.49787
2008 12 12 54812 0.0245 0.1568 -0.49823
2008 12 13 54813 0.0213 0.1575 -0.49850
2008 12 14 54814 0.0182 0.1582 -0.49887
2008 12 15 54815 0.0150 0.1589 -0.49944
2008 12 16 54816 0.0119 0.1597 -0.50027
2008 12 17 54817 0.0088 0.1605 -0.50131
2008 12 18 54818 0.0057 0.1614 -0.50246
2008 12 19 54819 0.0026 0.1623 -0.50358
2008 12 20 54820 -0.0005 0.1633 -0.50453
2008 12 21 54821 -0.0035 0.1643 -0.50534
2008 12 22 54822 -0.0065 0.1654 -0.50597
2008 12 23 54823 -0.0095 0.1665 -0.50642
2008 12 24 54824 -0.0125 0.1677 -0.50675
2008 12 25 54825 -0.0155 0.1689 -0.50694
2008 12 26 54826 -0.0184 0.1701 -0.50701
2008 12 27 54827 -0.0214 0.1714 -0.50710
2008 12 28 54828 -0.0242 0.1728 -0.50731
2008 12 29 54829 -0.0271 0.1742 -0.50767
2008 12 30 54830 -0.0300 0.1756 -0.50823
2008 12 31 54831 -0.0328 0.1771 -0.50899
2009 1 1 54832 -0.0356 0.1786 -0.50991
2009 1 2 54833 -0.0384 0.1802 -0.51098
2009 1 3 54834 -0.0411 0.1818 -0.51215
2009 1 4 54835 -0.0438 0.1834 -0.51328
2009 1 5 54836 -0.0465 0.1851 -0.51434
2009 1 6 54837 -0.0491 0.1868 -0.51526
2009 1 7 54838 -0.0518 0.1886 -0.51601
2009 1 8 54839 -0.0544 0.1904 -0.51653
2009 1 9 54840 -0.0569 0.1922 -0.51686
2009 1 10 54841 -0.0595 0.1941 -0.51718
2009 1 11 54842 -0.0619 0.1960 -0.51763
2009 1 12 54843 -0.0644 0.1980 -0.51835
2009 1 13 54844 -0.0668 0.2000 -0.51934
2009 1 14 54845 -0.0692 0.2020 -0.52056
2009 1 15 54846 -0.0716 0.2041 -0.52190
2009 1 16 54847 -0.0739 0.2062 -0.52320
2009 1 17 54848 -0.0762 0.2084 -0.52438
2009 1 18 54849 -0.0784 0.2106 -0.52539
2009 1 19 54850 -0.0806 0.2128 -0.52617
2009 1 20 54851 -0.0828 0.2150 -0.52675
2009 1 21 54852 -0.0849 0.2173 -0.52719
2009 1 22 54853 -0.0870 0.2196 -0.52754
2009 1 23 54854 -0.0891 0.2220 -0.52785
2009 1 24 54855 -0.0911 0.2243 -0.52827
2009 1 25 54856 -0.0930 0.2268 -0.52879
2009 1 26 54857 -0.0950 0.2292 -0.52947
2009 1 27 54858 -0.0969 0.2317 -0.53029
2009 1 28 54859 -0.0987 0.2342 -0.53135
2009 1 29 54860 -0.1005 0.2367 -0.53264
2009 1 30 54861 -0.1022 0.2392 -0.53404
2009 1 31 54862 -0.1040 0.2418 -0.53552
2009 2 1 54863 -0.1056 0.2444 -0.53698
2009 2 2 54864 -0.1072 0.2471 -0.53838
2009 2 3 54865 -0.1088 0.2497 -0.53960
2009 2 4 54866 -0.1103 0.2524 -0.54061
2009 2 5 54867 -0.1118 0.2551 -0.54149
2009 2 6 54868 -0.1133 0.2578 -0.54236
2009 2 7 54869 -0.1146 0.2606 -0.54329
2009 2 8 54870 -0.1160 0.2633 -0.54440
2009 2 9 54871 -0.1173 0.2661 -0.54575
2009 2 10 54872 -0.1185 0.2689 -0.54731
2009 2 11 54873 -0.1197 0.2718 -0.54902
2009 2 12 54874 -0.1209 0.2746 -0.55073
2009 2 13 54875 -0.1220 0.2775 -0.55235
2009 2 14 54876 -0.1230 0.2804 -0.55381
2009 2 15 54877 -0.1240 0.2833 -0.55502
2009 2 16 54878 -0.1250 0.2862 -0.55598
2009 2 17 54879 -0.1259 0.2891 -0.55673
2009 2 18 54880 -0.1267 0.2921 -0.55734
2009 2 19 54881 -0.1275 0.2950 -0.55785
2009 2 20 54882 -0.1283 0.2980 -0.55834
2009 2 21 54883 -0.1290 0.3010 -0.55893
2009 2 22 54884 -0.1297 0.3039 -0.55969
2009 2 23 54885 -0.1303 0.3069 -0.56066
2009 2 24 54886 -0.1308 0.3100 -0.56184
2009 2 25 54887 -0.1313 0.3130 -0.56317
2009 2 26 54888 -0.1318 0.3160 -0.56458
2009 2 27 54889 -0.1322 0.3190 -0.56605
2009 2 28 54890 -0.1325 0.3221 -0.56756
2009 3 1 54891 -0.1328 0.3251 -0.56900
2009 3 2 54892 -0.1330 0.3282 -0.57032
2009 3 3 54893 -0.1332 0.3312 -0.57144
2009 3 4 54894 -0.1334 0.3343 -0.57235
2009 3 5 54895 -0.1335 0.3373 -0.57321
2009 3 6 54896 -0.1335 0.3404 -0.57417
2009 3 7 54897 -0.1335 0.3435 -0.57536
2009 3 8 54898 -0.1334 0.3465 -0.57681
2009 3 9 54899 -0.1333 0.3496 -0.57850
2009 3 10 54900 -0.1331 0.3526 -0.58040
2009 3 11 54901 -0.1329 0.3557 -0.58240
2009 3 12 54902 -0.1326 0.3587 -0.58438
2009 3 13 54903 -0.1323 0.3618 -0.58624
2009 3 14 54904 -0.1319 0.3648 -0.58789
2009 3 15 54905 -0.1315 0.3679 -0.58923
2009 3 16 54906 -0.1310 0.3709 -0.59034
2009 3 17 54907 -0.1305 0.3739 -0.59140
2009 3 18 54908 -0.1299 0.3769 -0.59252
2009 3 19 54909 -0.1293 0.3799 -0.59364
2009 3 20 54910 -0.1286 0.3829 -0.59482
2009 3 21 54911 -0.1279 0.3859 -0.59611
2009 3 22 54912 -0.1271 0.3889 -0.59754
2009 3 23 54913 -0.1263 0.3919 -0.59921
2009 3 24 54914 -0.1254 0.3948 -0.60110
2009 3 25 54915 -0.1245 0.3978 -0.60314
2009 3 26 54916 -0.1235 0.4007 -0.60523
2009 3 27 54917 -0.1225 0.4036 -0.60719
2009 3 28 54918 -0.1214 0.4065 -0.60898
2009 3 29 54919 -0.1203 0.4094 -0.61056
2009 3 30 54920 -0.1191 0.4122 -0.61188
2009 3 31 54921 -0.1179 0.4151 -0.61296
2009 4 1 54922 -0.1166 0.4179 -0.61392
2009 4 2 54923 -0.1153 0.4207 -0.61492
2009 4 3 54924 -0.1140 0.4235 -0.61603
2009 4 4 54925 -0.1126 0.4262 -0.61732
2009 4 5 54926 -0.1111 0.4290 -0.61880
2009 4 6 54927 -0.1096 0.4317 -0.62043
2009 4 7 54928 -0.1081 0.4344 -0.62218
2009 4 8 54929 -0.1065 0.4371 -0.62393
2009 4 9 54930 -0.1049 0.4397 -0.62556
2009 4 10 54931 -0.1032 0.4423 -0.62700
These predictions are based on all announced leap seconds.

CELESTIAL POLE OFFSET SERIES:
NEOS Celestial Pole Offset Series
MJD dpsi error deps error
(msec. of arc)
54546 -57.14 .33 -6.06 .14
54547 -57.12 .45 -5.97 .19
54548 -57.09 .45 -6.08 .19
54549 -57.14 .74 -6.30 .13
54550 -57.28 .74 -6.48 .13
54551 -57.40 .74 -6.54 .34
54552 -57.34 .74 -6.59 .34
54553 -57.09 .74 -6.78 .34

IAU2000A Celestial Pole Offset Series
MJD dX error dY error
(msec. of arc)
54546 0.270 0.131 -0.000 0.140
54547 0.252 0.179 0.011 0.190
54548 0.246 0.179 0.004 0.190
54549 0.248 0.294 -0.013 0.130
54550 0.255 0.294 -0.043 0.130
54551 0.258 0.294 -0.067 0.340
54552 0.270 0.294 -0.085 0.340
54553 0.291 0.294 -0.098 0.340

The IERS Conventions recommended software to predict celestial pole
offsets is available at ftp://maia.usno.navy.mil/conventions/chapter5/
ceppred.f and software for the calculation of the dX and dY with
respect to IAU2000A Nutation/Precession Theory can be found at the
<http://maia.usno.navy.mil/conv2003.html> web site in section chapter 5.

And GAMIT keeps kidding me on GPSDC.

sh_gamit run by each day can obtain solutions for about 32 stations. However, if put it into batch run, it will only got solutions for some 8 sites. I am indeed mad at this!

The head drives me crazy.

The head utility has different behavior on h-files which I generated and SIO supplied. Its -n option works well with SIO h-files, but always returns one more line when dealing with h-files created by local GAMIT run (cmon/bjgp). What the hell is it doing? I don't think the format of those h-files are different. They are all the same, ASCII, unit style line ending. Did I miss someting?

How to get information from GAMIT H-file or SINEX files?

Such information as sites names, coordinates, etc.

Format of mk_pos_ps input map list file.

To plot time series,
mk_pos_ps map.list sites.list output_directory
Actually, mk_pos_ps shell script calls plot_tseries program to do the work.
plot_tseries map.list SITE_GPS 2
This will create a Postscript file SITE.ps in current directory.

Format of map lists file:
[tianyf@tyf ~/tmp]$ cat qoca_map.list
virtual_run.map
That is, just those file names, no total number of files as the first line.

Format of site list file:
[tianyf@tyf ~/tmp]$ cat map.sit
bjfs
wuhn
daej

There seems a error with mk_pos_ps script:
[tianyf@tyf ~/tmp]$ mk_pos_ps qoca_map.list map.sit ps
plot_tseries qoca_map.list bjfs 2
plot time series for site bjfs
minmax: No input data found!
minmax: No input data found!
minmax: No input data found!
minmax: No input data found!
plot_tseries qoca_map.list wuhn 2
plot time series for site wuhn
minmax: No input data found!
minmax: No input data found!
minmax: No input data found!
minmax: No input data found!
plot_tseries qoca_map.list daej 2
plot time series for site daej
minmax: No input data found!
minmax: No input data found!
minmax: No input data found!
minmax: No input data found!

Is this a QOCA bug? The QOCA is 1.33 version (Linux).

Netcdf building problem.

[tianyf@tyf netcdf-3.6.2]$ make
make all-recursive
make[1]: Entering directory `/usr/local/GMT4.2.1/netcdf-3.6.2'
Making all in fortran
make[2]: Entering directory `/usr/local/GMT4.2.1/netcdf-3.6.2/fortran'
make all-am
make[3]: Entering directory `/usr/local/GMT4.2.1/netcdf-3.6.2/fortran'
depbase=`echo fort-attio.lo | sed 's|[^/]*$|.deps/&|;s|\.lo$||'`; \
if /bin/sh ../libtool --tag=CC --mode=compile cc -DHAVE_CONFIG_H -I. -I. -I.. -I. -I../libsrc -g -O2 -MT fort-attio.lo -MD -MP -MF "$depbase.Tpo" -c -o fort-attio.lo fort-attio.c; \
then mv -f "$depbase.Tpo" "$depbase.Plo"; else rm -f "$depbase.Tpo"; exit 1; fi
libtool: compile: cc -DHAVE_CONFIG_H -I. -I. -I.. -I. -I../libsrc -g -O2 -MT fort-attio.lo -MD -MP -MF .deps/fort-attio.Tpo -c fort-attio.c -o fort-attio.o
In file included from ncfortran.h:12,
from fort-attio.c:7:
cfortran.h:138:3: #error "cfortran.h: Can't find your environment among: - MIPS cc and f77 2.0. (e.g. Silicon Graphics, DECstations, ...) - IBM AIX XL C and FORTRAN Compiler/6000 Version 01.01.0000.0000 - VAX VMS CC 3.1 and FORTRAN 5.4. - Alpha VMS DEC C 1.3 and DEC FORTRAN 6.0. - Alpha OSF DEC C and DEC Fortran for OSF/1 AXP Version 1.2 - Apollo DomainOS 10.2 (sys5.3) with f77 10.7 and cc 6.7. - CRAY - NEC SX-4 SUPER-UX - CONVEX - Sun - PowerStation Fortran with Visual C++ - HP9000s300/s700/s800 Latest test with: HP-UX A.08.07 A 9000/730 - LynxOS: cc or gcc with f2c. - VAXUltrix: vcc,cc or gcc with f2c. gcc or cc with f77. - f77 with vcc works; but missing link magic for f77 I/O. - NO fort. None of gcc, cc or vcc generate required names. - f2c : Use #define f2cFortran, or cc -Df2cFortran - NAG f90: Use #define NAGf90Fortran, or cc -DNAGf90Fortran - Absoft UNIX F77: Use #define AbsoftUNIXFortran or cc -DAbsoftUNIXFortran - Absoft Pro Fortran: Use #define AbsoftProFortran - Portland Group Fortran: Use #define pgiFortran - PathScale Fortran: Use #define PATHSCALE_COMPILER"
make[3]: *** [fort-attio.lo] Error 1
make[3]: Leaving directory `/usr/local/GMT4.2.1/netcdf-3.6.2/fortran'
make[2]: *** [all] Error 2
make[2]: Leaving directory `/usr/local/GMT4.2.1/netcdf-3.6.2/fortran'
make[1]: *** [all-recursive] Error 1
make[1]: Leaving directory `/usr/local/GMT4.2.1/netcdf-3.6.2'
make: *** [all] Error 2

The last blank line in qoca.drv file matters.

It is mandatory that there is a blank line at the end of the QOCA driver file. Otherwise, the qoca will stop at an error of fmt string. For example:


tianyf@tyf tmp]$ qoca q1.drv

--------------------------------------------------------------------------------
                 QOCA version:      1.32 12/28/2007
--------------------------------------------------------------------------------

  drive file name: q1.drv
  coordinate file name: itrf00.net
  input list file name: mrgqob/snxes.qob
  output file name: virtual_run.out
  mapping file name: virtual_run.map
  UTPM table file: pmu.bull_a
  leap-second table file: leap.sec
  screening file: tmp.screen
  reference frame: WGS84
  rtime,comode:  1986.465 2
  solution combination mode: 2
  filtering mode: forward
  network adjust parameter: xyzuvws 1
fmt: end of file
apparent state: unit 8 named q1.drv
last format: (a)
lately reading sequential formatted external IO
Aborted (core dumped)

This thing bothered me much. Ha..

Friday, April 11, 2008

QOCA: cov_eigen error.


[tianyf@tyf tmp]$ cov_eigen mrgqob/snxes.qob coveigen.out 1
    Begin to check <<< qobs/sioigs13650.qob
  trouble in dpttrf of dpteqr: 456 456  0.
 problem in dpteqr
  the  456-th term is not positive defined
    Begin to check <<< qobs/sioigs13651.qob
  trouble in dpttrf of dpteqr: 474 474  0.
 problem in dpteqr
  the  474-th term is not positive defined
    Begin to check <<< qobs/sioigs13652.qob
  trouble in dpttrf of dpteqr: 480 480  0.
 problem in dpteqr
  the  480-th term is not positive defined
    Begin to check <<< qobs/sioigs13653.qob
  trouble in dpttrf of dpteqr: 483 483  0.
 problem in dpteqr
  the  483-th term is not positive defined
    Begin to check <<< qobs/sioigs13654.qob
  trouble in dpttrf of dpteqr: 465 465  0.
 problem in dpteqr
  the  465-th term is not positive defined
    Begin to check <<< qobs/sioigs13655.qob
  trouble in dpttrf of dpteqr: 462 462  0.
 problem in dpteqr
  the  462-th term is not positive defined
    Begin to check <<< qobs/sioigs13656.qob
  trouble in dpttrf of dpteqr: 477 477  0.
 problem in dpteqr
  the  477-th term is not positive defined

The output file:


              INF  4.254973E+02  0.000000E+00  453  456  456   qobs/sioigs13650.qob
   adding term:  0.1000E-05
     4.694180E+08  4.694180E+02  1.000000E-06    1  456  456   qobs/sioigs13650.qob
              INF  4.146400E+02  0.000000E+00  471  474  474   qobs/sioigs13651.qob
   adding term:  0.1000E-05
     4.531920E+08  4.531920E+02  1.000000E-06    1  474  474   qobs/sioigs13651.qob
              INF  4.084898E+02  0.000000E+00  477  480  480   qobs/sioigs13652.qob
   adding term:  0.1000E-05
     4.496876E+08  4.496876E+02  1.000000E-06    1  480  480   qobs/sioigs13652.qob
              INF  4.067726E+02  0.000000E+00  480  483  483   qobs/sioigs13653.qob
   adding term:  0.1000E-05
     4.452042E+08  4.452042E+02  1.000000E-06    1  483  483   qobs/sioigs13653.qob
              INF  4.198978E+02  0.000000E+00  462  465  465   qobs/sioigs13654.qob
   adding term:  0.1000E-05
     4.699089E+08  4.699089E+02  1.000000E-06    1  465  465   qobs/sioigs13654.qob
              INF  4.236193E+02  0.000000E+00  459  462  462   qobs/sioigs13655.qob
   adding term:  0.1000E-05
     4.703942E+08  4.703942E+02  1.000000E-06    1  462  462   qobs/sioigs13655.qob
              INF  4.191456E+02  0.000000E+00  474  477  477   qobs/sioigs13656.qob
   adding term:  0.1000E-05
     4.567178E+08  4.567178E+02  1.000000E-06    1  477  477   qobs/sioigs13656.qob

QOCA failed to convert GAMIT H-files to QOB format.


[tianyf@tyf tmp]$ mk_qob mk_qob.drv
  --------------------------------------------------
               MK_QOB version: 1.06 beta
  --------------------------------------------------
  input file list: mk_qob.list
  input data file style: gamith 6
  unused site:  0
    Processing  <<<  hfiles/higs3a.08007
  Begining to extract site adjusts from h-files ....
read start: end of file
apparent state: unit 21 named tmp_fil2
last format: (a)
lately reading sequential formatted external IO
Aborted (core dumped)

The same problem with either .Z or uncompressed files.

Install CATS under Cygwin.

It's a great challenge to install CATS under Cygwin@Windows. Actually, it is very simple. Just follow the below general steps.

a). Install Cygwin with GCC compiler;
b). Compile CLAPACK with gcc;
c). Compile CATS with gcc.

However, there are some tricks to make it work.
For step a), please refer to those instructions on the internet.
Step b).
1. Download the latest clapack.tgz from netbig site. Unpack it.
tar zxvf clapack.tgz
This will create a directory named CLPACK-3.1.1.
2. Change to the CLAPACK-3.1.1 directory, type make to compile, and get a cup of coffee.
cd CLAPACK-3.1.1
make
This will create several static libraries:
CLAPACK-3.1.1/blas_LINUX.a
CLAPACK-3.1.1/F2CLIBS/libf2c.a
CLAPACK-3.1.1/lapack_LINUX.a
CLAPACK-3.1.1/tmglib_LINUX.a
It's done.

Step c). Compile CATS
1. Unpack CATS distribution.
tar zxvf cats_v3.1.2.tar.gz
This will create a cats_v3.1.2 directory.
2. Compile CATS libraries.
First, modify the make.inc files to suit your needs.
make.inc:
####################################################################
# CATS make include file. #
# CATS, Version 2.0 #
# August 16, 2002 #
####################################################################

BINDIR = /home/tianyf/gpsf/cgps/bin
CLAPACK = /home/tianyf/CLAPACK-3.1.1

LAPACKLIB = /home/tianyf/CLAPACK-3.1.1/lapack_LINUX.a
BLASLIB = /home/tianyf/CLAPACK-3.1.1/blas_LINUX.a
F77LIB = /home/tianyf/CLAPACK-3.1.1/F2CLIBS/libf2c.a
#-lF77_LINUX
I77LIB = /home/tianyf/CLAPACK-3.1.1/tmglib_LINUX.a
#-lI77_LINUX

COMPILER = gcc
AR = ar
INSTALL = install
#CFLAGS = -static -g3 -pg -O -Q
CFLAGS = -static -O
#CFLAGS = -O

Second, cd lib/time sub-directories and type make.
It is recommended that first clean the existing files by "make cleaner".

Finally, change to analysis directory, type make to compile the cats program.

That's all.

Thursday, April 10, 2008

Delete outliers when use multibase.

I finally got an idea why the three components time series from PANGA are in different length. Those guys at PANGA group may use a tool behaving like multibase to eliminate outliers. To my surprise, multibase only delete outliers in separate component time series, leaving it difficult to concatenate the N/E/U series into an integral one--the number of epochs of each component for one site may differ.

multibase example:
multibase val -y -max_signe .05 -max_sigu .1

Do you know how ensum handle command-line parameters?

After dig into its source code (~/gg/kf/blsum/ensum.f), I got that this is a trick related to system shell. For a command like "ensum 2 sum val *.org", when call the program ensum, shell (say, csh) just pass all the file names matching "*.org" into the ensum as command-line parameters. Thus, the final calling syntax may look like:
ensum 2 sum val 1.org 2.org 3.org ...

Then, how to know how many files are there?
Just check count the number of non-blank parameters. (For example, it could be supposed that there maybe maximum 10000 *.org files.)

Script to download brdc/nav files from SOPAC data center.

Examples:
sh_wget_nav -yrs 2005 -yre 2006 -doys 360 -doye 3

sh_wget_nav:

#!/bin/sh

local_root=/cygdrive/e/igs/pub/nav
yearStart=2005
yearStart=1999
#yearStart=2001
yearEnd=2007

doys=001
doye=`ymd2ydoy ${yearStart} 12 31 | awk '{print $2}'`


# process command-line parameters
#####################################################
#COMMAND LINE PARAMETERS
#####################################################
while [ "$1" != "" ]
do
 #echo $1
 case $1 in
     -dir)
   cd $2
   local_root=$2
   ;;
     -yrs)
   yearStart=$2
   yr=`echo $yearStart | awk '{print substr($0,3,2)}'`
   y=`echo $yearStart | awk '{print substr($0,4,1)}'`
   ndays=`ymd2ydoy ${yearStart} 12 31 | awk '{print $2}'`
   ;;
     -yre)
   yearEnd=$2
   ;;
     -doys)
   doys=$2
   ;;
     -doye)
   doye=$2
   ;;
     -sites)
   sites=$2
   sites=`echo $sites | awk -F, '{for (i=1;i<=NF;i++) print $i}'`
   ;;
     -sitefile)
   sfile=$2
  
   sites=`sh_site_read $sfile`
   #echo $sites
   ;;
     *)
   echo "invalid options: $1"
   exit 1
   ;;
 esac
 shift 2
done
#echo $sites


#echo $sites

#exit

year=${yearStart}
while [ ${year} -le ${yearEnd} ]; do

   if [ ${year} -eq ${yearStart} ]; then
 doy=${doys}
   else
      doy=001
   fi

   if [ ${year} -eq ${yearEnd} ]; then
 ndays=${doye}
   else
 ndays=`ymd2ydoy ${year} 12 31 | awk '{print $2}'`
   fi

   yy=`echo $year | awk '{print substr($0,3,2)}'`
   while [ ${doy} -le ${ndays} ]; do
       doy=`echo $doy | awk '{printf("%03d",$0)'}`
       day_dir=${local_root}/${year}/${doy}
       mkdir -p ${day_dir}
       cd ${day_dir}
       pwd

 
 echo "wget -c --user=anonymous --password=tianyf@gmail.com http://garner.ucsd.edu/pub/rinex/${year}/${doy}/auto${doy}0.${yy}n.Z -O brdc${doy}0.${yy}n.Z"
       #from HTTP Server
 wget -c --user=anonymous --password=tianyf@gmail.com http://garner.ucsd.edu/pub/rinex/${year}/${doy}/auto${doy}0.${yy}n.Z -O brdc${doy}0.${yy}n.Z
 #from FTP Server
 #wget -c ftp://garner.ucsd.edu/pub/rinex/${year}/${doy}/auto${doy}0.${yy}n.Z -O brdc${doy}0.${yy}n.Z
 doy=`expr $doy + 1`
   done
   year=`expr $year + 1`
done

Sad at CMONOC Rerun V2 result.

Platform:

GPSDC - Dell Precision 390, 1G RAM, Solaris 10 X86, GAMIT/GLOBK 10.32 (????) with Sun Studio 12;
GPSAC - Dell Optiplex 755, 4G RAM, RHEL-5_X64, GAMIT/GLOBK 10.32 (????) with Intel Fortran/C Compiler 10.

*Note: The version of GAMIT/GLOBK is different in two systems. The parameters may vary in different sessions.

The data are split into four sessions.

2000-2001 (gpsac);
2002-2003 (gpsac);
2004-2005 (gpsdc);
2006-2007 (gpsdc).

It took those tow computers one week to complete the solutions.

The resulting time series (combined with SIO igs1/igs2/igs3) have large offsets at the year boundaries of 2002.0, 2004.0, and 2006.0 epochs, especially for North component.

Why those jumps?

Strange uncompress problem.

The rerun script failed to uncompress the igs14523.sp3.Z file. The error is:
igs14523.sp3.Z is not in compressed format.
The same error occurred with file downloaded from CDDIS server. After unzip/zip it in Cygwin, the Solaris 10 X86 uncompress utility worked.
And there is no problem with SOL10X86 in VMWare.

How to combine files column by column?

It's simple--paste.
e.g.:
paste mb_WHJF_GPS.dat1 mb_WHJF_GPS.dat12 mb_WHJF_GPS.dat13 > whjf.neu
Extra work is needed to remove superfluous lines at the beginning of the file, and delete duplicate column. Also, you need add YEAR DOY column to conform to the SIO/NEU convenction.

Wednesday, April 09, 2008

How to extract certain kinds of files from tared file?

First, here is a example:

$ tar -xf foo.tar --wildcards '*.txt' --no-wildcards '[remarks]'

instructs tar to extract from ‘foo.tar’ all files whose names end in ‘.txt’ and the file named ‘[remarks]’.

Another example:
> gtar zxvf gsoln.v2.tar.gz --wildcards "*.org"
will extract all glorg output files *.org to gsoln directory.

Power outage in the office.

There was a sudden power outage in the office this afternoon. The GPS-DC and GPS-AC servers boot normally after the power was resumed.
GPS-DC: Solaris 10 X86_64
GPS-AC: RHEL5_X86_64
I think they both are great OS.

Windows XP crashed two time recently.

Both cases are related to high load of CPU when VMWare workstation Guest OS was running. Is this a bug of VMware? Or the result of a fatal flaw of XP itself? Or the hardware temperature is too high?

Tuesday, April 08, 2008

Defnition of Besselian Year

Besselian year.
A quantity introduced by F. W. Bessel in the nineteenth century that has been used into the twentieth century. Bessel introduced a system whereby it would be convenient to identify any instant of time by giving the year and the decimal fraction of the year to a few places, but the starting time of the year was not convenient for dynamical studies that utilize Julian dates (see definition for Julian date), differing by 0.5 day, and the Besselian year varies slowly. The recent change to Julian year usage in dynamical astronomy (and the J2000.0 equinox) took effect in solar-system ephemerides of the Minor Planet Center and Central Bureau for Astronomical Telegrams on Jan. 1, 1992. (See Julian year.)

Julian date (JD).
The interval of time in days (and fraction of a day) since Greenwich noon on Jan. 1, 4713 BC. The JD is always half a day off from Universal Time, because the current definition of JD was introduced when the astronomical day was defined to start at noon (prior to 1925) instead of midnight. Thus, 1995 Oct. 10.0 UT = JD 2450000.5.

Julian year.
Exactly 365.25 days, in which a century (100 years) is exactly 36525 days and in which 1900.0 corresponds exactly to 1900 January 0.5 (from the Julian-date system, which is half a day different from civil time or UT). The standard epoch J2000.0, now used for new star-position catalogues and in solar-system-orbital calculations, means 2000 Jan. 1.5 Barycentric Dynamical Time (TDB) = Julian Date 2451545.0 TDB. When this dynamical, artificial "Julian year" is employed, a letter "J" prefixes the year.

Intranet IP at ICD

IP: 10.4.21.23
Gateway: 10.4.21.254
DNS: 202.106.0.20

Binded to the MAC address of ethernet card.

Monday, April 07, 2008

autcln.out SO large!

It took up all left space of the hard disk.
56G 016/2005 CMONOC@GPSDC
46G 136/2006 BJGPS@GPSAC

What happened?

NOT only project, but ALSO dissertation.

Otherwise, I think I might failed to obtain the degree forever.

Summary paper submitted finally.

O.M.G! It took me so long time to send it out. Waiting for responses. Wish good news.

Why Suse Linux makes it looks like Windows?

And also, set Windows style theme as the default. Is Microsoft Windows the best one? Can't understand! Be just yourself, man.

Sunday, April 06, 2008

Setup iPod Touch Win32 Cygwin development environment.

Follow the directions in http://wiki.iphonegameover.com/Windows_Cygwin_Binary_Toolchain_Installation.

To avoid downloading the huge XCode 2.5 image, you could get it directly from rapidshare: http://rapidshare.com/files/66764258/MacOSX10.4u.sdk.rar. Make sure that your iTouch has SSHD installed.

It goes easily if you follow the above directions.

HellWorld test program runs well. Haha... So did the classic helloWorld.c program.

When first start cygwin terminal, it will create a archive of iTouch filesystem:
/private/var/iphone-root-fs.tar
Then download and extract it in cygwin.

Wednesday, April 02, 2008

GSAC Help Documents


     USAGE: gsac-client --f[unction]  {options}
     ABOUT: This is a the front-end (client) interface for the GSAC.
            It essentially provides a GSAC user with a "basic"
            means of acquiring either information or actual data files
            from a GSAC Retailer (specified via -r{etailer} )
            matching the user's criteria.
 FUNCTIONS: status {acquire basic information about given retailer}
            count  {display number of GSAC data records matching request}
            info   {retrieve and display GSAC data records matching request}
            sites  {retrieve and display GSAC monument records matching request}
            data   {retrieve actual data files matching request}
            help   {display this usage}
            version {version of this executable}
   OPTIONS:
     --Operational Configuration Options--
            --v[erbose] {Invoke verbose mode; Synonymous flags: --v or -v}
            --q[uiet] {No stderr output whatsoever; Synonymous flags: --q or -q}
            --debug {Turn on detailed debugging}
            --r[etailer] 
               {Specify which GSAC Retailer Server to use.  Can also be set
                in $GSAC_RETAILER environment variable, where multiple
                retailer hosts are separated by colons;  Defaults
                to :
                    #1 choice :  gsac.ucsd.edu
                    #2 choice :  gsac01.unavco.org
                Synonymous flags: --r or -r}
            --e[mail] 
               {Used for compatibility/identification purposes at most
                public ftp servers; No default - attempts to collect
                data from some ftp servers may fail without this option;
                Can also be set in $GSAC_USER_EMAIL_ADDRESS
                environment variable; Synonymous flags: --e or -e}
            --un[compress] {Valid only with the "data" function;  Indicates an attempt
                will be made to uncompress fetched files before writing to filesystem;
                Assumes "gzip" utility is available and in users $PATH environment
                variable; Synonymous flags: --un or -un}
            --unhat {Valid only with the "data" function; Indicates an attempt
                will be made to remove "Hatanaka" rinex compression from rinex observation
                files (whereever applicable) before writing to local filesystem;
                Assumes "crx2rnx" utility is available and in users $PATH environment
                variable; Synonymouc flags: --uh or -uh}
            --du[plicates] {Valid only with the "data" function; Indicates duplicate
                files from different wholesalers should be collected and written with
                wholesaler prefixes; Synonymous flags: --du or -du}
            --delim[iter]  {Specify the character you wish to have output
                  fields delimited by;  Default is "," (comma)}
            --l[og] {Log transactions to file specified by $GSAC_LOG environment variable;
                  No default; Synonymous flags: --l or -l}
            --offline-log-file 
                  {Path to file containing references to offline data;
                  No default; Synonymous flags: --ol or -ol}
            --output-directory 
                  {Directory to write fetched files; Can also be set using
                  $GSAC_OUTPUT_DIRECTORY environment variable;
                  {Default: .; Synonymous flags: --od or -od}
            --max-files <# files>
                  {Maximum number of files to fetch; Default: 10;
                  Synonymous flags: --mf or -mf}
            --max-mb 
                  {Maximum cumulative size of files to retrieve, in megabytes;
                  Default: 10; Synonymous flags: --mm or -mm}
            --max-parallel-hosts <# hosts>
                  {Max number of hosts to collect files from in parallel;
                  Default: 10; Synonymous flags: --mph or -mph}
            --max-parallel-requests <# requests>
                  {Maximum number of parallel file requests per host; Default: 3;
                  Synonymous flags: --mpr or -mpr}

     --Query Constraint Options--
            --d[ay] <3-digit day> {Day of data to consider; No default;
                  Synonymous flags: --d or -d}
            --y[ear] <4-digit year> {Year of data to consider; No default;
                  Synonymous flags: --y or -y}
            --start-year <4-digit startyear> {Earliest year of data to consider;
                  No default; Synonymous flags: --sy or -sy}
            --start-day <3-digit startday> {Earliest day of data to consider;
                  No default; Synonymous flags: --sd or -sd}
            --end-year <4-digit endyear> {Most recent year of data to consider;
                  No default; Synonymous flags: --ey or -ey}
            --end-day <3-digit endday> {Most recent day of data to consider;
                  No default; Synonymous flags: --ed or -ed}
            --data-type  {Data type to consider; Default is all data types;
                  Can be used more than once; Synonymous flags: --t or -t}
            --exclude-data-type  {Data type to explicitly ignore;
                  No default; Can be used more than once; Synonymous flags: --xt or -xt}
            --w[holesaler]  {GSAC wholesaler to consider; Default
                  is all wholesalers; Can be used more than once;
                  Synonymous flags: --w or -w}
            --exclude-wholesaler  {GSAC wholesaler to explicitly ignore;
                  No default; Synonymous flags: --xw or -xw}
            --site-id  {Unique site identifier to consider;
                  Default is all sites/monuments; Can be used more than once;
                  Synonymous flags: --s or -s}
            --exclude-site-id  {Unique site identifier to
                  explicitly ignore; No default; Can be used more than once;
                  Synonymous flags: --xs or -xs}
            --char-id  {Site character identifier to consider;
                  Default is all sites/monuments; Can be used more than once;
                  Synonymous flags: --c or -c}
            --exclude-char-id  {Site character identifier to
                  explicitly ignore; No default; Can be used more than once;
                  Synonymous flags: --xc or -xc}
            --b[ounding-box]  {Bounding box in longitude latitude
                  decimal degrees; No default; Synonymous flags: --b or -b or --bb or -bb}
            --pr[oximity]  {Longitude,latitude in decimal degrees plus radius
                  in kilometers; No default; Synonymous flags: --pr or -pr}
DATA TYPES: rinex_obs, rinex_met, rinex_nav, raw_gps, igs_site_log, orbit_sp3, sinex

GSAC cannot handle the last day of leap years!

> gsac-client -y 2000 -d 366 -s bjfs0000 -f info
Day option argument (-d) formatted incorrectly...ignoring
When specifying a year and day (via -y and -d) you need both parameteters (see -
-function help)

> doy 2000 366
Date 2000/12/31 0:00 hrs, DOY 366 JD 2451909.5000 MJD 51909.0000
GPS Week 1095 Day of week 0, GPS Seconds 0 Day of Week Sun
Decimal Year 2000.997268
>

It is confirmed that gsac cannot handle the last day of leap years, whether in Windows or Linux/Unix environment.

If it is true, then this is a serious bug. How could this happen?

I can browse GBlog in office.

It's a good news. The internet network in office is much faster sometimes than that in my home.

Tuesday, April 01, 2008

Why the GSAC error?

tian@jiang ~
$ gsac-client -q -f info -d 366 -y 2004 -t rinex_obs -s auck0000 -s ban201 -s bili0000 -s bjfs0000 -s brmu0000 -s cagl0000 -s cedu0000 -s chat0000 -s chur0000 -s coco0000 -s daej0000 -s dubo0000 -s elat0000 -s fair0000 -s glps0000 -s gode0000 -s gol20000 -s graz0000 -s guam0000 -s guao01 -s hers0000 -s hnlc0000 -s hob20000 -s hyde0000 -s iisc0000 -s invk0000 -s irkt0000 -s ispa01 -s kely0000 -s kunm0000 -s lpgs0000 -s mac10000 -s mad20000 -s madr0000 -s mali0000 -s mas10000 -s mate0000 -s mcm40000 -s mdo10000 -s mdvj01 -s mkea0000 -s mobn0000 -s nklg0000 -s nlib0000 -s nrc10000 -s ntus0000 -s nvsk0000 -s nya10000 -s nyal0000 -s onsa0000 -s pert0000 -s petp0000 -s pets01 -s pie10000 -s pimo0000 -s pol20000 -s pots0000 -s qaq10000 -s reyk0000 -s sant0000 -s sele0000 -s sey10000 -s shao0000 -s ssia0000 -s stjo0000 -s suwn0000 -s syog0000 -s thti0000 -s tidb0000 -s tixi0000 -s tnml01 -s tow20000 -s tro10000 -s trom0000 -s ulab0000 -s urum0000 -s usno0000 -s vill0000 -s wes20000 -s whit0000 -s wsrt0000 -s wtzr0000 -s wuhn0000 -s yakt0000 -s yell0000 -s yssk0000 -s zwe201
When specifying a year and day (via -y and -d) you need both parameteters (see --help)

tian@jiang ~
$ gsac-client -q -f info -d 365 -y 2004 -t rinex_obs -s ban201 -s brmu0000 -s cagl0000 -s chur0000 -s elat0000 -s fair0000 -s gol20000 -s guam0000 -s guao01 -s hnlc0000 -s hyde0000 -s iisc0000 -s invk0000 -s irkt0000 -s ispa01 -s kely0000 -s kunm0000 -s lpgs0000 -s mac10000 -s mad20000 -s madr0000 -s mali0000 -s mas10000 -s mate0000 -s mcm40000 -s mdo10000 -s mdvj01 -s mkea0000 -s mobn0000 -s nklg0000 -s nlib0000 -s nrc10000 -s ntus0000 -s nvsk0000 -s nya10000 -s nyal0000 -s onsa0000 -s pert0000 -s petp0000 -s pets01 -s pie10000 -s pimo0000 -s pol20000 -s pots0000 -s qaq10000 -s reyk0000 -s sant0000 -s sele0000 -s sey10000 -s shao0000 -s ssia0000 -s stjo0000 -s suwn0000 -s syog0000 -s thti0000 -s tidb0000 -s tixi0000 -s tnml01 -s tow20000 -s ulab0000 -s urum0000 -s usno0000 -s vill0000 -s wes20000 -s whit0000 -s wsrt0000 -s wuhn0000 -s yakt0000 -s yell0000 -s yssk0000 -s zwe201
ftp://garner.ucsd.edu/pub/rinex/2004/365/syog3650.04d.Z,293877,rinex_obs,5073665,sopac
ftp://garner.ucsd.edu/pub/rinex/2004/365/nvsk3650.04d.Z,295887,rinex_obs,5096030,sopac
ftp://garner.ucsd.edu/pub/rinex/2004/365/ulab3650.04d.Z,441577,rinex_obs,5075728,sopac
ftp://garner.ucsd.edu/pub/rinex/2004/365/pots3650.04d.Z,415042,rinex_obs,5075731,sopac
ftp://garner.ucsd.edu/pub/rinex/2004/365/pets3650.04d.Z,247809,rinex_obs,5071845,sopac
ftp://garner.ucsd.edu/pub/rinex/2004/365/stjo3650.04d.Z,301586,rinex_obs,5071846,sopac
ftp://garner.ucsd.edu/pub/rinex/2004/365/sant3650.04d.Z,316549,rinex_obs,5071848,sopac
ftp://garner.ucsd.edu/pub/rinex/2004/365/guam3650.04d.Z,315443,rinex_obs,5085008,sopac
ftp://garner.ucsd.edu/pub/rinex/2004/365/hyde3650.04d.Z,397493,rinex_obs,5072896,sopac
ftp://garner.ucsd.edu/pub/rinex/2004/365/mobn3650.04d.Z,200569,rinex_obs,5071722,sopac
ftp://garner.ucsd.edu/pub/rinex/2004/365/kely3650.04d.Z,364103,rinex_obs,5071728,sopac
ftp://garner.ucsd.edu/pub/rinex/2004/365/mad23650.04d.Z,295741,rinex_obs,5071729,sopac
ftp://garner.ucsd.edu/pub/rinex/2004/365/whit3650.04d.Z,369047,rinex_obs,5071730,sopac
ftp://garner.ucsd.edu/pub/rinex/2004/365/pol23650.04d.Z,293203,rinex_obs,5071732,sopac
ftp://garner.ucsd.edu/pub/rinex/2004/365/madr3650.04d.Z,297455,rinex_obs,5071733,sopac
ftp://garner.ucsd.edu/pub/rinex/2004/365/nrc13650.04d.Z,303177,rinex_obs,5071734,sopac
ftp://garner.ucsd.edu/pub/rinex/2004/365/tidb3650.04d.Z,290621,rinex_obs,5071735,sopac
ftp://garner.ucsd.edu/pub/rinex/2004/365/petp3650.04d.Z,259289,rinex_obs,5072190,sopac
...