sfc_climo_gen 1.14.0
Loading...
Searching...
No Matches
output.f90
Go to the documentation of this file.
1
4
18 subroutine output(data_one_tile, lat_one_tile, lon_one_tile, i_mdl, j_mdl, &
19 tile, record, time, field_idx)
20
21 use mpi
22 use esmf
23 use netcdf
24 use utils
25 use source_grid, only : field_names, source, num_fields, &
26 num_time_recs, num_records, day_of_rec
27 use model_grid, only : missing, grid_tiles
28 use program_setup, only : halo
29
30 implicit none
31
32 integer, intent(in) :: i_mdl, j_mdl, tile
33 integer, intent(in) :: record, time, field_idx
34
35 real(esmf_kind_r4), intent(in) :: data_one_tile(i_mdl,j_mdl)
36 real(esmf_kind_r4), intent(in) :: lat_one_tile(i_mdl,j_mdl)
37 real(esmf_kind_r4), intent(in) :: lon_one_tile(i_mdl,j_mdl)
38
39 character(len=200) :: out_file
40 character(len=200) :: out_file_with_halo
41
42 integer :: initialsiz, fsize, error, j
43 integer :: dim_x, dim_y, id_data
44 integer :: dim_time, id_times, ierr
45 integer :: header_buffer_val = 16384
46 integer :: i_out, j_out, id_lat, id_lon
47 integer :: i_start, i_end, j_start, j_end
48 integer, save :: ncid(6), ncid_with_halo
49
50 select case (field_names(field_idx))
51 case ('substrate_temperature')
52 out_file = "./substrate_temperature." // grid_tiles(tile) // ".nc"
53 out_file_with_halo = "./substrate_temperature." // grid_tiles(tile) // ".halo.nc"
54 case ('vegetation_greenness')
55 out_file = "./vegetation_greenness." // grid_tiles(tile) // ".nc"
56 out_file_with_halo = "./vegetation_greenness." // grid_tiles(tile) // ".halo.nc"
57 case ('maximum_snow_albedo')
58 out_file = "./maximum_snow_albedo." // grid_tiles(tile) // ".nc"
59 out_file_with_halo = "./maximum_snow_albedo." // grid_tiles(tile) // ".halo.nc"
60 case ('leaf_area_index')
61 out_file = "./leaf_area_index." // grid_tiles(tile) // ".nc"
62 out_file_with_halo = "./leaf_area_index." // grid_tiles(tile) // ".halo.nc"
63 case ('visible_black_sky_albedo', 'visible_white_sky_albedo', 'near_IR_black_sky_albedo', 'near_IR_white_sky_albedo')
64 out_file = "./snowfree_albedo." // grid_tiles(tile) // ".nc"
65 out_file_with_halo = "./snowfree_albedo." // grid_tiles(tile) // ".halo.nc"
66 case ('facsf')
67 out_file = "./facsf." // grid_tiles(tile) // ".nc"
68 out_file_with_halo = "./facsf." // grid_tiles(tile) // ".halo.nc"
69 case ('slope_type')
70 out_file = "./slope_type." // grid_tiles(tile) // ".nc"
71 out_file_with_halo = "./slope_type." // grid_tiles(tile) // ".halo.nc"
72 case ('soil_type')
73 out_file = "./soil_type." // grid_tiles(tile) // ".nc"
74 out_file_with_halo = "./soil_type." // grid_tiles(tile) // ".halo.nc"
75 case ('soil_color')
76 out_file = "./soil_color." // grid_tiles(tile) // ".nc"
77 out_file_with_halo = "./soil_color." // grid_tiles(tile) // ".halo.nc"
78 case ('vegetation_type')
79 out_file = "./vegetation_type." // grid_tiles(tile) // ".nc"
80 out_file_with_halo = "./vegetation_type." // grid_tiles(tile) // ".halo.nc"
81 case default
82 print*,'- FATAL ERROR IN ROUTINE OUTPUT. UNIDENTIFIED FIELD : ', field_names(field_idx)
83 call mpi_abort(mpi_comm_world, 67, ierr)
84 end select
85
86!----------------------------------------------------------------------
87! If user specified a halo (for running stand-alone regional grid),
88! remove it.
89!----------------------------------------------------------------------
90
91 if (halo > 0) then
92 print*,"- WILL REMOVE HALO REGION OF ", halo, " ROWS/COLS."
93 i_start = 1 + halo
94 i_end = i_mdl - halo
95 j_start = 1 + halo
96 j_end = j_mdl - halo
97 else
98 i_start = 1
99 i_end = i_mdl
100 j_start = 1
101 j_end = j_mdl
102 endif
103
104 i_out = i_end - i_start + 1
105 j_out = j_end - j_start + 1
106
107 if (record == 1) then
108
109 initialsiz = 0
110 fsize = 65536
111 error = nf90_create(out_file, ior(nf90_netcdf4,nf90_classic_model), &
112 ncid(tile), initialsize=initialsiz, chunksize=fsize)
113 call netcdf_err(error, 'ERROR IN NF90_CREATE' )
114 error = nf90_def_dim(ncid(tile), 'nx', i_out, dim_x)
115 call netcdf_err(error, 'DEFINING NX DIMENSION' )
116 error = nf90_def_dim(ncid(tile), 'ny', j_out, dim_y)
117 call netcdf_err(error, 'DEFINING NY DIMENSION' )
118 error = nf90_def_dim(ncid(tile), 'time', num_time_recs, dim_time)
119 call netcdf_err(error, 'DEFINING TIME DIMENSION' )
120 error = nf90_def_var(ncid(tile), 'time', nf90_float, dim_time, id_times)
121 call netcdf_err(error, 'DEFINING TIME VARIABLE' )
122 error = nf90_put_att(ncid(tile), id_times, "units", "days since 2015-1-1")
123 call netcdf_err(error, 'DEFINING TIME ATTRIBUTE' )
124 if (len_trim(source) > 0) then
125 error = nf90_put_att(ncid(tile), nf90_global, 'source', source)
126 call netcdf_err(error, 'DEFINING GLOBAL SOURCE ATTRIBUTE' )
127 endif
128
129 error = nf90_def_var(ncid(tile), 'geolat', nf90_float, (/dim_x,dim_y/), id_lat)
130 call netcdf_err(error, 'DEFINING GEOLAT FIELD' )
131 error = nf90_put_att(ncid(tile), id_lat, "long_name", "Latitude")
132 call netcdf_err(error, 'DEFINING GEOLAT NAME ATTRIBUTE' )
133 error = nf90_put_att(ncid(tile), id_lat, "units", "degrees_north")
134 call netcdf_err(error, 'DEFINING GEOLAT UNIT ATTRIBUTE' )
135 error = nf90_def_var(ncid(tile), 'geolon', nf90_float, (/dim_x,dim_y/), id_lon)
136 call netcdf_err(error, 'DEFINING GEOLON FIELD' )
137 error = nf90_put_att(ncid(tile), id_lon, "long_name", "Longitude")
138 call netcdf_err(error, 'DEFINING GEOLON NAME ATTRIBUTE' )
139 error = nf90_put_att(ncid(tile), id_lon, "units", "degrees_east")
140 call netcdf_err(error, 'DEFINING GEOLON UNIT ATTRIBUTE' )
141
142 do j = 1, num_fields
143 error = nf90_def_var(ncid(tile), trim(field_names(j)), nf90_float, (/dim_x,dim_y,dim_time/), id_data)
144 call netcdf_err(error, 'DEFINING FIELD' )
145 error = nf90_put_att(ncid(tile), id_data, "missing_value", missing)
146 call netcdf_err(error, 'DEFINING FIELD ATTRIBUTE' )
147 error = nf90_put_att(ncid(tile), id_data, "coordinates", "geolon geolat")
148 call netcdf_err(error, 'DEFINING COORD ATTRIBUTE' )
149 enddo
150
151 error = nf90_enddef(ncid(tile), header_buffer_val,4,0,4)
152 call netcdf_err(error, 'IN NF90_ENDDEF' )
153
154 error = nf90_put_var( ncid(tile), id_times, day_of_rec)
155 call netcdf_err(error, 'WRITING TIME FIELD' )
156
157 error = nf90_put_var( ncid(tile), id_lat, lat_one_tile(i_start:i_end,j_start:j_end), &
158 start=(/1,1/), count=(/i_out,j_out/))
159 call netcdf_err(error, 'IN NF90_PUT_VAR FOR GEOLAT' )
160
161 error = nf90_put_var( ncid(tile), id_lon, lon_one_tile(i_start:i_end,j_start:j_end), &
162 start=(/1,1/), count=(/i_out,j_out/))
163 call netcdf_err(error, 'IN NF90_PUT_VAR FOR GEOLON' )
164
165 endif
166
167 print*,'- WRITE DATA FOR RECORD: ',record
168 error = nf90_inq_varid( ncid(tile), field_names(field_idx), id_data)
169 call netcdf_err(error, 'IN NF90_INQ_VARID' )
170 error = nf90_put_var( ncid(tile), id_data, data_one_tile(i_start:i_end,j_start:j_end), &
171 start=(/1,1,time/), count=(/i_out,j_out,1/))
172 call netcdf_err(error, 'IN NF90_PUT_VAR' )
173
174 if (record == num_records) then
175 error = nf90_close(ncid(tile))
176 endif
177
178!----------------------------------------------------------------------
179! For regional nests, also output files including the halo
180!----------------------------------------------------------------------
181
182 if (halo == 0) return
183
184 print*,"- WRITE OUT FILES THAT INCLUDE HALO REGION."
185
186 if (record == 1) then
187
188 initialsiz = 0
189 fsize = 65536
190 error = nf90_create(out_file_with_halo, ior(nf90_netcdf4,nf90_classic_model), &
191 ncid_with_halo, initialsize=initialsiz, chunksize=fsize)
192 call netcdf_err(error, 'IN NF90_CREATE' )
193 error = nf90_def_dim(ncid_with_halo, 'nx', i_mdl, dim_x)
194 call netcdf_err(error, 'DEFINING NX DIMENSION' )
195 error = nf90_def_dim(ncid_with_halo, 'ny', j_mdl, dim_y)
196 call netcdf_err(error, 'DEFINING NY DIMENSION' )
197 error = nf90_def_dim(ncid_with_halo, 'time', num_time_recs, dim_time)
198 call netcdf_err(error, 'DEFINING TIME DIMENSION' )
199 error = nf90_def_var(ncid_with_halo, 'time', nf90_float, dim_time, id_times)
200 call netcdf_err(error, 'DEFINING TIME VARIABLE' )
201 error = nf90_put_att(ncid_with_halo, id_times, "units", "days since 2015-1-1")
202 call netcdf_err(error, 'DEFINING TIME ATTRIBUTE' )
203 if (len_trim(source) > 0) then
204 error = nf90_put_att(ncid_with_halo, nf90_global, 'source', source)
205 call netcdf_err(error, 'DEFINING GLOBAL SOURCE ATTRIBUTE' )
206 endif
207
208 error = nf90_def_var(ncid_with_halo, 'geolat', nf90_float, (/dim_x,dim_y/), id_lat)
209 call netcdf_err(error, 'DEFINING GEOLAT FIELD' )
210 error = nf90_put_att(ncid_with_halo, id_lat, "long_name", "Latitude")
211 call netcdf_err(error, 'DEFINING GEOLAT NAME ATTRIBUTE' )
212 error = nf90_put_att(ncid_with_halo, id_lat, "units", "degrees_north")
213 call netcdf_err(error, 'DEFINING GEOLAT UNIT ATTRIBUTE' )
214 error = nf90_def_var(ncid_with_halo, 'geolon', nf90_float, (/dim_x,dim_y/), id_lon)
215 call netcdf_err(error, 'DEFINING GEOLON FIELD' )
216 error = nf90_put_att(ncid_with_halo, id_lon, "long_name", "Longitude")
217 call netcdf_err(error, 'DEFINING GEOLON NAME ATTRIBUTE' )
218 error = nf90_put_att(ncid_with_halo, id_lon, "units", "degrees_east")
219 call netcdf_err(error, 'DEFINING GEOLON UNIT ATTRIBUTE' )
220
221 do j = 1, num_fields
222 error = nf90_def_var(ncid_with_halo, field_names(j), nf90_float, (/dim_x,dim_y,dim_time/), id_data)
223 call netcdf_err(error, 'DEFINING FIELD VARIABLE' )
224 error = nf90_put_att(ncid_with_halo, id_data, "missing_value", missing)
225 call netcdf_err(error, 'DEFINING FIELD ATTRIBUTE' )
226 error = nf90_put_att(ncid_with_halo, id_data, "coordinates", "geolon geolat")
227 call netcdf_err(error, 'DEFINING COORD ATTRIBUTE' )
228 enddo
229
230 error = nf90_enddef(ncid_with_halo, header_buffer_val,4,0,4)
231 call netcdf_err(error, 'WRITING HEADER ENDDEF' )
232
233 error = nf90_put_var(ncid_with_halo, id_times, day_of_rec)
234 call netcdf_err(error, 'WRITING TIME VARIABLE' )
235
236 error = nf90_put_var( ncid_with_halo, id_lat, lat_one_tile, &
237 start=(/1,1/), count=(/i_mdl,j_mdl/))
238 call netcdf_err(error, 'IN NF90_PUT_VAR FOR GEOLAT' )
239
240 error = nf90_put_var( ncid_with_halo, id_lon, lon_one_tile, &
241 start=(/1,1/), count=(/i_mdl,j_mdl/))
242 call netcdf_err(error, 'IN NF90_PUT_VAR FOR GEOLON' )
243
244 endif
245
246 print*,'- WRITE DATA FOR RECORD: ',record
247 error = nf90_inq_varid(ncid_with_halo, field_names(field_idx), id_data)
248 call netcdf_err(error, 'IN NF90_INQ_VARID' )
249
250 error = nf90_put_var(ncid_with_halo, id_data, data_one_tile, &
251 start=(/1,1,time/), count=(/i_mdl,j_mdl,1/))
252 call netcdf_err(error, 'IN NF90_PUT_VAR' )
253
254 if (record == num_records) then
255 error = nf90_close(ncid_with_halo)
256 endif
257
258 return
259
260 end subroutine output
model_grid
This module defines the model grid.
Definition model_grid.F90:10
model_grid::missing
real(kind=4), public missing
Value assigned to undefined points (i.e., ocean points for a land field).
Definition model_grid.F90:25
model_grid::grid_tiles
character(len=5), dimension(:), allocatable, public grid_tiles
Array of model grid tile names.
Definition model_grid.F90:18
program_setup
Set up program execution.
Definition program_setup.f90:17
program_setup::halo
integer, public halo
Number of row/cols defining the lateral boundary halo.
Definition program_setup.f90:48
source_grid
Read grid specs, date information and land/sea mask for the source data that will be interpolated to ...
Definition source_grid.F90:12
source_grid::source
character(len=75), public source
Original source of the data.
Definition source_grid.F90:22
source_grid::num_time_recs
integer, public num_time_recs
Number of time records.
Definition source_grid.F90:27
source_grid::num_records
integer, public num_records
Number of fields times time records.
Definition source_grid.F90:26
source_grid::field_names
character(len=50), dimension(:), allocatable, public field_names
Names of fields to be processed.
Definition source_grid.F90:21
source_grid::num_fields
integer, public num_fields
Number of fields in the file.
Definition source_grid.F90:28
source_grid::day_of_rec
integer, dimension(:), allocatable, public day_of_rec
Day of each time record with respect to Jan 1.
Definition source_grid.F90:31
utils
Utilities.
Definition utils.f90:8
utils::netcdf_err
subroutine, public netcdf_err(err, string)
Handle netCDF error codes.
Definition utils.f90:23
output
subroutine output(data_one_tile, lat_one_tile, lon_one_tile, i_mdl, j_mdl, tile, record, time, field_idx)
Output model data for a single tile and a single record in netcdf format.
Definition output.f90:20