Context Navigation

output_wrf.F @ 2749

Last change on this file since 2749 was 2747, checked in by aslmd, 3 years ago
adding unmodified code from WRFV3.0.1.1, expurged from useless data +1M size
File size: 32.8 KB

Line
1	!WRF:MEDIATION:IO
2	! ---principal wrf output routine (called from routines in module_io_domain )
3	SUBROUTINE output_wrf ( fid , grid , config_flags, switch , ierr )
4	USE module_io
5	USE module_wrf_error
6	USE module_io_wrf
7	USE module_domain
8	USE module_state_description
9	USE module_configure
10	! USE module_date_time
11	USE module_model_constants
12	USE module_utility
13	IMPLICIT NONE
14	#include <wrf_io_flags.h>
15	#include <wrf_status_codes.h>
16	TYPE(domain) :: grid
17	TYPE(grid_config_rec_type), INTENT(INOUT) :: config_flags
18	INTEGER, INTENT(IN) :: fid, switch
19	INTEGER, INTENT(INOUT) :: ierr
20
21	! Local data
22	INTEGER ids , ide , jds , jde , kds , kde , &
23	ims , ime , jms , jme , kms , kme , &
24	ips , ipe , jps , jpe , kps , kpe
25
26	INTEGER , DIMENSION(3) :: domain_start , domain_end
27	INTEGER , DIMENSION(3) :: memory_start , memory_end
28	INTEGER , DIMENSION(3) :: patch_start , patch_end
29	INTEGER i,j
30	INTEGER julyr, julday, idt, iswater , map_proj
31	INTEGER filestate
32	LOGICAL dryrun
33	REAL gmt, cen_lat, cen_lon, bdyfrq , truelat1 , truelat2 , moad_cen_lat , stand_lon
34	INTEGER km_opt, diff_opt, damp_opt, &
35	mp_physics, ra_lw_physics, ra_sw_physics, sf_sfclay_physics, &
36	sf_surface_physics, bl_pbl_physics, cu_physics
37	REAL khdif, kvdif, swrad_scat
38	INTEGER ucmcall, w_damping, smooth_option, feedback, surface_input_source, sst_update
39	#if (EM_CORE == 1)
40	INTEGER grid_id , parent_id , i_parent_start , j_parent_start , parent_grid_ratio
41	INTEGER diff_6th_opt
42	REAL diff_6th_factor
43	INTEGER grid_fdda, gfdda_interval_m, gfdda_end_h, if_ramping, &
44	obs_nudge_opt, obs_nudge_wind, obs_nudge_temp, obs_nudge_mois, obs_nudge_pstr, obs_idynin, obs_ionf
45	REAL fgdt, guv, gt, gq, dtramp_min, &
46	obs_coef_wind, obs_coef_temp, obs_coef_mois, obs_coef_pstr, obs_dtramp, fdda_end
47	LOGICAL pd_moist, pd_scalar, pd_tke
48	#endif
49	CHARACTER (len=19) simulation_start_date
50	CHARACTER (len=len_current_date) current_date_save
51	INTEGER simulation_start_year , &
52	simulation_start_month , &
53	simulation_start_day , &
54	simulation_start_hour , &
55	simulation_start_minute , &
56	simulation_start_second
57	INTEGER rc
58	INTEGER :: io_form
59	LOGICAL, EXTERNAL :: multi_files
60	INTEGER, EXTERNAL :: use_package
61	INTEGER p_hr, p_min, p_sec, p_ms
62
63	CHARACTER*256 message
64	CHARACTER*80 fname
65	CHARACTER*80 char_junk
66	INTEGER ibuf(1)
67	REAL rbuf(1)
68	TYPE(WRFU_TimeInterval) :: bdy_increment
69	TYPE(WRFU_Time) :: next_time, currentTime, startTime
70	CHARACTER*40 :: next_datestr
71	INTEGER :: start_year , start_month , start_day , start_hour , start_minute , start_second
72	LOGICAL :: adjust
73
74	WRITE(wrf_err_message,*)'output_wrf: begin, fid = ',fid
75	CALL wrf_debug( 300 , wrf_err_message )
76
77	CALL wrf_inquire_filename ( fid , fname , filestate , ierr )
78	IF ( ierr /= 0 ) THEN
79	WRITE(wrf_err_message,*)'module_io_wrf: output_wrf: wrf_inquire_filename Status = ',ierr
80	CALL wrf_error_fatal( wrf_err_message )
81	ENDIF
82
83	WRITE(wrf_err_message,*)'output_wrf: fid,filestate = ',fid,filestate
84	CALL wrf_debug( 300 , wrf_err_message )
85
86	! io_form is used to determine if multi-file I/O is enabled and to
87	! control writing of format-specific time-independent metadata
88	IF ( switch .EQ. model_input_only ) THEN
89	CALL nl_get_io_form_input( 1, io_form )
90	ELSE IF ( switch .EQ. aux_model_input1_only ) THEN
91	CALL nl_get_io_form_auxinput1( 1, io_form )
92	ELSE IF ( switch .EQ. aux_model_input2_only ) THEN
93	CALL nl_get_io_form_auxinput2( 1, io_form )
94	ELSE IF ( switch .EQ. aux_model_input3_only ) THEN
95	CALL nl_get_io_form_auxinput3( 1, io_form )
96	ELSE IF ( switch .EQ. aux_model_input4_only ) THEN
97	CALL nl_get_io_form_auxinput4( 1, io_form )
98	ELSE IF ( switch .EQ. aux_model_input5_only ) THEN
99	CALL nl_get_io_form_auxinput5( 1, io_form )
100	ELSE IF ( switch .EQ. aux_model_input6_only ) THEN
101	CALL nl_get_io_form_auxinput6( 1, io_form )
102	ELSE IF ( switch .EQ. aux_model_input7_only ) THEN
103	CALL nl_get_io_form_auxinput7( 1, io_form )
104	ELSE IF ( switch .EQ. aux_model_input8_only ) THEN
105	CALL nl_get_io_form_auxinput8( 1, io_form )
106	ELSE IF ( switch .EQ. aux_model_input9_only ) THEN
107	CALL nl_get_io_form_auxinput9( 1, io_form )
108	ELSE IF ( switch .EQ. aux_model_input10_only ) THEN
109	CALL nl_get_io_form_gfdda( 1, io_form )
110	ELSE IF ( switch .EQ. aux_model_input11_only ) THEN
111	CALL nl_get_io_form_auxinput11( 1, io_form )
112
113	ELSE IF ( switch .EQ. history_only ) THEN
114	CALL nl_get_io_form_history( 1, io_form )
115	ELSE IF ( switch .EQ. aux_hist1_only ) THEN
116	CALL nl_get_io_form_auxhist1( 1, io_form )
117	ELSE IF ( switch .EQ. aux_hist2_only ) THEN
118	CALL nl_get_io_form_auxhist2( 1, io_form )
119	ELSE IF ( switch .EQ. aux_hist3_only ) THEN
120	CALL nl_get_io_form_auxhist3( 1, io_form )
121	ELSE IF ( switch .EQ. aux_hist4_only ) THEN
122	CALL nl_get_io_form_auxhist4( 1, io_form )
123	ELSE IF ( switch .EQ. aux_hist5_only ) THEN
124	CALL nl_get_io_form_auxhist5( 1, io_form )
125	ELSE IF ( switch .EQ. aux_hist6_only ) THEN
126	CALL nl_get_io_form_auxhist6( 1, io_form )
127	ELSE IF ( switch .EQ. aux_hist7_only ) THEN
128	CALL nl_get_io_form_auxhist7( 1, io_form )
129	ELSE IF ( switch .EQ. aux_hist8_only ) THEN
130	CALL nl_get_io_form_auxhist8( 1, io_form )
131	ELSE IF ( switch .EQ. aux_hist9_only ) THEN
132	CALL nl_get_io_form_auxhist9( 1, io_form )
133	ELSE IF ( switch .EQ. aux_hist10_only ) THEN
134	CALL nl_get_io_form_auxhist10( 1, io_form )
135	ELSE IF ( switch .EQ. aux_hist11_only ) THEN
136	CALL nl_get_io_form_auxhist11( 1, io_form )
137
138	ELSE IF ( switch .EQ. restart_only ) THEN
139	CALL nl_get_io_form_restart( 1, io_form )
140	ELSE IF ( switch .EQ. boundary_only ) THEN
141	CALL nl_get_io_form_boundary( 1, io_form )
142	ELSE ! default: use history
143	CALL nl_get_io_form_history( 1, io_form )
144	ENDIF
145
146	dryrun = ( filestate .EQ. WRF_FILE_OPENED_NOT_COMMITTED )
147
148	WRITE(wrf_err_message,*)'output_wrf: dryrun = ',dryrun
149	CALL wrf_debug( 300 , wrf_err_message )
150
151	CALL get_ijk_from_grid ( grid , &
152	ids, ide, jds, jde, kds, kde, &
153	ims, ime, jms, jme, kms, kme, &
154	ips, ipe, jps, jpe, kps, kpe )
155
156	call nl_get_diff_opt ( 1, diff_opt )
157	call nl_get_km_opt ( 1, km_opt )
158	call nl_get_damp_opt ( 1, damp_opt )
159	call nl_get_khdif ( grid%id, khdif )
160	call nl_get_kvdif ( grid%id, kvdif )
161	call nl_get_mp_physics ( grid%id, mp_physics )
162	call nl_get_ra_lw_physics ( grid%id, ra_lw_physics )
163	call nl_get_ra_sw_physics ( grid%id, ra_sw_physics )
164	call nl_get_sf_sfclay_physics ( grid%id, sf_sfclay_physics )
165	call nl_get_sf_surface_physics ( grid%id, sf_surface_physics )
166	call nl_get_bl_pbl_physics ( grid%id, bl_pbl_physics )
167	call nl_get_cu_physics ( grid%id, cu_physics )
168
169	! add nml variables in 2.2
170	call nl_get_surface_input_source ( 1 , surface_input_source )
171	call nl_get_sst_update ( 1 , sst_update )
172	call nl_get_feedback ( 1 , feedback )
173	call nl_get_smooth_option ( 1 , smooth_option )
174	call nl_get_swrad_scat ( 1 , swrad_scat )
175	call nl_get_ucmcall ( 1 , ucmcall )
176	call nl_get_w_damping ( 1 , w_damping )
177
178	#if (EM_CORE == 1)
179	CALL nl_get_pd_moist ( grid%id , pd_moist )
180	CALL nl_get_pd_scalar ( grid%id , pd_scalar )
181	CALL nl_get_pd_tke ( grid%id , pd_tke )
182	CALL nl_get_diff_6th_opt ( grid%id , diff_6th_opt )
183	CALL nl_get_diff_6th_factor ( grid%id , diff_6th_factor )
184	CALL nl_get_grid_fdda ( grid%id , grid_fdda )
185	CALL nl_get_gfdda_end_h( grid%id , gfdda_end_h )
186	CALL nl_get_gfdda_interval_m ( grid%id , gfdda_interval_m )
187
188	IF ( grid_fdda == 1 ) THEN
189	CALL nl_get_fgdt ( grid%id , fgdt )
190	CALL nl_get_guv ( grid%id , guv )
191	CALL nl_get_gt ( grid%id , gt )
192	CALL nl_get_gq ( grid%id , gq )
193	CALL nl_get_if_ramping ( 1 , if_ramping )
194	CALL nl_get_dtramp_min ( 1 , dtramp_min )
195	ENDIF
196
197	CALL nl_get_obs_nudge_opt ( grid%id , obs_nudge_opt )
198	IF ( obs_nudge_opt == 1 ) THEN
199	CALL nl_get_fdda_end ( grid%id , fdda_end )
200	CALL nl_get_obs_nudge_wind ( grid%id , obs_nudge_wind )
201	CALL nl_get_obs_coef_wind ( grid%id , obs_coef_wind )
202	CALL nl_get_obs_nudge_temp ( grid%id , obs_nudge_temp )
203	CALL nl_get_obs_coef_temp ( grid%id , obs_coef_temp )
204	CALL nl_get_obs_nudge_mois ( grid%id , obs_nudge_mois )
205	CALL nl_get_obs_coef_mois ( grid%id , obs_coef_mois )
206	CALL nl_get_obs_nudge_pstr ( grid%id , obs_nudge_pstr )
207	CALL nl_get_obs_coef_pstr ( grid%id , obs_coef_pstr )
208	CALL nl_get_obs_ionf ( 1 , obs_ionf )
209	CALL nl_get_obs_idynin ( 1 , obs_idynin )
210	CALL nl_get_obs_dtramp ( 1 , obs_dtramp )
211	ENDIF
212	#endif
213
214	! julday and gmt can be set in namelist_03 for ideal.exe run
215	CALL nl_get_gmt (grid%id, gmt)
216	CALL nl_get_julyr (grid%id, julyr)
217	CALL nl_get_julday (grid%id, julday)
218	CALL nl_get_mminlu ( 1, char_junk(1:4) )
219	CALL nl_get_iswater (grid%id, iswater )
220	CALL nl_get_cen_lat ( grid%id , cen_lat )
221	CALL nl_get_cen_lon ( grid%id , cen_lon )
222	CALL nl_get_truelat1 ( grid%id , truelat1 )
223	CALL nl_get_truelat2 ( grid%id , truelat2 )
224	CALL nl_get_moad_cen_lat ( grid%id , moad_cen_lat )
225	CALL nl_get_stand_lon ( grid%id , stand_lon )
226	CALL nl_get_map_proj ( grid%id , map_proj )
227
228	#if (EM_CORE == 1)
229	CALL nl_get_parent_id ( grid%id , parent_id )
230	CALL nl_get_i_parent_start ( grid%id , i_parent_start )
231	CALL nl_get_j_parent_start ( grid%id , j_parent_start )
232	CALL nl_get_parent_grid_ratio ( grid%id , parent_grid_ratio )
233	#endif
234
235	CALL domain_clockprint(150, grid, &
236	'DEBUG output_wrf(): before call to domain_clock_get,')
237	CALL domain_clock_get( grid, current_time=currentTime, &
238	start_time=startTime, &
239	current_timestr=current_date )
240	WRITE ( wrf_err_message , * ) 'output_wrf: begin, current_date=',current_date
241	CALL wrf_debug ( 300 , wrf_err_message )
242
243	WRITE( message , * ) "OUTPUT FROM " , TRIM(program_name)
244	CALL wrf_put_dom_ti_char ( fid , 'TITLE' , TRIM(message) , ierr )
245	! added grib-specific metadata: Todd Hutchinson 8/21/2005
246	IF ( ( use_package( io_form ) == IO_GRIB1 ) .OR. &
247	( use_package( io_form ) == IO_GRIB2 ) ) THEN
248	CALL wrf_put_dom_ti_char ( fid, 'PROGRAM_NAME', TRIM(program_name) , ierr )
249	ENDIF
250	CALL nl_get_start_year(grid%id,start_year)
251	CALL nl_get_start_month(grid%id,start_month)
252	CALL nl_get_start_day(grid%id,start_day)
253	CALL nl_get_start_hour(grid%id,start_hour)
254	CALL nl_get_start_minute(grid%id,start_minute)
255	CALL nl_get_start_second(grid%id,start_second)
256	#ifdef PLANET
257	WRITE ( start_date , FMT = '(I4.4,"-",I5.5,"_",I2.2,":",I2.2,":",I2.2)' ) &
258	start_year,start_day,start_hour,start_minute,start_second
259	#else
260	WRITE ( start_date , FMT = '(I4.4,"-",I2.2,"-",I2.2,"_",I2.2,":",I2.2,":",I2.2)' ) &
261	start_year,start_month,start_day,start_hour,start_minute,start_second
262	#endif
263	CALL wrf_put_dom_ti_char ( fid , 'START_DATE', TRIM(start_date) , ierr )
264	IF ( switch .EQ. model_input_only) THEN
265	CALL wrf_put_dom_ti_char ( fid , 'SIMULATION_START_DATE', TRIM(start_date) , ierr )
266	ELSE IF ( ( switch .EQ. restart_only ) .OR. ( switch .EQ. history_only ) ) THEN
267	CALL nl_get_simulation_start_year ( 1, simulation_start_year )
268	CALL nl_get_simulation_start_month ( 1, simulation_start_month )
269	CALL nl_get_simulation_start_day ( 1, simulation_start_day )
270	CALL nl_get_simulation_start_hour ( 1, simulation_start_hour )
271	CALL nl_get_simulation_start_minute ( 1, simulation_start_minute )
272	CALL nl_get_simulation_start_second ( 1, simulation_start_second )
273	#ifdef PLANET
274	WRITE ( simulation_start_date , FMT = '(I4.4,"-",I5.5,"_",I2.2,":",I2.2,":",I2.2)' ) &
275	simulation_start_year,simulation_start_day,&
276	simulation_start_hour,simulation_start_minute,simulation_start_second
277	#else
278	WRITE ( simulation_start_date , FMT = '(I4.4,"-",I2.2,"-",I2.2,"_",I2.2,":",I2.2,":",I2.2)' ) &
279	simulation_start_year,simulation_start_month,simulation_start_day,&
280	simulation_start_hour,simulation_start_minute,simulation_start_second
281	#endif
282	CALL wrf_put_dom_ti_char ( fid , 'SIMULATION_START_DATE', TRIM(simulation_start_date) , ierr )
283	END IF
284
285	ibuf(1) = config_flags%e_we - config_flags%s_we + 1
286	CALL wrf_put_dom_ti_integer ( fid , 'WEST-EAST_GRID_DIMENSION' , ibuf , 1 , ierr )
287
288	ibuf(1) = config_flags%e_sn - config_flags%s_sn + 1
289	CALL wrf_put_dom_ti_integer ( fid , 'SOUTH-NORTH_GRID_DIMENSION' , ibuf , 1 , ierr )
290
291	ibuf(1) = config_flags%e_vert - config_flags%s_vert + 1
292	CALL wrf_put_dom_ti_integer ( fid , 'BOTTOM-TOP_GRID_DIMENSION' , ibuf , 1 , ierr )
293
294	#if (EM_CORE == 1)
295	IF (grid%map_proj == 6) THEN
296	! Global ... dx determined automatically
297	! Don't use value from namelist ... used derived value instead
298	CALL wrf_put_dom_ti_real ( fid , 'DX' , grid%dx , 1 , ierr )
299	CALL wrf_put_dom_ti_real ( fid , 'DY' , grid%dy , 1 , ierr )
300	ELSE
301	CALL wrf_put_dom_ti_real ( fid , 'DX' , config_flags%dx , 1 , ierr )
302	CALL wrf_put_dom_ti_real ( fid , 'DY' , config_flags%dy , 1 , ierr )
303	END IF
304	#endif
305
306	! added this metadatum for H. Chuan, NCEP, 030417, JM
307	#if (NMM_CORE == 1)
308	CALL wrf_put_dom_ti_char ( fid , 'GRIDTYPE', 'E' , ierr )
309	#endif
310	#if (EM_CORE == 1)
311	CALL wrf_put_dom_ti_char ( fid , 'GRIDTYPE', 'C' , ierr )
312	#endif
313	#if (COAMPS_CORE == 1 )
314	CALL wrf_put_dom_ti_char ( fid , 'GRIDTYPE', 'B' , ierr )
315	#endif
316
317	! added these fields for W. Skamarock, 020402, JM
318	ibuf(1) = diff_opt
319	CALL wrf_put_dom_ti_integer ( fid , 'DIFF_OPT' , ibuf , 1 , ierr )
320	ibuf(1) = km_opt
321	CALL wrf_put_dom_ti_integer ( fid , 'KM_OPT' , ibuf , 1 , ierr )
322	ibuf(1) = damp_opt
323	CALL wrf_put_dom_ti_integer ( fid , 'DAMP_OPT' , ibuf , 1 , ierr )
324	rbuf(1) = khdif
325	CALL wrf_put_dom_ti_real ( fid , 'KHDIF' , rbuf , 1 , ierr )
326	rbuf(1) = kvdif
327	CALL wrf_put_dom_ti_real ( fid , 'KVDIF' , rbuf , 1 , ierr )
328	ibuf(1) = mp_physics
329	CALL wrf_put_dom_ti_integer ( fid , 'MP_PHYSICS' , ibuf , 1 , ierr )
330	ibuf(1) = ra_lw_physics
331	CALL wrf_put_dom_ti_integer ( fid , 'RA_LW_PHYSICS' , ibuf , 1 , ierr )
332	ibuf(1) = ra_sw_physics
333	CALL wrf_put_dom_ti_integer ( fid , 'RA_SW_PHYSICS' , ibuf , 1 , ierr )
334	ibuf(1) = sf_sfclay_physics
335	CALL wrf_put_dom_ti_integer ( fid , 'SF_SFCLAY_PHYSICS' , ibuf , 1 , ierr )
336	ibuf(1) = sf_surface_physics
337	CALL wrf_put_dom_ti_integer ( fid , 'SF_SURFACE_PHYSICS' , ibuf , 1 , ierr )
338	ibuf(1) = bl_pbl_physics
339	CALL wrf_put_dom_ti_integer ( fid , 'BL_PBL_PHYSICS' , ibuf , 1 , ierr )
340	ibuf(1) = cu_physics
341	CALL wrf_put_dom_ti_integer ( fid , 'CU_PHYSICS' , ibuf , 1 , ierr )
342
343	! added netcdf-specific metadata:
344	IF ( ( use_package( io_form ) == IO_NETCDF ) .OR. &
345	( use_package( io_form ) == IO_PHDF5 ) .OR. &
346	( use_package( io_form ) == IO_PNETCDF ) ) THEN
347	CALL wrf_put_dom_ti_integer ( fid, 'SURFACE_INPUT_SOURCE', surface_input_source , 1 , ierr )
348	CALL wrf_put_dom_ti_integer ( fid, 'SST_UPDATE', sst_update , 1 , ierr )
349	#if (EM_CORE == 1)
350	CALL wrf_put_dom_ti_integer ( fid, 'GRID_FDDA', grid_fdda , 1 , ierr )
351	CALL wrf_put_dom_ti_integer ( fid, 'GFDDA_INTERVAL_M', gfdda_interval_m , 1 , ierr )
352	CALL wrf_put_dom_ti_integer ( fid, 'GFDDA_END_H', gfdda_end_h , 1 , ierr )
353	#endif
354
355	IF ( switch .EQ. history_only ) THEN
356	CALL wrf_put_dom_ti_integer ( fid, 'UCMCALL', ucmcall , 1 , ierr )
357	CALL wrf_put_dom_ti_integer ( fid, 'FEEDBACK', feedback , 1 , ierr )
358	CALL wrf_put_dom_ti_integer ( fid, 'SMOOTH_OPTION', smooth_option , 1 , ierr )
359	CALL wrf_put_dom_ti_real ( fid, 'SWRAD_SCAT', swrad_scat , 1 , ierr )
360	CALL wrf_put_dom_ti_integer ( fid, 'W_DAMPING', w_damping , 1 , ierr )
361
362	#if (EM_CORE == 1)
363	CALL wrf_put_dom_ti_logical ( fid, 'PD_MOIST', pd_moist , 1 , ierr )
364	CALL wrf_put_dom_ti_logical ( fid, 'PD_SCALAR', pd_scalar , 1 , ierr )
365	CALL wrf_put_dom_ti_logical ( fid, 'PD_TKE', pd_tke , 1 , ierr )
366	CALL wrf_put_dom_ti_integer ( fid, 'DIFF_6TH_OPT', diff_6th_opt , 1 , ierr )
367	CALL wrf_put_dom_ti_real ( fid, 'DIFF_6TH_FACTOR', diff_6th_factor , 1 , ierr )
368
369	IF ( grid_fdda == 1 ) THEN
370	CALL wrf_put_dom_ti_real ( fid, 'FGDT', fgdt , 1 , ierr )
371	CALL wrf_put_dom_ti_real ( fid, 'GUV', guv , 1 , ierr )
372	CALL wrf_put_dom_ti_real ( fid, 'GT', gt , 1 , ierr )
373	CALL wrf_put_dom_ti_real ( fid, 'GQ', gq , 1 , ierr )
374	CALL wrf_put_dom_ti_integer ( fid, 'IF_RAMPING', if_ramping , 1 , ierr )
375	CALL wrf_put_dom_ti_real ( fid, 'DTRAMP_MIN', dtramp_min , 1 , ierr )
376	ENDIF
377
378	CALL wrf_put_dom_ti_integer ( fid, 'OBS_NUDGE_OPT', obs_nudge_opt , 1 , ierr )
379	IF ( obs_nudge_opt == 1 ) THEN
380	CALL wrf_put_dom_ti_real ( fid, 'FDDA_END', fdda_end , 1 , ierr )
381	CALL wrf_put_dom_ti_integer ( fid, 'OBS_NUDGE_WIND', obs_nudge_wind , 1 , ierr )
382	CALL wrf_put_dom_ti_real ( fid, 'OBS_COEF_WIND', obs_coef_wind , 1 , ierr )
383	CALL wrf_put_dom_ti_integer ( fid, 'OBS_NUDGE_TEMP', obs_nudge_temp , 1 , ierr )
384	CALL wrf_put_dom_ti_real ( fid, 'OBS_COEF_TEMP', obs_coef_temp , 1 , ierr )
385	CALL wrf_put_dom_ti_integer ( fid, 'OBS_NUDGE_MOIS', obs_nudge_mois , 1 , ierr )
386	CALL wrf_put_dom_ti_real ( fid, 'OBS_COEF_MOIS', obs_coef_mois , 1 , ierr )
387	CALL wrf_put_dom_ti_integer ( fid, 'OBS_NUDGE_PSTR', obs_nudge_pstr , 1 , ierr )
388	CALL wrf_put_dom_ti_real ( fid, 'OBS_COEF_PSTR', obs_coef_pstr , 1 , ierr )
389	CALL wrf_put_dom_ti_integer ( fid, 'OBS_IONF', obs_ionf , 1 , ierr )
390	CALL wrf_put_dom_ti_integer ( fid, 'OBS_IDYNIN', obs_idynin , 1 , ierr )
391	CALL wrf_put_dom_ti_real ( fid, 'OBS_DTRAMP', obs_dtramp , 1 , ierr )
392	ENDIF
393	#endif
394	ENDIF ! history_only
395	ENDIF
396
397	! added these fields for use by reassembly programs , 010831, JM
398	! modified these fields so "patch" == "domain" when multi-file output
399	! formats are not used. 051018, TBH
400
401	ibuf(1) = MAX(ips,ids)
402	IF ( .NOT. multi_files ( io_form ) ) ibuf(1) = ids
403	CALL wrf_put_dom_ti_integer ( fid , 'WEST-EAST_PATCH_START_UNSTAG' , ibuf , 1 , ierr )
404	ibuf(1) = MIN(ipe,ide-1)
405	IF ( .NOT. multi_files ( io_form ) ) ibuf(1) = ide - 1
406	CALL wrf_put_dom_ti_integer ( fid , 'WEST-EAST_PATCH_END_UNSTAG' , ibuf , 1 , ierr )
407	ibuf(1) = MAX(ips,ids)
408	IF ( .NOT. multi_files ( io_form ) ) ibuf(1) = ids
409	CALL wrf_put_dom_ti_integer ( fid , 'WEST-EAST_PATCH_START_STAG' , ibuf , 1 , ierr )
410	ibuf(1) = MIN(ipe,ide)
411	IF ( .NOT. multi_files ( io_form ) ) ibuf(1) = ide
412	CALL wrf_put_dom_ti_integer ( fid , 'WEST-EAST_PATCH_END_STAG' , ibuf , 1 , ierr )
413	ibuf(1) = MAX(jps,jds)
414	IF ( .NOT. multi_files ( io_form ) ) ibuf(1) = jds
415	CALL wrf_put_dom_ti_integer ( fid , 'SOUTH-NORTH_PATCH_START_UNSTAG' , ibuf , 1 , ierr )
416	ibuf(1) = MIN(jpe,jde-1)
417	IF ( .NOT. multi_files ( io_form ) ) ibuf(1) = jde - 1
418	CALL wrf_put_dom_ti_integer ( fid , 'SOUTH-NORTH_PATCH_END_UNSTAG' , ibuf , 1 , ierr )
419	ibuf(1) = MAX(jps,jds)
420	IF ( .NOT. multi_files ( io_form ) ) ibuf(1) = jds
421	CALL wrf_put_dom_ti_integer ( fid , 'SOUTH-NORTH_PATCH_START_STAG' , ibuf , 1 , ierr )
422	ibuf(1) = MIN(jpe,jde)
423	IF ( .NOT. multi_files ( io_form ) ) ibuf(1) = jde
424	CALL wrf_put_dom_ti_integer ( fid , 'SOUTH-NORTH_PATCH_END_STAG' , ibuf , 1 , ierr )
425
426	ibuf(1) = MAX(kps,kds)
427	IF ( .NOT. multi_files ( io_form ) ) ibuf(1) = kds
428	CALL wrf_put_dom_ti_integer ( fid , 'BOTTOM-TOP_PATCH_START_UNSTAG' , ibuf , 1 , ierr )
429	ibuf(1) = MIN(kpe,kde-1)
430	IF ( .NOT. multi_files ( io_form ) ) ibuf(1) = kde - 1
431	CALL wrf_put_dom_ti_integer ( fid , 'BOTTOM-TOP_PATCH_END_UNSTAG' , ibuf , 1 , ierr )
432	ibuf(1) = MAX(kps,kds)
433	IF ( .NOT. multi_files ( io_form ) ) ibuf(1) = kds
434	CALL wrf_put_dom_ti_integer ( fid , 'BOTTOM-TOP_PATCH_START_STAG' , ibuf , 1 , ierr )
435	ibuf(1) = MIN(kpe,kde)
436	IF ( .NOT. multi_files ( io_form ) ) ibuf(1) = kde
437	CALL wrf_put_dom_ti_integer ( fid , 'BOTTOM-TOP_PATCH_END_STAG' , ibuf , 1 , ierr )
438	#if (EM_CORE == 1)
439	ibuf(1) = grid%id
440	CALL wrf_put_dom_ti_integer ( fid , 'GRID_ID' , ibuf , 1 , ierr )
441	ibuf(1) = parent_id
442	CALL wrf_put_dom_ti_integer ( fid , 'PARENT_ID' , ibuf , 1 , ierr )
443	ibuf(1) = i_parent_start
444	CALL wrf_put_dom_ti_integer ( fid , 'I_PARENT_START' , ibuf , 1 , ierr )
445	ibuf(1) = j_parent_start
446	CALL wrf_put_dom_ti_integer ( fid , 'J_PARENT_START' , ibuf , 1 , ierr )
447	ibuf(1) = parent_grid_ratio
448	CALL wrf_put_dom_ti_integer ( fid , 'PARENT_GRID_RATIO' , ibuf , 1 , ierr )
449	#endif
450
451	! end add 010831 JM
452
453	#if (EM_CORE != 1)
454	CALL wrf_put_dom_ti_real ( fid , 'DX' , config_flags%dx , 1 , ierr )
455	CALL wrf_put_dom_ti_real ( fid , 'DY' , config_flags%dy , 1 , ierr )
456	#endif
457	! Updated by T. Hutchinson to use adaptive time step.
458	CALL wrf_put_dom_ti_real ( fid , 'DT' , grid%dt , 1 , ierr )
459	! CALL wrf_put_dom_ti_real ( fid , 'DT' , config%dt , 1 , ierr )
460	CALL wrf_put_dom_ti_real ( fid , 'CEN_LAT' , config_flags%cen_lat , 1 , ierr )
461	CALL wrf_put_dom_ti_real ( fid , 'CEN_LON' , config_flags%cen_lon , 1 , ierr )
462	CALL wrf_put_dom_ti_real ( fid , 'TRUELAT1', config_flags%truelat1, 1 , ierr )
463	CALL wrf_put_dom_ti_real ( fid , 'TRUELAT2', config_flags%truelat2, 1 , ierr )
464	CALL wrf_put_dom_ti_real ( fid , 'MOAD_CEN_LAT', config_flags%moad_cen_lat, 1 , ierr )
465	CALL wrf_put_dom_ti_real ( fid , 'STAND_LON', config_flags%stand_lon, 1 , ierr )
466	IF ( switch .NE. boundary_only ) THEN
467	#ifdef PLANET
468	! When writing to restart files, use the values of the instantaneous
469	! time for determining the values of JULYR, JULDAY, and GMT. If the
470	! original values in config_flags are used, this assumes that the
471	! restart simulation will start with an itimestep NE 0. If we use
472	! the instantaneous time, we can start a restart simulation with a
473	! different value of delta-t for timestep and still get the clocks
474	! calendars (and orbital information!) correct.
475	!
476	! Current time is still defined from above call to WRF_UTIL_ClockGet
477	CALL WRFU_TimeGet( currentTime, YY=julyr, dayOfYear=julday, H=p_hr, M=p_min, S=p_sec, MS=p_ms, rc=rc)
478	WRITE(wrf_err_message,*)'output_wrf: julyr,julday,hr,min,sec,ms = ',julyr,julday,p_hr,p_min,p_sec,p_ms
479	CALL wrf_debug( 100 , wrf_err_message )
480	gmt = REAL(p_hr)+REAL(p_min)/60.+REAL(p_sec)/3600.+REAL(p_ms)/3600000.
481	CALL wrf_put_dom_ti_real ( fid , 'GMT' , gmt , 1 , ierr )
482	CALL wrf_put_dom_ti_integer ( fid , 'JULYR' , julyr , 1 , ierr )
483	CALL wrf_put_dom_ti_integer ( fid , 'JULDAY' , julday , 1 , ierr )
484	#else
485	CALL wrf_put_dom_ti_real ( fid , 'GMT' , config_flags%gmt , 1 , ierr )
486	CALL wrf_put_dom_ti_integer ( fid , 'JULYR' , config_flags%julyr , 1 , ierr )
487	CALL wrf_put_dom_ti_integer ( fid , 'JULDAY' , config_flags%julday , 1 , ierr )
488	#endif
489	ENDIF
490	#if (NMM_CORE == 1)
491	write(0,*) 'MMINLU would be: ', MMINLU(1:4)
492	MMINLU(1:4)='USGS'
493	write(0,*) 'MMINLU now: ', MMINLU(1:4)
494	#endif
495	CALL wrf_put_dom_ti_integer ( fid , 'MAP_PROJ' , config_flags%map_proj , 1 , ierr )
496	CALL wrf_put_dom_ti_char ( fid , 'MMINLU', mminlu(1:4) , ierr )
497	CALL wrf_put_dom_ti_integer ( fid , 'ISWATER' , iswater , 1 , ierr )
498	! CALL wrf_put_dom_ti_integer ( fid , 'ISWATER' , config_flags%iswater , 1 , ierr )
499	CALL wrf_put_dom_ti_integer ( fid , 'ISICE' , config_flags%isice , 1 , ierr )
500	CALL wrf_put_dom_ti_integer ( fid , 'ISURBAN' , config_flags%isurban , 1 , ierr )
501	CALL wrf_put_dom_ti_integer ( fid , 'ISOILWATER' , config_flags%isoilwater , 1 , ierr )
502	! added these fields for restarting of moving nests, JM
503	CALL wrf_put_dom_ti_integer ( fid , 'I_PARENT_START' , config_flags%i_parent_start , 1 , ierr )
504	CALL wrf_put_dom_ti_integer ( fid , 'J_PARENT_START' , config_flags%j_parent_start , 1 , ierr )
505
506
507	IF ( switch .EQ. boundary_only ) THEN
508	CALL WRFU_TimeIntervalSet( bdy_increment, S=NINT(config_flags%bdyfrq),rc=rc)
509	next_time = currentTime + bdy_increment
510	CALL wrf_timetoa ( next_time, next_datestr )
511	CALL wrf_put_dom_td_char ( fid , 'THISBDYTIME' , current_date(1:19), current_date(1:19), ierr )
512	CALL wrf_put_dom_td_char ( fid , 'NEXTBDYTIME' , current_date(1:19), next_datestr(1:19), ierr )
513	ENDIF
514
515	! added grib2-specific metadata: Todd Hutchinson 8/21/2005
516	IF ( use_package( io_form ) == IO_GRIB2 ) THEN
517	CALL wrf_put_dom_ti_integer ( fid , 'BACKGROUND_PROC_ID' , config_flags%background_proc_id , 1 , ierr )
518	CALL wrf_put_dom_ti_integer ( fid , 'FORECAST_PROC_ID' , config_flags%forecast_proc_id , 1 , ierr )
519	CALL wrf_put_dom_ti_integer ( fid , 'PRODUCTION_STATUS' , config_flags%production_status , 1 , ierr )
520	CALL wrf_put_dom_ti_integer ( fid , 'COMPRESSION' , config_flags%compression , 1 , ierr )
521	ENDIF
522
523	CALL nl_get_adjust_output_times( grid%id, adjust )
524	current_date_save = current_date
525	#if 1
526	IF ( switch .EQ. model_input_only ) THEN
527	CALL wrf_debug ( 300 , 'output_wrf: calling code in wrf_inputout.inc' )
528	CALL wrf_inputout( fid , grid , config_flags, switch, dryrun, ierr )
529	ELSE IF ( switch .EQ. aux_model_input1_only ) THEN
530	CALL wrf_debug ( 300 , 'output_wrf: calling code in wrf_auxinput1out.inc' )
531	CALL wrf_auxinput1out( fid , grid , config_flags, switch, dryrun, ierr )
532	ELSE IF ( switch .EQ. aux_model_input2_only ) THEN
533	CALL wrf_debug ( 300 , 'output_wrf: calling code in wrf_auxinput2out.inc' )
534	CALL wrf_auxinput2out( fid , grid , config_flags, switch, dryrun, ierr )
535	ELSE IF ( switch .EQ. aux_model_input3_only ) THEN
536	CALL wrf_debug ( 300 , 'output_wrf: calling code in wrf_auxinput3out.inc' )
537	CALL wrf_auxinput3out( fid , grid , config_flags, switch, dryrun, ierr )
538	ELSE IF ( switch .EQ. aux_model_input4_only ) THEN
539	CALL wrf_debug ( 300 , 'output_wrf: calling code in wrf_auxinput4out.inc' )
540	CALL wrf_auxinput4out( fid , grid , config_flags, switch, dryrun, ierr )
541	ELSE IF ( switch .EQ. aux_model_input5_only ) THEN
542	CALL wrf_debug ( 300 , 'output_wrf: calling code in wrf_auxinput5out.inc' )
543	CALL wrf_auxinput5out( fid , grid , config_flags, switch, dryrun, ierr )
544	ELSE IF ( switch .EQ. aux_model_input6_only ) THEN
545	CALL wrf_debug ( 300 , 'output_wrf: calling code in wrf_auxinput6out.inc' )
546	CALL wrf_auxinput6out( fid , grid , config_flags, switch, dryrun, ierr )
547	ELSE IF ( switch .EQ. aux_model_input7_only ) THEN
548	CALL wrf_debug ( 300 , 'output_wrf: calling code in wrf_auxinput7out.inc' )
549	CALL wrf_auxinput7out( fid , grid , config_flags, switch, dryrun, ierr )
550	ELSE IF ( switch .EQ. aux_model_input8_only ) THEN
551	CALL wrf_debug ( 300 , 'output_wrf: calling code in wrf_auxinput8out.inc' )
552	CALL wrf_auxinput8out( fid , grid , config_flags, switch, dryrun, ierr )
553	ELSE IF ( switch .EQ. aux_model_input9_only ) THEN
554	CALL wrf_debug ( 300 , 'output_wrf: calling code in wrf_auxinput9out.inc' )
555	CALL wrf_auxinput9out( fid , grid , config_flags, switch, dryrun, ierr )
556	ELSE IF ( switch .EQ. aux_model_input10_only ) THEN
557	CALL wrf_debug ( 300 , 'output_wrf: calling code in wrf_auxinput10out.inc' )
558	CALL wrf_auxinput10out( fid , grid , config_flags, switch, dryrun, ierr )
559	ELSE IF ( switch .EQ. aux_model_input11_only ) THEN
560	CALL wrf_debug ( 300 , 'output_wrf: calling code in wrf_auxinput11out.inc' )
561	CALL wrf_auxinput11out( fid , grid , config_flags, switch, dryrun, ierr )
562	ELSE IF ( switch .EQ. history_only ) THEN
563	CALL wrf_debug ( 300 , 'output_wrf: calling code in wrf_histout.inc' )
564	IF ( adjust ) CALL adjust_io_timestr( grid%io_intervals( HISTORY_ALARM ), currentTime, startTime, current_date )
565	CALL wrf_histout( fid , grid , config_flags, switch, dryrun, ierr )
566	ELSE IF ( switch .EQ. aux_hist1_only ) THEN
567	CALL wrf_debug ( 300 , 'output_wrf: calling code in wrf_auxhist1out' )
568	IF ( adjust ) CALL adjust_io_timestr( grid%io_intervals( AUXHIST1_ALARM ), currentTime, startTime, current_date )
569	CALL wrf_auxhist1out( fid , grid , config_flags, switch, dryrun, ierr )
570	ELSE IF ( switch .EQ. aux_hist2_only ) THEN
571	CALL wrf_debug ( 300 , 'output_wrf: calling code in wrf_auxhist2out.inc' )
572	IF ( adjust ) CALL adjust_io_timestr( grid%io_intervals( AUXHIST2_ALARM ), currentTime, startTime, current_date )
573	CALL wrf_auxhist2out( fid , grid , config_flags, switch, dryrun, ierr )
574	ELSE IF ( switch .EQ. aux_hist3_only ) THEN
575	CALL wrf_debug ( 300 , 'output_wrf: calling code in wrf_auxhist3out.inc' )
576	IF ( adjust ) CALL adjust_io_timestr( grid%io_intervals( AUXHIST3_ALARM ), currentTime, startTime, current_date )
577	CALL wrf_auxhist3out( fid , grid , config_flags, switch, dryrun, ierr )
578	ELSE IF ( switch .EQ. aux_hist4_only ) THEN
579	CALL wrf_debug ( 300 , 'output_wrf: calling code in wrf_auxhist4out.inc' )
580	IF ( adjust ) CALL adjust_io_timestr( grid%io_intervals( AUXHIST4_ALARM ), currentTime, startTime, current_date )
581	CALL wrf_auxhist4out( fid , grid , config_flags, switch, dryrun, ierr )
582	ELSE IF ( switch .EQ. aux_hist5_only ) THEN
583	CALL wrf_debug ( 300 , 'output_wrf: calling code in wrf_auxhist5out.inc' )
584	IF ( adjust ) CALL adjust_io_timestr( grid%io_intervals( AUXHIST5_ALARM ), currentTime, startTime, current_date )
585	CALL wrf_auxhist5out( fid , grid , config_flags, switch, dryrun, ierr )
586	ELSE IF ( switch .EQ. aux_hist6_only ) THEN
587	CALL wrf_debug ( 300 , 'output_wrf: calling code in wrf_auxhist6out.inc' )
588	IF ( adjust ) CALL adjust_io_timestr( grid%io_intervals( AUXHIST6_ALARM ), currentTime, startTime, current_date )
589	CALL wrf_auxhist6out( fid , grid , config_flags, switch, dryrun, ierr )
590	ELSE IF ( switch .EQ. aux_hist7_only ) THEN
591	CALL wrf_debug ( 300 , 'output_wrf: calling code in wrf_auxhist7out.inc' )
592	IF ( adjust ) CALL adjust_io_timestr( grid%io_intervals( AUXHIST7_ALARM ), currentTime, startTime, current_date )
593	CALL wrf_auxhist7out( fid , grid , config_flags, switch, dryrun, ierr )
594	ELSE IF ( switch .EQ. aux_hist8_only ) THEN
595	CALL wrf_debug ( 300 , 'output_wrf: calling code in wrf_auxhist8out.inc' )
596	IF ( adjust ) CALL adjust_io_timestr( grid%io_intervals( AUXHIST8_ALARM ), currentTime, startTime, current_date )
597	CALL wrf_auxhist8out( fid , grid , config_flags, switch, dryrun, ierr )
598	ELSE IF ( switch .EQ. aux_hist9_only ) THEN
599	CALL wrf_debug ( 300 , 'output_wrf: calling code in wrf_auxhist9out.inc' )
600	IF ( adjust ) CALL adjust_io_timestr( grid%io_intervals( AUXHIST9_ALARM ), currentTime, startTime, current_date )
601	CALL wrf_auxhist9out( fid , grid , config_flags, switch, dryrun, ierr )
602	ELSE IF ( switch .EQ. aux_hist10_only ) THEN
603	CALL wrf_debug ( 300 , 'output_wrf: calling code in wrf_auxhist10out.inc' )
604	IF ( adjust ) CALL adjust_io_timestr( grid%io_intervals( AUXHIST10_ALARM ), currentTime, startTime, current_date )
605	CALL wrf_auxhist10out( fid , grid , config_flags, switch, dryrun, ierr )
606	ELSE IF ( switch .EQ. aux_hist11_only ) THEN
607	CALL wrf_debug ( 300 , 'output_wrf: calling code in wrf_auxhist11out.inc' )
608	IF ( adjust ) CALL adjust_io_timestr( grid%io_intervals( AUXHIST11_ALARM ), currentTime, startTime, current_date )
609	CALL wrf_auxhist11out( fid , grid , config_flags, switch, dryrun, ierr )
610	ELSE IF ( switch .EQ. restart_only ) THEN
611	CALL wrf_debug ( 300 , 'output_wrf: calling code in wrf_restartout.inc' )
612	CALL wrf_restartout( fid , grid , config_flags, switch, dryrun, ierr )
613	ELSE IF ( switch .EQ. boundary_only ) THEN
614	CALL wrf_debug ( 300 , 'output_wrf: calling code in wrf_bdyout.inc' )
615	CALL wrf_bdyout( fid , grid , config_flags, switch, dryrun, ierr )
616	ENDIF
617	#else
618	CALL wrf_message ( "ALL I/O DISABLED IN share/module_io_wrf.F")
619	#endif
620	current_date = current_date_save
621
622	IF ( .NOT. dryrun ) THEN
623	CALL wrf_debug ( 300 , 'output_wrf: calling wrf_iosync ' )
624	CALL wrf_iosync ( fid , ierr )
625	CALL wrf_debug ( 300 , 'output_wrf: back from wrf_iosync ' )
626	ENDIF
627
628	WRITE(wrf_err_message,*)'output_wrf: end, fid = ',fid
629	CALL wrf_debug( 300 , wrf_err_message )
630
631	RETURN
632	END SUBROUTINE output_wrf

Note: See TracBrowser for help on using the repository browser.

Download in other formats:

Original Format