tag:blogger.com,1999:blog-1020595799365545552012-02-21T14:15:42.040+08:00nearly almost everything about my lifeEnodhttp://www.blogger.com/profile/08433335265518211485noreply@blogger.comBlogger2721100tag:blogger.com,1999:blog-102059579936554555.post-60355696578665057022010-04-11T21:27:00.002+08:002010-04-11T21:29:44.292+08:00I can't stand any more today. Then I bought one VPN account. It's to connect the outside world now. How are you, my old friends?<br />+ Google Group<br />+ www.blogger.com<br />+ ...<div class="blogger-post-footer"><img width='1' height='1' src='https://blogger.googleusercontent.com/tracker/102059579936554555-6035569657866505702?l=yftian.blogspot.com' alt='' /></div>Enodhttp://www.blogger.com/profile/08433335265518211485noreply@blogger.com4tag:blogger.com,1999:blog-102059579936554555.post-19220957723162325772009-05-13T16:05:00.002+08:002009-05-13T16:10:58.528+08:00szip 2.1<br />hdf 4.2r4<br /><br />OS: RHEL5_x64/PCSol10_x64<br /><br />To compile szlib 2.1: <br /> ./configure --prefix=/usr/local/szlib CC=icc F77=ifort CXX=icpc<br /> make<br /> make install<br /><br />To compile HDF4.2r4:<br /> ./configure --prefix=/usr/local/hdf4 --with-szlib=/usr/local/szlib CC=icc CXX=icpc F77=ifort<br /><br />The error message is:<br /><span style="font-style:italic;">checking szlib.h usability... yes<br />checking szlib.h presence... yes<br />checking for szlib.h... yes<br />checking for SZ_BufftoBuffCompress in -lsz... yes<br />checking for szlib encoder... checking for xdr library support... yes<br />checking for ANSI C header files... (cached) yes<br />checking whether cc understands -c and -o together... yes<br />checking how to run the C++ preprocessor... CC -E<br />checking whether f77 understands -c and -o together... yes<br />checking for int*... yes<br />checking size of int*... configure: error: cannot compute sizeof (int*)<br />See `config.log' for more details.</span><br /><br />I cannot figure out what happened.<div class="blogger-post-footer"><img width='1' height='1' src='https://blogger.googleusercontent.com/tracker/102059579936554555-1922095772316232577?l=yftian.blogspot.com' alt='' /></div>Enodhttp://www.blogger.com/profile/08433335265518211485noreply@blogger.com24tag:blogger.com,1999:blog-102059579936554555.post-52121110010141477652009-04-22T15:52:00.000+08:002009-04-22T15:53:04.103+08:00<pre><br /> ALBH AMC2 CAGL CHUR DRAO GODE HRAO JPLM MAD2 NRC1 PIE1 REYK SCH2 STJO SUTH TSKB WES2 WHIT WTZR #<br />1999.16300 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 0 0 13<br />1999.16580 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 1 0 14<br />1999.16850 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 0 0 13<br />1999.17120 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 0 0 13<br />1999.17400 1 1 0 1 1 1 1 1 0 0 1 0 1 1 1 0 1 1 0 13<br />1999.17670 1 1 0 1 1 1 1 1 0 0 1 1 1 1 1 0 1 1 0 14<br />1999.17950 1 1 0 1 1 1 1 1 0 0 1 0 1 1 1 0 1 1 0 13<br />1999.18220 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 0 0 13<br />1999.18490 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 1 0 14<br />1999.18770 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 1 0 14<br />1999.19040 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 1 0 14<br />1999.19320 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 1 0 14<br />1999.19590 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 0 0 13<br />1999.19860 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 1 0 14<br />1999.20140 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 0 0 13<br />1999.20410 1 1 0 1 1 1 1 1 0 1 1 1 1 1 1 0 1 1 0 15<br />1999.20680 1 1 0 1 1 1 1 1 0 1 1 1 1 1 1 0 1 1 0 15<br />1999.20960 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 1 0 14<br />1999.21230 1 1 0 1 1 1 1 1 0 1 1 1 1 1 1 0 1 1 0 15<br />1999.21510 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 1 0 14<br />1999.21780 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 1 0 14<br />1999.22050 1 1 0 1 1 1 1 1 0 0 1 0 1 1 1 0 1 1 0 13<br />1999.22330 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 1 0 14<br />1999.22600 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 1 0 14<br />1999.22880 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 1 0 14<br />1999.23150 1 1 0 1 1 1 1 0 0 1 1 0 1 1 1 0 1 1 0 13<br />1999.23420 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 1 0 14<br />1999.23700 1 1 0 1 1 1 1 1 0 1 1 1 1 1 1 0 1 1 0 15<br />1999.23970 1 1 0 1 1 1 1 1 0 1 1 1 0 1 1 0 1 1 0 14<br />1999.24250 0 1 0 1 1 1 1 1 0 1 1 0 0 1 1 0 1 0 0 11<br />1999.24520 1 0 0 1 1 1 1 1 0 1 1 0 0 1 1 0 1 1 0 12<br />1999.24790 1 1 1 1 1 1 1 1 0 1 1 1 0 1 1 0 1 1 1 16<br />1999.25070 1 1 0 1 1 1 1 1 0 1 1 1 0 1 1 0 1 1 0 14<br />1999.25340 1 1 1 1 1 1 0 1 0 1 1 1 0 1 0 0 1 1 1 14<br />1999.25620 1 1 1 1 1 1 0 1 0 1 1 1 0 1 0 0 1 1 1 14<br />1999.25890 1 1 0 1 1 1 1 1 0 1 1 1 0 1 1 0 1 1 0 14<br />1999.26160 1 1 0 1 1 1 1 1 0 1 1 1 0 1 1 0 1 1 1 15<br />1999.26440 1 1 0 1 1 1 1 1 0 1 1 1 0 1 1 0 1 1 0 14<br />1999.26710 1 1 0 1 1 1 1 1 0 1 1 1 0 1 1 0 1 1 0 14<br />1999.26990 1 0 0 1 1 1 1 1 0 1 1 1 0 1 1 0 1 1 1 14<br />1999.27260 1 0 0 1 1 1 1 1 0 1 1 0 0 1 1 0 1 1 0 12<br />1999.27530 1 0 0 1 1 1 1 1 0 1 1 1 0 1 1 0 1 1 0 13<br />1999.27810 1 0 0 1 1 1 1 1 0 1 1 1 0 1 1 0 1 1 0 13<br />1999.28080 1 1 0 1 1 1 1 1 0 1 1 0 0 1 1 0 1 1 0 13<br />1999.28360 1 0 0 1 1 1 1 1 0 1 0 0 1 1 1 0 1 0 0 11<br />1999.28630 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 1 0 14<br />1999.28900 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 1 0 14<br />1999.29180 1 1 0 0 1 1 1 1 0 1 1 0 1 1 1 0 1 0 0 12<br />1999.29450 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 1 0 14<br />1999.29730 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 1 0 14<br />1999.30000 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 0 0 13<br />1999.30270 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 1 0 14<br />1999.30550 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 1 0 14<br />1999.30820 1 0 0 1 1 1 1 1 0 1 1 1 1 1 1 0 1 1 0 14<br />1999.31100 1 1 0 1 1 1 1 1 0 1 1 1 1 1 1 0 1 1 0 15<br />1999.31370 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 0 0 13<br />1999.31640 1 1 1 1 1 1 1 1 0 1 1 1 1 1 1 0 1 1 1 17<br />1999.31920 1 1 0 1 1 1 1 1 0 1 1 1 1 1 1 0 1 1 0 15<br />1999.32190 1 1 0 1 1 1 1 1 0 1 1 1 1 1 1 0 1 1 1 16<br />1999.32470 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 1 0 14<br />1999.32740 1 1 0 0 1 1 1 1 0 1 1 0 1 1 1 0 1 1 0 13<br />1999.33010 0 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 0 0 12<br />1999.33290 0 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 1 0 13<br />1999.33560 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 1 0 14<br />1999.33840 1 0 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 1 0 13<br />1999.34110 1 1 0 0 1 1 1 1 0 1 1 1 1 1 1 0 0 1 0 13<br />1999.34380 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 1 0 14<br />1999.34660 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 1 0 14<br />1999.34930 0 1 0 1 1 1 1 1 0 1 1 1 1 1 1 0 1 0 0 13<br />1999.35210 1 1 0 1 1 1 1 1 0 1 1 1 1 1 1 0 1 1 0 15<br />1999.35480 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 1 0 14<br />1999.35750 1 1 0 1 1 1 1 1 0 1 1 1 1 1 1 0 1 1 0 15<br />1999.36030 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 1 0 14<br />1999.36300 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 0 0 13<br />1999.36580 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 1 0 14<br />1999.36850 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 1 0 14<br />1999.37120 1 1 0 1 1 1 1 1 0 1 1 1 1 1 1 0 1 1 0 15<br />1999.37400 1 1 0 1 1 1 1 1 0 1 1 1 1 1 1 0 1 1 0 15<br />1999.37670 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 1 0 14<br />1999.37950 1 1 1 1 1 0 1 1 0 1 1 1 1 1 1 0 1 1 1 16<br />1999.38220 0 1 0 1 1 0 1 1 0 1 1 0 1 1 1 0 1 1 0 12<br />1999.38490 1 1 0 1 1 0 1 1 0 1 1 1 1 1 1 0 1 1 0 14<br />1999.38770 1 1 0 1 1 0 1 1 0 1 1 1 1 1 1 0 1 1 1 15<br />1999.39040 1 0 0 1 1 0 1 1 0 1 1 1 1 1 1 0 1 1 0 13<br />1999.39320 0 1 1 1 1 0 1 1 0 1 1 1 1 1 1 0 1 1 1 15<br />1999.39590 1 1 1 1 1 0 0 1 0 1 1 0 1 1 0 0 1 0 1 12<br />1999.39860 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 1 0 14<br />1999.40140 1 1 0 1 1 1 1 1 0 1 1 1 1 1 1 0 1 1 0 15<br />1999.40410 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 1 0 14<br />1999.40680 1 0 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 1 0 13<br />1999.40960 1 0 0 1 1 1 1 1 0 1 1 1 1 1 1 0 1 1 0 14<br />1999.41230 1 0 0 1 1 1 1 1 0 1 1 1 1 1 1 0 1 1 0 14<br />1999.41510 1 1 0 1 1 1 1 1 0 1 1 1 1 1 1 0 0 1 0 14<br />1999.41780 1 1 1 1 1 1 1 1 0 1 1 1 1 1 1 0 1 1 1 17<br />1999.42050 1 1 0 1 1 1 1 1 0 1 1 1 1 1 1 0 1 1 0 15<br />1999.42330 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 1 0 14<br />1999.42600 1 1 0 1 1 1 1 1 0 1 1 1 1 1 1 0 1 1 0 15<br />1999.42880 1 1 0 1 1 1 1 1 0 1 1 1 1 1 1 0 1 1 0 15<br />1999.43150 1 0 0 1 1 1 1 1 0 1 1 1 1 1 1 0 1 1 0 14<br />1999.43420 1 1 0 1 1 1 1 1 0 1 1 1 1 1 1 0 1 1 0 15<br />1999.43700 1 1 0 1 1 1 1 1 0 1 1 1 1 1 1 0 1 1 0 15<br />1999.43970 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 1 0 14<br />1999.44250 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 1 0 14<br />1999.44520 1 1 0 1 1 1 1 1 0 0 1 1 1 1 1 0 1 1 0 14<br />1999.44790 1 0 0 1 1 1 1 1 0 1 1 1 1 1 1 0 1 1 1 15<br />1999.45070 1 1 0 1 1 1 1 1 0 1 1 1 1 1 1 0 1 1 1 16<br />1999.45340 1 1 1 1 1 1 1 0 0 1 1 1 1 1 1 0 1 1 1 16<br />1999.45620 1 1 0 1 1 1 1 1 0 1 1 1 1 1 1 0 1 1 0 15<br />1999.45890 1 1 0 1 1 1 1 1 0 1 1 1 1 1 1 0 1 1 0 15<br />1999.46160 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 1 0 14<br />1999.46440 1 1 0 1 1 1 1 1 0 1 1 1 1 1 1 0 1 1 0 15<br />1999.46710 1 1 0 1 1 1 1 1 0 1 1 1 1 1 1 0 1 1 0 15<br />1999.46990 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 0 0 13<br />1999.47260 1 1 0 1 1 1 1 1 0 1 1 1 1 1 1 0 1 1 0 15<br />1999.47530 1 1 0 1 1 1 1 1 0 1 1 1 1 1 1 0 1 1 0 15<br />1999.47810 1 0 0 1 1 1 1 1 0 1 1 1 1 1 1 0 1 1 0 14<br />1999.48080 1 1 0 1 1 1 1 1 0 1 1 1 1 1 1 0 1 1 0 15<br />1999.48360 0 0 1 1 1 1 1 1 0 1 1 1 1 1 1 0 1 1 1 15<br />1999.48630 1 0 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 1 0 13<br />1999.48900 1 0 0 1 1 1 1 1 0 1 1 1 1 1 1 0 1 1 0 14<br />1999.49180 1 1 0 1 1 1 1 1 0 1 1 1 1 1 1 0 1 1 1 16<br />1999.49450 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 1 0 14<br />1999.49730 1 0 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 1 0 13<br />1999.50000 1 1 1 1 1 1 1 1 0 1 1 1 1 1 1 0 1 1 1 17<br />1999.50270 1 1 1 1 1 1 1 1 0 1 1 1 1 1 1 0 1 1 1 17<br />1999.50550 1 1 0 1 1 1 1 1 0 1 1 1 1 1 1 0 1 1 0 15<br />1999.50820 1 1 0 1 1 1 1 1 0 1 1 1 1 1 1 0 1 1 0 15<br />1999.51100 1 1 1 1 1 1 1 1 0 1 1 1 1 1 1 0 1 1 1 17<br />1999.51640 1 1 1 1 1 1 1 1 0 1 1 1 1 1 1 0 1 1 1 17<br />1999.51920 1 1 0 1 1 1 1 1 0 1 1 1 1 1 1 0 1 1 0 15<br />1999.52190 1 1 0 1 1 1 1 1 0 1 1 1 1 1 1 0 0 1 0 14<br />1999.52470 1 1 1 1 1 1 0 1 0 1 1 1 1 1 0 0 1 1 1 15<br />1999.52740 1 0 1 1 1 0 1 1 0 1 1 1 1 1 1 0 1 1 1 15<br />1999.53010 1 1 0 1 1 1 1 0 0 1 1 0 1 1 1 0 1 1 0 13<br />1999.53290 1 1 0 1 1 1 1 1 0 1 1 1 1 1 1 0 1 1 0 15<br />1999.53560 1 1 1 1 1 1 1 1 0 1 1 1 1 1 1 0 1 1 1 17<br />1999.53840 1 1 0 1 1 1 1 1 0 1 1 1 1 1 1 0 1 1 1 16<br />1999.54110 1 1 1 1 1 1 1 1 0 1 1 1 1 1 1 0 1 1 1 17<br />1999.54380 1 1 1 1 1 1 1 0 0 1 1 1 1 1 1 0 1 1 1 16<br />1999.54660 1 1 1 1 1 1 1 0 0 1 1 1 1 1 1 0 1 1 1 16<br />1999.54930 1 1 0 1 1 1 1 1 0 1 1 1 1 1 1 0 1 1 0 15<br />1999.55210 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 1 0 14<br />1999.55480 1 1 0 1 1 1 1 0 0 1 1 0 1 1 1 0 1 1 0 13<br />1999.55750 1 1 0 1 1 1 1 0 0 1 1 1 1 1 1 0 1 1 0 14<br />1999.56030 1 1 1 1 1 1 1 1 0 0 1 1 1 1 1 0 1 1 1 16<br />1999.56300 1 1 1 1 1 1 1 1 0 0 1 1 1 1 1 0 1 1 1 16<br />1999.56580 1 1 1 1 1 1 1 1 0 1 1 1 1 1 1 0 1 1 1 17<br />1999.56850 1 1 1 1 1 1 1 1 0 1 1 1 1 1 1 0 1 1 1 17<br />1999.57120 1 1 1 1 1 1 1 0 0 1 1 1 1 1 1 0 1 1 1 16<br />1999.57400 1 1 0 1 1 1 1 1 0 1 1 1 1 1 1 0 1 1 0 15<br />1999.57670 1 1 1 1 1 1 1 1 0 1 1 1 1 1 1 0 1 0 1 16<br />1999.57950 1 1 1 1 1 1 1 1 0 1 1 1 1 1 1 0 1 1 1 17<br />1999.58220 1 1 0 1 1 1 1 1 0 1 1 1 1 1 1 0 1 1 1 16<br />1999.58490 1 1 0 1 1 1 1 1 0 1 1 1 1 1 1 0 1 1 0 15<br />1999.58770 1 1 1 1 1 1 1 1 0 1 1 1 1 1 1 0 1 1 1 17<br />1999.59040 1 1 1 1 1 1 1 1 0 1 1 1 1 1 1 0 1 1 1 17<br />1999.59320 1 1 0 1 1 1 1 1 0 1 1 1 1 1 1 0 1 1 0 15<br />1999.59590 1 1 1 1 1 1 1 1 0 1 1 1 1 1 1 0 1 1 0 16<br />1999.59860 1 1 1 1 1 1 1 1 0 1 1 1 1 1 1 0 1 1 1 17<br />1999.60140 0 1 1 1 1 1 1 1 0 1 1 1 1 1 1 0 1 1 1 16<br />1999.60410 1 1 0 1 1 1 1 1 0 1 1 1 1 1 1 0 1 1 0 15<br />1999.60680 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 19<br />1999.60960 1 0 0 1 1 1 1 1 0 1 1 1 1 1 1 0 1 1 0 14<br />1999.61230 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 0 0 13<br />1999.61510 1 1 0 1 1 1 1 1 0 1 1 1 1 1 1 0 1 1 0 15<br />1999.61780 1 1 1 1 1 0 1 1 0 1 1 0 1 1 1 0 1 1 1 15<br />1999.62050 1 1 0 1 1 1 1 1 0 1 1 1 1 1 1 0 1 1 1 16<br />1999.62330 1 1 1 1 1 1 1 1 0 1 1 1 1 1 1 0 1 1 1 17<br />1999.62600 0 1 1 1 1 1 1 1 0 1 1 1 1 1 1 0 1 1 1 16<br />1999.62880 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 1 0 14<br />1999.63150 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 1 0 14<br />1999.63420 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 1 0 14<br />1999.63700 1 0 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 1 0 13<br />1999.63970 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 1 0 14<br />1999.64250 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 0 0 13<br />1999.64520 1 1 0 1 1 1 1 1 0 1 1 1 1 1 1 0 1 1 0 15<br />1999.64790 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 0 1 0 13<br />1999.65070 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 1 0 14<br />1999.65340 1 1 0 1 1 1 0 1 0 1 1 0 1 1 0 0 1 1 0 12<br />1999.65620 1 1 1 1 1 1 0 1 0 1 1 0 1 1 0 0 1 0 0 12<br />1999.65890 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 1 0 14<br />1999.66160 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 1 0 14<br />1999.66440 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 1 0 14<br />1999.66710 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 1 0 14<br />1999.66990 1 1 0 1 1 1 1 1 0 1 1 1 1 1 1 0 1 1 0 15<br />1999.67260 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 0 0 13<br />1999.67530 1 1 0 1 1 0 1 1 0 1 1 0 1 1 1 0 1 1 0 13<br />1999.67810 1 1 1 1 1 0 1 1 0 1 1 1 1 1 1 0 1 1 1 16<br />1999.68080 1 1 0 1 1 0 1 1 0 1 1 1 1 1 1 0 1 1 0 14<br />1999.68360 1 1 0 1 1 0 1 1 0 1 1 1 1 1 1 0 1 1 0 14<br />1999.68630 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 1 0 14<br />1999.68900 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 1 0 14<br />1999.69180 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 1 0 14<br />1999.69450 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 1 0 14<br />1999.69730 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 0 0 13<br />1999.70000 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 0 0 13<br />1999.70270 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 1 0 14<br />1999.70550 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 1 0 14<br />1999.70820 1 1 0 1 1 0 1 1 0 1 1 0 1 1 1 0 1 1 0 13<br />1999.71100 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 1 0 14<br />1999.71370 1 1 0 1 1 1 1 0 0 1 1 0 1 0 1 0 0 0 0 10<br />1999.71640 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 0 1 0 13<br />1999.71920 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 1 0 14<br />1999.72190 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 1 0 14<br />1999.72470 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 0 0 13<br />1999.72740 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 1 0 14<br />1999.73010 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 1 0 14<br />1999.73290 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 1 0 14<br />1999.73560 1 1 0 1 1 1 1 1 0 1 1 0 1 0 1 0 1 1 0 13<br />1999.73840 1 1 0 1 1 1 0 1 0 0 1 0 1 1 0 0 1 0 0 10<br />1999.74110 1 1 0 1 1 1 1 0 0 1 0 0 1 1 1 0 1 0 0 11<br />1999.74380 1 1 1 1 1 1 0 1 0 1 1 0 1 0 0 0 1 0 1 12<br />1999.74660 1 1 0 1 1 0 1 1 0 0 1 0 0 0 1 0 0 1 0 9<br />1999.74930 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 1 0 14<br />1999.75210 1 1 0 1 1 1 1 0 0 1 1 0 1 1 1 0 1 1 0 13<br />1999.75480 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 1 0 14<br />1999.75750 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 1 0 14<br />1999.76030 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 1 0 14<br />1999.76300 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 1 0 14<br />1999.76580 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 1 0 14<br />1999.76850 1 1 0 1 1 0 1 1 0 1 1 0 1 1 1 0 1 1 0 13<br />1999.77120 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 1 0 14<br />1999.77400 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 1 0 14<br />1999.77670 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 1 0 14<br />1999.77950 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 0 0 13<br />1999.78220 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 1 0 14<br />1999.78490 1 1 0 1 1 0 1 1 0 1 1 0 1 1 1 0 1 1 0 13<br />1999.78770 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 0 0 13<br />1999.79040 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 1 0 14<br />1999.79320 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 0 1 0 13<br />1999.79590 1 1 0 1 1 0 1 1 0 1 1 0 1 1 1 0 1 1 0 13<br />1999.79860 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 0 0 13<br />1999.80140 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 1 0 14<br />1999.80410 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 1 0 14<br />1999.80680 1 1 0 1 1 1 1 1 0 0 1 0 1 1 1 0 0 0 0 11<br />1999.80960 1 1 0 1 1 1 0 1 0 1 1 0 1 1 0 0 1 1 0 12<br />1999.81230 1 0 1 1 1 1 0 1 0 1 1 0 1 0 0 0 1 0 1 11<br />1999.81510 1 0 0 1 1 1 0 1 0 1 1 0 1 1 0 0 1 1 0 11<br />1999.81780 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 0 0 13<br />1999.82050 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 1 0 14<br />1999.82330 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 0 0 13<br />1999.82600 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 0 0 13<br />1999.82880 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 0 0 13<br />1999.83150 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 0 0 13<br />1999.83420 1 1 0 1 1 0 1 1 0 1 1 0 1 1 1 0 1 0 0 12<br />1999.83700 1 1 0 1 1 0 1 1 0 1 1 0 1 1 1 0 1 0 0 12<br />1999.83970 1 1 0 1 1 0 1 1 0 1 1 1 1 1 1 0 0 1 0 13<br />1999.84250 1 1 0 1 1 1 1 1 0 1 1 0 0 1 1 0 1 1 0 13<br />1999.84520 1 1 0 1 1 1 1 0 0 1 1 0 1 1 1 0 1 1 0 13<br />1999.84790 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 1 0 14<br />1999.85070 1 1 0 1 1 0 1 1 0 1 1 0 1 1 1 0 1 0 0 12<br />1999.85340 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 0 0 13<br />1999.85620 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 0 0 13<br />1999.85890 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 1 0 14<br />1999.86160 1 1 0 1 1 0 1 1 0 1 1 0 1 1 1 0 1 0 0 12<br />1999.86440 1 1 0 1 1 0 1 1 0 1 1 0 1 1 1 0 1 0 0 12<br />1999.86710 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 0 0 13<br />1999.86990 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 1 0 14<br />1999.87260 1 1 0 1 1 0 1 1 0 1 1 0 1 0 1 0 1 1 0 12<br />1999.87530 1 1 0 1 1 0 1 1 0 1 1 0 1 0 1 0 1 0 0 11<br />1999.87810 0 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 0 0 12<br />1999.88080 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 0 0 13<br />1999.88360 1 1 0 1 1 0 1 1 0 1 1 0 1 1 1 0 1 0 0 12<br />1999.88630 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 0 0 13<br />1999.88900 1 1 0 1 1 1 0 1 0 1 1 0 1 1 0 0 1 0 0 11<br />1999.89180 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 1 0 14<br />1999.89450 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 0 0 13<br />1999.89730 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 1 0 14<br />1999.90000 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 1 0 14<br />1999.90270 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 0 0 13<br />1999.90550 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 1 0 14<br />1999.90820 1 1 0 1 1 0 1 1 0 1 1 0 1 0 1 0 1 0 0 11<br />1999.91100 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 1 0 14<br />1999.91370 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 0 0 13<br />1999.91640 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 0 0 13<br />1999.91920 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 0 0 13<br />1999.92190 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 1 0 14<br />1999.92470 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 1 0 14<br />1999.92740 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 1 0 14<br />1999.93010 1 1 0 1 1 0 1 1 0 1 1 0 1 1 1 0 1 1 0 13<br />1999.93290 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 1 0 14<br />1999.93560 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 1 0 14<br />1999.93840 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 1 0 14<br />1999.94110 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 1 0 14<br />1999.94380 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 1 0 14<br />1999.94660 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 1 0 14<br />1999.94930 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 0 0 13<br />1999.95210 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 1 0 14<br />1999.95480 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 1 0 14<br />1999.95750 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 1 0 14<br />1999.96030 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 1 0 14<br />1999.96300 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 1 0 14<br />1999.96580 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 0 0 13<br />1999.96850 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 1 0 14<br />1999.97120 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 1 0 14<br />1999.97400 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 1 0 14<br />1999.97670 1 1 0 0 1 1 1 1 0 1 1 0 1 1 1 0 1 0 0 12<br />1999.97950 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 0 0 13<br />1999.98220 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 1 0 14<br />1999.98490 1 1 0 0 1 1 1 1 0 1 1 0 0 1 1 0 1 0 0 11<br />1999.98770 1 1 0 0 1 1 1 1 0 1 1 0 1 1 1 0 1 0 0 12<br />1999.99040 1 1 0 0 1 1 1 1 0 1 1 0 1 1 1 0 1 0 0 12<br />1999.99320 1 1 0 0 1 1 1 1 0 1 1 0 1 1 1 0 1 0 0 12<br />1999.99590 1 1 0 1 1 0 1 1 0 1 1 0 1 0 1 0 1 0 0 11<br />1999.99860 1 1 0 1 1 1 1 1 0 1 1 0 1 1 1 0 1 0 0 13<br /> % 96 92 13 97 100 90 96 96 0 96 99 34 93 97 96 0 96 77 15<br /></pre><div class="blogger-post-footer"><img width='1' height='1' src='https://blogger.googleusercontent.com/tracker/102059579936554555-5212111001014147765?l=yftian.blogspot.com' alt='' /></div>Enodhttp://www.blogger.com/profile/08433335265518211485noreply@blogger.com1tag:blogger.com,1999:blog-102059579936554555.post-26243794795727987662009-04-19T07:43:00.002+08:002009-04-19T07:46:31.375+08:00<a onblur="try {parent.deselectBloggerImageGracefully();} catch(e) {}" href="http://2.bp.blogspot.com/_2ATje7JNtMw/Sepl4yzSQsI/AAAAAAAAAK0/V60HQtcL698/s1600-h/sites.around.china.jpg"><img style="display:block; margin:0px auto 10px; text-align:center;cursor:pointer; cursor:hand;width: 320px; height: 194px;" src="http://2.bp.blogspot.com/_2ATje7JNtMw/Sepl4yzSQsI/AAAAAAAAAK0/V60HQtcL698/s320/sites.around.china.jpg" border="0" alt=""id="BLOGGER_PHOTO_ID_5326181535702991554" /></a><br /><br />It turns out the longitude of COSA was wrong in the log file.<br /><br /><pre><br />2. Site Location Information<br /><br /> City or Town : Scottsdale<br /> State or Province : Arizona<br /> Country : USA<br /> Tectonic Plate : North America<br /> Approximate Position (ITRF)<br /> X coordinate (m) : -1982826.826<br /> Y coordinate (m) : -4936878.058<br /> Z coordinate (m) : 3506904.185<br /> Latitude (N is +) : 333407.32<br /> Longitude (E is +) : 1115255.92<br /> Elevation (m,ellips.) : 391.349 <br /> Additional Information : ARP ITRF96 POSITION (EPOCH 1997.0) computed in<br /> Oct. 1998 using 12 days of data.<br /><br /></pre><div class="blogger-post-footer"><img width='1' height='1' src='https://blogger.googleusercontent.com/tracker/102059579936554555-2624379479572798766?l=yftian.blogspot.com' alt='' /></div>Enodhttp://www.blogger.com/profile/08433335265518211485noreply@blogger.com0tag:blogger.com,1999:blog-102059579936554555.post-16605780777083730552009-04-11T23:59:00.002+08:002009-04-12T00:02:02.268+08:00<a onblur="try {parent.deselectBloggerImageGracefully();} catch(e) {}" href="http://3.bp.blogspot.com/_2ATje7JNtMw/SeC-vESw_qI/AAAAAAAAAKk/pB8z5tQaZyw/s1600-h/idl-icon-text.JPG"><img style="margin: 0px auto 10px; display: block; text-align: center; cursor: pointer; width: 320px; height: 226px;" src="http://3.bp.blogspot.com/_2ATje7JNtMw/SeC-vESw_qI/AAAAAAAAAKk/pB8z5tQaZyw/s320/idl-icon-text.JPG" alt="" id="BLOGGER_PHOTO_ID_5323464475367505570" border="0" /></a><br />It's strange! How could this happen? It never happened before.<div class="blogger-post-footer"><img width='1' height='1' src='https://blogger.googleusercontent.com/tracker/102059579936554555-1660578077708373055?l=yftian.blogspot.com' alt='' /></div>Enodhttp://www.blogger.com/profile/08433335265518211485noreply@blogger.com0tag:blogger.com,1999:blog-102059579936554555.post-84050400706195190162009-04-10T00:22:00.002+08:002009-04-18T21:19:59.757+08:00<a onblur="try {parent.deselectBloggerImageGracefully();} catch(e) {}" href="http://2.bp.blogspot.com/_2ATje7JNtMw/SenTbCz-nuI/AAAAAAAAAKs/bMjqlT-uYIs/s1600-h/psbase_ensum.AQUI.jpg"><img style="margin: 0px auto 10px; display: block; text-align: center; cursor: pointer; width: 226px; height: 320px;" src="http://2.bp.blogspot.com/_2ATje7JNtMw/SenTbCz-nuI/AAAAAAAAAKs/bMjqlT-uYIs/s320/psbase_ensum.AQUI.jpg" alt="" id="BLOGGER_PHOTO_ID_5326020495907135202" border="0" /></a> <br /> <br />Coseismic displacements of AQUI for Mw 6.3 central Italy earthquake: <br /><meta equiv="Content-Type" content="text/html; charset=utf-8"><meta name="ProgId" content="Word.Document"><meta name="Generator" content="Microsoft Word 11"><meta name="Originator" content="Microsoft Word 11"><link rel="File-List" href="file:///G:%5CDOCUME%7E1%5Ctianyf%5CLOCALS%7E1%5CTemp%5Cmsohtml1%5C01%5Cclip_filelist.xml"><!--[if gte mso 9]><xml> <w:worddocument> <w:view>Normal</w:View> <w:zoom>0</w:Zoom> <w:punctuationkerning/> <w:validateagainstschemas/> <w:saveifxmlinvalid>false</w:SaveIfXMLInvalid> <w:ignoremixedcontent>false</w:IgnoreMixedContent> <w:alwaysshowplaceholdertext>false</w:AlwaysShowPlaceholderText> <w:compatibility> <w:breakwrappedtables/> <w:snaptogridincell/> <w:wraptextwithpunct/> <w:useasianbreakrules/> <w:dontgrowautofit/> <w:usefelayout/> </w:Compatibility> <w:browserlevel>MicrosoftInternetExplorer4</w:BrowserLevel> </w:WordDocument> </xml><![endif]--><!--[if gte mso 9]><xml> <w:latentstyles deflockedstate="false" latentstylecount="156"> </w:LatentStyles> </xml><![endif]--><style> <!-- /* Font Definitions */ @font-face {font-family:Courier; panose-1:2 7 4 9 2 2 5 2 4 4; mso-font-charset:0; mso-generic-font-family:modern; mso-font-format:other; mso-font-pitch:fixed; mso-font-signature:3 0 0 0 1 0;} @font-face {font-family:宋体; panose-1:2 1 6 0 3 1 1 1 1 1; mso-font-alt:SimSun; mso-font-charset:134; mso-generic-font-family:auto; mso-font-pitch:variable; mso-font-signature:3 135135232 16 0 262145 0;} @font-face {font-family:"\@宋体"; panose-1:2 1 6 0 3 1 1 1 1 1; mso-font-charset:134; mso-generic-font-family:auto; mso-font-pitch:variable; mso-font-signature:3 135135232 16 0 262145 0;} /* Style Definitions */ p.MsoNormal, li.MsoNormal, div.MsoNormal {mso-style-parent:""; margin:0in; margin-bottom:.0001pt; mso-pagination:widow-orphan; font-size:12.0pt; font-family:"Times New Roman"; mso-fareast-font-family:宋体;} @page Section1 {size:8.5in 11.0in; margin:1.0in 1.25in 1.0in 1.25in; mso-header-margin:.5in; mso-footer-margin:.5in; mso-paper-source:0;} div.Section1 {page:Section1;} --> </style><!--[if gte mso 10]> <style> /* Style Definitions */ table.MsoNormalTable {mso-style-name:"Table Normal"; mso-tstyle-rowband-size:0; mso-tstyle-colband-size:0; mso-style-noshow:yes; mso-style-parent:""; mso-padding-alt:0in 5.4pt 0in 5.4pt; mso-para-margin:0in; mso-para-margin-bottom:.0001pt; mso-pagination:widow-orphan; font-size:10.0pt; font-family:"Times New Roman"; mso-fareast-font-family:"Times New Roman"; mso-ansi-language:#0400; mso-fareast-language:#0400; mso-bidi-language:#0400;} </style> <![endif]--> <p class="MsoNormal" style="font-family:courier new;">N offset 1:<span style=""> </span>-0.04025 +-<span style=""> </span>0.00263 (2009.26164)<o:p></o:p></p> <p class="MsoNormal" style="font-family:courier new;">E offset 1:<span style=""> </span>0.00961 +-<span style=""> </span>0.00233 (2009.26164)<o:p></o:p></p> <p class="MsoNormal"><span style="font-family:Courier;"><span style="font-family:courier new;">U offset 1: -0.05367 +- 0.00525 (2009.26164)</span></span></p><p class="MsoNormal"> <br /></p><p class="MsoNormal">Acknowledgement: Raw data from EUREF. <br /><span style="font-family:Courier;"><span style="font-family:courier new;"></span><o:p></o:p></span></p> <br /><div class="blogger-post-footer"><img width='1' height='1' src='https://blogger.googleusercontent.com/tracker/102059579936554555-8405040070619519016?l=yftian.blogspot.com' alt='' /></div>Enodhttp://www.blogger.com/profile/08433335265518211485noreply@blogger.com0tag:blogger.com,1999:blog-102059579936554555.post-41445293569175092402009-03-20T20:18:00.001+08:002009-03-20T20:20:21.311+08:00Big volcanoes and strong earthquake have been occurring in Tonga region. Were they related?<div class="blogger-post-footer"><img width='1' height='1' src='https://blogger.googleusercontent.com/tracker/102059579936554555-4144529356917509240?l=yftian.blogspot.com' alt='' /></div>Enodhttp://www.blogger.com/profile/08433335265518211485noreply@blogger.com0tag:blogger.com,1999:blog-102059579936554555.post-78466913080852267222009-03-19T20:46:00.000+08:002009-03-19T20:47:52.362+08:00It is said that sp3 is much more faster than sp2.<br />It feels good as of now.<div class="blogger-post-footer"><img width='1' height='1' src='https://blogger.googleusercontent.com/tracker/102059579936554555-7846691308085226722?l=yftian.blogspot.com' alt='' /></div>Enodhttp://www.blogger.com/profile/08433335265518211485noreply@blogger.com0tag:blogger.com,1999:blog-102059579936554555.post-67648949416543129482009-03-19T20:44:00.001+08:002009-03-19T20:46:39.142+08:00Run /usr/bin/cygserver-config first!<br /><br />Install proftpd as a XP service:<br />cygrunsrv --install proftpd --path /usr/sbin/proftpd --args --nodaemon --termsig TERM --shutdown <br /><br />Configure /etc/proftpd.conf.<br />Configure /etc/xinetd.d/ftpd as the directions.<div class="blogger-post-footer"><img width='1' height='1' src='https://blogger.googleusercontent.com/tracker/102059579936554555-6764894941654312948?l=yftian.blogspot.com' alt='' /></div>Enodhttp://www.blogger.com/profile/08433335265518211485noreply@blogger.com0tag:blogger.com,1999:blog-102059579936554555.post-40894348400147014612009-03-16T11:14:00.000+08:002009-03-16T11:15:18.300+08:00use bchunk.<div class="blogger-post-footer"><img width='1' height='1' src='https://blogger.googleusercontent.com/tracker/102059579936554555-4089434840014701461?l=yftian.blogspot.com' alt='' /></div>Enodhttp://www.blogger.com/profile/08433335265518211485noreply@blogger.com0tag:blogger.com,1999:blog-102059579936554555.post-50177679479662874902009-03-09T10:37:00.002+08:002009-03-09T10:40:06.527+08:00The links does not work in Vsftpd or ProFtpd servers. You must *mount* the shared directories.<br /><br />In Redhat 5, use<br /> mount --bind src_dir target_dir<br />and in Solaris, use<br /> mount -F lofs src_dir target_dir<br />to mount local directories at other points. <br /><br />This is crucial for construction of a FTP file server.<div class="blogger-post-footer"><img width='1' height='1' src='https://blogger.googleusercontent.com/tracker/102059579936554555-5017767947966287490?l=yftian.blogspot.com' alt='' /></div>Enodhttp://www.blogger.com/profile/08433335265518211485noreply@blogger.com0tag:blogger.com,1999:blog-102059579936554555.post-78898330701826167562009-03-06T10:40:00.003+08:002009-03-06T10:44:53.605+08:00There are errors in my files. I'm distressed.<br /><br />[tianyf@gpsac geoxyz]$ geoxyz 388042.6123 -740382.4073 6302001.9040 2<br /> 1.43977935246514 5.19506102217517 79.0583605368302<br /><br />The correct values are:<br />XYZ coordinates (m): 388042.6102 -740382.4131 6302001.911<br />[Ref. epoch: 2009.0726] [Ref. frame: ITRF2005]<br />Lat/Lon/Ellip. Ht (deg/deg/m): 82.49429410 -62.34046905 78.1523<br />[Ref. epoch: 2009.1767] [Geodetic datum: WGS84]<div class="blogger-post-footer"><img width='1' height='1' src='https://blogger.googleusercontent.com/tracker/102059579936554555-7889833070182616756?l=yftian.blogspot.com' alt='' /></div>Enodhttp://www.blogger.com/profile/08433335265518211485noreply@blogger.com0tag:blogger.com,1999:blog-102059579936554555.post-30133513413122638152009-03-05T08:58:00.001+08:002009-03-05T09:00:00.716+08:00http://www.giss.nasa.gov/tools/panoply/download_gen.html<br /><br />I like it. It benefit me a lot when deal with NetCDF files. Thanks those guys who developed it. :)<div class="blogger-post-footer"><img width='1' height='1' src='https://blogger.googleusercontent.com/tracker/102059579936554555-3013351341312263815?l=yftian.blogspot.com' alt='' /></div>Enodhttp://www.blogger.com/profile/08433335265518211485noreply@blogger.com0tag:blogger.com,1999:blog-102059579936554555.post-40293547177505525382009-02-17T16:12:00.002+08:002009-02-17T16:14:40.800+08:00In the popup Firefox window, the URL http://127.0.0.1:55647/help/index.jsp loads with errors:<br /><br />HTTP ERROR: 500<br /><br />Unable to compile class for JSP<br /><br />Generated servlet error:<br />The type java.lang.Object cannot be resolved. It is indirectly referenced from required .class files<br /><br /><br />RequestURI=/help/index.jsp<br /><br />Powered by Jetty://<div class="blogger-post-footer"><img width='1' height='1' src='https://blogger.googleusercontent.com/tracker/102059579936554555-4029354717750552538?l=yftian.blogspot.com' alt='' /></div>Enodhttp://www.blogger.com/profile/08433335265518211485noreply@blogger.com0tag:blogger.com,1999:blog-102059579936554555.post-48282713523530910832009-02-08T22:28:00.002+08:002009-02-08T22:33:25.668+08:00I found I cannot switch on the wireless radio today. It turns out to be the result of accidentally deleting of HP Wireless Assistant software (which I though useless) several days ago. Furthermore, I found that the hardware switch on the keyboard should be first turned on, then one could turn on/off the wireless in HP Wireless Assistant. Thus, it's a good idea to always keep the hardware radio switch on and use the software switch to enable or disable the wireless signal.<div class="blogger-post-footer"><img width='1' height='1' src='https://blogger.googleusercontent.com/tracker/102059579936554555-4828271352353091083?l=yftian.blogspot.com' alt='' /></div>Enodhttp://www.blogger.com/profile/08433335265518211485noreply@blogger.com0tag:blogger.com,1999:blog-102059579936554555.post-80694238785356009032009-02-08T15:43:00.002+08:002009-02-08T15:44:38.989+08:00According to one newest published paper in GRL, the Tibet Plateau is uplifting with a rate of about 1.14 mm/yr.<br /><br />"4. GPS Data Reveal Uplift of the Tibetan Plateau<br /><br />[7] We use GPS data at the three stations to obtain the vertical displacement rate. Since the early 1990s, several GPS networks for tectonic studies have been established in China and neighboring regions [Wang et al., 2001]. Especially, 25 permanent GPS stations have been established and continuous observations have been made to the present day [Gan et al., 2007] The Lhasa, Kunming, and Dali stations belong to parts of the permanent network; high-quality data are obtained. Here we use GIPSY software to reduce phase and pseudorange data into a site position. The GPS data collected at ca. 20 IGS stations throughout Asia and 18 continuous sites within China are organized into discrete 24 h segments and are processed simultaneously to yield a daily network solution. The precise satellite orbit and clock are fixed to the non-fiducial products provided by the JPL. We solve for 3D coordinates, clock biases, and atmospheric refraction per site with loose constraints. Phase ambiguities are also estimated, but no ambiguity resolution was done (phase ambiguities are estimated as real number and not fixed to integers). We apply azimuth and elevation-dependent antenna phase center models, following the tables recommended by the IGS. We select a subset of IGS sites with a longer time span for definition of reference frame, and rotate the loosely constrained daily solutions onto the ITRF2000 by seven-parameter similarity transformation. We then obtain a time series of vertical components of each site. The time series is then fit in a weighted least-squares sense to a linear function that includes the initial coordinate plus the secular motion rate. Furthermore, assuming a time series characterized by a white noise process, we derive a post-fit RMS of 1.1–1.5 mm for the vertical component of site position and standard deviation of <0.1 mm/yr (1 sigma) for the vertical rate. The rate uncertainty might be overoptimistic because of neglect of the time-correlated noise inherited in the geodetic time series. The rate error might be underestimated by a factor of 5–11 if a pure white noise model is assumed for the coordinate time series [Mao et al., 1999]. Zhang et al. [1997] found that the rate uncertainty might be 2–4 time larger than that for fractional white noise and 3–6 for white plus flicker noises. We adopt a moderate scale factor of five because the time-correlation noise averaged down gradually over time. Multiplying the standard deviation of 0.1 mm/yr by five engenders the rate uncertainty of ca. 0.5 mm/yr, which is comparable with vertical rate precision of 0.4–0.6 mm/yr assigned by Prawirodirdjo and Bock [2004] and Calais et al. [2006]. Finally, we obtain the aggregate station coordinates of the three stations and plot them in Figure 2 (right). Results show that all three stations exhibit positive displacements vertically, implying that the Tibetan Plateau is continuing its uplift. The uplift rates are, respectively, +0.8 ± 0.5, +2.3 ± 0.5 and +0.5 ± 0.5 mm/yr for Lhasa, Kunming, and Dali. Such a mean vertical displacement rate of 1.2 mm/yr agrees well with the conclusion (1.14 mm/yr) obtained on a geological scale [Rowley and Currie, 2006]."<br /><br />Is this really the case?<div class="blogger-post-footer"><img width='1' height='1' src='https://blogger.googleusercontent.com/tracker/102059579936554555-8069423878535600903?l=yftian.blogspot.com' alt='' /></div>Enodhttp://www.blogger.com/profile/08433335265518211485noreply@blogger.com0tag:blogger.com,1999:blog-102059579936554555.post-87734535725649321892009-01-13T01:33:00.002+08:002009-01-13T01:39:12.205+08:00gldas data is in grib1 format. It contains data types that cannot be recognized by ncl_convert2nc. Thus, we can only extract variables using wgrib.<br /><br /><p><br />The gribtab table can be found in the Readme document.<div class="blogger-post-footer"><img width='1' height='1' src='https://blogger.googleusercontent.com/tracker/102059579936554555-8773453572564932189?l=yftian.blogspot.com' alt='' /></div>Enodhttp://www.blogger.com/profile/08433335265518211485noreply@blogger.com144tag:blogger.com,1999:blog-102059579936554555.post-83100922607932137092009-01-13T01:30:00.000+08:002009-01-13T01:32:57.412+08:00large differences. Why?<div class="blogger-post-footer"><img width='1' height='1' src='https://blogger.googleusercontent.com/tracker/102059579936554555-8310092260793213709?l=yftian.blogspot.com' alt='' /></div>Enodhttp://www.blogger.com/profile/08433335265518211485noreply@blogger.com0tag:blogger.com,1999:blog-102059579936554555.post-61911772430577029792009-01-10T23:06:00.001+08:002009-01-10T23:07:54.285+08:00See my origin post <a href="http://groups.google.com/group/comp.lang.fortran/browse_thread/thread/a168181d1bd0481a?hl=en#">http://groups.google.com/group/comp.lang.fortran/browse_thread/thread/a168181d1bd0481a?hl=en#</a>.<br /><br />I found these in Intel website.<br /><br /><h4 style="font-style: italic;">Segmentation Violation with static link to libpthreads</h4> <p style="font-style: italic;">Applications built with libpthreads.a statically linked, (libpthreads.a is used by default when -static is used), may fail with a segmentation violation on some versions of Linux when the applications use more than 2MB of stack space. This is a known limitation of Linuxthreads. If you encounter this problem, link libpthreads dynamically. As an alternative, on Red Hat Linux 9 and Red Hat Enterprise Linux 3.0, you can install the <code>nptl-devel</code> package and pass "<code>-I/usr/include/nptl -L/usr/lib/nptl</code>" on the ifort command line. This will create a statically-linked binary which will run under nptl only, but which does not have the stack size limitation.</p><div class="blogger-post-footer"><img width='1' height='1' src='https://blogger.googleusercontent.com/tracker/102059579936554555-6191177243057702979?l=yftian.blogspot.com' alt='' /></div>Enodhttp://www.blogger.com/profile/08433335265518211485noreply@blogger.com1tag:blogger.com,1999:blog-102059579936554555.post-56909173990740827712009-01-08T23:27:00.003+08:002009-01-15T01:02:52.588+08:00wgrib gives:<br /><br />tianyf@tyf ~/gldas]$ wgrib -s GLDAS_NOAH025_M.A200701.001.grb<br />1:0:d=07010100:NSWRS:sfc:0-3hr ave:NAve=12<br />2:563720:d=07010100:NLWRS:sfc:0-3hr ave:NAve=12<br />3:1097066:d=07010100:LHTFL:sfc:0-3hr ave:NAve=12<br />4:1660786:d=07010100:SHTFL:sfc:0-3hr ave:NAve=12<br /><br />Undefined parameter table (center 221-221 table 1), using NCEP-opn<br />5:2224506:d=07010100:GFLUX:sfc:0-3hr ave:NAve=12<br />6:2727476:d=07010100:LFTX:sfc:anl:NAve=12<br />7:3291196:d=07010100:4LFTX:sfc:anl:NAve=12<br />8:3885292:d=07010100:EVP:sfc:anl:NAve=12<br />9:4388262:d=07010100:SSRUN:sfc:anl:NAve=12<br />10:4951982:d=07010100:BGRUN:sfc:anl:NAve=12<br />11:5546078:d=07010100:SNOM:sfc:0-3hr ave:NAve=12<br />12:5988298:d=07010100:BVF2:sfc:anl:NAve=12<br />13:6521644:d=07010100:WEASD:sfc:anl:NAve=12<br />14:7115740:d=07010100:TSOIL:0-4 cm down:anl:NAve=12<br />15:7649086:d=07010100:TSOIL:0-3 cm down:anl:NAve=12<br />16:8182432:d=07010100:TSOIL:0-2 cm down:anl:NAve=12<br />17:8715778:d=07010100:TSOIL:0-1 cm down:anl:NAve=12<br />18:9218748:d=07010100:SOILM:0-4 cm down:anl:NAve=12<br />19:9812844:d=07010100:SOILM:0-3 cm down:anl:NAve=12<br />20:10406940:d=07010100:SOILM:0-2 cm down:anl:NAve=12<br />21:11001036:d=07010100:SOILM:0-1 cm down:anl:NAve=12<br />22:11595132:d=07010100:TCDC:sfc:anl:NAve=12<br />23:12098102:d=07010100:WIND:sfc:anl:NAve=12<br />24:12540322:d=07010100:TMP:sfc:anl:NAve=12<br />25:13043292:d=07010100:SPFH:sfc:anl:NAve=12<br />26:13303258:d=07010100:PRES:sfc:anl:NAve=12<br />27:13897354:d=07010100:DSWRF:sfc:0-3hr ave:NAve=12<br />28:14491450:d=07010100:DLWRF:sfc:0-3hr ave:NAve=12<br /><br />However, when convert it to .nc format using ncl,<br /><br />[tianyf@tyf ~/gldas]$ ncl_convert2nc GLDAS_NOAH025_M.A200701.001.grb -v weasd<br />Processing file: GLDAS_NOAH025_M.A200701.001.grb...<br />Copyright (C) 1995-2007 - All Rights Reserved<br />University Corporation for Atmospheric Research<br />NCAR Command Language Version 5.0.0<br />The use of this software is governed by a License Agreement.<br />See http://www.ncl.ucar.edu/ for more details.<br />warning:NclGRIB: Unknown grib parameter number detected (155, center 221, table version 1 grib record 5), using default variable name (VAR_155)<br />warning:NclGRIB: Unknown grib parameter number detected (131, center 221, table version 1 grib record 6), using default variable name (VAR_131)<br />warning:NclGRIB: Unknown grib parameter number detected (132, center 221, table version 1 grib record 7), using default variable name (VAR_132)<br />warning:NclGRIB: Unknown grib parameter number detected (235, center 221, table version 1 grib record 9), using default variable name (VAR_235)<br />warning:NclGRIB: Unknown grib parameter number detected (234, center 221, table version 1 grib record 10), using default variable name (VAR_234)<br />warning:NclGRIB: Unknown grib parameter number detected (138, center 221, table version 1 grib record 12), using default variable name (VAR_138)<br />warning:NclGRIB: Unknown grib parameter number detected (204, center 221, table version 1 grib record 27), using default variable name (VAR_204)<br />warning:NclGRIB: Unknown grib parameter number detected (205, center 221, table version 1 grib record 28), using default variable name (VAR_205)<br />warning:getfilevartypes: variable (weasd) does not exist in file (GLDAS_NOAH025_M.A200701.001)<br />fatal:The result of the conditional expression yields a missing value. NCL can not determine branch, see ismissing function<br />fatal:Execute: Error occurred at or near line 351 in file /tmp/tmp8895.ncl<br /><br /><br />ncl does handle other grib file correctly,<br /><br />[tianyf@tyf ~/ds607]$ ncl_convert2nc 19981201.NOAH_d.grb<br />Processing file: 19981201.NOAH_d.grb...<br />Copyright (C) 1995-2007 - All Rights Reserved<br />University Corporation for Atmospheric Research<br />NCAR Command Language Version 5.0.0<br />The use of this software is governed by a License Agreement.<br />See http://www.ncl.ucar.edu/ for more details.<br /><br />Thus, I can only use wgrib to extract variables from grib files.<br /><br />wgrib -s GLDAS_NOAH025_M.A200701.001.grb | grep WEASD | wgrib -i -s -text GLDAS_NOAH025_M.A200701.001.grb -o o<div class="blogger-post-footer"><img width='1' height='1' src='https://blogger.googleusercontent.com/tracker/102059579936554555-5690917399074082771?l=yftian.blogspot.com' alt='' /></div>Enodhttp://www.blogger.com/profile/08433335265518211485noreply@blogger.com0tag:blogger.com,1999:blog-102059579936554555.post-5961950818553497402009-01-08T13:31:00.001+08:002009-01-08T13:34:23.711+08:00Error "the application is being terminated because of memory corruption" when you exit Acrobat 7 (on Windows)<br /><br />Issue<br /><br />When you exit Acrobat 7, Acrobat returns the error message,"the application is being terminated because of memory corruption."<br /><br />Solutions<br /><br />Do one of the following solutions:<br /><br />Solution 1: Install the Acrobat 7.0.8 update.<br /><br />If you have Acrobat 7.0, choose Help > Check For Updates Now to install the update. Alternately, you can download the update from the Adobe website at www.adobe.com/downloads/ .<br /><br />Solution 2: Remove the files from theAcrobat Organizer folder.<br /><br />When you remove files from the Acrobat Organizer folder, you force Acrobat to rebuild the Organizer database. This process replaces any damaged files in the Organizer that may be causing the error message.<br /><br />To remove the files from the Acrobat Organizer folder:<br /><br />1.Choose Start > Run.<br /><br />2.In the Open text box, type:%UserProfile%\Application Data\Adobe\Acrobat\7.0\Organizer70.<br /><br />3.Delete or move all filesfrom this directory.<br /><br />Solution3: Remove the files from the Updaterfolder.<br /><br />You may receive this error message after a failed update, indicating that the files in the Updater folder are corrupted.<br /><br />To remove the files from the Updater folder:<br /><br />1. Choose Start > Run.<br /><br />2. In the Open text box, type:%UserProfile%\Application Data\Adobe\Acrobat\7.0\Updater.<br /><br />3. Delete or move all the files fromthis directory.<br /><br />Background information<br /><br />In rare instances, the Acrobat 7.0.0 or 7.0.1 Organizer database is corrupted on exit if the update manager has been run. This issue is resolved in the Acrobat 7.0.5 update.<div class="blogger-post-footer"><img width='1' height='1' src='https://blogger.googleusercontent.com/tracker/102059579936554555-596195081855349740?l=yftian.blogspot.com' alt='' /></div>Enodhttp://www.blogger.com/profile/08433335265518211485noreply@blogger.com2tag:blogger.com,1999:blog-102059579936554555.post-33509147926812877882009-01-07T15:49:00.001+08:002009-01-07T15:51:22.527+08:00<pre><br /><br />--------[ Summary ]-----------------------------------------------------------------------------------------------------<br /><br /> Computer:<br /> Computer Type ACPI Multiprocessor PC (Mobile)<br /> Operating System Microsoft Windows XP Professional<br /> OS Service Pack Service Pack 2<br /> Internet Explorer 6.0.2900.2180 (IE 6.0 SP2)<br /> DirectX 4.09.00.0904 (DirectX 9.0c)<br /> Motherboard:<br /> CPU Type Mobile DualCore Intel Core 2 Duo T5500, 1666 MHz (10 x 167)<br /> Motherboard Name Hewlett-Packard HP Compaq nx6330 (RQ708PA#AB2)<br /> Motherboard Chipset Mobile Intel Calistoga-PM i945PM<br /> System Memory 3072 MB (DDR2-667 DDR2 SDRAM)<br /> DIMM1: Kingston 99U5295-011.A00LF 2 GB DDR2-667 DDR2 SDRAM (5-5-5-15 @ 333 MHz) (4-4-4-12 @ 266 MHz) (3-3-3-9 @ 200 MHz)<br /> DIMM3: Kingston 1 GB DDR2-667 DDR2 SDRAM (5-5-5-15 @ 333 MHz) (4-4-4-12 @ 266 MHz) (3-3-3-9 @ 200 MHz)<br /> BIOS Type Compaq (07/31/07)<br /><br /> Display:<br /> Video Adapter ATI Mobility Radeon X1300 (512 MB)<br /> Video Adapter ATI Mobility Radeon X1300 (512 MB)<br /> 3D Accelerator ATI Mobility Radeon X1300 (M52)<br /> Monitor Generic Television<br /> Monitor LG Philips LP141WX1-TLA1 [14.1" LCD]<br /><br /> Multimedia:<br /> Audio Adapter Analog Devices AD1981HD @ Intel 82801GBM ICH7-M - High Definition Audio Controller [A-1]<br /><br /> Storage:<br /> IDE Controller Intel(R) 82801G (ICH7 Family) Ultra ATA Storage Controllers - 27DF<br /> IDE Controller Intel(R) 82801GBM SATA AHCI Controller<br /> Storage Controller ALNVHAHV IDE Controller<br /> Disk Drive TOSHIBA MK8034GSX (60 GB, 5400 RPM, SATA)<br /> Optical Drive JAPWFQZ X6F89QR0 SCSI CdRom Device<br /> Optical Drive MATSHITA UJDA770 DVD/CDRW (DVD:8x, CD:24x/24x/24x DVD-ROM/CD-RW)<br /> SMART Hard Disks Status OK<br /><br /> Partitions:<br /> C: (NTFS) 15001 MB (1293 MB free)<br /> D: (NTFS) 20487 MB (2341 MB free)<br /> E: (NTFS) 9996 MB (413 MB free)<br /> F: (NTFS) 3750 MB (1039 MB free)<br /> H: (NTFS) 9501 MB (2823 MB free)<br /> Total Size 58735 MB (7909 MB free)<br /><br /> Input:<br /> Keyboard HID Keyboard Device<br /> Keyboard Standard 101/102-Key or Microsoft Natural PS/2 Keyboard<br /> Mouse HID-compliant mouse<br /> Mouse Synaptics PS/2 Port TouchPad<br /><br /> Network:<br /> Primary IP Address 192.168.11.23<br /> Primary MAC Address 00-17-A4-DA-C9-61<br /> Network Adapter Broadcom NetLink (TM) Gigabit Ethernet (192.168.11.23)<br /> Network Adapter Intel(R) PRO/Wireless 3945ABG Network Connection<br /> Network Adapter VMware Virtual Ethernet Adapter for VMnet8 (192.168.213.1)<br /> Modem Agere Systems HDA Modem<br /><br /> Peripherals:<br /> Printer HP LaserJet 1020<br /> Printer LaserJet 1010 Series Driver<br /> Printer Adobe PDF<br /> Printer Auto HP LaserJet 1020 on JIANG<br /> Printer HP LaserJet 1020<br /> Printer HP Mobile Printing PS<br /> Printer Microsoft Office Document Image Writer<br /> Printer Microsoft XPS Document Writer<br /> Printer PDFConverter<br /> USB1 Controller Intel 82801GBM ICH7-M - USB Universal Host Controller [A-1]<br /> USB1 Controller Intel 82801GBM ICH7-M - USB Universal Host Controller [A-1]<br /> USB1 Controller Intel 82801GBM ICH7-M - USB Universal Host Controller [A-1]<br /> USB1 Controller Intel 82801GBM ICH7-M - USB Universal Host Controller [A-1]<br /> USB2 Controller Intel 82801GBM ICH7-M - Enhanced USB2 Controller [A-1]<br /> USB Device USB Human Interface Device<br /> Battery Microsoft AC Adapter<br /> Battery Microsoft ACPI-Compliant Control Method Battery<br /> Battery Microsoft ACPI-Compliant Control Method Battery<br /><br /> DMI:<br /> DMI BIOS Vendor Hewlett-Packard<br /> DMI BIOS Version 68YTS Ver. F.08<br /> DMI System Manufacturer Hewlett-Packard<br /> DMI System Product HP Compaq nx6330 (RQ708PA#AB2)<br /> DMI System Version F.08<br /> DMI System Serial Number <br /> DMI System UUID 13F32941-97D7DD11-2B906D99-0E144529<br /> DMI Motherboard Manufacturer Hewlett-Packard<br /> DMI Motherboard Product 0A4C<br /> DMI Motherboard Version KBC Version 63.09<br /> DMI Motherboard Serial Number <br /> DMI Chassis Manufacturer Hewlett-Packard<br /> DMI Chassis Version <br /> DMI Chassis Serial Number <br /> DMI Chassis Asset Tag <br /> DMI Chassis Type Notebook<br /><br /></pre><div class="blogger-post-footer"><img width='1' height='1' src='https://blogger.googleusercontent.com/tracker/102059579936554555-3350914792681287788?l=yftian.blogspot.com' alt='' /></div>Enodhttp://www.blogger.com/profile/08433335265518211485noreply@blogger.com14tag:blogger.com,1999:blog-102059579936554555.post-21304106966689890262009-01-06T01:27:00.000+08:002009-01-06T01:29:04.015+08:00much much smaller that of ncep rearead-2.<div class="blogger-post-footer"><img width='1' height='1' src='https://blogger.googleusercontent.com/tracker/102059579936554555-2130410696668989026?l=yftian.blogspot.com' alt='' /></div>Enodhttp://www.blogger.com/profile/08433335265518211485noreply@blogger.com0tag:blogger.com,1999:blog-102059579936554555.post-14662412465861696832009-01-05T19:34:00.001+08:002009-01-05T19:34:50.702+08:00Try set this environment variable:<br /><pre>export CPPFLAGS=-DpgiFortran<br /><br /></pre><div class="blogger-post-footer"><img width='1' height='1' src='https://blogger.googleusercontent.com/tracker/102059579936554555-1466241246586169683?l=yftian.blogspot.com' alt='' /></div>Enodhttp://www.blogger.com/profile/08433335265518211485noreply@blogger.com0tag:blogger.com,1999:blog-102059579936554555.post-89516297313316541342008-12-20T22:43:00.000+08:002008-12-20T22:44:04.510+08:00<a href="http://hpiers.obspm.fr/eop-pc/">http://hpiers.obspm.fr/eop-pc/</a><br /><br /><pre>A leap second in 2008<br />Daniel Gambis<br />Observatoire de Paris<br />December 2008<br /><br /><br />On December 31, 2008 a "leap second" will be added to the world's clocks at<br />23 hours, 59 minutes and 59 seconds Coordinated Universal Time (UTC).<br />In France, due to the +1 hour time offset with respect to the meridian origin,<br />the talking clock will announce the affected "fourth beep" at one<br />o'clock legal time.<br /><br />The Coordinated Universal Time (UTC, replacing GMT) is the reference time<br />scale derived from the Temps Atomique International (TAI) calculated by the<br />Bureau International des Poids et Mesures (BIPM) using a worldwide network of<br />atomic clocks. UTC differs from TAI by an integer number of seconds. It is the<br />basis of all activities in the world. UT1 is the time scale based on the<br />observation of the Earth's rotation. It is now derived mostly from the<br />observation of extragalactic radio-sources by Very Long Baseline<br />Interferometry techniques (VLBI).<br /><br />The various irregular fluctuations progressively detected in the rotation rate<br />of the Earth lead in 1972 to the replacement of UT1 by UTC as the reference<br />time scale. However, it was desired by the scientific community to maintain<br />the difference UT1-UTC smaller than 0.9 second to ensure agreement between<br />the physical and astronomical time scales.<br /><br />Why this extra second? It exists because the rotation of the Earth on its axis,<br />which determines the passing of days and nights, slows down over a long period,<br />mainly as a consequence of Moon-Sun attraction effects. In addition, the Earth<br />is affected by its internal (core, mantle) and external (atmosphere, oceans)<br />constituents. Nowadays, though, time is measured by procedures impervious to<br />our planet’s moods, thanks to around 250 atomic clocks belonging to several<br />countries. Together they are used to calculate UTC. In addition, we have to<br />consider that the length of the day is nowadays 2 ms longer than in the year 1820.<br /><br />Not surprisingly then, the Earth's rotation slowly gets out of synchronization<br />with UTC.<br /><br />In view of a 1972 international agreement stipulating that the difference<br />between the two should never exceed one second, it is necessary from time<br />to time to add intercalated or leap seconds to UTC. On 1st January 2009, the<br />difference will be 34 seconds. Since 1972, leap seconds have been added with a<br />rate interval varying from six months to seven years, with the last being<br />inserted on December 31, 2005.<br /><br />Since the adoption of this system in 1972, partly due to the initial choice of<br />the value of the second and secondly to the general slowing down of the Earth's<br />rotation, it has been necessary to add 23 s to UTC.<br /><br />The decision to introduce a leap second in UTC is the responsibility of the<br />Earth Orientation Center of the International Earth Rotation and reference<br />System Service (IERS).<br /><br />This center is located at the Observatoire de Paris at SYRTE department.<br /><br /></pre><div class="blogger-post-footer"><img width='1' height='1' src='https://blogger.googleusercontent.com/tracker/102059579936554555-8951629731331654134?l=yftian.blogspot.com' alt='' /></div>Enodhttp://www.blogger.com/profile/08433335265518211485noreply@blogger.com0tag:blogger.com,1999:blog-102059579936554555.post-9340056394391651862008-12-11T09:44:00.001+08:002008-12-11T09:46:42.258+08:00Google will terminate your account in accordance with the terms of service if you fail to login to your account for a period of <span style="font-weight: bold;">nine months.<br /><br /></span>Yahoo<span style="font-weight: bold;">: 4 months.<br /></span><div class="blogger-post-footer"><img width='1' height='1' src='https://blogger.googleusercontent.com/tracker/102059579936554555-934005639439165186?l=yftian.blogspot.com' alt='' /></div>Enodhttp://www.blogger.com/profile/08433335265518211485noreply@blogger.com0tag:blogger.com,1999:blog-102059579936554555.post-79167104147129002822008-12-05T23:01:00.000+08:002008-12-05T23:02:29.806+08:00<a href="http://sail.ucsd.edu/%7Ewalker/grdraster.html">http://sail.ucsd.edu/~walker/grdraster.html</a><br /><br /><h2>Getting grdraster to work</h2> I am referring to the GMT program called "grdraster". Basically, I wanted a global topography database with 30-sec resolution, so I downloaded all the GTOPO30 database DEM files, and wanted to piece them together and form one big NetCDF grd file, which is the gridding format GMT uses. Well, I had problems with seg faults, so here is what I found. First, you need to make sure that the filename in the grdraster.info file includes the full path. This critical file is located in the $GMTHOME/share/dbase directory as a default for GMT3.3.6. To save you time if you don't already have it, here is an example of what you need in your grdraster.info file: <p> </p><pre>8 "GTOPO30 E020N40" "m" -R20/60/-10/40 -I0.5m P i 1 0 -9999 /janeway/usr/local/GMT3.3.6/share/dbase/E020N40.DEM<br />9 "GTOPO30 E020N90" "m" -R20/60/40/90 -I0.5m P i 1 0 -9999 /janeway/usr/local/GMT3.3.6/share/dbase/E020N90.DEM<br />10 "GTOPO30 E020S10" "m" -R20/60/-60/-10 -I0.5m P i 1 0 -9999 /janeway/usr/local/GMT3.3.6/share/dbase/E020S10.DEM<br />11 "GTOPO30 E060N40" "m" -R60/100/-10/40 -I0.5m P i 1 0 -9999 /janeway/usr/local/GMT3.3.6/share/dbase/E060N40.DEM<br />12 "GTOPO30 E060N90" "m" -R60/100/40/90 -I0.5m P i 1 0 -9999 /janeway/usr/local/GMT3.3.6/share/dbase/E060N90.DEM<br />13 "GTOPO30 E060S10" "m" -R60/100/-60/-10 -I0.5m P i 1 0 -9999 /janeway/usr/local/GMT3.3.6/share/dbase/E060S10.DEM<br />14 "GTOPO30 E060S60" "m" -R60/100/-90/-60 -I0.5m P i 1 0 -9999 /janeway/usr/local/GMT3.3.6/share/dbase/E060S60.DEM<br />15 "GTOPO30 E100N40" "m" -R100/140/-10/40 -I0.5m P i 1 0 -9999 /janeway/usr/local/GMT3.3.6/share/dbase/E100N40.DEM<br />16 "GTOPO30 E100N90" "m" -R100/140/40/90 -I0.5m P i 1 0 -9999 /janeway/usr/local/GMT3.3.6/share/dbase/E100N90.DEM<br />17 "GTOPO30 E100S10" "m" -R100/140/-60/-10 -I0.5m P i 1 0 -9999 /janeway/usr/local/GMT3.3.6/share/dbase/E100S10.DEM<br />18 "GTOPO30 E120S60" "m" -R120/160/-90/-60 -I0.5m P i 1 0 -9999 /janeway/usr/local/GMT3.3.6/share/dbase/E120S60.DEM<br />19 "GTOPO30 E140N40" "m" -R140/180/-10/40 -I0.5m P i 1 0 -9999 /janeway/usr/local/GMT3.3.6/share/dbase/E140N40.DEM<br />20 "GTOPO30 E140N90" "m" -R140/180/40/90 -I0.5m P i 1 0 -9999 /janeway/usr/local/GMT3.3.6/share/dbase/E140N90.DEM<br />22 "GTOPO30 E140S10" "m" -R140/180/-60/-10 -I0.5m P i 1 0 -9999 /janeway/usr/local/GMT3.3.6/share/dbase/E140S10.DEM<br />23 "GTOPO30 W000S60" "m" -R0/60/-90/-60 -I0.5m P i 1 0 -9999 /janeway/usr/local/GMT3.3.6/share/dbase/W000S60.DEM<br />24 "GTOPO30 W020N40" "m" -R-20/20/-10/40 -I0.5m P i 1 0 -9999 /janeway/usr/local/GMT3.3.6/share/dbase/W020N40.DEM<br />25 "GTOPO30 W020N90" "m" -R-20/20/40/90 -I0.5m P i 1 0 -9999 /janeway/usr/local/GMT3.3.6/share/dbase/W020N90.DEM<br />26 "GTOPO30 W020S10" "m" -R-20/20/-60/-10 -I0.5m P i 1 0 -9999 /janeway/usr/local/GMT3.3.6/share/dbase/W020S10.DEM<br />27 "GTOPO30 W060N40" "m" -R-60/-20/-10/40 -I0.5m P i 1 0 -9999 /janeway/usr/local/GMT3.3.6/share/dbase/W060N40.DEM<br />28 "GTOPO30 W060N90" "m" -R-60/-20/40/90 -I0.5m P i 1 0 -9999 /janeway/usr/local/GMT3.3.6/share/dbase/W060N90.DEM<br />29 "GTOPO30 W060S10" "m" -R-60/-20/-60/-10 -I0.5m P i 1 0 -9999 /janeway/usr/local/GMT3.3.6/share/dbase/W060S10.DEM<br />30 "GTOPO30 W060S60" "m" -R-60/-20/-90/-60 -I0.5m P i 1 0 -9999 /janeway/usr/local/GMT3.3.6/share/dbase/W060S60.DEM<br />31 "GTOPO30 W100N40" "m" -R-100/-60/-10/40 -I0.5m P i 1 0 -9999 /janeway/usr/local/GMT3.3.6/share/dbase/W100N40.DEM<br />32 "GTOPO30 W100N90" "m" -R-100/-60/40/90 -I0.5m P i 1 0 -9999 /janeway/usr/local/GMT3.3.6/share/dbase/W100N90.DEM<br />33 "GTOPO30 W100S10" "m" -R-100/-60/-60/-10 -I0.5m P i 1 0 -9999 /janeway/usr/local/GMT3.3.6/share/dbase/W100S10.DEM<br />34 "GTOPO30 W120S60" "m" -R-120/-80/-90/-60 -I0.5m P i 1 0 -9999 /janeway/usr/local/GMT3.3.6/share/dbase/W120S60.DEM<br />35 "GTOPO30 W140N40" "m" -R-140/-100/-10/40 -I0.5m P i 1 0 -9999 /janeway/usr/local/GMT3.3.6/share/dbase/W140N40.DEM<br />36 "GTOPO30 W140N90" "m" -R-140/-100/40/90 -I0.5m P i 1 0 -9999 /janeway/usr/local/GMT3.3.6/share/dbase/W140N90.DEM<br />37 "GTOPO30 W140S10" "m" -R-140/-100/-60/-10 -I0.5m P i 1 0 -9999 /janeway/usr/local/GMT3.3.6/share/dbase/W140S10.DEM<br />38 "GTOPO30 W180N40" "m" -R-180/-140/-10/40 -I0.5m P i 1 0 -9999 /janeway/usr/local/GMT3.3.6/share/dbase/W180N40.DEM<br />39 "GTOPO30 W180N90" "m" -R-180/-140/40/90 -I0.5m P i 1 0 -9999 /janeway/usr/local/GMT3.3.6/share/dbase/W180N90.DEM<br />40 "GTOPO30 W180S10" "m" -R-180/-140/-60/-10 -I0.5m P i 1 0 -9999 /janeway/usr/local/GMT3.3.6/share/dbase/W180S10.DEM<br />41 "GTOPO30 W180S60" "m" -R-180/-140/-90/-60 -I0.5m P i 1 0 -9999 /janeway/usr/local/GMT3.3.6/share/dbase/W180S60.DEM<br /></pre> <p> After that is in place, then you should be able to run the following: </p><pre>> grdraster 11 -V -R60/100/-10/40 -I0.5m -Gtestoutput.grd<br />grdraster: Reading from raster /janeway/usr/local/GMT3.3.6/share/dbase/E060N40.DEM<br />grdraster: Finished reading from /janeway/usr/local/GMT3.3.6/share/dbase/E060N40.DEM<br />grdraster: min max and # NaN found: -32768 32539 0<br /></pre> <p> Make sure you use the -V option to be verbose if you are having problems. Also make sure you don't forget the -G in front of the filename, and that there are no spaces between the "G" and the first letter of the filename. </p><p> To check that everything went fine, try grdinfo on your newly created grd file: </p><pre>> grdinfo testoutput.grd<br />testoutput.grd: Title: GTOPO30 E060N40<br />testoutput.grd: Command: grdraster 4 -V -R60/100/-10/40 -I0.5m -Gtest.grd<br />testoutput.grd: Remark: /janeway/usr/local/GMT3.3.6/share/dbase/E060N40.DEM<br />testoutput.grd: Pixel node registration used<br />testoutput.grd: grdfile format # 0<br />testoutput.grd: x_min: 60 x_max: 100 x_inc: 0.00833333 units: Longitude in degrees nx: 4800<br />testoutput.grd: y_min: -10 y_max: 40 y_inc: 0.00833333 units: Latitude in degrees ny: 6000<br />testoutput.grd: z_min: -32768 z_max: 32539 units: m<br />testoutput.grd: scale_factor: 1 add_offset: 0<br /></pre> <h3>Additional Notes</h3> If you are troubleshooting, here are some notes that may help: <p> </p><li>grdraster.info only needs accurate paths/filenames of those files being called on </li><li>grdraster.info file numbers can be out of order or can have gaps in the sequence </li><li>compiling grdraster.info needs the netcdf.h include file or link to in the $GMTHOME/include directory <h3>Converting and pasting the DEM files to one GRD file</h3> This will come in the near future. Sorry for the inconvenience. </li><div class="blogger-post-footer"><img width='1' height='1' src='https://blogger.googleusercontent.com/tracker/102059579936554555-7916710414712900282?l=yftian.blogspot.com' alt='' /></div>Enodhttp://www.blogger.com/profile/08433335265518211485noreply@blogger.com96tag:blogger.com,1999:blog-102059579936554555.post-12277887433521077142008-12-04T16:01:00.000+08:002008-12-04T16:02:19.311+08:00ftp://cddis.gsfc.nasa.gov/gps/products/wwww/repro#/cc#wwwwd.typ.Z<br /># Reprocessing campaign number: 1, 2, …<br />cc first two characters of the AC code<br />ig# for combined IGS products<br />wwww GPS week<br />d Day of week<br />typ File type: snx, ssc, sp3, clk, erp, sum<div class="blogger-post-footer"><img width='1' height='1' src='https://blogger.googleusercontent.com/tracker/102059579936554555-1227788743352107714?l=yftian.blogspot.com' alt='' /></div>Enodhttp://www.blogger.com/profile/08433335265518211485noreply@blogger.com1tag:blogger.com,1999:blog-102059579936554555.post-41279990935185800982008-12-03T21:56:00.000+08:002008-12-03T21:58:19.457+08:00the resulting phase is relative to the first observation.<br /><br />Omg<div class="blogger-post-footer"><img width='1' height='1' src='https://blogger.googleusercontent.com/tracker/102059579936554555-4127999093518580098?l=yftian.blogspot.com' alt='' /></div>Enodhttp://www.blogger.com/profile/08433335265518211485noreply@blogger.com0tag:blogger.com,1999:blog-102059579936554555.post-63583202659986097262008-12-03T21:17:00.001+08:002008-12-03T21:54:36.243+08:00even several days<div class="blogger-post-footer"><img width='1' height='1' src='https://blogger.googleusercontent.com/tracker/102059579936554555-6358320265998609726?l=yftian.blogspot.com' alt='' /></div>Enodhttp://www.blogger.com/profile/08433335265518211485noreply@blogger.com0tag:blogger.com,1999:blog-102059579936554555.post-77595987465144963572008-12-03T21:01:00.000+08:002008-12-03T21:16:44.351+08:00the wefit and  icosta cannot bebew supported.<div class="blogger-post-footer"><img width='1' height='1' src='https://blogger.googleusercontent.com/tracker/102059579936554555-7759598746514496357?l=yftian.blogspot.com' alt='' /></div>Enodhttp://www.blogger.com/profile/08433335265518211485noreply@blogger.com0tag:blogger.com,1999:blog-102059579936554555.post-28668103600968848452008-12-03T20:46:00.001+08:002008-12-03T20:46:27.866+08:00bash-3.2# su mobile<br />$ bash<br />bash-3.2$ /var/Applications/Installer.app/Installer --launchedFromSB<br />2008-12-03 20:43:10.289 Installer[305:d03] ATInstaller: Initializing...<br />2008-12-03 20:43:10.298 Installer[305:d03] ATPackageManager: Initializing...<br />DataFile: file open error: /var/mobile/Library/Keyboard/zh_CN-dynamic-text.dat,<br />(Permission denied)<br />^C<div class="blogger-post-footer"><img width='1' height='1' src='https://blogger.googleusercontent.com/tracker/102059579936554555-2866810360096884845?l=yftian.blogspot.com' alt='' /></div>Enodhttp://www.blogger.com/profile/08433335265518211485noreply@blogger.com0tag:blogger.com,1999:blog-102059579936554555.post-31882006378436158982008-12-02T19:30:00.001+08:002008-12-02T19:30:41.272+08:00<table border="0" cellpadding="0" cellspacing="0" width="358"><tbody><tr> <td><a id="Free access for the next 30 days, read and download all articles from these journals" name="Free access for the next 30 days, read and download all articles from these journals"><span originaltag="yes" style="font-family:Arial, Helvetica, sans-serif;font-size:130%;color:#ee7900;">Free access for the next 30 days, read and download all articles from these journals</span></a></td></tr> <tr> <td height="4"><br /></td> </tr><tr> <td bgcolor="#ee7900" height="1"><br /></td> </tr><tr> <td height="10"><br /></td> </tr><tr> <td height="10"><br /></td> </tr><tr> <td> <table border="0" cellpadding="0" cellspacing="0" width="358"> <tbody> <tr> <td width="10"><br /></td> <td valign="top" width="100"> <table cellpading="0" border="0" cellspacing="0"> <tbody> <tr> <td><a href="http://springer.r.delivery.net/r/r?2.1.Ee.2Tp.1iXXnb.BxM2Yw..H.GP%5fG.34q8.EaUEcM00"><img alt="Journal of Geodesy" src="http://a676.g.akamaitech.net/f/676/773/60m/images.delivery.net/cm50content/1364/55776/springerCA9HIQ33.jpg" border="0" width="95" height="126" /></a></td></tr> <tr> <td height="6"><br /></td> </tr><tr> <td><a href="http://springer.r.delivery.net/r/r?2.1.Ee.2Tp.1iXXnb.BxM2Yw..H.GP%5fG.34q8.EaUEcM00"><span originaltag="yes" style="font-family:Arial, Helvetica, sans-serif;font-size:85%;"><b><span originaltag="yes" style="font-family:Arial, Helvetica, sans-serif;font-size:85%;color:#0000cc;">Journal of Geodesy</span></b></span></a></td></tr></tbody></table></td> <td width="19"><br /></td> <td valign="top" width="100"> <table cellpading="0" border="0" cellspacing="0"> <tbody> <tr> <td><a href="http://springer.r.delivery.net/r/r?2.1.Ee.2Tp.1iXXnb.BxM2Yw..H.GP%5fI.34q8.FAIEcO00"><img alt="GPS Solutions" src="http://a676.g.akamaitech.net/f/676/773/60m/images.delivery.net/cm50content/1364/55776/springerCA25VCQ0.jpg" border="0" width="95" height="126" /></a></td></tr> <tr> <td height="6"><br /></td> </tr><tr> <td><a href="http://springer.r.delivery.net/r/r?2.1.Ee.2Tp.1iXXnb.BxM2Yw..H.GP%5fI.34q8.FAIEcO00"><span originaltag="yes" style="font-family:Arial, Helvetica, sans-serif;font-size:85%;"><b><span originaltag="yes" style="font-family:Arial, Helvetica, sans-serif;font-size:85%;color:#0000cc;">GPS Solutions</span></b></span></a></td></tr></tbody></table></td> <td width="19"><br /></td> <td valign="top" width="100"> <table cellpading="0" border="0" cellspacing="0"> <tbody> <tr> <td><a href="http://springer.r.delivery.net/r/r?2.1.Ee.2Tp.1iXXnb.BxM2Yw..H.GP%5fK.34q8.FFcEcQ00"><img alt="Bulletin of Earthquake Engineering" src="http://a676.g.akamaitech.net/f/676/773/60m/images.delivery.net/cm50content/1364/55776/springerCAJN26AH.jpg" border="0" width="95" height="142" /></a></td></tr> <tr> <td height="6"><br /></td> </tr><tr> <td><a href="http://springer.r.delivery.net/r/r?2.1.Ee.2Tp.1iXXnb.BxM2Yw..H.GP%5fK.34q8.FFcEcQ00"><span originaltag="yes" style="font-family:Arial, Helvetica, sans-serif;font-size:85%;"><b><span originaltag="yes" style="font-family:Arial, Helvetica, sans-serif;font-size:85%;color:#0000cc;">Bulletin of Earthquake Engineering</span></b></span></a></td></tr></tbody></table></td></tr></tbody></table></td></tr></tbody></table><div class="blogger-post-footer"><img width='1' height='1' src='https://blogger.googleusercontent.com/tracker/102059579936554555-3188200637843615898?l=yftian.blogspot.com' alt='' /></div>Enodhttp://www.blogger.com/profile/08433335265518211485noreply@blogger.com0tag:blogger.com,1999:blog-102059579936554555.post-87522054039910129632008-12-02T15:51:00.002+08:002008-12-02T16:03:05.241+08:00By default, est_noise6ac will take up about 1.8GB RAM to run. Although RAM is very cheap; however, not everyone had upgraded their hardware configuration. Fortunately, the memory usage can be greatly cut by small modifications to the raw source code.<br /><br />The below is the original variables declaration:<br /><span style="font-size:85%;"><span style="font-family: courier new;"> character*3 net,instname(30),netx</span><br /><span style="font-family: courier new;"> character*132 string,filename</span><br /><span style="font-family: courier new;"> character*8 bsta(82), line(82), header(9)</span><br /><span style="font-family: courier new;"> character*1 ans</span><br /><span style="font-family: courier new;"> integer instnum(30),insol(30),intot(30),infunc(30)</span><br /><span style="font-family: courier new;"> integer instart(30),instop(30)</span><br /><span style="font-family: courier new;"> real inparam(30,10)</span><br /><span style="font-family: courier new;"> character*7 choice(10),exp_choice(10),exp_type(10)</span><br /><span style="font-family: courier new;"> dimension ran1(65536),ran2(65536),ran3(65536),ran4(65536)</span><br /><span style="font-family: courier new;"> double precision t_start,t_stop,rat_chg1(82),per(82),off(82),</span><br /><span style="font-family: courier new;"> &amp; dest(82),timex,t_year(8103),texp(10),tpress,dtime(8103),</span><br /><span style="font-family: courier new;"> &amp; dtime_orig(8103),rat_chg2(82), dec_timed,fjul,fjul1,fjul2,</span><br /><span style="font-family: courier new;"> &amp; t(8103),rat_chng_norm(82)</span><br /><span style="font-family: courier new;"> dimension az(82),rc(20),d(8103),A(8103,82)</span><br /><span style="font-family: courier new;"> dimension d_orig(8103),A_orig(8103,82),t_orig(8103)</span><br /><span style="font-family: courier new;"> dimension res(8103),covinv(8103,8103),covar(8103,8103),x(82),</span><br /><span style="font-family: courier new;"> &amp; e(82),dith(7),xx(7),bexp(10)</span><br /><span style="font-family: courier new;"> dimension ambp(8,7),amby(8),acov(8,8),cov(8,8),wz(23)</span><br /><span style="font-family: courier new;"> dimension vsig(8),vamp1(8),vexp1(8),vamp_bp(8),valpha(8),vamp2(8),</span><br /><span style="font-family: courier new;"> &amp; vexp2(8)</span><br /><span style="font-family: courier new;"> dimension covarpl1(8103,8103),covarpl2(8103,8103),</span><br /><span style="font-family: courier new;"> &amp; covarbp(8103,8103)</span><br /><span style="font-family: courier new;"> dimension press(8103),aux_norm(82)</span></span><br /><br />Note that, the default size for most arrays is 8103. This threshold has exceeded the observation length of most CGPS site around the globe. For example, the site WUHN (1996.0697 - 2008.8374), a value of 5000 is suffice. Thus, we can modify the source to as follows:<br /><span style="font-size:85%;"><span style="font-family: courier new;"> integer*4 maxn</span><br /><span style="font-family: courier new;">c parameter(maxn=3999)</span><br /><span style="font-family: courier new;"> parameter(maxn=5000)</span><br /><span style="font-family: courier new;"> character*3 net,instname(30),netx</span><br /><span style="font-family: courier new;"> character*132 string,filename</span><br /><span style="font-family: courier new;"> character*8 bsta(82), line(82), header(9)</span><br /><span style="font-family: courier new;"> character*1 ans</span><br /><span style="font-family: courier new;"> integer instnum(30),insol(30),intot(30),infunc(30)</span><br /><span style="font-family: courier new;"> integer instart(30),instop(30)</span><br /><span style="font-family: courier new;"> real inparam(30,10)</span><br /><span style="font-family: courier new;"> character*7 choice(10),exp_choice(10),exp_type(10)</span><br /><span style="font-family: courier new;"> dimension ran1(65536),ran2(65536),ran3(65536),ran4(65536)</span><br /><span style="font-family: courier new;"> double precision t_start,t_stop,rat_chg1(82),per(82),off(82),</span><br /><span style="font-family: courier new;"> &amp; dest(82),timex,t_year(maxn),texp(10),tpress,dtime(maxn),</span><br /><span style="font-family: courier new;"> &amp; dtime_orig(maxn),rat_chg2(82), dec_timed,fjul,fjul1,fjul2,</span><br /><span style="font-family: courier new;"> &amp; t(maxn),rat_chng_norm(82)</span><br /><span style="font-family: courier new;"> dimension az(82),rc(20),d(maxn),A(maxn,82)</span><br /><span style="font-family: courier new;"> dimension d_orig(maxn),A_orig(maxn,82),t_orig(maxn)</span><br /><span style="font-family: courier new;"> dimension res(maxn),covinv(maxn,maxn),covar(maxn,maxn),x(82),</span><br /><span style="font-family: courier new;"> &amp; e(82),dith(7),xx(7),bexp(10)</span><br /><span style="font-family: courier new;"> dimension ambp(8,7),amby(8),acov(8,8),cov(8,8),wz(23)</span><br /><span style="font-family: courier new;"> dimension vsig(8),vamp1(8),vexp1(8),vamp_bp(8),valpha(8),vamp2(8),</span><br /><span style="font-family: courier new;"> &amp; vexp2(8)</span><br /><span style="font-family: courier new;"> dimension covarpl1(maxn,maxn),covarpl2(maxn,maxn),</span><br /><span style="font-family: courier new;"> &amp; covarbp(maxn,maxn)</span><br /><span style="font-family: courier new;"> dimension press(maxn),aux_norm(82)</span></span><br />This will reduce the memory required to run est_noise from 1.8GB to about 632MB.<br /><br />Warning: if you want to process sites with a long history, say from 1992 to present, you may need to set the threshold <span style="font-family: courier new;">(maxn)</span> larger.<div class="blogger-post-footer"><img width='1' height='1' src='https://blogger.googleusercontent.com/tracker/102059579936554555-8752205403991012963?l=yftian.blogspot.com' alt='' /></div>Enodhttp://www.blogger.com/profile/08433335265518211485noreply@blogger.com1tag:blogger.com,1999:blog-102059579936554555.post-48636325997614760902008-12-02T15:40:00.000+08:002008-12-02T15:41:43.328+08:00<span style="font-family: courier new;">f90 -f77 -ftrap=%none -O3 -xlic_lib=sunperf est_noise6ac.f -o est_noise6ac</span><br /><br />Sun Studio with Performance Libarary<br />Sun Solaris X86 10<div class="blogger-post-footer"><img width='1' height='1' src='https://blogger.googleusercontent.com/tracker/102059579936554555-4863632599761476090?l=yftian.blogspot.com' alt='' /></div>Enodhttp://www.blogger.com/profile/08433335265518211485noreply@blogger.com0tag:blogger.com,1999:blog-102059579936554555.post-29648138511074358072008-11-29T21:38:00.000+08:002008-11-29T21:39:07.923+08:00Program estimates the power spectral density (PSD)<br /> of data. PSD functions can be a combination of:<br /> 1) white noise<br /> 2) Simple power law noise P/f^n<br /> 3) a Gauss Markov version P/(fa^n + f^2)<br /> where fa=alpha/2*pi<br /> 4) A second power law<br /> 5) Band-pass filtered noise<br /> <br /> While estimating the PSD function, program also estimates<br /> various parameters that describe the time series including<br /> 1) DC term<br /> 2) rate<br /> 3) sinusoidal amplitudes of specified frequencies<br /> 4) rate changes<br /> 5) offset<br /> <br /> Program can handle data in various formats including<br /> 2-color EDM data<br /> <br /> Input the data type for automatic processing<br /> most users can will use otr otd or otx<br /> otr=data, such as inferred monument displacements<br /> format of yr, julian day, data, error bar<br /> otd=data, such as inferred monument displacements<br /> format of yrmoda, data, error bar<br /> otx=data, such as inferred monument displacements<br /> format of yr mo da, data, error bar<br /> input number of time series<br /> program will estimate only one set of PSD functions<br /> input the period of interest<br /> start and stop times; year julian_day year julian_day<br /> "julian day" may be decimal day<br /> Day numbers from 1960 13175.0000000000 17839.0000000000 <br /> <br /> Input the parameters of time series to be estimated<br /> <br /> Will rate be estimated? y/n<br /> Rate changes:<br /> time is the "hinge point"<br /> Input the number of rate changes in data<br /> <br /> Input the number of periodicities in the data<br /> Input the period (in days) of 1 period<br /> Input the period (in days) of 2 period<br /> <br /> Offset: specify first day after offset<br /> Input the number of offsets<br /> Input the time (year jul_day) or (year month day) of 1 offset<br /> Offset at 15375.6161000000 <br /> Input the time (year jul_day) or (year month day) of 2 offset<br /> Offset at 15656.9400000000 <br /> <br /> Input the number of exponentials in time series<br /> The number of exponential time constants that are fixed: 0<br /> rate renomalization is 6.384668 <br /> Input the format style of baseline data (otr, otx, otd, tmr, pkf)<br /> Input name of file for baseline number 1<br /> column of A matrix prior to input pressure data 8<br /> Number of files of Auxillary data (pressure)<br /> Number of data read is 4304<br /> Number of model parameters is 8<br /> <br /> Average Sampling interval: 1.083895 days<br /> Standard deviation of Ave. interval: 0.8538907 days<br /> Shortest Sampling interval: 1.000000 days<br /> Longest Sampling Interval: 29.00000 days<br /> <br /> Input the minimum sampling interval in days to use<br /> Number of points in time series for analysis is 4665<br /> <br /> Do you want to get rid of "redundent" data? y/n<br /> Coupled with minimum sampling interval, this can decimate data.<br /> Or, in the case of two baselines with measurements at same time,<br /> it can screw-up design matrix AND make the covariance singular for some cases<br /> After resampling data and getting rid of redundent data<br /> Number of observations is 4304<br /> <br /> time, data, and A matrix are in prob1(2).out<br /> prob1.out has stuff after redundencies are removed<br /> prob2.out has stuff in cronological order before redundancies removed<br /> <br /> Do you want to substitute real data with random numbers? y/n<br /> Estimate of white noise component of data is: 7.206384 <br /> methods to decimate data<br /> 0 = no decimation<br /> 1 = keep 2, skip 1, keep 2, skip 1<br /> 2 = keep 2, skip 1, keep 1, skip 2, keep 2, skip 1...<br /> 3 = keep 2, skip 1, keep 1, skip 2, keep 1, skip 3, keep 2 skip 1...<br /> 4 or more...make it option 3<br /> input choice<br /> Eigenvalues 1 127.200546298917 2<br /> 153.657986228781 3 1221.91656232437 4<br /> 2113.71630746652 5 2130.43663134457 6<br /> 2170.10470602075 7 2200.15505332048 8<br /> 12557.3695800425 <br /> Using 8 out of 8 eigenvalues<br /> RMS fit using white noise model is 9.792912 <br /> Log MLE for white noise error model is -15923.37 <br /> <br /> Nomimal value for baseline 1 -28.64 +/- 0.04<br /> Rate in units per year -0.1596 +/- 0.0095<br /> Period of 365.250 days, cos amp= -3.49 +/- 0.02 sin amp= 3.26 +/- 0.02 magnitude= 4.78 +/- 0.02<br /> Period of 182.625 days, cos amp= -0.44 +/- 0.02 sin amp= -0.40 +/- 0.02 magnitude= 0.59 +/- 0.02<br /> Offset number 1 at 2002 34.616 is 40.23 +/- 0.07<br /> Offset number 2 at 2002 315.940 is 6.90 +/- 0.07<br /> sampling interval in yrs 2.737851E-03<br /> <br /> Input the initial parameters of the PSD and whether the item is "fix" or "f<br /> loat"<br /> if "fix", then the item is not estimated<br /> if "float", then the item is estimated<br /> <br /> Input the white noise "instrument precision" and fix/float<br /> Input the amplitude first Power law function and fix/float<br /> Input the exponent 1 < n < 3 and fix/float<br /> Input the time constant alpha in c/yr and fix/float<br /> <br /> Input the parameters for band-passed filtered noise<br /> Input low and high freq stop band in c/yr<br /> low frequency stop band is 0.5000000 c/yr<br /> high frequency stop band is 2.000000 c/yr<br /> number of poles between 1 and 4<br /> Input the amplitude and fix/float<br /> Input the exponent of second Power Law function fix/float<br /> Input the amplitude of second PL function fix/float<br /> <br /> Sometimes, it may be necessary to add white noise to data so<br /> that a better estimate of long period PSD parameters can be made.<br /> This is especially true for data that is predominantly power noise<br /> Enter value of white noise to be added (nominal it should be 0)<br /> Calculating power law covariance for first set<br /> Calculating power law covariance for second set<br /> <br /> list of trial covariance parameters<br /> white noise PL_1 amp PL_1 exp GM freq BP amp PL_2 amp PL_2 exp determinant chi^2 MLE cpu<br /> 1.0000 1.0000 1.0000 0.0000 0.0000 0.0000 2.0000 0.193E+03 292980.438 -150638.297 72.680<br /> 0.7000 1.0000 1.0000 0.0000 0.0000 0.0000 2.0000 -0.119E+04 531153.438 -268338.031 65.760<br /> 1.0000 1.9500 1.0000 0.0000 0.0000 0.0000 2.0000 0.551E+03 228566.516 -118789.180 65.550<br /> Initial solutions for Amoeba<br /> 150638.3 1.000000 1.000000 1.000000 0.000000E+00<br /> 0.000000E+00 0.000000E+00 2.000000 <br /> 268338.0 0.7000000 1.000000 1.000000 0.000000E+00<br /> 0.000000E+00 0.000000E+00 2.000000 <br /> 118789.2 1.000000 1.950000 1.000000 0.000000E+00<br /> 0.000000E+00 0.000000E+00 2.000000 <br /> 150638.3 1.000000 1.000000 1.000000 0.000000E+00<br /> 0.000000E+00 0.000000E+00 2.000000 <br /> 150638.3 1.000000 1.000000 1.000000 0.000000E+00<br /> 0.000000E+00 0.000000E+00 2.000000 <br /> 150638.3 1.000000 1.000000 1.000000 0.000000E+00<br /> 0.000000E+00 0.000000E+00 2.000000 <br /> 150638.3 1.000000 1.000000 1.000000 0.000000E+00<br /> 0.000000E+00 0.000000E+00 2.000000 <br /> 150638.3 1.000000 1.000000 1.000000 0.000000E+00<br /> 0.000000E+00 0.000000E+00 2.000000 <br /> 1.3000 1.2714 1.0000 0.0000 0.0000 0.0000 2.0000 0.132E+04 174504.484 -92522.672 65.520<br /> 1.6000 1.4071 1.0000 0.0000 0.0000 0.0000 2.0000 0.218E+04 118746.250 -65507.523 65.440<br /> 1.1714 1.3878 1.0000 0.0000 0.0000 0.0000 2.0000 0.935E+03 202482.875 -106131.211 65.520<br /> 1.2204 1.4985 1.0000 0.0000 0.0000 0.0000 2.0000 0.113E+04 184345.734 -97253.906 65.590<br /> 1.2834 1.6410 1.0000 0.0000 0.0000 0.0000 2.0000 0.136E+04 164415.156 -87523.047 65.570<br /> 1.3643 1.8241 1.0000 0.0000 0.0000 0.0000 2.0000 0.164E+04 143224.766 -77212.164 65.400<br /> 1.4684 2.0596 1.0000 0.0000 0.0000 0.0000 2.0000 0.199E+04 121519.352 -66700.039 65.410<br /> 1.6023 2.3623 1.0000 0.0000 0.0000 0.0000 2.0000 0.239E+04 100181.875 -56433.719 65.510<br /> 1.9034 3.0435 1.0000 0.0000 0.0000 0.0000 2.0000 0.318E+04 68773.023 -41518.789 65.630<br /> 1.8604 1.7247 1.0000 0.0000 0.0000 0.0000 2.0000 0.284E+04 86536.969 -50066.039 65.500<br /> 1.8858 2.3833 1.0000 0.0000 0.0000 0.0000 2.0000 0.301E+04 76452.273 -45192.828 65.400<br /> 2.0270 2.5253 1.0000 0.0000 0.0000 0.0000 2.0000 0.332E+04 66499.344 -40520.672 65.510<br /> 2.4303 3.0386 1.0000 0.0000 0.0000 0.0000 2.0000 0.410E+04 46203.660 -31155.395 65.520<br /> 2.2917 2.7822 1.0000 0.0000 0.0000 0.0000 2.0000 0.383E+04 52479.074 -34027.707 65.560<br /> 2.4757 2.8727 1.0000 0.0000 0.0000 0.0000 2.0000 0.415E+04 45642.617 -30923.447 65.500<br /> 3.0313 3.3970 1.0000 0.0000 0.0000 0.0000 2.0000 0.501E+04 30785.910 -24353.248 65.520<br /> 2.8181 3.0194 1.0000 0.0000 0.0000 0.0000 2.0000 0.468E+04 36126.281 -26693.301 65.440<br /> 3.0346 4.1325 1.0000 0.0000 0.0000 0.0000 2.0000 0.509E+04 28761.488 -23430.217 65.430<br /> 3.7519 5.4952 1.0000 0.0000 0.0000 0.0000 2.0000 0.605E+04 18316.992 -19159.568 65.450<br /> 3.3146 4.8922 1.0000 0.0000 0.0000 0.0000 2.0000 0.552E+04 23400.576 -21172.289 65.440<br /> 3.6974 4.9505 1.0000 0.0000 0.0000 0.0000 2.0000 0.594E+04 19491.828 -19637.287 65.400<br /> 4.1924 4.8351 1.0000 0.0000 0.0000 0.0000 2.0000 0.641E+04 15947.724 -18340.754 65.420<br /> 5.3369 5.7309 1.0000 0.0000 0.0000 0.0000 2.0000 0.742E+04 10066.285 -16412.479 65.450<br /> 4.6742 5.9389 1.0000 0.0000 0.0000 0.0000 2.0000 0.692E+04 12421.931 -17086.725 65.450<br /> 5.1767 6.5111 1.0000 0.0000 0.0000 0.0000 2.0000 0.736E+04 10162.461 -16391.967 65.400<br /> 6.5499 8.2473 1.0000 0.0000 0.0000 0.0000 2.0000 0.837E+04 6345.566 -15496.672 65.450<br /> 5.8551 8.0240 1.0000 0.0000 0.0000 0.0000 2.0000 0.793E+04 7708.284 -15735.510 65.890<br /> 6.4487 8.9684 1.0000 0.0000 0.0000 0.0000 2.0000 0.835E+04 6320.756 -15464.832 65.400<br /> 8.1573 11.7540 1.0000 0.0000 0.0000 0.0000 2.0000 0.938E+04 3898.782 -15284.209 65.370<br /> 7.5490 9.4338 1.0000 0.0000 0.0000 0.0000 2.0000 0.898E+04 4789.413 -15326.256 65.400<br /> 8.2667 10.6565 1.0000 0.0000 0.0000 0.0000 2.0000 0.938E+04 3951.573 -15312.186 65.420<br /> 9.5021 11.5864 1.0000 0.0000 0.0000 0.0000 2.0000 0.996E+04 3048.049 -15435.300 65.390<br /> 9.9592 12.7562 1.0000 0.0000 0.0000 0.0000 2.0000 0.102E+05 2728.618 -15499.472 65.380<br /> 10.6172 14.9715 1.0000 0.0000 0.0000 0.0000 2.0000 0.105E+05 2319.967 -15617.645 65.410<br /> 11.4596 14.6633 1.0000 0.0000 0.0000 0.0000 2.0000 0.108E+05 2061.600 -15769.482 65.460<br /> 7.2562 9.6839 1.0000 0.0000 0.0000 0.0000 2.0000 0.884E+04 5066.731 -15325.032 65.530<br /> 5.7372 6.2051 1.0000 0.0000 0.0000 0.0000 2.0000 0.774E+04 8691.115 -16038.730 65.430<br /> 9.3972 12.7799 1.0000 0.0000 0.0000 0.0000 2.0000 0.996E+04 3000.721 -15414.119 65.510<br /> 6.2346 8.4272 1.0000 0.0000 0.0000 0.0000 2.0000 0.819E+04 6832.184 -15560.886 65.510<br /> 9.0281 11.6739 1.0000 0.0000 0.0000 0.0000 2.0000 0.976E+04 3309.592 -15371.841 65.550<br /> 10.3520 13.9151 1.0000 0.0000 0.0000 0.0000 2.0000 0.104E+05 2483.092 -15566.000 65.460<br /> 7.5004 9.6643 1.0000 0.0000 0.0000 0.0000 2.0000 0.896E+04 4800.946 -15318.011 65.430<br /> 6.8279 10.0268 1.0000 0.0000 0.0000 0.0000 2.0000 0.862E+04 5525.248 -15342.189 65.450<br /> 6.1987 8.0467 1.0000 0.0000 0.0000 0.0000 2.0000 0.815E+04 7010.972 -15606.014 65.460<br /> 8.5976 11.5966 1.0000 0.0000 0.0000 0.0000 2.0000 0.957E+04 3595.442 -15324.635 65.480<br /> 6.4448 9.1306 1.0000 0.0000 0.0000 0.0000 2.0000 0.836E+04 6285.427 -15454.269 65.460<br /> 8.3823 11.0381 1.0000 0.0000 0.0000 0.0000 2.0000 0.945E+04 3814.846 -15313.561 65.430<br /> 9.0891 11.0667 1.0000 0.0000 0.0000 0.0000 2.0000 0.976E+04 3332.963 -15385.887 65.420<br /> 7.3932 10.2868 1.0000 0.0000 0.0000 0.0000 2.0000 0.894E+04 4808.054 -15297.027 65.460<br /> 8.3235 11.9034 1.0000 0.0000 0.0000 0.0000 2.0000 0.946E+04 3755.213 -15295.147 65.370<br /> 8.9212 12.2875 1.0000 0.0000 0.0000 0.0000 2.0000 0.974E+04 3314.293 -15353.544 65.540<br /> 7.6725 10.3348 1.0000 0.0000 0.0000 0.0000 2.0000 0.908E+04 4516.957 -15294.742 65.410<br /> 7.3155 10.0142 1.0000 0.0000 0.0000 0.0000 2.0000 0.888E+04 4939.837 -15309.153 65.460<br /> 8.3598 12.0465 1.0000 0.0000 0.0000 0.0000 2.0000 0.949E+04 3712.112 -15296.443 65.510<br /> Amoeba exceeding maximum iterations.<br /> <br /> Best fitting solutions<br /> MLE= -15284.21 <br /> white noise= 8.157310 <br /> Bandpass filter amplitude= 0.000000E+00<br /> power law noise 1<br /> amplitude= 11.75405 <br /> exponent= 1.000000 <br /> G-M freq= 0.000000E+00<br /> power law noise 2<br /> amplitude= 0.000000E+00<br /> exponent= 2.000000 <br /> <br /> 9.7888 11.7540 1.0000 0.0000 0.0000 0.0000 2.0000 0.101E+05 2886.873 -15476.197 65.460<br /> 8.1573 22.9204 1.0000 0.0000 0.0000 0.0000 2.0000 0.997E+04 2808.724 -15330.854 65.510<br /> Initial solutions for Amoeba<br /> 15284.21 8.157310 11.75405 1.000000 0.000000E+00<br /> 0.000000E+00 0.000000E+00 2.000000 <br /> 15476.20 9.788773 11.75405 1.000000 0.000000E+00<br /> 0.000000E+00 0.000000E+00 2.000000 <br /> 15330.85 8.157310 22.92039 1.000000 0.000000E+00<br /> 0.000000E+00 0.000000E+00 2.000000 <br /> 15284.21 8.157310 11.75405 1.000000 0.000000E+00<br /> 0.000000E+00 0.000000E+00 2.000000 <br /> 15284.21 8.157310 11.75405 1.000000 0.000000E+00<br /> 0.000000E+00 0.000000E+00 2.000000 <br /> 15284.21 8.157310 11.75405 1.000000 0.000000E+00<br /> 0.000000E+00 0.000000E+00 2.000000 <br /> 15284.21 8.157310 11.75405 1.000000 0.000000E+00<br /> 0.000000E+00 0.000000E+00 2.000000 <br /> 15284.21 8.157310 11.75405 1.000000 0.000000E+00<br /> 0.000000E+00 0.000000E+00 2.000000 <br /> 6.5258 14.9444 1.0000 0.0000 0.0000 0.0000 2.0000 0.877E+04 4950.469 -15202.870 65.590<br /> 4.8944 16.5396 1.0000 0.0000 0.0000 0.0000 2.0000 0.804E+04 6856.128 -15426.552 65.530<br /> 7.6912 1.4992 1.0000 0.0000 0.0000 0.0000 2.0000 0.879E+04 6659.046 -16076.798 65.520<br /> 8.0408 17.5651 1.0000 0.0000 0.0000 0.0000 2.0000 0.962E+04 3342.892 -15250.721 65.570<br /> 7.6579 14.3259 1.0000 0.0000 0.0000 0.0000 2.0000 0.928E+04 3973.027 -15220.808 65.480<br /> 7.5152 15.0607 1.0000 0.0000 0.0000 0.0000 2.0000 0.925E+04 3994.780 -15207.222 65.440<br /> 7.3317 16.0055 1.0000 0.0000 0.0000 0.0000 2.0000 0.923E+04 4022.129 -15192.372 65.470<br /> 6.9189 18.1312 1.0000 0.0000 0.0000 0.0000 2.0000 0.917E+04 4077.885 -15168.667 65.450<br /> 6.9779 17.8275 1.0000 0.0000 0.0000 0.0000 2.0000 0.918E+04 4070.602 -15171.293 65.450<br /> 6.6409 19.5628 1.0000 0.0000 0.0000 0.0000 2.0000 0.915E+04 4108.121 -15159.361 65.420<br /> 5.8827 23.4671 1.0000 0.0000 0.0000 0.0000 2.0000 0.913E+04 4143.565 -15155.591 65.450<br /> 5.9911 22.9093 1.0000 0.0000 0.0000 0.0000 2.0000 0.913E+04 4143.586 -15154.493 65.590<br /> 4.9079 28.4870 1.0000 0.0000 0.0000 0.0000 2.0000 0.920E+04 4049.646 -15183.117 65.510<br /> 5.5220 18.6252 1.0000 0.0000 0.0000 0.0000 2.0000 0.856E+04 5394.006 -15211.683 65.460<br /> 5.2946 23.0928 1.0000 0.0000 0.0000 0.0000 2.0000 0.885E+04 4715.053 -15166.986 65.440<br /> 7.3655 20.0699 1.0000 0.0000 0.0000 0.0000 2.0000 0.949E+04 3512.877 -15203.794 65.450<br /> 5.3296 25.0657 1.0000 0.0000 0.0000 0.0000 2.0000 0.905E+04 4322.794 -15162.197 65.420<br /> 4.8975 21.4838 1.0000 0.0000 0.0000 0.0000 2.0000 0.853E+04 5480.706 -15226.291 65.470<br /> 6.7485 20.4234 1.0000 0.0000 0.0000 0.0000 2.0000 0.926E+04 3905.840 -15165.524 65.450<br /> 5.8008 28.1747 1.0000 0.0000 0.0000 0.0000 2.0000 0.948E+04 3531.001 -15203.898 65.530<br /> 6.3446 18.2520 1.0000 0.0000 0.0000 0.0000 2.0000 0.892E+04 4572.325 -15162.458 65.520<br /> 5.1678 25.4129 1.0000 0.0000 0.0000 0.0000 2.0000 0.901E+04 4392.139 -15165.720 65.530<br /> 4.7265 27.1898 1.0000 0.0000 0.0000 0.0000 2.0000 0.902E+04 4412.730 -15179.054 65.560<br /> 6.3708 20.3958 1.0000 0.0000 0.0000 0.0000 2.0000 0.909E+04 4212.109 -15154.090 65.540<br /> 6.6583 21.4575 1.0000 0.0000 0.0000 0.0000 2.0000 0.929E+04 3837.766 -15167.026 65.510<br /> 5.6355 22.6840 1.0000 0.0000 0.0000 0.0000 2.0000 0.896E+04 4479.419 -15156.088 65.520<br /> 6.9187 18.3577 1.0000 0.0000 0.0000 0.0000 2.0000 0.919E+04 4047.370 -15168.460 65.440<br /> 5.6055 23.6491 1.0000 0.0000 0.0000 0.0000 2.0000 0.903E+04 4339.257 -15156.240 65.490<br /> 5.0114 24.2689 1.0000 0.0000 0.0000 0.0000 2.0000 0.885E+04 4745.451 -15174.772 65.490<br /> 6.3143 21.3848 1.0000 0.0000 0.0000 0.0000 2.0000 0.914E+04 4112.492 -15154.622 65.500<br /> 5.4067 27.3354 1.0000 0.0000 0.0000 0.0000 2.0000 0.927E+04 3896.867 -15177.643 65.450<br /> 6.1101 20.5228 1.0000 0.0000 0.0000 0.0000 2.0000 0.899E+04 4419.531 -15153.302 65.400<br /> 6.6447 19.2238 1.0000 0.0000 0.0000 0.0000 2.0000 0.913E+04 4153.084 -15159.452 65.400<br /> 6.3159 20.6842 1.0000 0.0000 0.0000 0.0000 2.0000 0.909E+04 4215.089 -15153.476 65.460<br /> 6.5717 19.7932 1.0000 0.0000 0.0000 0.0000 2.0000 0.914E+04 4132.414 -15157.688 65.460<br /> 5.8471 22.6851 1.0000 0.0000 0.0000 0.0000 2.0000 0.905E+04 4296.753 -15153.289 65.420<br /> 6.6022 20.7586 1.0000 0.0000 0.0000 0.0000 2.0000 0.922E+04 3973.179 -15160.868 65.420<br /> 5.8772 22.2027 1.0000 0.0000 0.0000 0.0000 2.0000 0.902E+04 4349.986 -15152.851 65.400<br /> 6.3534 19.6143 1.0000 0.0000 0.0000 0.0000 2.0000 0.903E+04 4346.274 -15154.950 65.420<br /> 6.2357 20.5775 1.0000 0.0000 0.0000 0.0000 2.0000 0.905E+04 4299.920 -15152.851 65.410<br /> 5.8994 21.4659 1.0000 0.0000 0.0000 0.0000 2.0000 0.897E+04 4452.893 -15153.516 65.500<br /> 6.1965 19.5290 1.0000 0.0000 0.0000 0.0000 2.0000 0.895E+04 4502.954 -15156.316 65.580<br /> 6.0424 22.0642 1.0000 0.0000 0.0000 0.0000 2.0000 0.908E+04 4232.460 -15152.586 65.530<br /> 5.7229 22.5192 1.0000 0.0000 0.0000 0.0000 2.0000 0.898E+04 4430.932 -15154.572 65.550<br /> 6.2088 20.9267 1.0000 0.0000 0.0000 0.0000 2.0000 0.906E+04 4268.928 -15152.562 65.570<br /> 6.2826 21.2950 1.0000 0.0000 0.0000 0.0000 2.0000 0.912E+04 4151.414 -15153.802 65.560<br /> 5.9952 21.4232 1.0000 0.0000 0.0000 0.0000 2.0000 0.901E+04 4374.979 -15152.454 65.550<br /> 5.7745 22.2878 1.0000 0.0000 0.0000 0.0000 2.0000 0.899E+04 4425.007 -15153.998 65.530<br /> 6.1806 21.0851 1.0000 0.0000 0.0000 0.0000 2.0000 0.906E+04 4268.485 -15152.417 65.500<br /> 6.0005 22.6099 1.0000 0.0000 0.0000 0.0000 2.0000 0.911E+04 4182.114 -15153.563 65.540<br /> 6.0827 21.0446 1.0000 0.0000 0.0000 0.0000 2.0000 0.902E+04 4359.525 -15152.402 65.550<br /> 6.3308 19.9788 1.0000 0.0000 0.0000 0.0000 2.0000 0.904E+04 4309.770 -15154.043 65.520<br /> 5.9680 22.0086 1.0000 0.0000 0.0000 0.0000 2.0000 0.905E+04 4303.811 -15152.383 65.620<br /> 6.3267 20.4059 1.0000 0.0000 0.0000 0.0000 2.0000 0.907E+04 4248.056 -15153.568 65.560<br /> 5.9896 21.7535 1.0000 0.0000 0.0000 0.0000 2.0000 0.903E+04 4326.403 -15152.272 65.540<br /> Amoeba exceeding maximum iterations.<br /> <br /> Best fitting solutions<br /> MLE= -15152.27 <br /> white noise= 5.989558 <br /> Bandpass filter amplitude= 0.000000E+00<br /> power law noise 1<br /> amplitude= 21.75346 <br /> exponent= 1.000000 <br /> G-M freq= 0.000000E+00<br /> power law noise 2<br /> amplitude= 0.000000E+00<br /> exponent= 2.000000 <br /> <br /> 5.9896 21.7535 1.0000 0.0000 0.0000 0.0000 2.0000 0.903E+04 4326.403 -15152.272 65.600<br /> Start the covariance calculations for noise model<br /> best estimate 5.989558 dither 1.000000E-02<br /> best estimate 21.75346 dither 1.000000E-02<br /> best estimate 1.000000 dither 5.000000E-02<br /> best estimate 0.000000E+00 dither 5.000000E-02<br /> best estimate 0.000000E+00 dither 5.000000E-02<br /> best estimate 0.000000E+00 dither 5.000000E-02<br /> best estimate 2.000000 dither 0.1000000 <br /> 5.9297 21.7535 1.0000 0.0000 0.0000 0.0000 2.0000 0.901E+04 4378.312 -15152.646 65.550<br /> Dither changed to 6.000000E-03<br /> 5.9536 21.7535 1.0000 0.0000 0.0000 0.0000 2.0000 0.902E+04 4357.449 -15152.455 65.590<br /> Dither changed to 3.600000E-03<br /> 5.9680 21.7535 1.0000 0.0000 0.0000 0.0000 2.0000 0.902E+04 4344.996 -15152.368 65.550<br /> Dither changed to 2.160000E-03<br /> 5.9766 21.7535 1.0000 0.0000 0.0000 0.0000 2.0000 0.903E+04 4337.546 -15152.324 65.570<br /> Dither changed to 1.296000E-03<br /> 5.9818 21.7535 1.0000 0.0000 0.0000 0.0000 2.0000 0.903E+04 4333.084 -15152.302 65.580<br /> Dither changed to 7.776001E-04<br /> 5.9849 21.7535 1.0000 0.0000 0.0000 0.0000 2.0000 0.903E+04 4330.411 -15152.289 65.560<br /> 5.9942 21.7535 1.0000 0.0000 0.0000 0.0000 2.0000 0.904E+04 4322.401 -15152.258 65.490<br /> 5.9896 21.5359 1.0000 0.0000 0.0000 0.0000 2.0000 0.902E+04 4361.599 -15152.377 65.610<br /> Dither changed to 6.000000E-03<br /> 5.9896 21.6229 1.0000 0.0000 0.0000 0.0000 2.0000 0.902E+04 4347.481 -15152.315 65.570<br /> Dither changed to 3.600000E-03<br /> 5.9896 21.6751 1.0000 0.0000 0.0000 0.0000 2.0000 0.903E+04 4339.036 -15152.291 65.550<br /> 5.9896 21.8318 1.0000 0.0000 0.0000 0.0000 2.0000 0.904E+04 4313.813 -15152.276 65.570<br /> 5.9942 21.8318 1.0000 0.0000 0.0000 0.0000 2.0000 0.904E+04 4309.832 -15152.267 65.570<br /> 5.9849 21.8318 1.0000 0.0000 0.0000 0.0000 2.0000 0.904E+04 4317.799 -15152.288 65.590<br /> 5.9942 21.6751 1.0000 0.0000 0.0000 0.0000 2.0000 0.903E+04 4335.014 -15152.271 65.590<br /> 5.9849 21.6751 1.0000 0.0000 0.0000 0.0000 2.0000 0.903E+04 4343.063 -15152.313 65.550<br /> <br /> Inverse covariance matrix<br /> 1 acov= 90.03841 14.05651 0.000000E+00 0.000000E+00<br /> 0.000000E+00 0.000000E+00 0.000000E+00<br /> 2 acov= 14.05651 3.662404 0.000000E+00 0.000000E+00<br /> 0.000000E+00 0.000000E+00 0.000000E+00<br /> 3 acov= 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00<br /> 0.000000E+00 0.000000E+00 0.000000E+00<br /> 4 acov= 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00<br /> 0.000000E+00 0.000000E+00 0.000000E+00<br /> 5 acov= 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00<br /> 0.000000E+00 0.000000E+00 0.000000E+00<br /> 6 acov= 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00<br /> 0.000000E+00 0.000000E+00 0.000000E+00<br /> 7 acov= 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00<br /> 0.000000E+00 0.000000E+00 0.000000E+00<br /> The eigenvalues 1.432470 92.26835 <br /> ier= 0<br /> <br /> the covariance matrix<br /> 1 cov= 2.770945E-02 -0.1063504 0.000000E+00 0.000000E+00<br /> 0.000000E+00 0.000000E+00 0.000000E+00<br /> 2 cov= -0.1063504 0.6812234 0.000000E+00 0.000000E+00<br /> 0.000000E+00 0.000000E+00 0.000000E+00<br /> 3 cov= 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00<br /> 0.000000E+00 0.000000E+00 0.000000E+00<br /> 4 cov= 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00<br /> 0.000000E+00 0.000000E+00 0.000000E+00<br /> 5 cov= 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00<br /> 0.000000E+00 0.000000E+00 0.000000E+00<br /> 6 cov= 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00<br /> 0.000000E+00 0.000000E+00 0.000000E+00<br /> 7 cov= 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00<br /> 0.000000E+00 0.000000E+00 0.000000E+00<br /> <br /> Cross correlation matrix<br /> 1 cross correlation 1.000000 -0.7740703 0.000000E+00<br /> 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00<br /> 2 cross correlation -0.7740703 1.000000 0.000000E+00<br /> 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00<br /> 3 cross correlation 0.000000E+00 0.000000E+00 0.000000E+00<br /> 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00<br /> 4 cross correlation 0.000000E+00 0.000000E+00 0.000000E+00<br /> 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00<br /> 5 cross correlation 0.000000E+00 0.000000E+00 0.000000E+00<br /> 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00<br /> 15152.27 <br /> <br /> <br /> row number for optimal solution is 115<br /> nmod= 8<br /> -32.07200 0.000000E+00 9.846000 0.000000E+00 -3.540000 <br /> 0.000000E+00 3.686000 0.000000E+00 -0.4670000 0.000000E+00<br /> -0.3010000 0.000000E+00 26.66000 0.000000E+00 5.183000 <br /> 0.000000E+00<br /> Eigenvalues 1 3.74611910120969D-002 2<br /> 7.17193826658182D-002 3 0.127016565108053 4<br /> 0.251911574073773 5 1.57312529521336 6<br /> 1.62081302617026 7 3.01568244130233 8<br /> 3.10585976417084 <br /> Using 8 out of 8 eigenvalues<br /> 5.9896 21.7535 1.0000 0.0000 0.0000 0.0000 2.0000 0.903E+04 4326.403 -15152.272 68.960<br /> number of rows is 117<br /> -32.07200 2.842000 9.846000 4.446000 -3.540000 <br /> 0.8020000 3.686000 0.8040000 -0.4670000 0.5710000 <br /> -0.3010000 0.5790000 26.66000 3.511000 5.183000 <br /> 3.502000 <br /> Residual, decimated data in resid_dec.out<br /> col 1 &amp; 2, time; 3 is residual,4 is calculated, 5 is data<br /> Residual, data in resid.out<br /> col 1 &amp; 2, time; 3 is residual,4 is calculated, 5 is data<br /> <br /> <br /> Nomimal value for baseline 1 -32.07 +/- 2.84<br /><span style="font-weight: bold;"> Rate in units per year 1.5421 +/- 0.6964</span><br /><span style="font-weight: bold;"> Period of 365.250 days, cos amp= -3.54 +/- 0.80 sin amp= 3.69 +/- 0.80 magnitude= 5.11 +/- 0.80</span><br /> Period of 182.625 days, cos amp= -0.47 +/- 0.57 sin amp= -0.30 +/- 0.58 magnitude= 0.56 +/- 0.57<br /><span style="font-weight: bold;"> Offset number 1 at 2002 34.616 is 26.66 +/- 3.51</span><br /> <span style="font-weight: bold;">Offset number 2 at 2002 315.940 is 5.18 +/- 3.50</span><br /> <br /> <br /> Best fitting solutions<br /> MLE= -15152.27 <br /> <span style="font-weight: bold;">white noise= 5.989558 +/- 0.1664616 </span><br /> Bandpass filter amplitude= 0.000000E+00 +/- 0.000000E+00<br /> power law noise 1<br /> <span style="font-weight: bold;">amplitude= 21.75346 +/- 0.8253626 </span> <br /> exponent= 1.000000 +/- 0.000000E+00<br /> G-M freq= 0.000000E+00 +/- 0.000000E+00<br /> power law noise 2<br /> amplitude= 0.000000E+00 +/- 0.000000E+00<br /> exponent= 2.000000 +/- 0.000000E+00<div class="blogger-post-footer"><img width='1' height='1' src='https://blogger.googleusercontent.com/tracker/102059579936554555-2964813851107435807?l=yftian.blogspot.com' alt='' /></div>Enodhttp://www.blogger.com/profile/08433335265518211485noreply@blogger.com2tag:blogger.com,1999:blog-102059579936554555.post-58191780471626427422008-11-29T21:28:00.002+08:002008-11-29T21:32:46.550+08:00I finally found what caused an error in time conversion in my code. When calculate factorial day I miss used seconds in one hour (3600) to divide the total seconds past since the begin of the day, instead of using 24*3600. Thus, the epochs definition for offsets are all erroneous. I have to rerun the code from the beginning.<br />I lost two weeks.<div class="blogger-post-footer"><img width='1' height='1' src='https://blogger.googleusercontent.com/tracker/102059579936554555-5819178047162642742?l=yftian.blogspot.com' alt='' /></div>Enodhttp://www.blogger.com/profile/08433335265518211485noreply@blogger.com0tag:blogger.com,1999:blog-102059579936554555.post-58733953512617012122008-11-22T11:28:00.000+08:002008-11-22T11:29:51.743+08:00?<div class="blogger-post-footer"><img width='1' height='1' src='https://blogger.googleusercontent.com/tracker/102059579936554555-5873395351261701212?l=yftian.blogspot.com' alt='' /></div>Enodhttp://www.blogger.com/profile/08433335265518211485noreply@blogger.com0tag:blogger.com,1999:blog-102059579936554555.post-42171958013230547592008-11-17T21:23:00.001+08:002008-11-17T21:25:01.684+08:00<span style="font-family: courier new;">[tianyf@tyf ~/tmp]$ ifort -cm -w95 -O3 -ip -ipo_obj -mp1 -static -o est_noise6ab est_noise6ac.f -L/opt/intel/mkl/10.1.0.015/lib/32 -lmkl_lapack -lmkl_ia32 -lguide -lpthread</span><br /><span style="font-family: courier new;">ifort: command line remark #10148: option '-ipo_obj' not supported</span><br /><span style="font-family: courier new;">/opt/intel/mkl/10.1.0.015/lib/32/libguide.a(z_Linux_util.o)(.text+0xa64): In function `__kmp_register_atfork':</span><br /><span style="font-family: courier new;">: undefined reference to `pthread_atfork'</span><br /><br />Where is the problem?<div class="blogger-post-footer"><img width='1' height='1' src='https://blogger.googleusercontent.com/tracker/102059579936554555-4217195801323054759?l=yftian.blogspot.com' alt='' /></div>Enodhttp://www.blogger.com/profile/08433335265518211485noreply@blogger.com0tag:blogger.com,1999:blog-102059579936554555.post-60143635892278063132008-11-17T20:16:00.002+08:002008-11-17T20:20:15.869+08:00<a href="http://www.intel.com/support/performancetools/sb/CS-027812.htm">http://www.intel.com/support/performancetools/sb/CS-027812.htm</a><br /><br />Switch to bash, and set the environment variables as follows.<br /><span style="font-size:85%;"><br /><span style="font-family: courier new;">export CC=icc</span><br /><span style="font-family: courier new;">export CXX=icpc</span><br /><span style="font-family: courier new;">export CFLAGS='-O3 -xT -ip -no-prec-div -static'</span><br /><span style="font-family: courier new;">export CXXFLAGS='-O3 -xT -ip -no-prec-div -static'</span><br /><span style="font-family: courier new;">export F77=ifort</span><br /><span style="font-family: courier new;">export FFLAGS='-O3 -xT -ip -no-prec-div -static'</span><br /><span style="font-family: courier new;">export CPP='icc -E'</span><br /><span style="font-family: courier new;">export CXXCPP='icpc -E'</span></span><br /><br />Configure the source, and build the library.<br /><span style="font-size:85%;"><span style="font-family: courier new;">./configure --prefix=/usr/local/netcdfi</span><br /><span style="font-family: courier new;">make</span><br /><span style="font-family: courier new;">make check</span><br /><span style="font-family: courier new;">make install</span></span><br /><br />You must be root to install NetCDF to /usr/local directory.<div class="blogger-post-footer"><img width='1' height='1' src='https://blogger.googleusercontent.com/tracker/102059579936554555-6014363589227806313?l=yftian.blogspot.com' alt='' /></div>Enodhttp://www.blogger.com/profile/08433335265518211485noreply@blogger.com1tag:blogger.com,1999:blog-102059579936554555.post-88070669370241230712008-11-16T19:02:00.002+08:002008-11-16T19:03:19.949+08:00I didn't set the man path for my GPSF programs, but it was appended to the system MANPATH variable automatically. Why?<br /><br />[tianyf@tyf ~]$ manpath<br />/usr/local/GMT4.3.1/GMT4.3.1/man:/opt/intel/mkl/10.1.0.015/man:/opt/intel/cc/10.1.021/man:/opt/intel/fc/10.1.021/man:/home/tianyf/gpsf/cgps/man:/usr/share/man:/usr/X11R6/man:/usr/local/share/man<div class="blogger-post-footer"><img width='1' height='1' src='https://blogger.googleusercontent.com/tracker/102059579936554555-8807066937024123071?l=yftian.blogspot.com' alt='' /></div>Enodhttp://www.blogger.com/profile/08433335265518211485noreply@blogger.com0tag:blogger.com,1999:blog-102059579936554555.post-26206101927061451162008-11-16T13:10:00.001+08:002008-11-16T13:11:17.397+08:00Because QOCA cannot run in Solaris x86, and RH9 is much smaller.<div class="blogger-post-footer"><img width='1' height='1' src='https://blogger.googleusercontent.com/tracker/102059579936554555-2620610192706145116?l=yftian.blogspot.com' alt='' /></div>Enodhttp://www.blogger.com/profile/08433335265518211485noreply@blogger.com1tag:blogger.com,1999:blog-102059579936554555.post-38248216435044798892008-11-15T23:16:00.002+08:002008-11-15T23:20:46.384+08:00do not court for the newest thing, the best thing, ...<br />the BEST is the suitable one. <br />Enough is the BEST.<div class="blogger-post-footer"><img width='1' height='1' src='https://blogger.googleusercontent.com/tracker/102059579936554555-3824821643504479889?l=yftian.blogspot.com' alt='' /></div>Enodhttp://www.blogger.com/profile/08433335265518211485noreply@blogger.com0tag:blogger.com,1999:blog-102059579936554555.post-10769657948553580762008-11-15T23:12:00.001+08:002008-11-15T23:16:14.809+08:00I wanna quit the Internet now. Just necessary communication, papers, and data.<div class="blogger-post-footer"><img width='1' height='1' src='https://blogger.googleusercontent.com/tracker/102059579936554555-1076965794855358076?l=yftian.blogspot.com' alt='' /></div>Enodhttp://www.blogger.com/profile/08433335265518211485noreply@blogger.com0tag:blogger.com,1999:blog-102059579936554555.post-55670922717746820482008-11-14T17:03:00.002+08:002008-11-14T17:09:07.489+08:00It's a tragic day for me today. My 320GB Hitachi SATA USB harddisk failed to appear in Explorer of Windows XP. The testing program reports many many "unreadable sector" errors. About 320GB data (much of them are important backup data) were lost. It seems no way to get them back.<br /><br />So, I should make it obligatory to make a copy of my most important files in the future.<div class="blogger-post-footer"><img width='1' height='1' src='https://blogger.googleusercontent.com/tracker/102059579936554555-5567092271774682048?l=yftian.blogspot.com' alt='' /></div>Enodhttp://www.blogger.com/profile/08433335265518211485noreply@blogger.com2tag:blogger.com,1999:blog-102059579936554555.post-32547417421108945412008-11-03T20:45:00.000+08:002008-11-03T20:46:14.280+08:00write(*,'(I)') i <br /> 1<br />Error: Nonnegative width required in format string at (1)<br /><br /><br />Must I specify the width each time? Why ifort, f77, and Sun Studio f90 support this kind of formatting?<div class="blogger-post-footer"><img width='1' height='1' src='https://blogger.googleusercontent.com/tracker/102059579936554555-3254741742110894541?l=yftian.blogspot.com' alt='' /></div>Enodhttp://www.blogger.com/profile/08433335265518211485noreply@blogger.com0tag:blogger.com,1999:blog-102059579936554555.post-60061714321312535702008-11-03T20:32:00.002+08:002008-11-03T20:35:55.008+08:00The error message in Cygwin looks like:<br /><span style="font-family: courier new;"><br /> nan=0d10/0.<br /> ^<br />Division by 0 (zero) at (^) (IEEE not yet supported)</span><br /><br />I checked the GNU web site. It says GFortran supports isNaN() function. I think g77 does not support it.<br /><br />In Solaris, I can write a function similar to isNaN(), because IEEE Not-a-number is supported.<br /><br />The intel &amp; imsl also has this FORTRAN intrinsic.<br /><br />What should I do in Cygwin with g77?<div class="blogger-post-footer"><img width='1' height='1' src='https://blogger.googleusercontent.com/tracker/102059579936554555-6006171432131253570?l=yftian.blogspot.com' alt='' /></div>Enodhttp://www.blogger.com/profile/08433335265518211485noreply@blogger.com0tag:blogger.com,1999:blog-102059579936554555.post-18400944342524148022008-10-22T14:48:00.001+08:002008-10-22T14:49:50.045+08:00I cannot view my blog page now. But I can write blogs.<div class="blogger-post-footer"><img width='1' height='1' src='https://blogger.googleusercontent.com/tracker/102059579936554555-1840094434252414802?l=yftian.blogspot.com' alt='' /></div>Enodhttp://www.blogger.com/profile/08433335265518211485noreply@blogger.com0tag:blogger.com,1999:blog-102059579936554555.post-1851008954170488672008-10-22T14:46:00.001+08:002008-10-22T14:48:08.918+08:00In solaris, the find command gives that results sorted by letters.<br />However, in Linux, find command prints the results randomly.<div class="blogger-post-footer"><img width='1' height='1' src='https://blogger.googleusercontent.com/tracker/102059579936554555-185100895417048867?l=yftian.blogspot.com' alt='' /></div>Enodhttp://www.blogger.com/profile/08433335265518211485noreply@blogger.com0tag:blogger.com,1999:blog-102059579936554555.post-7686449929455519422008-10-22T14:44:00.001+08:002008-10-22T14:45:58.991+08:00Yes IDL 7 provides many new interesting features. However, it took up much of the CPU of my notebook. I had to revert to IDL 6.4 for best performance.<div class="blogger-post-footer"><img width='1' height='1' src='https://blogger.googleusercontent.com/tracker/102059579936554555-768644992945551942?l=yftian.blogspot.com' alt='' /></div>Enodhttp://www.blogger.com/profile/08433335265518211485noreply@blogger.com0tag:blogger.com,1999:blog-102059579936554555.post-10036984289307823822008-10-14T23:00:00.000+08:002008-10-14T23:01:42.888+08:00<pre> <br /><0.5 mm. The observed annual variations, however, disagree; this difference may be caused by some combination of local precipitation-induced site motion, unmodeled loading from other nearby sources, errors in the GSL model, and atmospheric errors." /> <br /> <br />GEOPHYSICAL RESEARCH LETTERS, VOL. 30, NO. 3, 1111, doi:10.1029/2002GL016579, 2003 <br /> <br />Crustal loading near Great Salt Lake, Utah <br /> <br />P. Elósegui and J. L. Davis <br /> <br />Harvard-Smithsonian Center for Astrophysics, <br />Cambridge, Massachusetts, USA <br /> <br />J. X. Mitrovica <br /> <br />Department of Physics, University of Toronto, <br />Toronto, Canada <br /> <br />R. A. Bennett <br /> <br />Harvard-Smithsonian Center for Astrophysics, <br />Cambridge, Massachusetts, USA <br /> <br />B. P. Wernicke <br /> <br />Division of Geological and Planetary Sciences, California Institute of Technology, <br />Pasadena, California, USA <br /> <br />Abstract <br /> <br />[1] Two sites of the BARGEN GPS network are located ~30 km south of Great Salt Lake (GSL). Lake-level records since mid-1996 indicate seasonal water elevation variations of ~0.3 m amplitude superimposed on a roughly “decadal” feature of amplitude ~0.6 m. Using an elastic Green's function and a simplified load geometry for GSL, we calculate that these variations translate into radial crustal loading signals of ±0.5 mm (seasonal) and ±1 mm (decadal). The horizontal loading signals are a factor of ~2 smaller. Despite the small size of the expected loading signals, we conclude that we can observe them using GPS time series for the coordinates of these two sites. The observed amplitudes of the variations agree with the predicted decadal variations to <0.5 mm. The observed annual variations, however, disagree; this difference may be caused by some combination of local precipitation-induced site motion, unmodeled loading from other nearby sources, errors in the GSL model, and atmospheric errors. <br /> <br />Received 7 November 2002; revised 4 December 2002; accepted 23 December 2002; published 5 February 2003. <br /> <br />Index Terms: 1208 Geodesy and Gravity: Crustal movements—intraplate (8110); 1243 Geodesy and Gravity: Space geodetic surveys; 1299 Geodesy and Gravity: General or miscellaneous; 8164 Tectonophysics: Evolution of the Earth: Stresses—crust and lithosphere. <br /></pre><div class="blogger-post-footer"><img width='1' height='1' src='https://blogger.googleusercontent.com/tracker/102059579936554555-1003698428930782382?l=yftian.blogspot.com' alt='' /></div>Enodhttp://www.blogger.com/profile/08433335265518211485noreply@blogger.com3tag:blogger.com,1999:blog-102059579936554555.post-10900017768412629752008-10-13T23:31:00.000+08:002008-10-13T23:32:25.227+08:00I hate it.<div class="blogger-post-footer"><img width='1' height='1' src='https://blogger.googleusercontent.com/tracker/102059579936554555-1090001776841262975?l=yftian.blogspot.com' alt='' /></div>Enodhttp://www.blogger.com/profile/08433335265518211485noreply@blogger.com1tag:blogger.com,1999:blog-102059579936554555.post-5605011033504279492008-10-10T11:19:00.001+08:002008-10-10T11:22:16.517+08:00Becker, J., and M. Bevis (2004), Love’s problem, Geophys. J. Int., 156,171– 178.<br /><br />Bevis, M., E. Kendrick, A. Cser, and R. Smalley (2004), Geodetic measurement of the local elastic response to the changing mass of water in Lago Laja, Chile, Phys. Earth Planet. Inter., 141, 71–78.<br /><br />Dziewonski, A., and D. Anderson (1981), Preliminary reference Earth model, Phys. Earth Planet. Inter., 25, 297– 356.<br /><br />Heki, K. (2001), Seasonal modulation of interseismic strain buildup in Northeastern Japan driven by snow loads, Science, 293, 89– 92.<div class="blogger-post-footer"><img width='1' height='1' src='https://blogger.googleusercontent.com/tracker/102059579936554555-560501103350427949?l=yftian.blogspot.com' alt='' /></div>Enodhttp://www.blogger.com/profile/08433335265518211485noreply@blogger.com0tag:blogger.com,1999:blog-102059579936554555.post-28868122639584768822008-10-07T21:00:00.001+08:002008-10-07T21:00:59.281+08:00<table border="0" cellpadding="0" cellspacing="0"><tbody><tr><td bgcolor="#339999" valign="top" width="18"><br /></td> <td align="left" bg valign="middle" height="20" style="color:#339999;"><span style="font-family:Helvetica;font-size:+1;color:#ffffff;"><b>Love Number</b></span></td> <td align="left" bgcolor="#339999" valign="top" width="14"><img src="http://scienceworld.wolfram.com/images/tealtab_topright.gif" alt="" border="0" /><br /></td> </tr> </tbody></table> <!-- End Title Tab --> <!-- Main Page --> <table border="0" cellpadding="0" cellspacing="0"><tbody><tr valign="top"> <td colspan="6"><img src="http://scienceworld.wolfram.com/images/gradient-teal.gif" alt="" height="2" /></td> </tr> <tr> <td bgcolor="#339999"><img src="http://scienceworld.wolfram.com/images/spacer.gif" alt="" width="2" height="10" /></td> </tr> <tr> <td bgcolor="#339999"><img src="http://scienceworld.wolfram.com/images/spacer.gif" alt="" width="2" /></td> <td colspan="3"> </td> <td> <p> <i>h</i> is the ratio of the height of a body tide to the static marine tide (introduced by A. E. H. Love). <i>k</i> is ratio of additional potential produced by the redistribution of mass to the deforming potential (introduced by A. E. H. Love). <i>l</i> is the ratio of horizontal displacement of the crust to that of the equilibrium fluid tide (introduced by T. Shida). For a rigid body, </p><div align="center"> <img src="http://scienceworld.wolfram.com/physics/limg328.gif" border="0" width="98" height="24" /><br /><br /><a href="http://scienceworld.wolfram.com/physics/LoveNumber.html">http://scienceworld.wolfram.com/physics/LoveNumber.html</a><br /></div></td></tr></tbody></table><div class="blogger-post-footer"><img width='1' height='1' src='https://blogger.googleusercontent.com/tracker/102059579936554555-2886812263958476882?l=yftian.blogspot.com' alt='' /></div>Enodhttp://www.blogger.com/profile/08433335265518211485noreply@blogger.com0tag:blogger.com,1999:blog-102059579936554555.post-79795965216595217672008-09-23T16:25:00.000+08:002008-09-23T16:26:30.390+08:00<a href="http://www.pgroup.com/userforum/viewtopic.php?t=929&amp;highlight=object+windows">http://www.pgroup.com/userforum/viewtopic.php?t=929&amp;highlight=object+windows</a><br /><br />export PGI_OBJSUFFIX=o<br /><br />Then run pgf77.<div class="blogger-post-footer"><img width='1' height='1' src='https://blogger.googleusercontent.com/tracker/102059579936554555-7979596521659521767?l=yftian.blogspot.com' alt='' /></div>Enodhttp://www.blogger.com/profile/08433335265518211485noreply@blogger.com0tag:blogger.com,1999:blog-102059579936554555.post-66531866392506114452008-09-10T00:09:00.000+08:002008-09-10T00:13:28.223+08:00For this<br />study, daily and seasonal irregularities in the global<br />reference frame were minimized using a regional stabilization<br />based on subtracting average common-mode<br />position variations at six sites (ALBH, CABL, DRAO,<br />GOBS, PABH, and QUIN, Plate 1). This technique<br />[Bock et al., 1997; Wdowinski et al., 1997] is particularly<br />useful in cases where the length of observations i s<br />short, and hence the relative effect of uncorrected seasonal<br />noise is potentially large.<div class="blogger-post-footer"><img width='1' height='1' src='https://blogger.googleusercontent.com/tracker/102059579936554555-6653186639250611445?l=yftian.blogspot.com' alt='' /></div>Enodhttp://www.blogger.com/profile/08433335265518211485noreply@blogger.com0tag:blogger.com,1999:blog-102059579936554555.post-31670326063880698822008-09-08T20:54:00.001+08:002008-09-08T20:55:47.930+08:00For the two global solutions we find that as expected,<br />the horizontal components are less noisy for white and<br />flicker noise magnitudes than the vertical components by a<br />factor of 2–3. <span style="font-weight: bold; font-style: italic; color: rgb(255, 0, 0);">The east components, also as expected, are<br />somewhat noisier than the north components because of no<br />(SOPAC) or incomplete (JPL) integer-cycle phase ambiguity<br />resolution in the global solutions.</span> The white noise amplitudes<br />are about 2 times larger in the JPL global solutions<br />compared to the SOPAC global solutions, in all three<br />components, while the flicker noise amplitudes are only<br />marginally larger. The white noise amplitude differences are<br />probably due to processing strategies (what these may be are<br />beyond the scope of this paper), while the comparable flicker<br />noise amplitudes may reflect a common physical basis, such<br />as seasonal atmospheric mass distributions [e.g<div class="blogger-post-footer"><img width='1' height='1' src='https://blogger.googleusercontent.com/tracker/102059579936554555-3167032606388069882?l=yftian.blogspot.com' alt='' /></div>Enodhttp://www.blogger.com/profile/08433335265518211485noreply@blogger.com2tag:blogger.com,1999:blog-102059579936554555.post-36045314777474571162008-09-07T10:21:00.001+08:002008-09-07T10:23:27.369+08:00When start a bash of X-window in Cygwin, the system CPU is highly utilized. Thus, it's good to start a bash without X-window support.<div class="blogger-post-footer"><img width='1' height='1' src='https://blogger.googleusercontent.com/tracker/102059579936554555-3604531477747457116?l=yftian.blogspot.com' alt='' /></div>Enodhttp://www.blogger.com/profile/08433335265518211485noreply@blogger.com0tag:blogger.com,1999:blog-102059579936554555.post-4097997709729333012008-09-01T23:39:00.001+08:002008-09-01T23:43:23.917+08:00It's frustrating. Urgly GUI, bad style, unfunctional buttons, ...<br /><br />It's just rubbish.<div class="blogger-post-footer"><img width='1' height='1' src='https://blogger.googleusercontent.com/tracker/102059579936554555-409799770972933301?l=yftian.blogspot.com' alt='' /></div>Enodhttp://www.blogger.com/profile/08433335265518211485noreply@blogger.com0tag:blogger.com,1999:blog-102059579936554555.post-28084046232330218002008-08-31T10:56:00.000+08:002008-08-31T10:57:00.080+08:00Just removed it.<div class="blogger-post-footer"><img width='1' height='1' src='https://blogger.googleusercontent.com/tracker/102059579936554555-2808404623233021800?l=yftian.blogspot.com' alt='' /></div>Enodhttp://www.blogger.com/profile/08433335265518211485noreply@blogger.com0tag:blogger.com,1999:blog-102059579936554555.post-37411650136440742802008-08-30T23:49:00.002+08:002008-08-30T23:55:46.792+08:00I mistakenly shift+deleted the <span style="font-style: italic;">Papers </span>directory which contains all the papers I wrote in the past few years. I thought I was deleting the directory in the backup drive; however, I was working in the source drive. It's like a disaster after I realized what I just did.<br /><br />(... anxious &amp; searching &amp; downloading...)<br /><br />(about one hour later)<br /><br />Oh, my god! I finally rescued myself. NO. It is <span style="font-weight: bold; font-style: italic;">Easy Recovery</span> saved me. I got back all the things I had just lose. Thanks god!<br /><br />BTW. <span style="font-style: italic;">Easy Recovery i</span>s far better than<span style="font-style: italic;"> Final Data 2</span>.<div class="blogger-post-footer"><img width='1' height='1' src='https://blogger.googleusercontent.com/tracker/102059579936554555-3741165013644074280?l=yftian.blogspot.com' alt='' /></div>Enodhttp://www.blogger.com/profile/08433335265518211485noreply@blogger.com0tag:blogger.com,1999:blog-102059579936554555.post-20745962593518733662008-08-24T10:54:00.000+08:002008-08-24T10:55:12.312+08:00/usr/lib/acct/accton /var/adm/pacct<br /><br />lastcomm [username]<div class="blogger-post-footer"><img width='1' height='1' src='https://blogger.googleusercontent.com/tracker/102059579936554555-2074596259351873366?l=yftian.blogspot.com' alt='' /></div>Enodhttp://www.blogger.com/profile/08433335265518211485noreply@blogger.com0tag:blogger.com,1999:blog-102059579936554555.post-92081068348313752292008-08-19T13:12:00.002+08:002008-08-19T13:14:52.215+08:00There is a program called "<span style="font-weight: bold; font-style: italic;">msconfig</span>" which can configure whether start the programs when the system starts up and user logins in the system.<div class="blogger-post-footer"><img width='1' height='1' src='https://blogger.googleusercontent.com/tracker/102059579936554555-9208106834831375229?l=yftian.blogspot.com' alt='' /></div>Enodhttp://www.blogger.com/profile/08433335265518211485noreply@blogger.com0tag:blogger.com,1999:blog-102059579936554555.post-16511170480865253542008-08-17T10:52:00.003+08:002008-08-17T10:55:12.017+08:00Install Sun's xVM VirtualBox software following the Readme.txt file.<br /><br />However, when running VitualBox, it gave an error:<br />> VirtualBox<br />ld.so.1: VirtualBox: fatal: libGL.so: open failed: No such file or directory<br />Killed<br /><br />The libGL.so was found in /usr/X11/lib/mesa/amd64/. Thus, add this to the library searching path will solve the problem:<br />setenv LD_LIBRARY_PATH /usr/X11/lib/mesa/amd64/:${LD_LIBRARY_PATH}<br /><br />Now, I can running Windows XP SP2 in Solaris.<div class="blogger-post-footer"><img width='1' height='1' src='https://blogger.googleusercontent.com/tracker/102059579936554555-1651117048086525354?l=yftian.blogspot.com' alt='' /></div>Enodhttp://www.blogger.com/profile/08433335265518211485noreply@blogger.com0tag:blogger.com,1999:blog-102059579936554555.post-64461353498230089422008-08-13T14:22:00.003+08:002008-08-13T14:26:35.966+08:00It is about the version of zcat.<br /><br />In my solaris 10 x86 system, the first zcat found was /usr/bin/zcat.<br /><span style="font-style: italic;font-size:85%;" >NAME<br /></span><span style="font-size:78%;"> compress, uncompress, zcat - compress, uncompress files or<br /> display expanded files<br /><br />SYNOPSIS<br /> compress [-fv] [-b bits] [file...]<br /><br /> compress [-cfv] [-b bits] [file]<br /><br /> uncompress [-cfv] [file...]<br /><br /> zcat [file...]</span><br /><br /><br />However, I think the dvipdfm program really wants GNU version of zcat. In Solaris, it is /usr/bin/gzcat.<br /><span style="font-size:78%;">NAME<br /> gzip, gunzip, gzcat - compress or expand files<br /><br />SYNOPSIS<br /> gzip [ -acdfhlLnNrtvV19 ] [-S suffix] [ name ... ]<br /> gunzip [ -acfhlLnNrtvV ] [-S suffix] [ name ... ]<br /> gzcat [ -fhLV ] [ name ... ]</span><br /><br />Thus, link gzcat to zcat and preceeds the /usr/bin/zcat will solves the problem.<div class="blogger-post-footer"><img width='1' height='1' src='https://blogger.googleusercontent.com/tracker/102059579936554555-6446135349823008942?l=yftian.blogspot.com' alt='' /></div>Enodhttp://www.blogger.com/profile/08433335265518211485noreply@blogger.com0tag:blogger.com,1999:blog-102059579936554555.post-86589805168658018212008-08-13T13:58:00.000+08:002008-08-13T13:59:11.176+08:00Q: 如何去掉WinEdt的自动换行？<br />A: WinEdt的编辑窗口最下面的状态条中有一项wrap，将其点击为灰色<br /><br />Q: WinEdt总是自动把几个单词缩到一行里面，如何解决？<br />A: 1. 关掉wrap功能<br /> 2. 在每行的首字前面加一个空格或一个"Tab"<div class="blogger-post-footer"><img width='1' height='1' src='https://blogger.googleusercontent.com/tracker/102059579936554555-8658980516865801821?l=yftian.blogspot.com' alt='' /></div>Enodhttp://www.blogger.com/profile/08433335265518211485noreply@blogger.com0tag:blogger.com,1999:blog-102059579936554555.post-69337183243888750102008-08-13T11:33:00.004+08:002008-08-13T11:45:15.776+08:00南南姑妈和姐姐要会宜昌，没法子搞到票，只好自己出面去西站碰碰运气。提前1个多小时到售票厅，听一售票员说晚7点以后只有5号窗口之后的才售票。此时6号窗口还没开，前面有2人排队，于是排第3。6点多的时候窗口开了，但要买5天之后的票要到7点才开始卖，于是排队等着。买4天之内的票可直接上前购买，不过很多人不知道，白排了半天队。而且售票员一直不停地在出4日内各车次卧铺票（只要有票就出），目的不明，不象是黄牛倒票。快到7点的时候，排第2的突然走了，于是前面又少了一个人。^_^。7点，放票了，运气不错，居然让我买了2两k49硬卧下铺。哈哈。。回家啰。。往外走的时候才发现：k，后面排队的人真多！！<br /><br />小结：去的早 + 运气不错。<div class="blogger-post-footer"><img width='1' height='1' src='https://blogger.googleusercontent.com/tracker/102059579936554555-6933718324388875010?l=yftian.blogspot.com' alt='' /></div>Enodhttp://www.blogger.com/profile/08433335265518211485noreply@blogger.com0tag:blogger.com,1999:blog-102059579936554555.post-57877432294965186502008-08-13T11:19:00.000+08:002008-08-13T11:20:41.657+08:00It seems that there is no way to force overwriting every files when downloading files using wget. However, use -N option can surely force downloading and overwriting newer files.<div class="blogger-post-footer"><img width='1' height='1' src='https://blogger.googleusercontent.com/tracker/102059579936554555-5787743229496518650?l=yftian.blogspot.com' alt='' /></div>Enodhttp://www.blogger.com/profile/08433335265518211485noreply@blogger.com4tag:blogger.com,1999:blog-102059579936554555.post-26799435800644114632008-08-13T10:45:00.000+08:002008-08-13T10:46:15.982+08:00How to ?<div class="blogger-post-footer"><img width='1' height='1' src='https://blogger.googleusercontent.com/tracker/102059579936554555-2679943580064411463?l=yftian.blogspot.com' alt='' /></div>Enodhttp://www.blogger.com/profile/08433335265518211485noreply@blogger.com0tag:blogger.com,1999:blog-102059579936554555.post-82575230590591782202008-08-13T10:42:00.001+08:002008-08-13T10:43:58.588+08:00Copy /etc/inet/ntp.server to /etc/inet/net.conf. And edit the ntp.conf,<br /><br />server asia.pool.ntp.org<br />server cn.pool.ntp.org<br />fudge cn.pool.ntp.org stratum 0<br /><br />broadcast 224.0.1.1 ttl 4<br /><br /><br />Then svcadm enable ntp<div class="blogger-post-footer"><img width='1' height='1' src='https://blogger.googleusercontent.com/tracker/102059579936554555-8257523059059178220?l=yftian.blogspot.com' alt='' /></div>Enodhttp://www.blogger.com/profile/08433335265518211485noreply@blogger.com0tag:blogger.com,1999:blog-102059579936554555.post-82892646595000454332008-08-11T12:14:00.002+08:002008-08-11T12:15:52.425+08:00It's terrific.<br />”（运动员的表现）真是太可怕了.“<br /><br />谁翻译的？！<div class="blogger-post-footer"><img width='1' height='1' src='https://blogger.googleusercontent.com/tracker/102059579936554555-8289264659500045433?l=yftian.blogspot.com' alt='' /></div>Enodhttp://www.blogger.com/profile/08433335265518211485noreply@blogger.com0tag:blogger.com,1999:blog-102059579936554555.post-51111564585714719102008-07-05T23:54:00.000+08:002008-07-05T23:55:12.092+08:00Which one is faster? Which one is better?<br /><br />Need your attention.<div class="blogger-post-footer"><img width='1' height='1' src='https://blogger.googleusercontent.com/tracker/102059579936554555-5111156458571471910?l=yftian.blogspot.com' alt='' /></div>Enodhttp://www.blogger.com/profile/08433335265518211485noreply@blogger.com0tag:blogger.com,1999:blog-102059579936554555.post-72463276379935654652008-07-05T23:51:00.001+08:002008-07-05T23:53:18.774+08:00Frequency Analysis Mapping On Unusual Sampling<br /><br><br /><a href='ftp://ftp.obs-nice.fr/pub/mignard/'>ftp://ftp.obs-nice.fr/pub/mignard/</a><br /><br><br />What does SNR of the results mean?<div class="blogger-post-footer"><img width='1' height='1' src='https://blogger.googleusercontent.com/tracker/102059579936554555-7246327637993565465?l=yftian.blogspot.com' alt='' /></div>Enodhttp://www.blogger.com/profile/08433335265518211485noreply@blogger.com0tag:blogger.com,1999:blog-102059579936554555.post-17628832375592871862008-07-05T22:51:00.003+08:002008-07-05T22:55:49.548+08:00Previously, est_noise by Langbein compiles successfully but quit when run it. The only error message was "Killed"--nothing else. This proves to be the problem of large expense of ram--est_noise requires 1.7+GB memory to run. I found that the maximum number of observations is hardwired into the fortran code (8103). Considering my CGPS time series are not so long. I changed this threshold to 5002. It works! Great.<div class="blogger-post-footer"><img width='1' height='1' src='https://blogger.googleusercontent.com/tracker/102059579936554555-1762883237559287186?l=yftian.blogspot.com' alt='' /></div>Enodhttp://www.blogger.com/profile/08433335265518211485noreply@blogger.com0tag:blogger.com,1999:blog-102059579936554555.post-74758743516868866642008-07-05T18:37:00.000+08:002008-07-05T18:39:43.364+08:00北京站, 为什么把票先打出来再慢慢地卖?难怪代售点没有票.<div class="blogger-post-footer"><img width='1' height='1' src='https://blogger.googleusercontent.com/tracker/102059579936554555-7475874351686886664?l=yftian.blogspot.com' alt='' /></div>Enodhttp://www.blogger.com/profile/08433335265518211485noreply@blogger.com0tag:blogger.com,1999:blog-102059579936554555.post-10712991693791995662008-07-03T12:13:00.000+08:002008-07-03T12:15:58.706+08:00What are those terms in model.dat file?<br /><br />I learned the meaning of max.dat from the source code. But the model.dat does not give me much information. I should dip more into the code.<div class="blogger-post-footer"><img width='1' height='1' src='https://blogger.googleusercontent.com/tracker/102059579936554555-1071299169379199566?l=yftian.blogspot.com' alt='' /></div>Enodhttp://www.blogger.com/profile/08433335265518211485noreply@blogger.com0tag:blogger.com,1999:blog-102059579936554555.post-21962445911103798492008-07-03T08:48:00.000+08:002008-07-03T08:49:31.693+08:00The firefox is good indeed. However, I wanna try more browsers.<div class="blogger-post-footer"><img width='1' height='1' src='https://blogger.googleusercontent.com/tracker/102059579936554555-2196244591110379849?l=yftian.blogspot.com' alt='' /></div>Enodhttp://www.blogger.com/profile/08433335265518211485noreply@blogger.com1tag:blogger.com,1999:blog-102059579936554555.post-22037226763384673722008-07-03T08:37:00.000+08:002008-07-03T08:38:59.054+08:00The hot key cannot call out the Chinese input method toolbar. Can I?<div class="blogger-post-footer"><img width='1' height='1' src='https://blogger.googleusercontent.com/tracker/102059579936554555-2203722676338467372?l=yftian.blogspot.com' alt='' /></div>Enodhttp://www.blogger.com/profile/08433335265518211485noreply@blogger.com0tag:blogger.com,1999:blog-102059579936554555.post-60647337039201085292008-07-03T08:32:00.000+08:002008-07-03T08:35:12.879+08:00If you choose Chinese GUI for solaris, the Firefox will start with Chinese menu and buttons. It is started with more parameters,<br /><br />> ps -ef | grep firefox<br /> tianyf 16040 15961 0 08:17:58 pts/14 0:00 /bin/sh /usr/bin/firefox<br /> tianyf 16054 16040 0 08:17:58 pts/14 0:00 /bin/sh /usr/lib/firefox/run-mozilla.sh /usr/lib/firefox/firefox-bin -UILocale <br /> tianyf 16060 16054 1 08:17:58 pts/14 0:31 /usr/lib/firefox/firefox-bin -UILocale zh-CN -contentLocale CN<br /> <br />What if we start only /usr/lib/firefox/firefox-bin ?<br />Ha.. It's English menu in Chinese GUI now. Good.<div class="blogger-post-footer"><img width='1' height='1' src='https://blogger.googleusercontent.com/tracker/102059579936554555-6064733703920108529?l=yftian.blogspot.com' alt='' /></div>Enodhttp://www.blogger.com/profile/08433335265518211485noreply@blogger.com0tag:blogger.com,1999:blog-102059579936554555.post-15518759242023552522008-07-02T22:34:00.001+08:002008-07-02T22:34:54.836+08:00<pre><br />#!/bin/sh<br /># sample script to run est_noise6ac<br /># test data is CARH north from Tom Herring's 'combined' SCIGN solutions<br /><br /><br />rm -f seed.dat<br />echo 82789427 > seed.dat # create a random number for doing hypothesis testing a la Langbein [2004] Figure 4<br /><br />rm -f junk.in<br />cat > junk.in <<EOF<br />otr # format style of inputting specifications (using julian day style; year, month and day could be input, too)<br />1 # Number of data sets (almost always 1)<br />1995 1 2005 1 # time span to be considered<br />y # rate # Estimate the secular rate; almost alway 'y' for yes<br />0 # num of rate change; If you think there are any rate changes in data, this should be the number of rate changes; if non-zero, then you'll need to supply the time of the rate change<br />1 # num of periods; I'm solving for a 1-year periodicity in the data; this is NOT a band-pass filter<br />365.25 # period in days<br />1 # Number of offsets in data<br />2003 357.0 # time of offset due to San Simeon EQ (M6+)<br />0 # number of exponential exp(-t/tau) or Log (1 + t/tau) or Log (tau + t ) terms for characterizing post-seismic responce. This seems to work okay; outputs searchs for optimal tau on file fort.89<br />otr # format style of input data<br />CARH.n <br />0 # Number of time-series of auxillary data --- This is useful for analyzing Dilatometer data -- this would be the presure time series<br />1 # mininum sampling in days (usually 1 day for CGPS)<br />n # This is usually n for 'no'. With 2-color EDM stuff, sometimes I'd have more than one measurement on a single day and I'd want to combine these data<br />n # n for 'no' unless you want to do simulations for hypothesis testing (Figures 2-4 in Langbein, 2004)<br />0 # decimation parmeter; choice of 0, 1, 2, or 3. Use 0 for 'best' results as it uses all of the available data; the other choices toss-out observations but are useful for 'quick and dirty' results-- more exploritory in nature<br />1 float #white noise Initial guess of white noise is 1 mm and you are asking the program to estimate the white noise amplitude<br />1 float #power law amplitude; initial guess of PL amplitude is 1 and you are asking the program to estimate a better amplitude<br />1 float # power law index; initial guess of PL index is 1 (flicker) and you are asking the program to estimate a better index<br />0 fix #Gauss Markov; fixing the GM time constant to 0; do not make any estimates<br />.5 2. # Pass-band limits of the band-pass filtered noise 0.5 cycles/year to 2 cycles/year<br />1 # Number of poles in the BP filter; (goes 1, 2, 3, 4, 5, and 6)<br />0 fix # Amplitude of BP noise; in this example, it is set to 0 and not estimated<br />2 fix # power law index; Index for a second power law; it is set to 2 (random walk) and not updated in this example<br />0.0 fix # power law amplitude; amplitude for a second power law; it is set to 0 and not updated in this example<br />0 # Additional white noise to add to data (data.in); usually set to 0 but for cases such as strain or creep data where the power-law index is > 2.5, it becomes critical to add some white noise to the data in order to keep the data covariance from becoming singular.<br />EOF<br /><br />time ~/proglib/Strain/Progs/est_noise6ac < junk.in<br /><br />exit<br />EOF<br /><br /></pre><div class="blogger-post-footer"><img width='1' height='1' src='https://blogger.googleusercontent.com/tracker/102059579936554555-1551875924202355252?l=yftian.blogspot.com' alt='' /></div>Enodhttp://www.blogger.com/profile/08433335265518211485noreply@blogger.com1tag:blogger.com,1999:blog-102059579936554555.post-89739391066675745452008-07-02T22:31:00.003+08:002008-07-02T22:32:49.648+08:00My HP notebook keeps running with high speed fanning. The network usage was usually 20-30%, and the CPU usage about 15%.<div class="blogger-post-footer"><img width='1' height='1' src='https://blogger.googleusercontent.com/tracker/102059579936554555-8973939106667574545?l=yftian.blogspot.com' alt='' /></div>Enodhttp://www.blogger.com/profile/08433335265518211485noreply@blogger.com0tag:blogger.com,1999:blog-102059579936554555.post-24722092098402020462008-07-02T22:29:00.002+08:002008-07-02T22:31:07.025+08:00It seems that cron doesn't load the login script (say, ~/.profile); thus, I should manually source it at the first part of cron script file.<div class="blogger-post-footer"><img width='1' height='1' src='https://blogger.googleusercontent.com/tracker/102059579936554555-2472209209840202046?l=yftian.blogspot.com' alt='' /></div>Enodhttp://www.blogger.com/profile/08433335265518211485noreply@blogger.com0tag:blogger.com,1999:blog-102059579936554555.post-81985880382241512602008-07-02T22:25:00.002+08:002008-07-02T22:28:57.804+08:00The est_noise6ac.f program failed to compile with f77 compiler of Sun Studio 12. It prints an error saying something about a missing function "cpu_time()" which may be a f90/gcc function. This program was finally compiled with f90 of Sun Studio 12, with the sun performance library (-xlic_lib=sunperf).<br /><br />The est_noise program took about 1700+MB ram when running.<div class="blogger-post-footer"><img width='1' height='1' src='https://blogger.googleusercontent.com/tracker/102059579936554555-8198588038224151260?l=yftian.blogspot.com' alt='' /></div>Enodhttp://www.blogger.com/profile/08433335265518211485noreply@blogger.com0tag:blogger.com,1999:blog-102059579936554555.post-23757203745265479952008-07-02T22:21:00.002+08:002008-07-02T22:24:53.125+08:00I found that the "crontab -e" command cannot work for the newly created user "datadm".<br />It prints "0" and "?", then seems hanging there without ending. It turns out to be the problem of default editor setting. Set it to "vi" will solve this problem. Now, put the following lines input your .profile if you are using "sh".<br /><br />EDITOR=vi<br />export EDITOR<br /><br />That's all. You can now use "crontal -e" to create scheduled tasks.<div class="blogger-post-footer"><img width='1' height='1' src='https://blogger.googleusercontent.com/tracker/102059579936554555-2375720374526547995?l=yftian.blogspot.com' alt='' /></div>Enodhttp://www.blogger.com/profile/08433335265518211485noreply@blogger.com0tag:blogger.com,1999:blog-102059579936554555.post-64846309158759057192008-07-02T22:14:00.002+08:002008-07-02T22:19:54.508+08:00The same input data and parameters, it took solaris (dell precision 390 with 2GB RAM) 6 min to complete. While rhel-5 (dell optiplex 755 with 4GB ram) obtained results within fifty-six seconds.<br /><br />Why so large difference in speed? Is the solaris not suitable for heavy load of numerical computations? Or the difference of RAM matters?<div class="blogger-post-footer"><img width='1' height='1' src='https://blogger.googleusercontent.com/tracker/102059579936554555-6484630915875905719?l=yftian.blogspot.com' alt='' /></div>Enodhttp://www.blogger.com/profile/08433335265518211485noreply@blogger.com0tag:blogger.com,1999:blog-102059579936554555.post-69211652780177423792008-07-02T22:11:00.000+08:002008-07-02T22:12:38.612+08:00算了，快完了，打酱油去了。<br />顺便做做俯卧撑。。<div class="blogger-post-footer"><img width='1' height='1' src='https://blogger.googleusercontent.com/tracker/102059579936554555-6921165278017742379?l=yftian.blogspot.com' alt='' /></div>Enodhttp://www.blogger.com/profile/08433335265518211485noreply@blogger.com0tag:blogger.com,1999:blog-102059579936554555.post-14120781536982865562008-07-02T08:58:00.001+08:002008-07-02T09:00:57.648+08:00<pre><br />$ . .shrc<br />GPSMET_DC_HOME=/export/home/tianyf/gpsf/external/gpsmet_dc: 不是一标识符<br /><br />$ which sh<br />/usr/bin/sh<br /><br />$ uname -a<br />SunOS gpsdc 5.10 Generic_120012-14 i86pc i386 i86pc<br /><br />$ more .shrc<br />GPSMET_DC_DATA=/data0/igs0/pub; export GPSMET_DC_DATA<br />export GPSMET_DC_HOME=/export/home/tianyf/gpsf/external/gpsmet_dc<br />export PATH=$PATH:$GPSMET_DC_HOME/sh<br /><br /></pre><br />Where is the problem? This script works well with bash.<div class="blogger-post-footer"><img width='1' height='1' src='https://blogger.googleusercontent.com/tracker/102059579936554555-1412078153698286556?l=yftian.blogspot.com' alt='' /></div>Enodhttp://www.blogger.com/profile/08433335265518211485noreply@blogger.com1tag:blogger.com,1999:blog-102059579936554555.post-22456787400816180612008-07-01T15:07:00.002+08:002008-07-01T15:14:01.503+08:00两天来，被贵州那个什么事件搅的无法安心写paper了。持续关注中。。。<br /><br />疑惑：<br />1 这种事情有什么好遮遮掩掩的？秉公办理不就行了。否则难逃存在腐败、失职之说。<br />2 把上万群众说成暴徒。老百姓没事干去政府那干啥？怎么能把责任推给“少数别有用心”、“黑社会”分子？<br />3 定调不对。说什么“稳定压倒一切”。我们要的是公平、正义。<br /><br />中国这一锅粥到底糊成什么样了？<div class="blogger-post-footer"><img width='1' height='1' src='https://blogger.googleusercontent.com/tracker/102059579936554555-2245678740081618061?l=yftian.blogspot.com' alt='' /></div>Enodhttp://www.blogger.com/profile/08433335265518211485noreply@blogger.com0tag:blogger.com,1999:blog-102059579936554555.post-58772039280384920042008-07-01T15:02:00.002+08:002008-07-01T15:07:32.812+08:00当通过xmanager连接solaris时，使用ctrl+space切换输入法，结果切换的确实windows的输入法。怎么切换xmanager里solaris的输入法？<br /><br />1。右键。<br />2。工具-输入法首选项<br />3。触发键-添加-（如ctrl+f1）<br /><br />Done。<br /><br />还可以设置用，和。翻页----很方便。不错。<div class="blogger-post-footer"><img width='1' height='1' src='https://blogger.googleusercontent.com/tracker/102059579936554555-5877203928038492004?l=yftian.blogspot.com' alt='' /></div>Enodhttp://www.blogger.com/profile/08433335265518211485noreply@blogger.com0tag:blogger.com,1999:blog-102059579936554555.post-55973042043500749032008-07-01T14:57:00.000+08:002008-07-01T15:01:51.826+08:00I got Xmanager to support Solaris Chinese GUI. Haha..<br /><br />1. Download /usr/openwin/lib/locale/zh.GBK/X11/fonts/75dpi/ to C:\Program Files\NetSarang\Xmanager Enterprise\Fonts. You'd better rename it to 75dpigbk because there is already a direcoty call 75dpi.<br />2. In Xconfig, configure the default profile to add the 75dpigkb to font searching paths.<br />3. Restart Xmanager to take effects.<br /><br />How, at the login screen of solaris, select zh_GBK will work.<br /><br />It's great. Isn't it?<br /><br />I switched my internet connection from winxp nootbook to Solaris X86 workstation today, it's important to type Chinese characters in solaris sometimes.<div class="blogger-post-footer"><img width='1' height='1' src='https://blogger.googleusercontent.com/tracker/102059579936554555-5597304204350074903?l=yftian.blogspot.com' alt='' /></div>Enodhttp://www.blogger.com/profile/08433335265518211485noreply@blogger.com0tag:blogger.com,1999:blog-102059579936554555.post-28603989882656129652008-06-30T08:51:00.002+08:002008-06-30T08:54:17.844+08:00In Cygwin, it seems that quit the XWin by right click its tray icon and select "Exit" doesn't affect the running bash.exe processes. Thus, It's better to exit all terminals before quit XWin.<div class="blogger-post-footer"><img width='1' height='1' src='https://blogger.googleusercontent.com/tracker/102059579936554555-2860398988265612965?l=yftian.blogspot.com' alt='' /></div>Enodhttp://www.blogger.com/profile/08433335265518211485noreply@blogger.com0tag:blogger.com,1999:blog-102059579936554555.post-82055687524377903232008-06-28T19:27:00.002+08:002008-06-28T19:32:11.887+08:00It took IDL tens of minutes to completed. However, the work was done in several seconds by native f77 program. It's really amazing!<br /><br />There are totally 857703 lines and 4 columns data. The task is to sort the first column, find those unique values and do summations of each unique value of the first column. The output file turns out to have 343352 lines. It involves sorting, finding unique values, searching and summation operations. I cann't believe f77 program can finished it so fast.<div class="blogger-post-footer"><img width='1' height='1' src='https://blogger.googleusercontent.com/tracker/102059579936554555-8205568752437790323?l=yftian.blogspot.com' alt='' /></div>Enodhttp://www.blogger.com/profile/08433335265518211485noreply@blogger.com0tag:blogger.com,1999:blog-102059579936554555.post-25820201024947367022008-06-28T19:21:00.001+08:002008-06-28T19:22:33.453+08:00g77 in Cygwin can find the error:<br /><br />psdsios.f: In program `psdsios':<br />psdsios.f:90:<br /> integer*4 i,j,k,l,m,n<br /> 1<br />psdsios.f:119: (continued):<br /> INTEGER*4 jmax,n,nout,np,nmax<br /> 2<br />Invalid declaration of or reference to symbol `n' at (2) [initially seen at (1)]<br />In file included from psdsios.f:81:<br /><br /><br />However, this program compiles well in Solaris using f77 of Sun Studio 12.<div class="blogger-post-footer"><img width='1' height='1' src='https://blogger.googleusercontent.com/tracker/102059579936554555-2582020102494736702?l=yftian.blogspot.com' alt='' /></div>Enodhttp://www.blogger.com/profile/08433335265518211485noreply@blogger.com1tag:blogger.com,1999:blog-102059579936554555.post-60458366939886729422008-06-27T15:58:00.003+08:002008-06-27T16:02:52.963+08:00<a onblur="try {parent.deselectBloggerImageGracefully();} catch(e) {}" href="http://4.bp.blogspot.com/_2ATje7JNtMw/SGSeH6biOyI/AAAAAAAAAGg/emKBzZ19bjY/s1600-h/ssh_problem.JPG"><img style="margin: 0px auto 10px; display: block; text-align: center; cursor: pointer;" src="http://4.bp.blogspot.com/_2ATje7JNtMw/SGSeH6biOyI/AAAAAAAAAGg/emKBzZ19bjY/s320/ssh_problem.JPG" alt="" id="BLOGGER_PHOTO_ID_5216468127182895906" border="0" /></a><br />Suddenly, the sshd process of Cygwin took up nearly 100% one of my core CPU. The fan of my notebook ran so fast.<div class="blogger-post-footer"><img width='1' height='1' src='https://blogger.googleusercontent.com/tracker/102059579936554555-6045836693988672942?l=yftian.blogspot.com' alt='' /></div>Enodhttp://www.blogger.com/profile/08433335265518211485noreply@blogger.com0tag:blogger.com,1999:blog-102059579936554555.post-62589333000216701792008-06-17T21:30:00.001+08:002008-06-17T21:50:44.513+08:00ENVI 4.5 is based upon the latest IDL 7.0.<br /><br /><br />It's a great news for geophysicists.<div class="blogger-post-footer"><img width='1' height='1' src='https://blogger.googleusercontent.com/tracker/102059579936554555-6258933300021670179?l=yftian.blogspot.com' alt='' /></div>Enodhttp://www.blogger.com/profile/08433335265518211485noreply@blogger.com0tag:blogger.com,1999:blog-102059579936554555.post-86323124074912015062008-06-15T21:10:00.000+08:002008-06-15T21:11:03.280+08:00<a href="http://code.google.com/p/iphone-dev/wiki/Building">http://code.google.com/p/iphone-dev/wiki/Building</a><div class="blogger-post-footer"><img width='1' height='1' src='https://blogger.googleusercontent.com/tracker/102059579936554555-8632312407491201506?l=yftian.blogspot.com' alt='' /></div>Enodhttp://www.blogger.com/profile/08433335265518211485noreply@blogger.com0tag:blogger.com,1999:blog-102059579936554555.post-62524681986668133652008-06-15T20:23:00.000+08:002008-06-15T20:24:13.962+08:00$ arm-apple-darwin-gcc -v<br />Using built-in specs.<br />Target: arm-apple-darwin<br />Configured with: ../llvm-gcc-4.0-iphone/configure --enable-llvm=/tmp/llvm-svn-42<br />498r/.build --enable-languages=c,c++,objc,obj-c++ --target=arm-apple-darwin --en<br />able-sjlj-exceptions --with-as=/usr/local/bin/arm-apple-darwin-as --with-ld=/usr<br />/local/bin/arm-apple-darwin-ld --with-heavenly=/usr/local/arm-apple-darwin/files<br />ystem<br />Thread model: posix<br />gcc version 4.0.1 (Apple Computer, Inc. build 5449)(LLVM build 2.0-svn-iphone-de<br />v-0.3-svn)<div class="blogger-post-footer"><img width='1' height='1' src='https://blogger.googleusercontent.com/tracker/102059579936554555-6252468198666813365?l=yftian.blogspot.com' alt='' /></div>Enodhttp://www.blogger.com/profile/08433335265518211485noreply@blogger.com0tag:blogger.com,1999:blog-102059579936554555.post-32482853500043073082008-06-15T09:49:00.000+08:002008-06-15T09:50:13.141+08:00It just display the binary content of MP3 files. Why?<div class="blogger-post-footer"><img width='1' height='1' src='https://blogger.googleusercontent.com/tracker/102059579936554555-3248285350004307308?l=yftian.blogspot.com' alt='' /></div>Enodhttp://www.blogger.com/profile/08433335265518211485noreply@blogger.com0tag:blogger.com,1999:blog-102059579936554555.post-33161958534723708762008-06-13T22:28:00.001+08:002008-06-13T22:29:54.152+08:00It can not only upload/download files, but also change file properties and modifying files. Without it, I cannot make my iTouch work.<div class="blogger-post-footer"><img width='1' height='1' src='https://blogger.googleusercontent.com/tracker/102059579936554555-3316195853472370876?l=yftian.blogspot.com' alt='' /></div>Enodhttp://www.blogger.com/profile/08433335265518211485noreply@blogger.com0tag:blogger.com,1999:blog-102059579936554555.post-36887729353436930812008-06-13T22:15:00.002+08:002008-06-13T22:28:18.399+08:00Old configuration: ipod touch 8G + 1.1.1 with jailbreak crack. It works like a charm. However, the system became unconformable with time goes by. Thus, I decided to upgrade.<br /><br />1. After downloaded the 1.1.4 firmware, open iTunes, press Shift and clike the restore button, select the file and restore the iTouch. Now, everything old are ruined. You will be happy if you did a little backup work before restore iTouch.<br /><br />2. Download Ziphone 3.0 to jailbreak the iTouch. Close iTunes, kill ipod...exe, itune..help.exe processes with task manager of Windows Xp. Then run Ziphone 3 GUI. Click advanced selection options. Be sure to uncheck the activation and ... checkboxes if you are not using iPhone. Otherwise, the cracked system is unstable. Now, click "Jailbreak [include ipod]" and wait for about one minute.<br /><br />3. The iTouch is now ready to install packages. Search for wireless network. You can login the iTouch by using ssh of Cygwin/WinXP. The password for root is "alpine". DO NOT change the password. Or, you will get continuous of crash of Summerboard and it looks like a piece of brick.<br /><br />The left steps is adding some sources, install softwares for your application. But, there are still a lot to go. And the mixed-usage of <span style="font-style: italic;">mobile </span>and <span style="font-style: italic;">root </span>for default user will make you programs stopping working. Googling for more details.<div class="blogger-post-footer"><img width='1' height='1' src='https://blogger.googleusercontent.com/tracker/102059579936554555-3688772935343693081?l=yftian.blogspot.com' alt='' /></div>Enodhttp://www.blogger.com/profile/08433335265518211485noreply@blogger.com0