File description
see below about the global data:
0.25 deg
grid binary data: monthly, daily binary,&
readme file
0.05 deg grid binary data: monthly, daily binary,&
readme file
- File names
(a) 1-channel PAR
data
MDS021KM_J20080201Avh_c121_2701_2601_PAR_le : PAR
binary data for the web image
(half-month)
MDS021KM_J20080201Avm_c121_2701_2601_PAR_le :
PAR binary data for the web image (monthly)
1-byte (8bit) binary data: see here.
(b) 20 channel variables (old
version)
MDS021KM_J20080201Avh_c121_2701_2601_par : par data
including subparameters (SW, DIF..)
(b') 32 channel variables (new
version: v601)
MDS021KM_J20080201Avh_v601_2701_2601_par :
par data including subparameters (SW,
DIF..)
ref(1-14),par,dpar,tipar,swr,uva,uvb,cie,taua1~4,alp(-1),cfr,tauc(-1,log10),chla(-2,log10),ptw,lst,ctt
(c) Scene-by-scene PAR in a one file
per
half-months
MDS021KM_J20080201Avh_c121_2701_2601_daily### :
PAR of each scene
(d) RGB image
binary
MDS021KM_J20080201Avh_c121_2701_2601_1Krgb : RGB
binary data for the web image (half-month)
MDS021KM_J20080201Avm_c121_2701_2601_1Krgb :
RGB binary data for the web image (monthly)
(e) Text file of prefecture averages
and at several sample
sites
MDS_PAR_JPN.txt :
text data for the web graph plot (revised each half-month)
file
labels:
MDS021KM_J20080201Avh_c121_2701_2601
YYYYMMDD pixel size, line
size
Avh: half-month average, Avm: monthly average,
daily###: store original scene-by-scene values
###: scene number in the half month
pixel size and line
size are variable (e.g., 7200 & 3600, 1440 & 721..)
- Data
size:
(a) *_PAR_le: 2701*2byte header + 2byte
x 2701(pixel) x 2601(line) * 1 channel
(b) *_par:
2701*2byte header + 2byte x 2701(pixel) x 2601(line) * 20
channels
(b') *_par: 2701*2byte header + 2byte x 2701(pixel)
x 2601(line) * 32 channels
(c) *_daily###: 2701*2byte header + 2byte
x 2701(pixel) x 2601(line) * 5 channels *
###scenes
(d) *_1Krgb: 2701*1byte header + 1byte
x 2701(pixel) x 2601(line) * 3 channels
-
Grid:
MODIS: upper-left grid location (grid
center): 50N, 123E; grid interbal: 0.01 deg
(a) * "_PAR_le"
file
- Read
header (size= pixel size *2byte):
character
head*1000
read(10,rec=1)
head
read(head,'(2i6,2f8.2,f8.4,2e12.5,a1,a8,a1,a40)')
&
npixel,nline,lon_min,lat_max,reso,slope,offset,',',
& para,',',outfile
- Read
data (e.g., fortran77)
parameter(nl=2701, ml=2701)
... open
file by "unformatted", "recl=nl*2(byte)" (,"bytereclen")
integer*2
i2buf(nl,ml)
do
m=1,ml
read(10,rec=1+m) (i2buf(n,m), n=1,nl)
do
n=1,nl
par=i2buf(n,m)*slope+offset
write(6,*) 'PAR[Ein/m^2/day]=',par
enddo
enddo
(b) * "_par" file ("c121": old
version)
1. ground surface reflectance at ch01
645.8nm
2 ground surface reflectance at ch02 856.9nm
3
ground surface reflectance at ch03 466.1nm
4 ground surface
reflectance at ch04 553.9nm
5 ground surface reflectance at ch05
1241.5nm
6 ground surface reflectance at ch06 1628.1nm
7
ground surface reflectance at ch07 2113.9nm
8 ground surface
reflectance at ch08 412.4nm
9 ground surface reflectance at
ch09 442.2nm
10 ground surface reflectance at ch11 530.1nm
11 ground surface reflectance at ch17 904.4nm
12 ground surface
reflectance at ch26 1382.3nm
13 brightness temp at ch20 3.789um
[K]
14 brightness temp at ch31 11.006um [K]
15 SST [K]
16 Aerosol
Optical thickness at 550nm [Tau]
17 Direct PAR ratio
18 Cloud Tau at 550nm
(averaged both fine and cloudy areas)
19 Daily Shortwave radiation
[W/m^2]
20 PAR [Ein/m^2/day]
header
format(2i6,2f8.2,f8.4,i3,20e12.5,20i3):
nl,ml,lon_min,lat_max,dd,15+5,(slope(k),k=1,15+5),(cho(k),k=1,15),(36+k,k=1,5)
- Parameter
conversion:
parameters(k) = DN * slope(k)
example:
2701 2601
123.00 50.00 0.0100 20
0.10000E-03 0.10000E-03
0.10000E-03 0.10000E-03 0.10000E-03 0.10000E-03 0.10000E-03 0.10000E-03
0.10000E-03 0.10000E-03 0.10000E-03 0.10000E-03 0.10000E-01 0.10000E-01
0.10000E-01 0.10000E-03 0.10000E-03 0.20000E-03 0.20000E-01
0.10000E-01
1 2 3 4 5 6 7
8 9 11 17 20 21 31 32 37 38 39 40 41
(b') * "_par" file (32channel v601)
- Parameter
conversion:
parameters(k) = DN * slope(k) (except for
26, 28, 29)
(c) * "_daily" file ("c121": old
version)
do i=1,scenes
do j=1,5
if j==1
Aerosol Optical thickness
at 550nm [Tau]
else if j==2
Direct PAR ratio
else if j==3
Cloud Tau at 550nm
else if
j==4
Daily Shortwave radiation [W/m^2]
else if
j==5
PAR [Ein/m^2/day]
endif
enddo
enddo
header
format(2i6,2f8.2,f8.4,i3,5e12.5,500i3):
nl,ml,lon_min,lat_max,dd,5,(slope(15+k),k=1,5),(day(k),k=1,nscenes)
example:
2701 2601
123.00 50.00 0.0100 5
0.10000E-03 0.10000E-03
0.20000E-03 0.10000E-01 0.10000E-01
1 1 1 2
2 2 3 3 3 3 4 4 4 5
5 5 6 6 6 7 7 7 8 8
8 9 9 9 10 10 10 11 11 11 12 12 12 13 13 14 14 14 15 15
15
(c') * "_daily" file (11 * Nscene channels; v601)
nb=1~11: par,dpar,tipar,swr,uva,uvb,cie,taua1~4
1-byte unsiged integer: 1-line header + 1byte * pixel * line *11channel
header
format(2i6,2f8.2,f8.4,i3,11e12.5,500i3):
nl,ml,lon_min,lat_max,dd,5,(slope(14+k),k=1,11),(day(k),k=1,nscenes)
(d) text data
format:
read(10,'(i4,3(1x,i2),80f7.2)') year,month,day_min,day_max,(par_average(k),k=1,49),(site_samples(k),k=1,nloc)
# pav: japan(1), prefecture ID+1:2~48, forign country(49),valsites(50~80)
# Prefectures (ID name)
#
00 All Japan
# 01 Hokkaido
#
02 Aomori
# 03 Iwate
# 04 Miyagi
#
05 Akita
# 06 Yamagata
#
07 Fukushima
# 08 Ibaraki
#
09 Tochigi
# 10 Gunma
# 11 Saitama
#
12 Chiba
# 13 Tokyo
# 14 Kanagawa
#
15 Niigata
# 16 Toyama
# 17 Ishikawa
#
18 Fukui
# 19 Yamanashi
# 20 Nagano
#
21 Gifu
# 22 Shizuoka
# 23 Aichi
#
24 Mie
# 25 Shiga
# 26 Kyoto
#
27 Osaka
# 28 Hyogo
# 29 Nara
#
30 Wakayama
# 31 Tottori
#
32 Shimane
# 33 Okayama
#
34 Hiroshima
# 35 Yamaguchi
#
36 Tokushima
# 37 Kagawa
# 38 Ehime
#
39 Kohchi
# 40 Fukuoka
# 41 Saga
#
42 Nagasaki
# 43 Kumamoto
# 44 Ooita
#
45 Miyazaki
# 46 Kagoshima
#
47 Okinawa
# 48 Forign country areas
# val sites:
# 45.0500,
142.1000, ! TSE
# 45.2786,
127.5784, ! LSH
# 44.3842,
142.3186, ! MBF
# 44.3219,
142.2614, ! MMF
# 42.9900,
141.3875, ! HEF
# 42.9868,
141.3853, ! SAP
# 42.7370,
141.5186, ! TMK
# 42.4333,
142.4833, ! SZN
# 40.0017,
140.9375, ! API
# 37.7500,
127.1500, ! GDK
# 37.0667,
126.0333, ! har
# 36.9167,
139.9667, ! NSS
# 36.5500,
140.1333, ! NFT
# 36.4000,
138.5833, ! KZW
# 36.2114,
139.9474, ! SJP
# 36.1462,
137.4231, ! TKY
# 36.1397,
137.3708, ! TKC
# 36.1138,
140.0944, ! TGF
# 36.0540,
140.0269, ! MSE
# 35.8693,
139.4903, ! KWG
# 35.4514,
138.7653, ! FJY
# 35.4333,
138.7500, ! FHK
# 35.2500,
137.0667, ! SMF
# 34.9619,
135.9958, ! KEW
# 34.7948,
135.8462, ! YMS
# 34.7349,
134.3743, ! AKO
# 34.5500,
126.5700, ! HFK
# 34.5345,
133.9302, ! yhm
# 33.1370,
130.7095, ! KHW
# 31.9667,
130.9333, ! KBY
# 31.7347,
131.0139/ ! kmtc