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)
 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 ch17  904.4nm
11. ground surface reflectance at ch26 1382.3nm
12. brightness temp at ch20  3.789um [K]
13. brightness temp at ch31 11um [K]
14. brightness temp at ch32 12um [K]
15. Daily mean PAR [Ein/m^2/day]
16. Daily mean of direct PAR [Ein/m^2/day]
17. Transmittance of instantaneous PAR at noon time []
18. Daily mean Shortwave radiation [W/m^2]
19. UV-A radiation  [W/m^2]
20. UV-B radiation [W/m^2]
21. CIE weighted UV radiation [W/m^2]
22. Aerosol Optical thickness at 466nm [Tau]
23. Aerosol Optical thickness at 554nm [Tau]
24. Aerosol Optical thickness at 646nm [Tau]
25. Aerosol Optical thickness at 857nm [Tau]
26. Aerosol angstrom exponent (parameter=DN*slope -1)
27. Cloud fraction
28. Cloud Tau at 550nm (average of cloudy areas; parameter=10^(DN*slope-1))
29. Ocean chlorophyll-a concentration [mg/m3] (parameter=10^(DN*slope-2))
30. Clear area precipitable water [mm]
31. Land/Ocean surface temperature [K]
32. Cloud top temperature  [K]

header format(2i6,2f8.2,f8.4,i3,32e12.5,32i3):
nl,ml,lon_min,lat_max,dd,nbl,(slope(k),k=1,nbl),(ich(k),k=1,nbl)
nbl=32

- 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