|
1 ;************************************************ |
|
2 ; Read half-hourly data and write monthly data |
|
3 ; input data is : half-hourly |
|
4 ; |
|
5 ; output data is : CO2_flux, RAD_FLUX, SH_FLUX, LH_FLUX |
|
6 ; date, lat, lon |
|
7 ;************************************************ |
|
8 load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl" |
|
9 load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl" |
|
10 load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl" |
|
11 ;************************************************ |
|
12 begin |
|
13 |
|
14 dir_root = "/fis/cgd/cseg/people/jeff/clamp_data/fluxnet/" |
|
15 station_name = "Hyytiala" |
|
16 file_name_add = "carboeurope" |
|
17 year_start = 1996 |
|
18 year_end = 2003 |
|
19 year_leap1 = 1996 |
|
20 year_leap2 = 2000 |
|
21 year_leap3 = 2004 |
|
22 |
|
23 lat_out = 61.8474 |
|
24 lon_out = 24.2948 |
|
25 |
|
26 ; final output |
|
27 diro = dir_root + station_name + "/" |
|
28 filo = station_name+"_"+year_start+"-"+year_end+"_monthly.nc" |
|
29 c = addfile(diro+filo,"c") |
|
30 print (filo) |
|
31 |
|
32 ; input dir |
|
33 diri = dir_root + station_name + "/" |
|
34 |
|
35 nyear = year_end - year_start + 1 |
|
36 nmon = 12 |
|
37 nlat = 1 |
|
38 nlon = 1 |
|
39 |
|
40 ; day_of_month = (/31.,28.,31.,30.,31.,30.,31.,31.,30.,31.,30.,31./) |
|
41 end_of_month1 = (/31.,59.,90.,120.,151.,181.,212.,243.,273.,304.,334.,365./) |
|
42 end_of_month2 = (/31.,60.,91.,121.,152.,182.,213.,244.,274.,305.,335.,366./) |
|
43 |
|
44 data_out1 = new((/nyear,nmon/),float) |
|
45 data_out2 = new((/nyear,nmon/),float) |
|
46 data_out3 = new((/nyear,nmon/),float) |
|
47 data_out4 = new((/nyear,nmon/),float) |
|
48 time_out = new((/nyear,nmon/),integer) |
|
49 |
|
50 do m = 0,nyear-1 |
|
51 year = year_start + m |
|
52 |
|
53 ; if (mod(year,4) .eq. 0) then |
|
54 if (year.eq.year_leap1.or.year.eq.year_leap2.or.year.eq.year_leap3) then |
|
55 hour_of_month = end_of_month2 * 24. |
|
56 else |
|
57 hour_of_month = end_of_month1 * 24. |
|
58 end if |
|
59 ; print (hour_of_month) |
|
60 |
|
61 ; input file |
|
62 fili = file_name_add+"."+station_name+"."+year+".nc" |
|
63 b = addfile(diri+fili,"r") |
|
64 print (fili) |
|
65 |
|
66 ; co2 flux unit: umol m-2 s-1 |
|
67 data1 = b->FCO2M |
|
68 ; net radiation flux unit: W m-2 |
|
69 data2 = b->RN |
|
70 ; sensible heat flux unit: W m-2 |
|
71 data3 = b->H |
|
72 ; latent heat flux unit: W m-2 |
|
73 data4 = b->LE |
|
74 |
|
75 ; time unit: hour since 01-01 (mm-dd) |
|
76 time = b->time |
|
77 |
|
78 do n= 0,nmon-1 |
|
79 if (n.eq.0) then |
|
80 timeL = 0. |
|
81 timeR = hour_of_month(0) |
|
82 else |
|
83 timeL = hour_of_month(n-1) |
|
84 timeR = hour_of_month(n) |
|
85 end if |
|
86 ; print (timeL) |
|
87 ; print (timeR) |
|
88 |
|
89 i = ind(time.ge.timeL .and. time.lt.timeR) |
|
90 ; print (i) |
|
91 |
|
92 data_out1(m,n) = avg(data1(i)) |
|
93 data_out2(m,n) = avg(data2(i)) |
|
94 data_out3(m,n) = avg(data3(i)) |
|
95 data_out4(m,n) = avg(data4(i)) |
|
96 |
|
97 time_out(m,n) = year*100 + n + 1 |
|
98 ; print (time_out(m,n)) |
|
99 |
|
100 delete (i) |
|
101 end do |
|
102 if (m.lt.nyear-1) then |
|
103 delete (data1) |
|
104 delete (data2) |
|
105 delete (data3) |
|
106 delete (data4) |
|
107 delete (time) |
|
108 end if |
|
109 end do |
|
110 |
|
111 data_out1!0 = "year" |
|
112 data_out1!1 = "month" |
|
113 data_out1@long_name = data1@long_name |
|
114 data_out1@units = data1@units |
|
115 data_out1@_FillValue = 1.e+36 |
|
116 ; print (data_out1) |
|
117 |
|
118 data_out2!0 = "year" |
|
119 data_out2!1 = "month" |
|
120 data_out2@long_name = data2@long_name |
|
121 data_out2@units = data2@units |
|
122 data_out2@_FillValue = 1.e+36 |
|
123 ; print (data_out2) |
|
124 |
|
125 data_out3!0 = "year" |
|
126 data_out3!1 = "month" |
|
127 data_out3@long_name = data3@long_name |
|
128 data_out3@units = data3@units |
|
129 data_out3@_FillValue = 1.e+36 |
|
130 ; print (data_out3) |
|
131 |
|
132 data_out4!0 = "year" |
|
133 data_out4!1 = "month" |
|
134 data_out4@long_name = data4@long_name |
|
135 data_out4@units = data4@units |
|
136 data_out4@_FillValue = 1.e+36 |
|
137 ; print (data_out4) |
|
138 |
|
139 time_out!0 = "year" |
|
140 time_out!1 = "month" |
|
141 time_out@long_name = "current date as yyyymm" |
|
142 time_out@units = "current date as yyyymm" |
|
143 print (time_out) |
|
144 |
|
145 lat_out!0 = "lat" |
|
146 lon_out!0 = "lon" |
|
147 lat_out@units = "degrees_north" |
|
148 lat_out@long_name = "Latitude" |
|
149 lon_out@units = "degrees_east" |
|
150 lon_out@long_name = "Longitude" |
|
151 print (lat_out) |
|
152 print (lon_out) |
|
153 |
|
154 c->lat = lat_out |
|
155 c->lon = lon_out |
|
156 c->date = time_out |
|
157 c->CO2_FLUX = data_out1 |
|
158 c->RAD_FLUX = data_out2 |
|
159 c->SH_FLUX = data_out3 |
|
160 c->LH_FLUX = data_out4 |
|
161 |
|
162 end |