GPS Analyzed Data and Relates

GPS Data

GPS PWV data files is separated yearly and in each site. Missing date are indicated by -999.0.
GPS PWV is estimated by Prof. Satomura, Shizuoka Univ. High temporal resolution data is available from Satomura Laboratory at Shizuoka Univercity.

Daily Mean GPS PWV
- 1998 from March
- 1999
- 2000
- 2001 Coming soon !
Data format
julian_day pwv(mm)
Pentad Mean GPS PWV
- 1998 from March
- 1999
- 2000
- 2001 Coming soon !
Data format
pentad_number pwv(mm)

Other Data

Hourly surface meteorological data
- description of original data
  This data was prepared for GPS Meteorology by Thai Meteorological Department.
  Some components are available only in three hourly temporal resolution.
  - Period : 1st/January to 30th/November of 1998
  - Available stations : six GPS observation stations at TMD observatories.
  - Temporal Resolution : 1 hour

3 hourly convection indices
- description of original data
  Please
  - GMS histogram data
  - Period : 16st/April to 15th/October of 1998
  - Spatial resolution 25km x 25km
  - Temporal Resolution : 3 hour
  - Coverage area : 60S-60N, 70E-170W
  - Type : Tbb, threshold values: 250, 260, 270, 280
    Please contact with us if you need the complete data. Because the size of data is some giga bytes, we supply them by tape media.
- Time series of convection indices Ic280, max(280-Tbb, 0)
  - Bangkok Premonsoon, Early monsoon, Late monsoon
  - Chiang Mai Premonsoon, Early monsoon, Late monsoon
    These files are binary files and you can use the same program below presented.
- Diurnal variation of convection indices Ic280, max(280-Tbb, 0)
  - Bangkok
  - Chaing Mai
    Data format
    
    hour Premonsoon, Earyl_monsoon, Late_monsoon
    where hour is Thai local standard time
Pentad mean meteorological variables over Indochina Penisula
- Common information
  - Period : 16th April. to 14th Oct. of 1998
  - Spatial averaging area : 13N-21N, 98E-105E
  - Temporal Resolution : pentad
    Data format
    
    pentad_number data_values
- Rainfall
  - description of original data
    - CPC Merged Analysis of Precipitation
    - Period : Jan. to Dec. of 1998
    - Spatial resolution 2.5deg x 2.5deg
    - Temporal Resolution : pentad
- Outergoing Longwaver Radiation
  - description of original data
    - NOAA Interpolated OLR
    - Period : Jan. to Dec. of 1998
    - Spatial resolution 2.5deg x 2.5deg
    - Temporal Resolution : daily
- Zonal wind at 850hPa
  - description of original data
    - GAME/Reanalysis
    - Period : Apri. to Oct. of 1998
    - Spatial resolution 1.25deg x 1.25deg
    - Temporal Resolution : 6 hourly

Data from GAME Reanalysis Ver.1.

description of original data
- Period : 1st/April to 31st/October of 1998
- Spatial Resolution : 1.25 x 1.25 degrees.
- Temporal Resolution : 6 hour

Precipitable water

averaged over Thailand (98E-106E,13N-21N)
Phuket
Bangkok
Ubon Ratchatani
Sri Samrong
Nong Khai

Chiang Mai

These files are binary direct files and a following sample program help you to read them. Data is stored as single precision real number in a big endian byte order.

Fortran

     program sample
c
c    Coded by Tosi Nakaegawa
c
      integer, parameter :: nmax = 856   ! Number of Data
      real :: tpw(nmax)
      open( 11, file='tpw.dr',form='unformatted'
     &    ,access='direct', recl=nmax*4 )
c      
      read( 11, rec=1 ) tpw
c
      stop
      end

Related Data Informations

The data listed here is essential to analyzing the GPS data and to compare the results.

GAME Reanalysis
(Meteorologica Research Institute:Yamazaki, Chief)
From here in Japanese only
- 詳しいことは、天気または水文水資源学会誌を参照のこと。天気の内容は上記リンク先に掲載されています。
- データ配付:2.5°はCD-ROM、1.25°、0.5°は8mmテープのみ。 filesvにもあるけど...。
- 8mmテープの読み方。月毎にデータが入っている。1.25°の2次元物理量を取り出したい時は、
```
4月は、readmeとsampleが入っているため、
tosi > mt -f /dev/rmt/2n fsf 6
tosi > tar xvf/dev/rmt/2n 
それ以外の月は
tosi > mt -f /dev/rmt/2n fsf 3
tosi > tar xvf/dev/rmt/2n 
```
- 物理量解析値は2.5°と1.25°のみ。
- 0.5°はGAME領域のみで、3次元基本量だけ。
- 全ファイル数 : 856
- データ在処 : wcl3e027:/work9/tnakaega/1.25
- タイ領域で可降水量は10mmほど小さい。
- Ver.2 作成作業中。(2001.11.10現在)
NCEP/NCAR Reanalysis-daily data
(NCEP/NOAA)
From here in Japanese only
- 1949年からの客観解析日データを取得することができます。
- データハンドリングなどについても、記述があります。
- 日々自動でデータを送ってもらう設定もできます。
- 但し、データは GRIB形式のため、気象業界デファクトソフト GrADSでも直接は読み込めません。
- 下のURLにあるデータセットも同一の再解析値です。
- データセット
  - 空間：2.5x2.5度、144x73グリッド、90N～90S、0E～357.5E、
  - 時間：1948 年1月1日～現在まで、6時間間隔
  - 鉛直レベル17 pressure levels (hPa): 1000, 925, 850, 700, 600, 500, 400, 300, 250, 200, 150, 100, 70, 50, 30, 20, 10
