2 ! ---principal wrf output routine (called from routines in module_io_domain )
3 SUBROUTINE output_wrf ( fid , grid , config_flags, switch , ierr )
8 USE module_domain_type, ONLY : fieldlist
9 USE module_state_description
11 ! USE module_date_time
12 USE module_model_constants
15 #include <wrf_io_flags.h>
16 #include <wrf_status_codes.h>
18 TYPE(grid_config_rec_type), INTENT(INOUT) :: config_flags
19 INTEGER, INTENT(IN) :: fid, switch
20 INTEGER, INTENT(INOUT) :: ierr
23 INTEGER ids , ide , jds , jde , kds , kde , &
24 ims , ime , jms , jme , kms , kme , &
25 ips , ipe , jps , jpe , kps , kpe
27 TYPE( fieldlist ), POINTER :: p
29 INTEGER newswitch, itrace
31 INTEGER , DIMENSION(3) :: domain_start , domain_end
32 INTEGER , DIMENSION(3) :: memory_start , memory_end
33 INTEGER , DIMENSION(3) :: patch_start , patch_end
35 INTEGER julyr, julday, idt, iswater , islake, map_proj
38 REAL gmt, cen_lat, cen_lon, bdyfrq , truelat1 , truelat2 , moad_cen_lat , stand_lon
39 INTEGER km_opt, diff_opt, damp_opt, &
40 mp_physics, ra_lw_physics, ra_sw_physics, sf_sfclay_physics, &
41 sf_surface_physics, bl_pbl_physics, cu_physics
42 REAL khdif, kvdif, swrad_scat, dampcoef
43 INTEGER sf_urban_physics, w_damping, smooth_option, feedback, surface_input_source, sst_update
45 INTEGER grid_id , parent_id , i_parent_start , j_parent_start , parent_grid_ratio
48 INTEGER grid_fdda, gfdda_interval_m, gfdda_end_h, if_ramping, &
49 obs_nudge_opt, obs_nudge_wind, obs_nudge_temp, obs_nudge_mois, obs_nudge_pstr, obs_idynin, obs_ionf
50 INTEGER grid_sfdda, sgfdda_interval_m, sgfdda_end_h
51 REAL fgdt, guv, gt, gq, gph, dtramp_min, &
52 obs_coef_wind, obs_coef_temp, obs_coef_mois, obs_coef_pstr, obs_dtramp, fdda_end
53 REAL guv_sfc, gt_sfc, gq_sfc, rinblw
54 INTEGER moist_adv_opt, scalar_adv_opt, tke_adv_opt
55 INTEGER save_topo_orig
57 CHARACTER (len=19) simulation_start_date
58 CHARACTER (len=len_current_date) current_date_save
59 INTEGER simulation_start_year , &
60 simulation_start_month , &
61 simulation_start_day , &
62 simulation_start_hour , &
63 simulation_start_minute , &
64 simulation_start_second
67 LOGICAL, EXTERNAL :: multi_files
68 INTEGER, EXTERNAL :: use_package
69 INTEGER p_hr, p_min, p_sec, p_ms
71 CHARACTER*80 dname, memord
74 CHARACTER*80 char_junk
75 CHARACTER(LEN=256) :: MMINLU
78 TYPE(WRFU_TimeInterval) :: bdy_increment
79 TYPE(WRFU_Time) :: next_time, currentTime, startTime
80 CHARACTER*40 :: next_datestr
81 INTEGER :: start_year , start_month , start_day , start_hour , start_minute , start_second
84 TYPE(WRFU_Time) :: ringTime, stopTime, curtime
85 TYPE(WRFU_TimeInterval) :: interval, tmpinterval
86 CHARACTER*80 alarmname, timestring, debuggal
87 INTEGER seconds, seconds2, iring
89 WRITE(wrf_err_message,*)'output_wrf: begin, fid = ',fid
90 CALL wrf_debug( 300 , wrf_err_message )
92 CALL modify_io_masks ( grid%id ) ! this adjusts the I/O masks according to the users run-time specs, if any
94 CALL wrf_inquire_filename ( fid , fname , filestate , ierr )
96 WRITE(wrf_err_message,*)'module_io_wrf: output_wrf: wrf_inquire_filename Status = ',ierr
97 CALL wrf_error_fatal( wrf_err_message )
100 WRITE(wrf_err_message,*)'output_wrf: fid,filestate = ',fid,filestate
101 CALL wrf_debug( 300 , wrf_err_message )
103 ! io_form is used to determine if multi-file I/O is enabled and to
104 ! control writing of format-specific time-independent metadata
105 io_form = io_form_for_stream( switch )
107 dryrun = ( filestate .EQ. WRF_FILE_OPENED_NOT_COMMITTED )
109 WRITE(wrf_err_message,*)'output_wrf: dryrun = ',dryrun
110 CALL wrf_debug( 300 , wrf_err_message )
112 CALL get_ijk_from_grid ( grid , &
113 ids, ide, jds, jde, kds, kde, &
114 ims, ime, jms, jme, kms, kme, &
115 ips, ipe, jps, jpe, kps, kpe )
117 call nl_get_diff_opt ( 1, diff_opt )
118 call nl_get_km_opt ( 1, km_opt )
119 call nl_get_damp_opt ( 1, damp_opt )
120 call nl_get_dampcoef ( grid%id, dampcoef )
121 call nl_get_khdif ( grid%id, khdif )
122 call nl_get_kvdif ( grid%id, kvdif )
123 call nl_get_mp_physics ( grid%id, mp_physics )
124 call nl_get_ra_lw_physics ( grid%id, ra_lw_physics )
125 call nl_get_ra_sw_physics ( grid%id, ra_sw_physics )
126 call nl_get_sf_sfclay_physics ( grid%id, sf_sfclay_physics )
127 call nl_get_sf_surface_physics ( grid%id, sf_surface_physics )
128 call nl_get_bl_pbl_physics ( grid%id, bl_pbl_physics )
129 call nl_get_cu_physics ( grid%id, cu_physics )
131 ! add nml variables in 2.2
132 call nl_get_surface_input_source ( 1 , surface_input_source )
133 call nl_get_sst_update ( 1 , sst_update )
134 call nl_get_feedback ( 1 , feedback )
135 call nl_get_smooth_option ( 1 , smooth_option )
136 call nl_get_swrad_scat ( 1 , swrad_scat )
137 call nl_get_sf_urban_physics ( 1 , sf_urban_physics )
138 call nl_get_w_damping ( 1 , w_damping )
141 CALL nl_get_moist_adv_opt ( grid%id , moist_adv_opt )
142 CALL nl_get_scalar_adv_opt ( grid%id , scalar_adv_opt )
143 CALL nl_get_tke_adv_opt ( grid%id , tke_adv_opt )
144 CALL nl_get_diff_6th_opt ( grid%id , diff_6th_opt )
145 CALL nl_get_diff_6th_factor ( grid%id , diff_6th_factor )
146 CALL nl_get_grid_fdda ( grid%id , grid_fdda )
147 CALL nl_get_auxinput10_end_h( grid%id , gfdda_end_h )
148 CALL nl_get_auxinput10_interval_m ( grid%id , gfdda_interval_m )
149 CALL nl_get_grid_sfdda ( grid%id , grid_sfdda )
150 CALL nl_get_auxinput9_end_h( grid%id , sgfdda_end_h )
151 CALL nl_get_auxinput9_interval_m ( grid%id , sgfdda_interval_m )
153 IF ( grid_fdda == 1 ) THEN
154 CALL nl_get_fgdt ( grid%id , fgdt )
155 CALL nl_get_guv ( grid%id , guv )
156 CALL nl_get_gt ( grid%id , gt )
157 CALL nl_get_gq ( grid%id , gq )
158 CALL nl_get_if_ramping ( 1 , if_ramping )
159 CALL nl_get_dtramp_min ( 1 , dtramp_min )
162 IF ( grid_fdda == 2 ) THEN
163 CALL nl_get_fgdt ( grid%id , fgdt )
164 CALL nl_get_guv ( grid%id , guv )
165 CALL nl_get_gt ( grid%id , gt )
166 CALL nl_get_gph ( grid%id , gph )
167 CALL nl_get_if_ramping ( 1 , if_ramping )
168 CALL nl_get_dtramp_min ( 1 , dtramp_min )
171 IF ( grid_sfdda == 1 ) THEN
172 CALL nl_get_guv_sfc ( grid%id , guv_sfc )
173 CALL nl_get_gt_sfc ( grid%id , gt_sfc )
174 CALL nl_get_gq_sfc ( grid%id , gq_sfc )
175 CALL nl_get_rinblw ( grid%id , rinblw )
178 CALL nl_get_obs_nudge_opt ( grid%id , obs_nudge_opt )
179 IF ( obs_nudge_opt == 1 ) THEN
180 CALL nl_get_fdda_end ( grid%id , fdda_end )
181 CALL nl_get_obs_nudge_wind ( grid%id , obs_nudge_wind )
182 CALL nl_get_obs_coef_wind ( grid%id , obs_coef_wind )
183 CALL nl_get_obs_nudge_temp ( grid%id , obs_nudge_temp )
184 CALL nl_get_obs_coef_temp ( grid%id , obs_coef_temp )
185 CALL nl_get_obs_nudge_mois ( grid%id , obs_nudge_mois )
186 CALL nl_get_obs_coef_mois ( grid%id , obs_coef_mois )
187 CALL nl_get_obs_nudge_pstr ( grid%id , obs_nudge_pstr )
188 CALL nl_get_obs_coef_pstr ( grid%id , obs_coef_pstr )
189 CALL nl_get_obs_ionf ( 1 , obs_ionf )
190 CALL nl_get_obs_idynin ( 1 , obs_idynin )
191 CALL nl_get_obs_dtramp ( 1 , obs_dtramp )
195 ! julday and gmt can be set in namelist_03 for ideal.exe run
196 CALL nl_get_gmt (grid%id, gmt)
197 CALL nl_get_julyr (grid%id, julyr)
198 CALL nl_get_julday (grid%id, julday)
199 CALL nl_get_mminlu ( grid%id, mminlu )
200 call wrf_debug(300,"OUTPUT_WRF: mminlu = " // mminlu )
201 CALL nl_get_iswater (grid%id, iswater )
202 CALL nl_get_islake (grid%id, islake )
203 CALL nl_get_cen_lat ( grid%id , cen_lat )
204 CALL nl_get_cen_lon ( grid%id , cen_lon )
205 CALL nl_get_truelat1 ( grid%id , truelat1 )
206 CALL nl_get_truelat2 ( grid%id , truelat2 )
207 CALL nl_get_moad_cen_lat ( grid%id , moad_cen_lat )
208 CALL nl_get_stand_lon ( grid%id , stand_lon )
209 CALL nl_get_map_proj ( grid%id , map_proj )
212 CALL nl_get_parent_id ( grid%id , parent_id )
213 CALL nl_get_i_parent_start ( grid%id , i_parent_start )
214 CALL nl_get_j_parent_start ( grid%id , j_parent_start )
215 CALL nl_get_parent_grid_ratio ( grid%id , parent_grid_ratio )
218 CALL domain_clockprint(150, grid, &
219 'DEBUG output_wrf(): before call to domain_clock_get,')
220 CALL domain_clock_get( grid, current_time=currentTime, &
221 start_time=startTime, &
222 current_timestr=current_date )
224 IF (switch .EQ. history_only) THEN
225 CALL nl_get_adjust_output_times( grid%id, adjust )
227 CALL adjust_io_timestr( grid%io_intervals( history_alarm ), currentTime, startTime, timestring )
228 current_date_save = current_date
229 current_date = timestring
233 WRITE ( wrf_err_message , * ) 'output_wrf: begin, current_date=',current_date
234 CALL wrf_debug ( 300 , wrf_err_message )
236 WRITE( message , * ) "OUTPUT FROM " , TRIM(program_name)
237 CALL wrf_put_dom_ti_char ( fid , 'TITLE' , TRIM(message) , ierr )
238 ! added grib-specific metadata: Todd Hutchinson 8/21/2005
239 IF ( ( use_package( io_form ) == IO_GRIB1 ) .OR. &
240 ( use_package( io_form ) == IO_GRIB2 ) ) THEN
241 CALL wrf_put_dom_ti_char ( fid, 'PROGRAM_NAME', TRIM(program_name) , ierr )
243 CALL nl_get_start_year(grid%id,start_year)
244 CALL nl_get_start_month(grid%id,start_month)
245 CALL nl_get_start_day(grid%id,start_day)
246 CALL nl_get_start_hour(grid%id,start_hour)
247 CALL nl_get_start_minute(grid%id,start_minute)
248 CALL nl_get_start_second(grid%id,start_second)
250 WRITE ( start_date , FMT = '(I4.4,"-",I5.5,"_",I2.2,":",I2.2,":",I2.2)' ) &
251 start_year,start_day,start_hour,start_minute,start_second
253 WRITE ( start_date , FMT = '(I4.4,"-",I2.2,"-",I2.2,"_",I2.2,":",I2.2,":",I2.2)' ) &
254 start_year,start_month,start_day,start_hour,start_minute,start_second
256 CALL wrf_put_dom_ti_char ( fid , 'START_DATE', TRIM(start_date) , ierr )
257 IF ( switch .EQ. input_only) THEN
258 CALL wrf_put_dom_ti_char ( fid , 'SIMULATION_START_DATE', TRIM(start_date) , ierr )
259 ELSE IF ( ( switch .EQ. restart_only ) .OR. ( switch .EQ. history_only ) ) THEN
260 CALL nl_get_simulation_start_year ( 1, simulation_start_year )
261 CALL nl_get_simulation_start_month ( 1, simulation_start_month )
262 CALL nl_get_simulation_start_day ( 1, simulation_start_day )
263 CALL nl_get_simulation_start_hour ( 1, simulation_start_hour )
264 CALL nl_get_simulation_start_minute ( 1, simulation_start_minute )
265 CALL nl_get_simulation_start_second ( 1, simulation_start_second )
267 WRITE ( simulation_start_date , FMT = '(I4.4,"-",I5.5,"_",I2.2,":",I2.2,":",I2.2)' ) &
268 simulation_start_year,simulation_start_day,&
269 simulation_start_hour,simulation_start_minute,simulation_start_second
271 WRITE ( simulation_start_date , FMT = '(I4.4,"-",I2.2,"-",I2.2,"_",I2.2,":",I2.2,":",I2.2)' ) &
272 simulation_start_year,simulation_start_month,simulation_start_day,&
273 simulation_start_hour,simulation_start_minute,simulation_start_second
275 CALL wrf_put_dom_ti_char ( fid , 'SIMULATION_START_DATE', TRIM(simulation_start_date) , ierr )
278 IF ( switch .EQ. restart_only ) THEN
279 ! add some delicious alarm metadata to the restart files
280 ibuf(1) = MAX_WRF_ALARMS
281 CALL wrf_put_dom_ti_integer( fid, 'MAX_WRF_ALARMS', ibuf, 1, ierr )
282 curtime = domain_get_current_time( grid )
283 DO i = 1, MAX_WRF_ALARMS
284 IF ( grid%alarms_created(i) ) THEN
285 IF ( i .LT. 10 ) THEN
286 write(alarmname,'("WRF_ALARM_ISRINGING_0",i1)')i
288 write(alarmname,'("WRF_ALARM_ISRINGING_",i2)')i
290 IF ( WRFU_AlarmIsRinging( grid%alarms( i ), rc=rc ) ) THEN
295 CALL wrf_put_dom_ti_integer( fid, TRIM(alarmname), iring, 1, ierr )
297 CALL WRFU_AlarmGet( grid%alarms(i),PrevRingTime=ringTime,RingInterval=interval,rc=rc)
300 !different, perhaps simpler way to do this, JM 20100319
301 IF ( i .LT. 10 ) THEN
302 write(alarmname,'("WRF_ALARM_PREVRINGTIME_0",i1)')i
304 write(alarmname,'("WRF_ALARM_PREVRINGTIME_",i2)')i
306 CALL WRFU_TimeGet( ringTime, timeString=prevringtime )
307 CALL wrf_put_dom_ti_char( fid, TRIM(alarmname), TRIM(prevringtime), 1, ierr )
310 ! compute time to next ring as interval minus time since last ring
311 tmpinterval = curtime - ringTime
313 IF ( i .LT. 10 ) THEN
314 write(alarmname,'("WRF_ALARM_SECS_TIL_NEXT_RING_0",i1)')i
316 write(alarmname,'("WRF_ALARM_SECS_TIL_NEXT_RING_",i2)')i
318 CALL WRFU_TimeIntervalGet(interval,S=seconds)
319 CALL WRFU_TimeIntervalGet(tmpinterval,S=seconds2)
320 IF ( seconds .GE. 1700000000 .OR. seconds .LE. -1700000000 ) THEN ! it is a forever value, do not change it
321 CALL wrf_put_dom_ti_integer( fid, TRIM(alarmname), seconds, 1, ierr )
323 CALL wrf_put_dom_ti_integer( fid, TRIM(alarmname), seconds-seconds2, 1, ierr )
326 if ( i .eq. history_alarm .or. i .eq. restart_alarm .or. i .eq. auxhist24_alarm) then
327 if ( i .eq. history_alarm ) write(*,*)'----------------- info for history alarm ',i
328 if ( i .eq. restart_alarm ) write(*,*)'----------------- info for restart alarm ',i
329 if ( i .eq. auxhist24_alarm ) write(*,*)'----------------- info for auxhist24 alarm ',i
330 CALL ESMF_TimeGet(curTime, TimeString=debuggal ) ; write(*,*)'curTime ',TRIM(debuggal)
331 CALL ESMF_TimeGet(ringTime, TimeString=debuggal ) ; write(*,*)'ringTime ',TRIM(debuggal)
332 CALL ESMF_TimeIntervalGet(tmpinterval, TimeString=debuggal ) ; write(*,*)'tmpinterval ',TRIM(debuggal)
333 CALL ESMF_TimeIntervalGet(interval, TimeString=debuggal ) ; write(*,*)'interval ',TRIM(debuggal)
334 write(*,*)'interval seconds ',seconds,' tmpinterval seconds ',seconds2,' seconds-seconds2 ',seconds-seconds2
335 write(*,*)'----------------- '
343 ibuf(1) = config_flags%e_we - config_flags%s_we + 1
344 CALL wrf_put_dom_ti_integer ( fid , 'WEST-EAST_GRID_DIMENSION' , ibuf , 1 , ierr )
346 ibuf(1) = config_flags%e_sn - config_flags%s_sn + 1
347 CALL wrf_put_dom_ti_integer ( fid , 'SOUTH-NORTH_GRID_DIMENSION' , ibuf , 1 , ierr )
349 ibuf(1) = config_flags%e_vert - config_flags%s_vert + 1
350 CALL wrf_put_dom_ti_integer ( fid , 'BOTTOM-TOP_GRID_DIMENSION' , ibuf , 1 , ierr )
353 IF (grid%map_proj == 6) THEN
354 ! Global ... dx determined automatically
355 ! Don't use value from namelist ... used derived value instead
356 CALL wrf_put_dom_ti_real ( fid , 'DX' , grid%dx , 1 , ierr )
357 CALL wrf_put_dom_ti_real ( fid , 'DY' , grid%dy , 1 , ierr )
359 CALL wrf_put_dom_ti_real ( fid , 'DX' , config_flags%dx , 1 , ierr )
360 CALL wrf_put_dom_ti_real ( fid , 'DY' , config_flags%dy , 1 , ierr )
365 !added the following for tc bogusing
366 if((config_flags%insert_bogus_storm) .or. (config_flags%remove_storm)) then
367 print *,"we have confirmed that insert or remove is true"
369 CALL wrf_put_dom_ti_integer ( fid , 'FLAG_METGRID' , ibuf , 1 , ierr )
370 if( grid%flag_snow .eq. 1) then
371 CALL wrf_put_dom_ti_integer ( fid , 'FLAG_SNOW' , ibuf , 1 , ierr )
373 if( grid%flag_mf_xy .eq. 1) then
374 CALL wrf_put_dom_ti_integer ( fid , 'FLAG_MF_XY' , ibuf , 1 , ierr )
376 if(grid%flag_psfc .eq. 1) then
377 CALL wrf_put_dom_ti_integer ( fid , 'FLAG_PSFC' , ibuf , 1 , ierr )
380 if(grid%flag_slp .eq. 1) then
381 CALL wrf_put_dom_ti_integer ( fid , 'FLAG_SLP' , ibuf , 1 , ierr )
384 if(grid%flag_sm000010 .eq. 1)then
385 CALL wrf_put_dom_ti_integer ( fid , 'FLAG_SM000010' , ibuf , 1 , ierr )
387 if(grid%flag_sm010040 .eq. 1)then
388 CALL wrf_put_dom_ti_integer ( fid , 'FLAG_SM010040' , ibuf , 1 , ierr )
390 if(grid%flag_sm040100 .eq. 1)then
391 CALL wrf_put_dom_ti_integer ( fid , 'FLAG_SM040100' , ibuf , 1 , ierr )
393 if(grid%flag_sm100200 .eq. 1)then
394 CALL wrf_put_dom_ti_integer ( fid , 'FLAG_SM100200' , ibuf , 1 , ierr )
397 if(grid%flag_st000010 .eq. 1)then
398 CALL wrf_put_dom_ti_integer ( fid , 'FLAG_ST000010' , ibuf , 1 , ierr )
400 if(grid%flag_st010040 .eq. 1)then
401 CALL wrf_put_dom_ti_integer ( fid , 'FLAG_ST010040' , ibuf , 1 , ierr )
403 if(grid%flag_st040100 .eq. 1)then
404 CALL wrf_put_dom_ti_integer ( fid , 'FLAG_ST040100' , ibuf , 1 , ierr )
406 if(grid%flag_st100200 .eq. 1)then
407 CALL wrf_put_dom_ti_integer ( fid , 'FLAG_ST100200' , ibuf , 1 , ierr )
410 ibuf(1) = grid%num_metgrid_levels
411 CALL wrf_put_dom_ti_integer ( fid , 'BOTTOM-TOP_GRID_DIMENSION' , ibuf , 1 , ierr )
413 CALL wrf_put_dom_ti_integer ( fid , 'num_metgrid_levels' , ibuf , 1 , ierr )
416 CALL wrf_put_dom_ti_char ( fid , 'SIMULATION_START_DATE', TRIM(start_date) , ierr )
418 WRITE( message , * ) "OUTPUT FROM TC BOGUS"
419 CALL wrf_put_dom_ti_char ( fid , 'TITLE' , TRIM(message) , ierr )
424 ! added this metadatum for H. Chuan, NCEP, 030417, JM
426 CALL wrf_put_dom_ti_char ( fid , 'GRIDTYPE', 'E' , ierr )
429 CALL wrf_put_dom_ti_char ( fid , 'GRIDTYPE', 'C' , ierr )
431 #if (COAMPS_CORE == 1 )
432 CALL wrf_put_dom_ti_char ( fid , 'GRIDTYPE', 'B' , ierr )
435 ! added these fields for W. Skamarock, 020402, JM
437 CALL wrf_put_dom_ti_integer ( fid , 'DIFF_OPT' , ibuf , 1 , ierr )
439 CALL wrf_put_dom_ti_integer ( fid , 'KM_OPT' , ibuf , 1 , ierr )
441 CALL wrf_put_dom_ti_integer ( fid , 'DAMP_OPT' , ibuf , 1 , ierr )
443 CALL wrf_put_dom_ti_real ( fid , 'DAMPCOEF' , rbuf , 1 , ierr )
445 CALL wrf_put_dom_ti_real ( fid , 'KHDIF' , rbuf , 1 , ierr )
447 CALL wrf_put_dom_ti_real ( fid , 'KVDIF' , rbuf , 1 , ierr )
449 CALL wrf_put_dom_ti_integer ( fid , 'MP_PHYSICS' , ibuf , 1 , ierr )
450 ibuf(1) = ra_lw_physics
451 CALL wrf_put_dom_ti_integer ( fid , 'RA_LW_PHYSICS' , ibuf , 1 , ierr )
452 ibuf(1) = ra_sw_physics
453 CALL wrf_put_dom_ti_integer ( fid , 'RA_SW_PHYSICS' , ibuf , 1 , ierr )
454 ibuf(1) = sf_sfclay_physics
455 CALL wrf_put_dom_ti_integer ( fid , 'SF_SFCLAY_PHYSICS' , ibuf , 1 , ierr )
456 ibuf(1) = sf_surface_physics
457 CALL wrf_put_dom_ti_integer ( fid , 'SF_SURFACE_PHYSICS' , ibuf , 1 , ierr )
458 ibuf(1) = bl_pbl_physics
459 CALL wrf_put_dom_ti_integer ( fid , 'BL_PBL_PHYSICS' , ibuf , 1 , ierr )
461 CALL wrf_put_dom_ti_integer ( fid , 'CU_PHYSICS' , ibuf , 1 , ierr )
463 ! added netcdf-specific metadata:
464 IF ( ( use_package( io_form ) == IO_NETCDF ) .OR. &
465 ( use_package( io_form ) == IO_PHDF5 ) .OR. &
466 ( use_package( io_form ) == IO_PNETCDF ) ) THEN
467 CALL wrf_put_dom_ti_integer ( fid, 'SURFACE_INPUT_SOURCE', surface_input_source , 1 , ierr )
468 CALL wrf_put_dom_ti_integer ( fid, 'SST_UPDATE', sst_update , 1 , ierr )
470 CALL wrf_put_dom_ti_integer ( fid, 'GRID_FDDA', grid_fdda , 1 , ierr )
471 CALL wrf_put_dom_ti_integer ( fid, 'GFDDA_INTERVAL_M', gfdda_interval_m , 1 , ierr )
472 CALL wrf_put_dom_ti_integer ( fid, 'GFDDA_END_H', gfdda_end_h , 1 , ierr )
473 CALL wrf_put_dom_ti_integer ( fid, 'GRID_SFDDA', grid_sfdda , 1 , ierr )
474 CALL wrf_put_dom_ti_integer ( fid, 'SGFDDA_INTERVAL_M', sgfdda_interval_m , 1 , ierr )
475 CALL wrf_put_dom_ti_integer ( fid, 'SGFDDA_END_H', sgfdda_end_h , 1 , ierr )
478 IF ( switch .EQ. history_only ) THEN
479 CALL wrf_put_dom_ti_integer ( fid, 'SF_URBAN_PHYSICS', sf_urban_physics , 1 , ierr )
480 CALL wrf_put_dom_ti_integer ( fid, 'FEEDBACK', feedback , 1 , ierr )
481 CALL wrf_put_dom_ti_integer ( fid, 'SMOOTH_OPTION', smooth_option , 1 , ierr )
482 CALL wrf_put_dom_ti_real ( fid, 'SWRAD_SCAT', swrad_scat , 1 , ierr )
483 CALL wrf_put_dom_ti_integer ( fid, 'W_DAMPING', w_damping , 1 , ierr )
486 CALL wrf_put_dom_ti_integer ( fid, 'MOIST_ADV_OPT', moist_adv_opt , 1 , ierr )
487 CALL wrf_put_dom_ti_integer ( fid, 'SCALAR_ADV_OPT', scalar_adv_opt , 1 , ierr )
488 CALL wrf_put_dom_ti_integer ( fid, 'TKE_ADV_OPT', tke_adv_opt , 1 , ierr )
489 CALL wrf_put_dom_ti_integer ( fid, 'DIFF_6TH_OPT', diff_6th_opt , 1 , ierr )
490 CALL wrf_put_dom_ti_real ( fid, 'DIFF_6TH_FACTOR', diff_6th_factor , 1 , ierr )
492 IF ( grid_fdda == 1 ) THEN
493 CALL wrf_put_dom_ti_real ( fid, 'FGDT', fgdt , 1 , ierr )
494 CALL wrf_put_dom_ti_real ( fid, 'GUV', guv , 1 , ierr )
495 CALL wrf_put_dom_ti_real ( fid, 'GT', gt , 1 , ierr )
496 CALL wrf_put_dom_ti_real ( fid, 'GQ', gq , 1 , ierr )
497 CALL wrf_put_dom_ti_integer ( fid, 'IF_RAMPING', if_ramping , 1 , ierr )
498 CALL wrf_put_dom_ti_real ( fid, 'DTRAMP_MIN', dtramp_min , 1 , ierr )
501 IF ( grid_fdda == 2 ) THEN
502 CALL wrf_put_dom_ti_real ( fid, 'FGDT', fgdt , 1 , ierr )
503 CALL wrf_put_dom_ti_real ( fid, 'GUV', guv , 1 , ierr )
504 CALL wrf_put_dom_ti_real ( fid, 'GT', gt , 1 , ierr )
505 CALL wrf_put_dom_ti_real ( fid, 'GPH', gph , 1 , ierr )
506 CALL wrf_put_dom_ti_integer ( fid, 'IF_RAMPING', if_ramping , 1 , ierr )
507 CALL wrf_put_dom_ti_real ( fid, 'DTRAMP_MIN', dtramp_min , 1 , ierr )
510 IF ( grid_sfdda == 1 ) THEN
511 CALL wrf_put_dom_ti_real ( fid, 'GUV_SFC', guv_sfc , 1 , ierr )
512 CALL wrf_put_dom_ti_real ( fid, 'GT_SFC', gt_sfc , 1 , ierr )
513 CALL wrf_put_dom_ti_real ( fid, 'GQ_SFC', gq_sfc , 1 , ierr )
514 CALL wrf_put_dom_ti_real ( fid, 'RINBLW', rinblw , 1 , ierr )
517 CALL wrf_put_dom_ti_integer ( fid, 'OBS_NUDGE_OPT', obs_nudge_opt , 1 , ierr )
518 IF ( obs_nudge_opt == 1 ) THEN
519 CALL wrf_put_dom_ti_real ( fid, 'FDDA_END', fdda_end , 1 , ierr )
520 CALL wrf_put_dom_ti_integer ( fid, 'OBS_NUDGE_WIND', obs_nudge_wind , 1 , ierr )
521 CALL wrf_put_dom_ti_real ( fid, 'OBS_COEF_WIND', obs_coef_wind , 1 , ierr )
522 CALL wrf_put_dom_ti_integer ( fid, 'OBS_NUDGE_TEMP', obs_nudge_temp , 1 , ierr )
523 CALL wrf_put_dom_ti_real ( fid, 'OBS_COEF_TEMP', obs_coef_temp , 1 , ierr )
524 CALL wrf_put_dom_ti_integer ( fid, 'OBS_NUDGE_MOIS', obs_nudge_mois , 1 , ierr )
525 CALL wrf_put_dom_ti_real ( fid, 'OBS_COEF_MOIS', obs_coef_mois , 1 , ierr )
526 CALL wrf_put_dom_ti_integer ( fid, 'OBS_NUDGE_PSTR', obs_nudge_pstr , 1 , ierr )
527 CALL wrf_put_dom_ti_real ( fid, 'OBS_COEF_PSTR', obs_coef_pstr , 1 , ierr )
528 CALL wrf_put_dom_ti_integer ( fid, 'OBS_IONF', obs_ionf , 1 , ierr )
529 CALL wrf_put_dom_ti_integer ( fid, 'OBS_IDYNIN', obs_idynin , 1 , ierr )
530 CALL wrf_put_dom_ti_real ( fid, 'OBS_DTRAMP', obs_dtramp , 1 , ierr )
533 CALL wrf_put_dom_ti_real ( fid, 'BUCKET_MM', config_flags%bucket_mm , 1 , ierr )
534 CALL wrf_put_dom_ti_real ( fid, 'BUCKET_J', config_flags%bucket_J , 1 , ierr )
535 CALL wrf_put_dom_ti_real ( fid, 'PREC_ACC_DT', config_flags%prec_acc_dt , 1 , ierr )
536 CALL wrf_put_dom_ti_integer ( fid, 'OMLCALL', config_flags%omlcall , 1 , ierr )
537 CALL wrf_put_dom_ti_integer ( fid, 'ISFTCFLX', config_flags%isftcflx , 1 , ierr )
538 CALL wrf_put_dom_ti_integer ( fid, 'ISHALLOW', config_flags%ishallow , 1 , ierr )
542 CALL wrf_put_dom_ti_integer ( fid, 'DFI_OPT', config_flags%dfi_opt , 1 , ierr )
543 CALL wrf_put_dom_ti_integer ( fid, 'SHCU_PHYSICS', config_flags%shcu_physics , 1 , ierr )
544 CALL wrf_put_dom_ti_integer ( fid, 'SF_URBAN_PHYSICS', config_flags%sf_urban_physics , 1 , ierr )
549 ! added these fields for use by reassembly programs , 010831, JM
550 ! modified these fields so "patch" == "domain" when multi-file output
551 ! formats are not used. 051018, TBH
553 ibuf(1) = MAX(ips,ids)
554 IF ( .NOT. multi_files ( io_form ) ) ibuf(1) = ids
555 CALL wrf_put_dom_ti_integer ( fid , 'WEST-EAST_PATCH_START_UNSTAG' , ibuf , 1 , ierr )
556 ibuf(1) = MIN(ipe,ide-1)
557 IF ( .NOT. multi_files ( io_form ) ) ibuf(1) = ide - 1
558 CALL wrf_put_dom_ti_integer ( fid , 'WEST-EAST_PATCH_END_UNSTAG' , ibuf , 1 , ierr )
559 ibuf(1) = MAX(ips,ids)
560 IF ( .NOT. multi_files ( io_form ) ) ibuf(1) = ids
561 CALL wrf_put_dom_ti_integer ( fid , 'WEST-EAST_PATCH_START_STAG' , ibuf , 1 , ierr )
562 ibuf(1) = MIN(ipe,ide)
563 IF ( .NOT. multi_files ( io_form ) ) ibuf(1) = ide
564 CALL wrf_put_dom_ti_integer ( fid , 'WEST-EAST_PATCH_END_STAG' , ibuf , 1 , ierr )
565 ibuf(1) = MAX(jps,jds)
566 IF ( .NOT. multi_files ( io_form ) ) ibuf(1) = jds
567 CALL wrf_put_dom_ti_integer ( fid , 'SOUTH-NORTH_PATCH_START_UNSTAG' , ibuf , 1 , ierr )
568 ibuf(1) = MIN(jpe,jde-1)
569 IF ( .NOT. multi_files ( io_form ) ) ibuf(1) = jde - 1
570 CALL wrf_put_dom_ti_integer ( fid , 'SOUTH-NORTH_PATCH_END_UNSTAG' , ibuf , 1 , ierr )
571 ibuf(1) = MAX(jps,jds)
572 IF ( .NOT. multi_files ( io_form ) ) ibuf(1) = jds
573 CALL wrf_put_dom_ti_integer ( fid , 'SOUTH-NORTH_PATCH_START_STAG' , ibuf , 1 , ierr )
574 ibuf(1) = MIN(jpe,jde)
575 IF ( .NOT. multi_files ( io_form ) ) ibuf(1) = jde
576 CALL wrf_put_dom_ti_integer ( fid , 'SOUTH-NORTH_PATCH_END_STAG' , ibuf , 1 , ierr )
578 ibuf(1) = MAX(kps,kds)
579 IF ( .NOT. multi_files ( io_form ) ) ibuf(1) = kds
580 CALL wrf_put_dom_ti_integer ( fid , 'BOTTOM-TOP_PATCH_START_UNSTAG' , ibuf , 1 , ierr )
581 ibuf(1) = MIN(kpe,kde-1)
582 IF ( .NOT. multi_files ( io_form ) ) ibuf(1) = kde - 1
583 CALL wrf_put_dom_ti_integer ( fid , 'BOTTOM-TOP_PATCH_END_UNSTAG' , ibuf , 1 , ierr )
584 ibuf(1) = MAX(kps,kds)
585 IF ( .NOT. multi_files ( io_form ) ) ibuf(1) = kds
586 CALL wrf_put_dom_ti_integer ( fid , 'BOTTOM-TOP_PATCH_START_STAG' , ibuf , 1 , ierr )
587 ibuf(1) = MIN(kpe,kde)
588 IF ( .NOT. multi_files ( io_form ) ) ibuf(1) = kde
589 CALL wrf_put_dom_ti_integer ( fid , 'BOTTOM-TOP_PATCH_END_STAG' , ibuf , 1 , ierr )
592 CALL wrf_put_dom_ti_integer ( fid , 'GRID_ID' , ibuf , 1 , ierr )
594 CALL wrf_put_dom_ti_integer ( fid , 'PARENT_ID' , ibuf , 1 , ierr )
595 ibuf(1) = i_parent_start
596 CALL wrf_put_dom_ti_integer ( fid , 'I_PARENT_START' , ibuf , 1 , ierr )
597 ibuf(1) = j_parent_start
598 CALL wrf_put_dom_ti_integer ( fid , 'J_PARENT_START' , ibuf , 1 , ierr )
599 ibuf(1) = parent_grid_ratio
600 CALL wrf_put_dom_ti_integer ( fid , 'PARENT_GRID_RATIO' , ibuf , 1 , ierr )
606 CALL wrf_put_dom_ti_real ( fid , 'DX' , config_flags%dx , 1 , ierr )
607 CALL wrf_put_dom_ti_real ( fid , 'DY' , config_flags%dy , 1 , ierr )
609 ! Updated by T. Hutchinson to use adaptive time step.
610 CALL wrf_put_dom_ti_real ( fid , 'DT' , grid%dt , 1 , ierr )
611 ! CALL wrf_put_dom_ti_real ( fid , 'DT' , config%dt , 1 , ierr )
612 CALL wrf_put_dom_ti_real ( fid , 'CEN_LAT' , config_flags%cen_lat , 1 , ierr )
613 CALL wrf_put_dom_ti_real ( fid , 'CEN_LON' , config_flags%cen_lon , 1 , ierr )
614 CALL wrf_put_dom_ti_real ( fid , 'TRUELAT1', config_flags%truelat1, 1 , ierr )
615 CALL wrf_put_dom_ti_real ( fid , 'TRUELAT2', config_flags%truelat2, 1 , ierr )
616 CALL wrf_put_dom_ti_real ( fid , 'MOAD_CEN_LAT', config_flags%moad_cen_lat, 1 , ierr )
617 CALL wrf_put_dom_ti_real ( fid , 'STAND_LON', config_flags%stand_lon, 1 , ierr )
619 CALL wrf_put_dom_ti_real ( fid , 'POLE_LAT', config_flags%pole_lat, 1 , ierr )
620 CALL wrf_put_dom_ti_real ( fid , 'POLE_LON', config_flags%pole_lon, 1 , ierr )
622 IF ( switch .NE. boundary_only .AND. switch .NE. auxinput9_only .AND. switch .NE. auxinput10_only ) THEN
624 ! When writing to restart files, use the values of the instantaneous
625 ! time for determining the values of JULYR, JULDAY, and GMT. If the
626 ! original values in config_flags are used, this assumes that the
627 ! restart simulation will start with an itimestep NE 0. If we use
628 ! the instantaneous time, we can start a restart simulation with a
629 ! different value of delta-t for timestep and still get the clocks
630 ! calendars (and orbital information!) correct.
632 ! Current time is still defined from above call to WRF_UTIL_ClockGet
633 CALL WRFU_TimeGet( currentTime, YY=julyr, dayOfYear=julday, H=p_hr, M=p_min, S=p_sec, MS=p_ms, rc=rc)
634 WRITE(wrf_err_message,*)'output_wrf: julyr,julday,hr,min,sec,ms = ',julyr,julday,p_hr,p_min,p_sec,p_ms
635 CALL wrf_debug( 100 , wrf_err_message )
636 gmt = REAL(p_hr)+REAL(p_min)/60.+REAL(p_sec)/3600.+REAL(p_ms)/3600000.
637 CALL wrf_put_dom_ti_real ( fid , 'GMT' , gmt , 1 , ierr )
638 CALL wrf_put_dom_ti_integer ( fid , 'JULYR' , julyr , 1 , ierr )
639 CALL wrf_put_dom_ti_integer ( fid , 'JULDAY' , julday , 1 , ierr )
641 CALL wrf_put_dom_ti_real ( fid , 'GMT' , config_flags%gmt , 1 , ierr )
642 CALL wrf_put_dom_ti_integer ( fid , 'JULYR' , config_flags%julyr , 1 , ierr )
643 CALL wrf_put_dom_ti_integer ( fid , 'JULDAY' , config_flags%julday , 1 , ierr )
649 CALL wrf_put_dom_ti_integer ( fid , 'MAP_PROJ' , config_flags%map_proj , 1 , ierr )
650 IF(MMINLU(1:1) .EQ. " ")THEN
651 CALL wrf_put_dom_ti_char ( fid , 'MMINLU', " " , ierr )
653 CALL wrf_put_dom_ti_char ( fid , 'MMINLU', TRIM(mminlu) , ierr )
655 call wrf_put_dom_ti_integer ( fid , 'NUM_LAND_CAT', config_flags%num_land_cat, 1, ierr)
656 CALL wrf_put_dom_ti_integer ( fid , 'ISWATER' , iswater , 1 , ierr )
657 CALL wrf_put_dom_ti_integer ( fid , 'ISLAKE' , islake , 1 , ierr )
658 ! CALL wrf_put_dom_ti_integer ( fid , 'ISWATER' , config_flags%iswater , 1 , ierr )
659 CALL wrf_put_dom_ti_integer ( fid , 'ISICE' , config_flags%isice , 1 , ierr )
660 CALL wrf_put_dom_ti_integer ( fid , 'ISURBAN' , config_flags%isurban , 1 , ierr )
661 CALL wrf_put_dom_ti_integer ( fid , 'ISOILWATER' , config_flags%isoilwater , 1 , ierr )
662 ! added these fields for restarting of moving nests, JM
665 CALL wrf_put_dom_ti_integer ( fid , 'I_PARENT_START' , grid%i_parent_start , 1 , ierr )
666 CALL wrf_put_dom_ti_integer ( fid , 'J_PARENT_START' , grid%j_parent_start , 1 , ierr )
668 CALL wrf_put_dom_ti_integer ( fid , 'I_PARENT_START' , config_flags%i_parent_start , 1 , ierr )
669 CALL wrf_put_dom_ti_integer ( fid , 'J_PARENT_START' , config_flags%j_parent_start , 1 , ierr )
672 IF ( switch .EQ. boundary_only ) THEN
673 CALL WRFU_TimeIntervalSet( bdy_increment, S=NINT(config_flags%bdyfrq),rc=rc)
674 next_time = currentTime + bdy_increment
675 CALL wrf_timetoa ( next_time, next_datestr )
676 CALL wrf_put_dom_td_char ( fid , 'THISBDYTIME' , current_date(1:19), current_date(1:19), ierr )
677 CALL wrf_put_dom_td_char ( fid , 'NEXTBDYTIME' , current_date(1:19), next_datestr(1:19), ierr )
680 ! added grib2-specific metadata: Todd Hutchinson 8/21/2005
681 IF ( use_package( io_form ) == IO_GRIB2 ) THEN
682 CALL wrf_put_dom_ti_integer ( fid , 'BACKGROUND_PROC_ID' , config_flags%background_proc_id , 1 , ierr )
683 CALL wrf_put_dom_ti_integer ( fid , 'FORECAST_PROC_ID' , config_flags%forecast_proc_id , 1 , ierr )
684 CALL wrf_put_dom_ti_integer ( fid , 'PRODUCTION_STATUS' , config_flags%production_status , 1 , ierr )
685 CALL wrf_put_dom_ti_integer ( fid , 'COMPRESSION' , config_flags%compression , 1 , ierr )
688 #if ( (EM_CORE == 1) && (DA_CORE != 1) )
689 save_topo_orig = grid%save_topo_from_real
692 IF ( (first_history .LE. switch .AND. switch .LE. last_history ) .OR. &
693 ( (switch .EQ. input_only) .AND. &
694 (program_name(1:7) .NE. 'REAL_EM') .AND. &
695 (grid%dfi_opt .EQ. DFI_NODFI ) ) .OR. &
696 ( switch .EQ. restart_only ) ) THEN
698 ! This flag sets the switch which defines the topography as the original
699 ! generated by real.exe. The "zero" value means that for the WRF model
700 ! output, the topography has been modified. All output from the ARW
701 ! WRF model has this flag set. However, to allow nests to still be
702 ! instantiated after a parent does IO, and to allow that nest domain
703 ! to have the topo adjusted, we save the incoming value of the save_topo
706 grid%save_topo_from_real=0
710 IF ( (first_input .LE. switch .AND. switch .LE. last_input) .OR. &
711 (first_history .LE. switch .AND. switch .LE. last_history ) .OR. &
712 switch .EQ. restart_only ) THEN
714 p => grid%head_statevars%next
715 DO WHILE ( ASSOCIATED( p ) )
716 IF ( p%ProcOrient .NE. 'X' .AND. p%ProcOrient .NE. 'Y' ) THEN ! no I/O for xposed variables
717 IF ( p%Ndim .EQ. 0 ) THEN
718 IF ((p%Restart.AND.switch.EQ.restart_only).OR.on_stream(p%streams,newswitch)) THEN
719 IF ( in_use_for_config(grid%id,TRIM(p%VarName)) ) THEN
721 IF (p%Ntl.GT.0.AND.switch.NE.restart_only)dname=dname(1:len(TRIM(dname))-2)
722 memord = p%MemoryOrder
723 IF ( p%Type .EQ. 'r' ) THEN
724 CALL wrf_ext_write_field ( &
726 current_date(1:19) , & ! DateStr
727 TRIM(p%DataName) , & ! Data Name
728 p%rfield_0d , & ! Field
729 WRF_FLOAT , & ! FieldType
730 grid%communicator , & ! Comm
731 grid%iocommunicator , & ! Comm
732 grid%domdesc , & ! Comm
733 grid%bdy_mask , & ! bdy_mask
735 '0' , & ! MemoryOrder
740 TRIM(p%Description) , & ! Desc
741 TRIM(p%Units) , & ! Units
742 __FILE__ // ' writing 0d real ' // TRIM(p%VarName) , & ! Debug message
743 1 , 1 , 1 , 1 , 1 , 1 , &
744 1 , 1 , 1 , 1 , 1 , 1 , &
745 1 , 1 , 1 , 1 , 1 , 1 , &
747 ELSE IF ( p%Type .EQ. 'd' ) THEN
748 CALL wrf_ext_write_field ( &
750 current_date(1:19) , & ! DateStr
751 TRIM(p%DataName) , & ! Data Name
752 p%dfield_0d , & ! Field
753 WRF_DOUBLE , & ! FieldType
754 grid%communicator , & ! Comm
755 grid%iocommunicator , & ! Comm
756 grid%domdesc , & ! Comm
757 grid%bdy_mask , & ! bdy_mask
759 '0' , & ! MemoryOrder
764 TRIM(p%Description) , & ! Desc
765 TRIM(p%Units) , & ! Units
766 __FILE__ // ' writing 0d double ' // TRIM(p%VarName) , & ! Debug message
767 1 , 1 , 1 , 1 , 1 , 1 , &
768 1 , 1 , 1 , 1 , 1 , 1 , &
769 1 , 1 , 1 , 1 , 1 , 1 , &
771 ELSE IF ( p%Type .EQ. 'i' ) THEN
772 CALL wrf_ext_write_field ( &
774 current_date(1:19) , & ! DateStr
775 TRIM(p%DataName) , & ! Data Name
776 p%ifield_0d , & ! Field
777 WRF_INTEGER , & ! FieldType
778 grid%communicator , & ! Comm
779 grid%iocommunicator , & ! Comm
780 grid%domdesc , & ! Comm
781 grid%bdy_mask , & ! bdy_mask
783 '0' , & ! MemoryOrder
788 TRIM(p%Description) , & ! Desc
789 TRIM(p%Units) , & ! Units
790 __FILE__ // ' writing 0d integer ' // TRIM(p%VarName) , & ! Debug message
791 1 , 1 , 1 , 1 , 1 , 1 , &
792 1 , 1 , 1 , 1 , 1 , 1 , &
793 1 , 1 , 1 , 1 , 1 , 1 , &
795 ELSE IF ( p%Type .EQ. 'l' ) THEN
796 CALL wrf_ext_write_field ( &
798 current_date(1:19) , & ! DateStr
799 TRIM(p%DataName) , & ! Data Name
800 p%lfield_0d , & ! Field
801 WRF_LOGICAL , & ! FieldType
802 grid%communicator , & ! Comm
803 grid%iocommunicator , & ! Comm
804 grid%domdesc , & ! Comm
805 grid%bdy_mask , & ! bdy_mask
807 '0' , & ! MemoryOrder
812 TRIM(p%Description) , & ! Desc
813 TRIM(p%Units) , & ! Units
814 __FILE__ // ' writing 0d logical ' // TRIM(p%VarName) , & ! Debug message
815 1 , 1 , 1 , 1 , 1 , 1 , &
816 1 , 1 , 1 , 1 , 1 , 1 , &
817 1 , 1 , 1 , 1 , 1 , 1 , &
822 ELSE IF ( p%Ndim .EQ. 1 ) THEN
823 IF ((p%Restart.AND.switch.EQ.restart_only).OR.on_stream(p%streams,newswitch)) THEN
824 IF ( in_use_for_config(grid%id,TRIM(p%VarName)) ) THEN
825 IF (switch.EQ.restart_only.OR.p%Ntl/100.EQ.mod(p%Ntl,100)) THEN
827 IF (p%Ntl.GT.0.AND.switch.NE.restart_only)dname=dname(1:len(TRIM(dname))-2)
828 memord = p%MemoryOrder
829 IF ( p%Type .EQ. 'r' ) THEN
830 CALL wrf_ext_write_field ( &
832 current_date(1:19) , & ! DateStr
833 TRIM(dname) , & ! Data Name
834 p%rfield_1d , & ! Field
835 WRF_FLOAT , & ! FieldType
836 grid%communicator , & ! Comm
837 grid%iocommunicator , & ! Comm
838 grid%domdesc , & ! Comm
839 grid%bdy_mask , & ! bdy_mask
841 TRIM(memord) , & ! MemoryOrder
842 TRIM(p%Stagger) , & ! Stagger
843 TRIM(p%dimname1) , & ! Dimname 1
844 TRIM(p%dimname2) , & ! Dimname 2
845 TRIM(p%dimname3) , & ! Dimname 3
846 TRIM(p%Description) , & ! Desc
847 TRIM(p%Units) , & ! Units
848 __FILE__ // ' writing 1d real ' // TRIM(p%VarName) , & ! Debug message
849 p%sd1 , p%ed1 , p%sd2 , p%ed2 , p%sd3 , p%ed3 , &
850 p%sm1 , p%em1 , p%sm2 , p%em2 , p%sm3 , p%em3 , &
851 p%sp1 , p%ep1 , p%sp2 , p%ep2 , p%sp3 , p%ep3 , &
853 ELSE IF ( p%Type .EQ. 'd' ) THEN
854 CALL wrf_ext_write_field ( &
856 current_date(1:19) , & ! DateStr
857 TRIM(dname) , & ! Data Name
858 p%dfield_1d , & ! Field
859 WRF_DOUBLE , & ! FieldType
860 grid%communicator , & ! Comm
861 grid%iocommunicator , & ! Comm
862 grid%domdesc , & ! Comm
863 grid%bdy_mask , & ! bdy_mask
865 TRIM(memord) , & ! MemoryOrder
866 TRIM(p%Stagger) , & ! Stagger
867 TRIM(p%dimname1) , & ! Dimname 1
868 TRIM(p%dimname2) , & ! Dimname 2
869 TRIM(p%dimname3) , & ! Dimname 3
870 TRIM(p%Description) , & ! Desc
871 TRIM(p%Units) , & ! Units
872 __FILE__ // ' writing 1d double ' // TRIM(p%VarName) , & ! Debug message
873 p%sd1 , p%ed1 , p%sd2 , p%ed2 , p%sd3 , p%ed3 , &
874 p%sm1 , p%em1 , p%sm2 , p%em2 , p%sm3 , p%em3 , &
875 p%sp1 , p%ep1 , p%sp2 , p%ep2 , p%sp3 , p%ep3 , &
877 ELSE IF ( p%Type .EQ. 'i' ) THEN
878 CALL wrf_ext_write_field ( &
880 current_date(1:19) , & ! DateStr
881 TRIM(dname) , & ! Data Name
882 p%ifield_1d , & ! Field
883 WRF_INTEGER , & ! FieldType
884 grid%communicator , & ! Comm
885 grid%iocommunicator , & ! Comm
886 grid%domdesc , & ! Comm
887 grid%bdy_mask , & ! bdy_mask
889 TRIM(memord) , & ! MemoryOrder
890 TRIM(p%Stagger) , & ! Stagger
891 TRIM(p%dimname1) , & ! Dimname 1
892 TRIM(p%dimname2) , & ! Dimname 2
893 TRIM(p%dimname3) , & ! Dimname 3
894 TRIM(p%Description) , & ! Desc
895 TRIM(p%Units) , & ! Units
896 __FILE__ // ' writing 1d integer ' // TRIM(p%VarName) , & ! Debug message
897 p%sd1 , p%ed1 , p%sd2 , p%ed2 , p%sd3 , p%ed3 , &
898 p%sm1 , p%em1 , p%sm2 , p%em2 , p%sm3 , p%em3 , &
899 p%sp1 , p%ep1 , p%sp2 , p%ep2 , p%sp3 , p%ep3 , &
901 ELSE IF ( p%Type .EQ. 'l' ) THEN
902 CALL wrf_ext_write_field ( &
904 current_date(1:19) , & ! DateStr
905 TRIM(dname) , & ! Data Name
906 p%lfield_1d , & ! Field
907 WRF_LOGICAL , & ! FieldType
908 grid%communicator , & ! Comm
909 grid%iocommunicator , & ! Comm
910 grid%domdesc , & ! Comm
911 grid%bdy_mask , & ! bdy_mask
913 TRIM(memord) , & ! MemoryOrder
914 TRIM(p%Stagger) , & ! Stagger
915 TRIM(p%dimname1) , & ! Dimname 1
916 TRIM(p%dimname2) , & ! Dimname 2
917 TRIM(p%dimname3) , & ! Dimname 3
918 TRIM(p%Description) , & ! Desc
919 TRIM(p%Units) , & ! Units
920 __FILE__ // ' writing 1d logical ' // TRIM(p%VarName) , & ! Debug message
921 p%sd1 , p%ed1 , p%sd2 , p%ed2 , p%sd3 , p%ed3 , &
922 p%sm1 , p%em1 , p%sm2 , p%em2 , p%sm3 , p%em3 , &
923 p%sp1 , p%ep1 , p%sp2 , p%ep2 , p%sp3 , p%ep3 , &
929 ELSE IF ( p%Ndim .EQ. 2 ) THEN
930 IF ((p%Restart.AND.switch.EQ.restart_only).OR.on_stream(p%streams,newswitch)) THEN
931 IF ( in_use_for_config(grid%id,TRIM(p%VarName)) .AND. &
932 ( .NOT. p%subgrid_x .OR. (p%subgrid_x .AND. grid%sr_x .GT. 0) ) .AND. &
933 ( .NOT. p%subgrid_y .OR. (p%subgrid_y .AND. grid%sr_y .GT. 0) ) &
935 IF (switch.EQ.restart_only.OR.p%Ntl/100.EQ.mod(p%Ntl,100)) THEN
937 IF (p%Ntl.GT.0.AND.switch.NE.restart_only)dname=dname(1:len(TRIM(dname))-2)
938 memord = p%MemoryOrder
939 IF ( p%Type .EQ. 'r' ) THEN
940 CALL wrf_ext_write_field ( &
942 current_date(1:19) , & ! DateStr
943 TRIM(dname) , & ! Data Name
944 p%rfield_2d , & ! Field
945 WRF_FLOAT , & ! FieldType
946 grid%communicator , & ! Comm
947 grid%iocommunicator , & ! Comm
948 grid%domdesc , & ! Comm
949 grid%bdy_mask , & ! bdy_mask
951 TRIM(memord) , & ! MemoryOrder
952 TRIM(p%Stagger) , & ! Stagger
953 TRIM(p%dimname1) , & ! Dimname 1
954 TRIM(p%dimname2) , & ! Dimname 2
955 TRIM(p%dimname3) , & ! Dimname 3
956 TRIM(p%Description) , & ! Desc
957 TRIM(p%Units) , & ! Units
958 __FILE__ // ' writing 2d real ' // TRIM(p%VarName) , & ! Debug message
959 p%sd1 , p%ed1 , p%sd2 , p%ed2 , p%sd3 , p%ed3 , &
960 p%sm1 , p%em1 , p%sm2 , p%em2 , p%sm3 , p%em3 , &
961 p%sp1 , p%ep1 , p%sp2 , p%ep2 , p%sp3 , p%ep3 , &
963 ELSE IF ( p%Type .EQ. 'd' ) THEN
964 CALL wrf_ext_write_field ( &
966 current_date(1:19) , & ! DateStr
967 TRIM(dname) , & ! Data Name
968 p%dfield_2d , & ! Field
969 WRF_DOUBLE , & ! FieldType
970 grid%communicator , & ! Comm
971 grid%iocommunicator , & ! Comm
972 grid%domdesc , & ! Comm
973 grid%bdy_mask , & ! bdy_mask
975 TRIM(memord) , & ! MemoryOrder
976 TRIM(p%Stagger) , & ! Stagger
977 TRIM(p%dimname1) , & ! Dimname 1
978 TRIM(p%dimname2) , & ! Dimname 2
979 TRIM(p%dimname3) , & ! Dimname 3
980 TRIM(p%Description) , & ! Desc
981 TRIM(p%Units) , & ! Units
982 __FILE__ // ' writing 2d double ' // TRIM(p%VarName) , & ! Debug message
983 p%sd1 , p%ed1 , p%sd2 , p%ed2 , p%sd3 , p%ed3 , &
984 p%sm1 , p%em1 , p%sm2 , p%em2 , p%sm3 , p%em3 , &
985 p%sp1 , p%ep1 , p%sp2 , p%ep2 , p%sp3 , p%ep3 , &
987 ELSE IF ( p%Type .EQ. 'i' ) THEN
988 CALL wrf_ext_write_field ( &
990 current_date(1:19) , & ! DateStr
991 TRIM(dname) , & ! Data Name
992 p%ifield_2d , & ! Field
993 WRF_INTEGER , & ! FieldType
994 grid%communicator , & ! Comm
995 grid%iocommunicator , & ! Comm
996 grid%domdesc , & ! Comm
997 grid%bdy_mask , & ! bdy_mask
999 TRIM(memord) , & ! MemoryOrder
1000 TRIM(p%Stagger) , & ! Stagger
1001 TRIM(p%dimname1) , & ! Dimname 1
1002 TRIM(p%dimname2) , & ! Dimname 2
1003 TRIM(p%dimname3) , & ! Dimname 3
1004 TRIM(p%Description) , & ! Desc
1005 TRIM(p%Units) , & ! Units
1006 __FILE__ // ' writing 2d integer ' // TRIM(p%VarName) , & ! Debug message
1007 p%sd1 , p%ed1 , p%sd2 , p%ed2 , p%sd3 , p%ed3 , &
1008 p%sm1 , p%em1 , p%sm2 , p%em2 , p%sm3 , p%em3 , &
1009 p%sp1 , p%ep1 , p%sp2 , p%ep2 , p%sp3 , p%ep3 , &
1011 ELSE IF ( p%Type .EQ. 'l' ) THEN
1012 CALL wrf_ext_write_field ( &
1013 fid , & ! DataHandle
1014 current_date(1:19) , & ! DateStr
1015 TRIM(dname) , & ! Data Name
1016 p%lfield_2d , & ! Field
1017 WRF_LOGICAL , & ! FieldType
1018 grid%communicator , & ! Comm
1019 grid%iocommunicator , & ! Comm
1020 grid%domdesc , & ! Comm
1021 grid%bdy_mask , & ! bdy_mask
1023 TRIM(memord) , & ! MemoryOrder
1024 TRIM(p%Stagger) , & ! Stagger
1025 TRIM(p%dimname1) , & ! Dimname 1
1026 TRIM(p%dimname2) , & ! Dimname 2
1027 TRIM(p%dimname3) , & ! Dimname 3
1028 TRIM(p%Description) , & ! Desc
1029 TRIM(p%Units) , & ! Units
1030 __FILE__ // ' writing 2d logical ' // TRIM(p%VarName) , & ! Debug message
1031 p%sd1 , p%ed1 , p%sd2 , p%ed2 , p%sd3 , p%ed3 , &
1032 p%sm1 , p%em1 , p%sm2 , p%em2 , p%sm3 , p%em3 , &
1033 p%sp1 , p%ep1 , p%sp2 , p%ep2 , p%sp3 , p%ep3 , &
1039 ELSE IF ( p%Ndim .EQ. 3 ) THEN
1040 IF ((p%Restart.AND.switch.EQ.restart_only).OR.on_stream(p%streams,newswitch)) THEN
1041 IF ( in_use_for_config(grid%id,TRIM(p%VarName)) .AND. &
1042 ( .NOT. p%subgrid_x .OR. (p%subgrid_x .AND. grid%sr_x .GT. 0) ) .AND. &
1043 ( .NOT. p%subgrid_y .OR. (p%subgrid_y .AND. grid%sr_y .GT. 0) ) &
1045 IF (switch.EQ.restart_only.OR.p%Ntl/100.EQ.mod(p%Ntl,100)) THEN
1047 IF (p%Ntl.GT.0.AND.switch.NE.restart_only)dname=dname(1:len(TRIM(dname))-2)
1048 memord = p%MemoryOrder
1049 IF ( p%Type .EQ. 'r' ) THEN
1050 CALL wrf_ext_write_field ( &
1051 fid , & ! DataHandle
1052 current_date(1:19) , & ! DateStr
1053 TRIM(dname) , & ! Data Name
1054 p%rfield_3d , & ! Field
1055 WRF_FLOAT , & ! FieldType
1056 grid%communicator , & ! Comm
1057 grid%iocommunicator , & ! Comm
1058 grid%domdesc , & ! Comm
1059 grid%bdy_mask , & ! bdy_mask
1061 TRIM(memord) , & ! MemoryOrder
1062 TRIM(p%Stagger) , & ! Stagger
1063 TRIM(p%dimname1) , & ! Dimname 1
1064 TRIM(p%dimname2) , & ! Dimname 2
1065 TRIM(p%dimname3) , & ! Dimname 3
1066 TRIM(p%Description) , & ! Desc
1067 TRIM(p%Units) , & ! Units
1068 __FILE__ // ' writing 3d real ' // TRIM(p%VarName) , & ! Debug message
1069 p%sd1 , p%ed1 , p%sd2 , p%ed2 , p%sd3 , p%ed3 , &
1070 p%sm1 , p%em1 , p%sm2 , p%em2 , p%sm3 , p%em3 , &
1071 p%sp1 , p%ep1 , p%sp2 , p%ep2 , p%sp3 , p%ep3 , &
1073 ELSE IF ( p%Type .EQ. 'd' ) THEN
1074 CALL wrf_ext_write_field ( &
1075 fid , & ! DataHandle
1076 current_date(1:19) , & ! DateStr
1077 TRIM(dname) , & ! Data Name
1078 p%dfield_3d , & ! Field
1079 WRF_DOUBLE , & ! FieldType
1080 grid%communicator , & ! Comm
1081 grid%iocommunicator , & ! Comm
1082 grid%domdesc , & ! Comm
1083 grid%bdy_mask , & ! bdy_mask
1085 TRIM(memord) , & ! MemoryOrder
1086 TRIM(p%Stagger) , & ! Stagger
1087 TRIM(p%dimname1) , & ! Dimname 1
1088 TRIM(p%dimname2) , & ! Dimname 2
1089 TRIM(p%dimname3) , & ! Dimname 3
1090 TRIM(p%Description) , & ! Desc
1091 TRIM(p%Units) , & ! Units
1092 __FILE__ // ' writing 3d double ' // TRIM(p%VarName) , & ! Debug message
1093 p%sd1 , p%ed1 , p%sd2 , p%ed2 , p%sd3 , p%ed3 , &
1094 p%sm1 , p%em1 , p%sm2 , p%em2 , p%sm3 , p%em3 , &
1095 p%sp1 , p%ep1 , p%sp2 , p%ep2 , p%sp3 , p%ep3 , &
1097 ELSE IF ( p%Type .EQ. 'i' ) THEN
1098 CALL wrf_ext_write_field ( &
1099 fid , & ! DataHandle
1100 current_date(1:19) , & ! DateStr
1101 TRIM(dname) , & ! Data Name
1102 p%ifield_3d , & ! Field
1103 WRF_INTEGER , & ! FieldType
1104 grid%communicator , & ! Comm
1105 grid%iocommunicator , & ! Comm
1106 grid%domdesc , & ! Comm
1107 grid%bdy_mask , & ! bdy_mask
1109 TRIM(memord) , & ! MemoryOrder
1110 TRIM(p%Stagger) , & ! Stagger
1111 TRIM(p%dimname1) , & ! Dimname 1
1112 TRIM(p%dimname2) , & ! Dimname 2
1113 TRIM(p%dimname3) , & ! Dimname 3
1114 TRIM(p%Description) , & ! Desc
1115 TRIM(p%Units) , & ! Units
1116 __FILE__ // ' writing 3d integer ' // TRIM(p%VarName) , & ! Debug message
1117 p%sd1 , p%ed1 , p%sd2 , p%ed2 , p%sd3 , p%ed3 , &
1118 p%sm1 , p%em1 , p%sm2 , p%em2 , p%sm3 , p%em3 , &
1119 p%sp1 , p%ep1 , p%sp2 , p%ep2 , p%sp3 , p%ep3 , &
1121 ! NOTE no io on logical arrays greater than 2d
1126 ELSE IF ( p%Ndim .EQ. 4 .AND. p%scalar_array ) THEN
1127 IF (switch.EQ.restart_only.OR.p%Ntl/100.EQ.mod(p%Ntl,100)) THEN
1128 ! Use a different write routine, wrf_ext_write_field_arr, and pass in the
1129 ! tracer indeces so that p%rfield_4d can be passsed in without arguments,
1130 ! avoiding the possiblity of a copy-in/copy-out problem for some compilers.
1132 DO itrace = PARAM_FIRST_SCALAR , p%num_table(grid%id)
1133 IF ((p%Restart.AND.switch.EQ.restart_only).OR.on_stream(p%streams_table(grid%id,itrace)%stream,newswitch)) THEN
1134 dname = p%dname_table( grid%id, itrace )
1135 IF (p%Ntl.GT.0.AND.switch.NE.restart_only)dname=dname(1:len(TRIM(dname))-2)
1136 memord = p%MemoryOrder
1137 IF ( p%Type .EQ. 'r' ) THEN
1138 CALL wrf_ext_write_field_arr ( &
1139 fid , & ! DataHandle
1140 current_date(1:19) , & ! DateStr
1141 TRIM(p%dname_table( grid%id, itrace )) , & ! Data Name
1142 p%rfield_4d , & ! Field
1143 itrace, 1, 1, 1 , & ! see comment above
1144 1, 1, 1 , & ! see comment above
1146 WRF_FLOAT , & ! FieldType
1147 grid%communicator , & ! Comm
1148 grid%iocommunicator , & ! Comm
1149 grid%domdesc , & ! Comm
1150 grid%bdy_mask , & ! bdy_mask
1152 TRIM(memord) , & ! MemoryOrder
1153 TRIM(p%Stagger) , & ! Stagger
1154 TRIM(p%dimname1) , & ! Dimname 1
1155 TRIM(p%dimname2) , & ! Dimname 2
1156 TRIM(p%dimname3) , & ! Dimname 3
1157 TRIM(p%desc_table( grid%id, itrace)) , & ! Desc
1158 TRIM(p%units_table( grid%id, itrace)) , & ! Units
1159 __FILE__ // ' writing 4d real ' // TRIM(p%dname_table(grid%id,itrace)) , & ! Debug message
1160 p%sd1 , p%ed1 , p%sd2 , p%ed2 , p%sd3 , p%ed3 , &
1161 p%sm1 , p%em1 , p%sm2 , p%em2 , p%sm3 , p%em3 , &
1162 p%sp1 , p%ep1 , p%sp2 , p%ep2 , p%sp3 , p%ep3 , &
1164 ELSE IF ( p%Type .EQ. 'd' ) THEN
1165 CALL wrf_ext_write_field_arr ( &
1166 fid , & ! DataHandle
1167 current_date(1:19) , & ! DateStr
1168 TRIM(p%dname_table( grid%id, itrace )) , & ! Data Name
1169 p%dfield_4d , & ! Field
1170 itrace, 1, 1, 1 , & ! see comment above
1171 1, 1, 1 , & ! see comment above
1173 WRF_DOUBLE , & ! FieldType
1174 grid%communicator , & ! Comm
1175 grid%iocommunicator , & ! Comm
1176 grid%domdesc , & ! Comm
1177 grid%bdy_mask , & ! bdy_mask
1179 TRIM(memord) , & ! MemoryOrder
1180 TRIM(p%Stagger) , & ! Stagger
1181 TRIM(p%dimname1) , & ! Dimname 1
1182 TRIM(p%dimname2) , & ! Dimname 2
1183 TRIM(p%dimname3) , & ! Dimname 3
1184 TRIM(p%desc_table( grid%id, itrace)) , & ! Desc
1185 TRIM(p%units_table( grid%id, itrace)) , & ! Units
1186 __FILE__ // ' writing 4d double ' // TRIM(p%dname_table(grid%id,itrace)) , & ! Debug message
1187 p%sd1 , p%ed1 , p%sd2 , p%ed2 , p%sd3 , p%ed3 , &
1188 p%sm1 , p%em1 , p%sm2 , p%em2 , p%sm3 , p%em3 , &
1189 p%sp1 , p%ep1 , p%sp2 , p%ep2 , p%sp3 , p%ep3 , &
1191 ELSE IF ( p%Type .EQ. 'i' ) THEN
1192 CALL wrf_ext_write_field_arr ( &
1193 fid , & ! DataHandle
1194 current_date(1:19) , & ! DateStr
1195 TRIM(p%dname_table( grid%id, itrace )) , & ! Data Name
1196 p%ifield_4d , & ! Field
1197 itrace, 1, 1, 1 , & ! see comment above
1198 1, 1, 1 , & ! see comment above
1200 WRF_INTEGER , & ! FieldType
1201 grid%communicator , & ! Comm
1202 grid%iocommunicator , & ! Comm
1203 grid%domdesc , & ! Comm
1204 grid%bdy_mask , & ! bdy_mask
1206 TRIM(memord) , & ! MemoryOrder
1207 TRIM(p%Stagger) , & ! Stagger
1208 TRIM(p%dimname1) , & ! Dimname 1
1209 TRIM(p%dimname2) , & ! Dimname 2
1210 TRIM(p%dimname3) , & ! Dimname 3
1211 TRIM(p%desc_table( grid%id, itrace)) , & ! Desc
1212 TRIM(p%units_table( grid%id, itrace)) , & ! Units
1213 __FILE__ // ' writing 4d integer ' // TRIM(p%dname_table(grid%id,itrace)) , & ! Debug message
1214 p%sd1 , p%ed1 , p%sd2 , p%ed2 , p%sd3 , p%ed3 , &
1215 p%sm1 , p%em1 , p%sm2 , p%em2 , p%sm3 , p%em3 , &
1216 p%sp1 , p%ep1 , p%sp2 , p%ep2 , p%sp3 , p%ep3 , &
1220 ENDDO ! loop over tracers
1221 ENDIF ! if-then-else over dim
1227 IF ( switch .EQ. boundary_only ) THEN
1228 CALL wrf_debug ( 300 , 'output_wrf: calling code in wrf_bdyout.inc' )
1229 CALL wrf_bdyout( fid , grid , config_flags, switch, dryrun, ierr )
1233 IF ( switch .EQ. history_only ) THEN
1235 current_date = current_date_save
1239 #if ( (EM_CORE == 1) && (DA_CORE != 1) )
1240 grid%save_topo_from_real = save_topo_orig
1243 IF ( .NOT. dryrun ) THEN
1244 CALL wrf_debug ( 300 , 'output_wrf: calling wrf_iosync ' )
1245 CALL wrf_iosync ( fid , ierr )
1246 CALL wrf_debug ( 300 , 'output_wrf: back from wrf_iosync ' )
1249 WRITE(wrf_err_message,*)'output_wrf: end, fid = ',fid
1250 CALL wrf_debug( 300 , wrf_err_message )
1253 END SUBROUTINE output_wrf
1255 SUBROUTINE traverse_statevars_debug (s,l)
1260 TYPE( fieldlist ), POINTER :: p
1261 p => head_grid%head_statevars%next
1262 ! write(0,*)'traverse_statevars_debug called ',TRIM(s),' ',l
1263 DO WHILE ( ASSOCIATED( p ) )
1264 ! IF ( p%Ndim .EQ. 4 .AND. p%scalar_array ) THEN
1265 ! DO itrace = PARAM_FIRST_SCALAR , p%num_table(1)
1266 !write(0,*)TRIM(s),l, itrace, p%streams_table(1,itrace)%stream
1269 ! if ( TRIM(p%VarName) .EQ. 'store_rand' .OR. TRIM(p%VarName) .EQ. 'STORE_RAND' ) then
1270 ! write(0,*)'traverse_statevars_debug sees ',TRIM(p%VarName),' >',p%Type,'<'
1275 END SUBROUTINE traverse_statevars_debug