NCEP/NCAR Reanalysis-6 hourly data
(CDC/NCEP)
When you do not have the surface meteorogical data to retrieve PWV from GPS signals, reanalysis data is very useful. The following web page by CDC/NCEP provides us with data essential to do it.
Microwave Radiometer at Sri Samrong
(CeReS, Univ. of Chiba, Prof. Takamura)
From here in Japanese only
```
TbSky23 23G 輝度温度
TbSky31 31G 輝度温度
TkBB    黒体温度
VapCM   可降水量cm
LiqCM   雲水量cm
DelCM   電波伝搬に於けるdelayの長さcm
AZact   Azimuth angle
ELact   Elevation angle（確か天頂が９０度ではなかったかと思います)
Tau23   23G 光学的厚さ
Tau31   31G 光学的厚さ
Tamb    気圧計を併設しているときの値
```
Tip curve 放射計の検定方法の一つで，大気輝度温度は，この波長では，大気輻射が絶対温度で数十度しかないために，液体窒素のような比較検定物を用いた検定が極めて難しく，放射伝達方程式から導かれる近似関係を用いて検定する方法です．簡単に言えば，横軸に大気の相対エアマス，縦軸に該当波長の光学的厚さをプロットしたとき，線形の関係があり，エアマスゼロの点が，2.7°Kの宇宙の背景放射になるというものです．下のELactの角度も多分それに使えるように選んであるのではないかと想像しています．従って，ELactの値が９０度のものが天頂の輝度温度と鉛直水蒸気量であり，それ以外は斜めに見たときのもので，輝度温度や水蒸気量も光学的厚さが増すのに対応して増えることになります．　これでよろしいでしょうか．radiometerの説明書は多分内山さんのところにあるので，それを見ていただければご理解頂けるのではと思います (千葉大学 CeReS 高村先生)．
Meteorological Surface observation data at SriSamrong
Raw data, namely in Voltage.

Meteorological Surface observation data at SriSamrong
Data of Physical Values

# column 1:  date                [yyyy/mm/dd]
# column 2:  hour of day         [o'clock]
# column 3:  solar zenith angle  [degree]
# column 4:  solar azimuth angle [degree]
# column 5:  relative air mass
# column 6:  sun-earth distance
# column 7:  downward visible irradiance (diffuse) [W/m^2]
# column 8:  upward visible irradiance [W/m^2]
# column 9:  direct solar irradiance [W/m^2]
# column 10: [3] + [5]*cos([2])   [W/m^2]
# column 11: surface albedo [10]/[8]
# column 12: PIR net  downward    [W/m^2]
# column 13: PIR body temp.       [K] 
# column 14: PIR dome temp.       [K]
# column 15: PIR net + body       [K]
# column 16: PIR net  upward      [W/m^2]
# column 17: PIR body temp.       [K]
# column 18: PIR dome temp.       [K]
# column 19: PIR net + body       [W/m^2]
# column 20: wind speed           [m/s]
# column 21: wind direction       [degree]
#            north(0deg)->east(90deg)->south(180deg)->west(270deg)
# column 22: pressure             [hPa]
# column 23: relative humidity    [%]
# column 24: temperature          [celsius]

# date,time
# CH.00 Pyranometer  (CM21)   diffuse
# CH.01 Pyranometer  (CM21)   reflect
# CH.02 Pyheliometer (MS-53A) direct
# CH.03 Pygeometer   (PIR)    downward
# CH.04 Pygeometer body temp
# CH.05 Pygeometer dome temp
# CH.06 Pygeometer   (PIR)    upward
# CH.07 Pygeometer body temp
# CH.08 Pygeometer dome temp
# CH.10 wind speed        [m/s]
# CH.11 wind direction    [degree]
#            north(0deg)->east(90deg)->south(180deg)->west(270deg)
# CH.12 pressure          [hPa]
# CH.13 relative humidity [%]
# CH.14 temperature       [celsius]
# CH.15 Reserved.
# CH.16 Reserved.
# CH.17 Reserved.
# CH.18 Reserved.