Context Navigation

← Previous Revision
Latest Revision
Next Revision →
Normal
Revision Log

mediation_integrate.F @ 7

Last change on this file since 7 was 1, checked in by lfita, 11 years ago

-- --- Opening of the WRF+LMDZ coupling repository --- -- -

WRF: version v3.3
LMDZ: version v1818

More details in:

File size: 96.8 KB

Rev	Line
[1]	1	!
	2	!WRF:MEDIATION_LAYER:IO
	3	!
	4	#if (DA_CORE != 1)
	5
	6	SUBROUTINE med_calc_model_time ( grid , config_flags )
	7	! Driver layer
	8	USE module_domain , ONLY : domain, domain_clock_get
	9	USE module_configure , ONLY : grid_config_rec_type
	10	! Model layer
	11	USE module_date_time
	12
	13	IMPLICIT NONE
	14
	15	! Arguments
	16	TYPE(domain) :: grid
	17	TYPE (grid_config_rec_type) , INTENT(IN) :: config_flags
	18
	19	! Local data
	20	REAL :: time
	21
	22	! this is now handled by with calls to time manager
	23	! time = head_grid%dt * head_grid%total_time_steps
	24	! CALL calc_current_date (grid%id, time)
	25
	26
	27	END SUBROUTINE med_calc_model_time
	28
	29	SUBROUTINE med_before_solve_io ( grid , config_flags )
	30	! Driver layer
	31	USE module_state_description
	32	USE module_domain , ONLY : domain, domain_clock_get
	33	USE module_configure , ONLY : grid_config_rec_type
	34	USE module_streams
	35	! Model layer
	36	USE module_utility
	37
	38	IMPLICIT NONE
	39
	40	! Arguments
	41	TYPE(domain) :: grid
	42	TYPE (grid_config_rec_type) , INTENT(IN) :: config_flags
	43	! Local
	44	INTEGER :: ialarm
	45	INTEGER :: rc
	46	TYPE(WRFU_Time) :: currTime, startTime
	47	#ifdef HWRF
	48	!zhang's doing
	49	! TYPE(WRFU_Time) :: CurrTime !zhang new
	50	INTEGER :: hr, min, sec, ms,julyr,julday
	51	REAL :: GMT
	52	!end of zhang's doing
	53	#endif
	54
	55	CHARACTER*256 :: message
	56
	57	! #if (EM_CORE == 1)
	58	CALL WRFU_ClockGet( grid%domain_clock, CurrTime=currTime, StartTime=startTime )
	59	IF( WRFU_AlarmIsRinging( grid%alarms( HISTORY_ALARM ), rc=rc ) .AND. &
	60	(grid%dfi_write_dfi_history .OR. grid%dfi_stage == DFI_FST .OR. grid%dfi_opt == DFI_NODFI) ) THEN
	61	! #else
	62	! IF( WRFU_AlarmIsRinging( grid%alarms( HISTORY_ALARM ), rc=rc )) THEN
	63	! #endif
	64	IF ( (config_flags%restart) .AND. ( currTime .EQ. startTime ) ) THEN
	65	! output history at beginning of restart if alarm is ringing
	66	CALL med_hist_out ( grid , HISTORY_ALARM, config_flags )
	67	ELSE
	68	CALL med_hist_out ( grid , HISTORY_ALARM, config_flags )
	69	END IF
	70	CALL WRFU_AlarmRingerOff( grid%alarms( HISTORY_ALARM ), rc=rc )
	71	#if (EM_CORE == 1)
	72	ELSE IF ( (config_flags%restart) .AND. ( currTime .EQ. startTime ) .AND. &
	73	( config_flags%write_hist_at_0h_rst ) ) THEN
	74	! output history at beginning of restart even if alarm is not ringing
	75	CALL med_hist_out ( grid , HISTORY_ALARM, config_flags )
	76	CALL WRFU_AlarmRingerOff( grid%alarms( HISTORY_ALARM ), rc=rc )
	77	#endif
	78	ENDIF
	79
	80	IF( WRFU_AlarmIsRinging( grid%alarms( INPUTOUT_ALARM ), rc=rc ) ) THEN
	81	CALL med_filter_out ( grid , config_flags )
	82	CALL WRFU_AlarmRingerOff( grid%alarms( INPUTOUT_ALARM ), rc=rc )
	83	ENDIF
	84
	85	DO ialarm = first_auxhist, last_auxhist
	86	IF ( .FALSE.) THEN
	87	rc = 1 ! dummy statement
	88	ELSE IF( WRFU_AlarmIsRinging( grid%alarms( ialarm ), rc=rc ) ) THEN
	89	CALL med_hist_out ( grid , ialarm, config_flags )
	90	CALL WRFU_AlarmRingerOff( grid%alarms( ialarm ), rc=rc )
	91	ENDIF
	92	ENDDO
	93
	94	DO ialarm = first_auxinput, last_auxinput
	95	IF ( .FALSE.) THEN
	96	rc = 1 ! dummy statement
	97	#ifdef WRF_CHEM
	98	! - Get chemistry data
	99	ELSE IF( ialarm .EQ. AUXINPUT5_ALARM .AND. config_flags%chem_opt > 0 ) THEN
	100	IF( config_flags%emiss_inpt_opt /= 0 ) THEN
	101	IF( WRFU_AlarmIsRinging( grid%alarms( ialarm ), rc=rc ) ) THEN
	102	call wrf_debug(15,' CALL med_read_wrf_chem_emiss ')
	103	CALL med_read_wrf_chem_emiss ( grid , config_flags )
	104	CALL WRFU_AlarmRingerOff( grid%alarms( ialarm ), rc=rc )
	105	call wrf_debug(15,' Back from CALL med_read_wrf_chem_emiss ')
	106	ENDIF
	107	ELSE
	108	IF( WRFU_AlarmIsRinging( grid%alarms( ialarm ), rc=rc ) ) THEN
	109	CALL med_auxinput_in ( grid, ialarm, config_flags )
	110	CALL WRFU_AlarmRingerOff( grid%alarms( ialarm ), rc=rc )
	111	ENDIF
	112	ENDIF
	113	#endif
	114	#if ( EM_CORE == 1 )
	115	ELSE IF( ialarm .EQ. AUXINPUT11_ALARM ) THEN
	116	IF( config_flags%obs_nudge_opt .EQ. 1) THEN
	117	CALL med_fddaobs_in ( grid , config_flags )
	118	ENDIF
	119	#endif
	120	ELSE IF( WRFU_AlarmIsRinging( grid%alarms( ialarm ), rc=rc ) ) THEN
	121	CALL med_auxinput_in ( grid, ialarm, config_flags )
	122	WRITE ( message , FMT='(A,i3,A,i3)' ) 'Input data processed for aux input ' , &
	123	ialarm - first_auxinput + 1, ' for domain ',grid%id
	124	CALL wrf_debug ( 0 , message )
	125	CALL WRFU_AlarmRingerOff( grid%alarms( ialarm ), rc=rc )
	126	ENDIF
	127	ENDDO
	128
	129	! - RESTART OUTPUT
	130	CALL WRFU_ClockGet( grid%domain_clock, CurrTime=currTime, StartTime=startTime )
	131	IF ( ( WRFU_AlarmIsRinging( grid%alarms( RESTART_ALARM ), rc=rc ) ) .AND. &
	132	( currTime .NE. startTime ) ) THEN
	133	#ifdef HWRF
	134	!zhang's doing
	135	CALL domain_clock_get( grid, current_time=CurrTime )
	136	CALL WRFU_TimeGet( CurrTime, YY=julyr, dayOfYear=julday, H=hr, M=min, S=sec, MS=ms, rc=rc)
	137	gmt=hr+real(min)/60.+real(sec)/3600.+real(ms)/(1000*3600)
	138	if (grid%id .eq. 2) call med_namelist_out ( grid , config_flags )
	139	!end of zhang's doing
	140	#endif
	141	IF ( grid%id .EQ. 1 ) THEN
	142	! Only the parent initiates the restart writing. Otherwise, different
	143	! domains may be written out at different times and with different
	144	! time stamps in the file names.
	145	CALL med_restart_out ( grid , config_flags )
	146	ENDIF
	147	CALL WRFU_AlarmRingerOff( grid%alarms( RESTART_ALARM ), rc=rc )
	148	ELSE
	149	CALL WRFU_AlarmRingerOff( grid%alarms( RESTART_ALARM ), rc=rc )
	150	ENDIF
	151
	152	! - Look for boundary data after writing out history and restart files
	153	CALL med_latbound_in ( grid , config_flags )
	154
	155	RETURN
	156	END SUBROUTINE med_before_solve_io
	157
	158	SUBROUTINE med_after_solve_io ( grid , config_flags )
	159	! Driver layer
	160	USE module_domain , ONLY : domain
	161	USE module_timing
	162	USE module_configure , ONLY : grid_config_rec_type
	163	! Model layer
	164
	165	IMPLICIT NONE
	166
	167	! Arguments
	168	TYPE(domain) :: grid
	169	TYPE (grid_config_rec_type) , INTENT(IN) :: config_flags
	170
	171	! Compute time series variables
	172	CALL calc_ts(grid)
	173
	174	RETURN
	175	END SUBROUTINE med_after_solve_io
	176
	177	SUBROUTINE med_pre_nest_initial ( parent , newid , config_flags )
	178	! Driver layer
	179	#ifdef MOVE_NESTS
	180	USE module_domain , ONLY : domain, domain_clock_get
	181	#else
	182	USE module_domain , ONLY : domain
	183	#endif
	184	USE module_utility , ONLY : WRFU_Time, WRFU_TimeEQ
	185	USE module_timing
	186	USE module_io_domain
	187	USE module_configure , ONLY : grid_config_rec_type
	188	! Model layer
	189
	190	IMPLICIT NONE
	191
	192	! Arguments
	193	TYPE(domain) , POINTER :: parent
	194	INTEGER, INTENT(IN) :: newid
	195	TYPE (grid_config_rec_type) , INTENT(INOUT) :: config_flags
	196	TYPE (grid_config_rec_type) :: nest_config_flags
	197
	198	! Local
	199	INTEGER :: itmp, fid, ierr, icnt
	200	CHARACTER*256 :: rstname, message, timestr
	201
	202	TYPE(WRFU_Time) :: strt_time, cur_time
	203
	204	#ifdef MOVE_NESTS
	205
	206	CALL domain_clock_get( parent, current_timestr=timestr, start_time=strt_time, current_time=cur_time )
	207	CALL construct_filename2a ( rstname , config_flags%rst_inname , newid , 2 , timestr )
	208
	209	IF ( config_flags%restart .AND. WRFU_TimeEQ(cur_time,strt_time) ) THEN
	210	WRITE(message,*)'RESTART: nest, opening ',TRIM(rstname),' for reading header information only'
	211	CALL wrf_message ( message )
	212	! note that the parent pointer is not strictly correct, but nest is not allocated yet and
	213	! only the i/o communicator fields are used from "parent" (and those are dummies in current
	214	! implementation.
	215	CALL open_r_dataset ( fid , TRIM(rstname) , parent , config_flags , "DATASET=RESTART", ierr )
	216	IF ( ierr .NE. 0 ) THEN
	217	WRITE( message , '("program wrf: error opening ",A32," for reading")') TRIM(rstname)
	218	CALL WRF_ERROR_FATAL ( message )
	219	ENDIF
	220
	221	! update the values of parent_start that were read in from the namelist (nest may have moved)
	222	CALL wrf_get_dom_ti_integer ( fid , 'I_PARENT_START' , itmp , 1 , icnt, ierr )
	223	IF ( ierr .EQ. 0 ) THEN
	224	config_flags%i_parent_start = itmp
	225	CALL nl_set_i_parent_start ( newid , config_flags%i_parent_start )
	226	ENDIF
	227	CALL wrf_get_dom_ti_integer ( fid , 'J_PARENT_START' , itmp , 1 , icnt, ierr )
	228	IF ( ierr .EQ. 0 ) THEN
	229	config_flags%j_parent_start = itmp
	230	CALL nl_set_j_parent_start ( newid , config_flags%j_parent_start )
	231	ENDIF
	232
	233	CALL close_dataset ( fid , config_flags , "DATASET=RESTART" )
	234	ENDIF
	235	#endif
	236
	237	END SUBROUTINE med_pre_nest_initial
	238
	239
	240	SUBROUTINE med_nest_initial ( parent , nest , config_flags )
	241	! Driver layer
	242	USE module_domain , ONLY : domain , domain_clock_get , get_ijk_from_grid
	243	USE module_timing
	244	USE module_io_domain
	245	USE module_configure , ONLY : grid_config_rec_type
	246	USE module_utility
	247	! Model layer
	248
	249	IMPLICIT NONE
	250
	251	! Arguments
	252	TYPE(domain) , POINTER :: parent, nest
	253	TYPE (grid_config_rec_type) , INTENT(IN) :: config_flags
	254	TYPE (grid_config_rec_type) :: nest_config_flags
	255
	256	! Local
	257	LOGICAL, EXTERNAL :: wrf_dm_on_monitor
	258	TYPE(WRFU_Time) :: strt_time, cur_time
	259	CHARACTER * 80 :: rstname , timestr
	260	CHARACTER * 256 :: message
	261	INTEGER :: fid
	262	INTEGER :: ierr
	263	INTEGER :: i , j, rc
	264	INTEGER :: ids , ide , jds , jde , kds , kde , &
	265	ims , ime , jms , jme , kms , kme , &
	266	ips , ipe , jps , jpe , kps , kpe
	267
	268	#if (EM_CORE == 1)
	269	#ifdef MOVE_NESTS
	270	TYPE (WRFU_TimeInterval) :: interval, TimeSinceStart
	271	INTEGER :: vortex_interval , n
	272	#endif
	273	INTEGER :: save_itimestep ! This is a kludge, correct fix will
	274	! involve integrating the time-step
	275	! counting into the time manager.
	276	! JM 20040604
	277	REAL, ALLOCATABLE, DIMENSION(:,:) :: save_acsnow &
	278	,save_acsnom &
	279	,save_cuppt &
	280	,save_rainc &
	281	,save_rainnc &
	282	,save_sfcevp &
	283	,save_sfcrunoff &
	284	,save_udrunoff
	285
	286
	287	INTERFACE
	288	SUBROUTINE med_interp_domain ( parent , nest )
	289	USE module_domain , ONLY : domain
	290	TYPE(domain) , POINTER :: parent , nest
	291	END SUBROUTINE med_interp_domain
	292
	293	SUBROUTINE med_initialdata_input_ptr( nest , config_flags )
	294	USE module_domain , ONLY : domain
	295	USE module_configure , ONLY : grid_config_rec_type
	296	TYPE (grid_config_rec_type), INTENT(IN) :: config_flags
	297	TYPE(domain) , POINTER :: nest
	298	END SUBROUTINE med_initialdata_input_ptr
	299
	300	SUBROUTINE med_nest_feedback ( parent , nest , config_flags )
	301	USE module_domain , ONLY : domain
	302	USE module_configure , ONLY : grid_config_rec_type
	303	TYPE (domain), POINTER :: nest , parent
	304	TYPE (grid_config_rec_type), INTENT(IN) :: config_flags
	305	END SUBROUTINE med_nest_feedback
	306
	307	SUBROUTINE start_domain ( grid , allowed_to_move )
	308	USE module_domain , ONLY : domain
	309	TYPE(domain) :: grid
	310	LOGICAL, INTENT(IN) :: allowed_to_move
	311	END SUBROUTINE start_domain
	312
	313	SUBROUTINE blend_terrain ( ter_interpolated , ter_input , &
	314	ids , ide , jds , jde , kds , kde , &
	315	ims , ime , jms , jme , kms , kme , &
	316	ips , ipe , jps , jpe , kps , kpe )
	317	INTEGER :: ids , ide , jds , jde , kds , kde , &
	318	ims , ime , jms , jme , kms , kme , &
	319	ips , ipe , jps , jpe , kps , kpe
	320	REAL , DIMENSION(ims:ime,jms:jme) :: ter_interpolated
	321	REAL , DIMENSION(ims:ime,jms:jme) :: ter_input
	322	END SUBROUTINE blend_terrain
	323
	324	SUBROUTINE copy_3d_field ( ter_interpolated , ter_input , &
	325	ids , ide , jds , jde , kds , kde , &
	326	ims , ime , jms , jme , kms , kme , &
	327	ips , ipe , jps , jpe , kps , kpe )
	328	INTEGER :: ids , ide , jds , jde , kds , kde , &
	329	ims , ime , jms , jme , kms , kme , &
	330	ips , ipe , jps , jpe , kps , kpe
	331	REAL , DIMENSION(ims:ime,jms:jme) :: ter_interpolated
	332	REAL , DIMENSION(ims:ime,jms:jme) :: ter_input
	333	END SUBROUTINE copy_3d_field
	334
	335	SUBROUTINE input_terrain_rsmas ( grid , &
	336	ids , ide , jds , jde , kds , kde , &
	337	ims , ime , jms , jme , kms , kme , &
	338	ips , ipe , jps , jpe , kps , kpe )
	339	USE module_domain , ONLY : domain
	340	TYPE ( domain ) :: grid
	341	INTEGER :: ids , ide , jds , jde , kds , kde , &
	342	ims , ime , jms , jme , kms , kme , &
	343	ips , ipe , jps , jpe , kps , kpe
	344	END SUBROUTINE input_terrain_rsmas
	345
	346	SUBROUTINE wrf_tsin ( grid , ierr )
	347	USE module_domain
	348	TYPE ( domain ), INTENT(INOUT) :: grid
	349	INTEGER, INTENT(INOUT) :: ierr
	350	END SUBROUTINE wrf_tsin
	351
	352	END INTERFACE
	353
	354	CALL domain_clock_get( parent, start_time=strt_time, current_time=cur_time )
	355
	356	IF ( .not. ( config_flags%restart .AND. strt_time .EQ. cur_time ) ) THEN
	357	nest%first_force = .true.
	358
	359	! initialize nest with interpolated data from the parent
	360	nest%imask_nostag = 1
	361	nest%imask_xstag = 1
	362	nest%imask_ystag = 1
	363	nest%imask_xystag = 1
	364
	365	#ifdef MOVE_NESTS
	366	parent%nest_pos = parent%ht
	367	where ( parent%nest_pos .gt. 0. ) parent%nest_pos = parent%nest_pos + 500. ! make a cliff
	368	#endif
	369
	370	! initialize some other constants (and 1d arrays in z)
	371	CALL init_domain_constants ( parent, nest )
	372
	373	! fill in entire fine grid domain with interpolated coarse grid data
	374	CALL med_interp_domain( parent, nest )
	375
	376	! De-reference dimension information stored in the grid data structure.
	377	CALL get_ijk_from_grid ( nest , &
	378	ids, ide, jds, jde, kds, kde, &
	379	ims, ime, jms, jme, kms, kme, &
	380	ips, ipe, jps, jpe, kps, kpe )
	381
	382	! get the nest config flags
	383	CALL model_to_grid_config_rec ( nest%id , model_config_rec , nest_config_flags )
	384
	385	IF ( nest_config_flags%input_from_file .OR. nest_config_flags%input_from_hires ) THEN
	386
	387	WRITE(message,FMT='(A,I2,A)') '*** Initializing nest domain #',nest%id,&
	388	' from an input file. ***'
	389	CALL wrf_debug ( 0 , message )
	390
	391	! Store horizontally interpolated terrain-based fields in temp location if the input
	392	! data is from a pristine, un-cycled model input file. For the original topo from
	393	! the real program, we will need to adjust the terrain (and a couple of other base-
	394	! state fields) so reflect the smoothing and matching between the parent and child
	395	! domains.
	396
	397	CALL copy_3d_field ( nest%ht_int , nest%ht , &
	398	ids , ide , jds , jde , 1 , 1 , &
	399	ims , ime , jms , jme , 1 , 1 , &
	400	ips , ipe , jps , jpe , 1 , 1 )
	401	CALL copy_3d_field ( nest%mub_fine , nest%mub , &
	402	ids , ide , jds , jde , 1 , 1 , &
	403	ims , ime , jms , jme , 1 , 1 , &
	404	ips , ipe , jps , jpe , 1 , 1 )
	405	CALL copy_3d_field ( nest%phb_fine , nest%phb , &
	406	ids , ide , jds , jde , kds , kde , &
	407	ims , ime , jms , jme , kms , kme , &
	408	ips , ipe , jps , jpe , kps , kpe )
	409
	410	IF ( nest_config_flags%input_from_file ) THEN
	411	! read input from dataset
	412	CALL med_initialdata_input_ptr( nest , nest_config_flags )
	413
	414	ELSE IF ( nest_config_flags%input_from_hires ) THEN
	415	! read in high res topography
	416	CALL input_terrain_rsmas ( nest, &
	417	ids , ide , jds , jde , 1 , 1 , &
	418	ims , ime , jms , jme , 1 , 1 , &
	419	ips , ipe , jps , jpe , 1 , 1 )
	420	ENDIF
	421
	422	! save elevation and mub for temp and qv adjustment
	423
	424	CALL copy_3d_field ( nest%ht_fine , nest%ht , &
	425	ids , ide , jds , jde , 1 , 1 , &
	426	ims , ime , jms , jme , 1 , 1 , &
	427	ips , ipe , jps , jpe , 1 , 1 )
	428	CALL copy_3d_field ( nest%mub_save , nest%mub , &
	429	ids , ide , jds , jde , 1 , 1 , &
	430	ims , ime , jms , jme , 1 , 1 , &
	431	ips , ipe , jps , jpe , 1 , 1 )
	432
	433	! blend parent and nest fields: terrain, mub, and phb. The ht, mub and phb are used in start_domain.
	434
	435	IF ( nest%save_topo_from_real == 1 ) THEN
	436	CALL blend_terrain ( nest%ht_int , nest%ht , &
	437	ids , ide , jds , jde , 1 , 1 , &
	438	ims , ime , jms , jme , 1 , 1 , &
	439	ips , ipe , jps , jpe , 1 , 1 )
	440	CALL blend_terrain ( nest%mub_fine , nest%mub , &
	441	ids , ide , jds , jde , 1 , 1 , &
	442	ims , ime , jms , jme , 1 , 1 , &
	443	ips , ipe , jps , jpe , 1 , 1 )
	444	CALL blend_terrain ( nest%phb_fine , nest%phb , &
	445	ids , ide , jds , jde , kds , kde , &
	446	ims , ime , jms , jme , kms , kme , &
	447	ips , ipe , jps , jpe , kps , kpe )
	448	ENDIF
	449
	450	! adjust temp and qv
	451
	452	CALL adjust_tempqv ( nest%mub , nest%mub_save , &
	453	nest%znw , nest%p_top , &
	454	nest%t_2 , nest%p , nest%moist(ims,kms,jms,P_QV) , &
	455	ids , ide , jds , jde , kds , kde , &
	456	ims , ime , jms , jme , kms , kme , &
	457	ips , ipe , jps , jpe , kps , kpe )
	458
	459	ELSE
	460	WRITE(message,FMT='(A,I2,A,I2,A)') '*** Initializing nest domain #',nest%id,&
	461	' by horizontally interpolating parent domain #' ,parent%id, &
	462	'. ***'
	463	CALL wrf_debug ( 0 , message )
	464
	465	#if (DA_CORE != 1)
	466	! For nests without an input file, we still need to read time series locations
	467	! from the tslist file
	468	CALL wrf_tsin( nest , ierr )
	469	#endif
	470	END IF
	471
	472
	473	! feedback, mostly for this new terrain, but it is the safe thing to do
	474	parent%ht_coarse = parent%ht
	475
	476	CALL med_nest_feedback ( parent , nest , config_flags )
	477
	478	! set some other initial fields, fill out halos, base fields; re-do parent due
	479	! to new terrain elevation from feedback
	480	nest%imask_nostag = 1
	481	nest%imask_xstag = 1
	482	nest%imask_ystag = 1
	483	nest%imask_xystag = 1
	484	nest%press_adj = .TRUE.
	485	CALL start_domain ( nest , .TRUE. )
	486	! kludge: 20040604
	487	CALL get_ijk_from_grid ( parent , &
	488	ids, ide, jds, jde, kds, kde, &
	489	ims, ime, jms, jme, kms, kme, &
	490	ips, ipe, jps, jpe, kps, kpe )
	491
	492	ALLOCATE( save_acsnow(ims:ime,jms:jme) )
	493	ALLOCATE( save_acsnom(ims:ime,jms:jme) )
	494	ALLOCATE( save_cuppt(ims:ime,jms:jme) )
	495	ALLOCATE( save_rainc(ims:ime,jms:jme) )
	496	ALLOCATE( save_rainnc(ims:ime,jms:jme) )
	497	ALLOCATE( save_sfcevp(ims:ime,jms:jme) )
	498	ALLOCATE( save_sfcrunoff(ims:ime,jms:jme) )
	499	ALLOCATE( save_udrunoff(ims:ime,jms:jme) )
	500	save_acsnow = parent%acsnow
	501	save_acsnom = parent%acsnom
	502	save_cuppt = parent%cuppt
	503	save_rainc = parent%rainc
	504	save_rainnc = parent%rainnc
	505	save_sfcevp = parent%sfcevp
	506	save_sfcrunoff = parent%sfcrunoff
	507	save_udrunoff = parent%udrunoff
	508	save_itimestep = parent%itimestep
	509	parent%imask_nostag = 1
	510	parent%imask_xstag = 1
	511	parent%imask_ystag = 1
	512	parent%imask_xystag = 1
	513
	514	parent%press_adj = .FALSE.
	515	CALL start_domain ( parent , .TRUE. )
	516
	517	parent%acsnow = save_acsnow
	518	parent%acsnom = save_acsnom
	519	parent%cuppt = save_cuppt
	520	parent%rainc = save_rainc
	521	parent%rainnc = save_rainnc
	522	parent%sfcevp = save_sfcevp
	523	parent%sfcrunoff = save_sfcrunoff
	524	parent%udrunoff = save_udrunoff
	525	parent%itimestep = save_itimestep
	526	DEALLOCATE( save_acsnow )
	527	DEALLOCATE( save_acsnom )
	528	DEALLOCATE( save_cuppt )
	529	DEALLOCATE( save_rainc )
	530	DEALLOCATE( save_rainnc )
	531	DEALLOCATE( save_sfcevp )
	532	DEALLOCATE( save_sfcrunoff )
	533	DEALLOCATE( save_udrunoff )
	534	! end of kludge: 20040604
	535
	536
	537	ELSE ! restart
	538
	539	IF ( wrf_dm_on_monitor() ) CALL start_timing
	540
	541	CALL domain_clock_get( nest, current_timestr=timestr )
	542	CALL construct_filename2a ( rstname , config_flags%rst_inname , nest%id , 2 , timestr )
	543
	544	WRITE(message,*)'RESTART: nest, opening ',TRIM(rstname),' for reading'
	545	CALL wrf_message ( message )
	546	CALL model_to_grid_config_rec ( nest%id , model_config_rec , nest_config_flags )
	547	CALL open_r_dataset ( fid , TRIM(rstname) , nest , nest_config_flags , "DATASET=RESTART", ierr )
	548	IF ( ierr .NE. 0 ) THEN
	549	WRITE( message , '("program wrf: error opening ",A32," for reading")') TRIM(rstname)
	550	CALL WRF_ERROR_FATAL ( message )
	551	ENDIF
	552	CALL input_restart ( fid, nest , nest_config_flags , ierr )
	553	CALL close_dataset ( fid , nest_config_flags , "DATASET=RESTART" )
	554
	555	IF ( wrf_dm_on_monitor() ) THEN
	556	WRITE ( message , FMT = '("processing restart file for domain ",I8)' ) nest%id
	557	CALL end_timing ( TRIM(message) )
	558	ENDIF
	559
	560	nest%imask_nostag = 1
	561	nest%imask_xstag = 1
	562	nest%imask_ystag = 1
	563	nest%imask_xystag = 1
	564	nest%press_adj = .FALSE.
	565	CALL start_domain ( nest , .TRUE. )
	566	#ifndef MOVE_NESTS
	567	! this doesn't need to be done for moving nests, since ht_coarse is part of the restart
	568	parent%ht_coarse = parent%ht
	569	#else
	570	# if 1
	571	! In case of a restart, assume that the movement has already occurred in the previous
	572	! run and turn off the alarm for the starting time. We must impose a requirement that the
	573	! run be restarted on-interval. Test for that and print a warning if it isn't.
	574	! Note, simulation_start, etc. should be available as metadata in the restart file, and
	575	! these will have gotten, set, and retrievable as rconfig data been set in share/input_wrf.F
	576	! using the nl_get routines below. JM 20060314
	577
	578	CALL nl_get_vortex_interval ( nest%id , vortex_interval )
	579	CALL WRFU_TimeIntervalSet( interval, M=vortex_interval, rc=rc )
	580
	581	CALL domain_clock_get( nest, timeSinceSimulationStart=TimeSinceStart )
	582	n = WRFU_TimeIntervalDIVQuot( TimeSinceStart , interval )
	583	IF ( ( interval * n ) .NE. TimeSinceStart ) THEN
	584	CALL wrf_message('WARNING: Restart is not on a vortex_interval time boundary.')
	585	CALL wrf_message('The code will work but results will not agree exactly with a ')
	586	CALL wrf_message('a run that was done straight-through, without a restart.')
	587	ENDIF
	588	!! In case of a restart, assume that the movement has already occurred in the previous
	589	!! run and turn off the alarm for the starting time. We must impose a requirement that the
	590	!! run be restarted on-interval. Test for that and print a warning if it isn't.
	591	!! Note, simulation_start, etc. should be available as metadata in the restart file, and
	592	!! these will have gotten, set, and retrievable as rconfig data been set in share/input_wrf.F
	593	!! using the nl_get routines below. JM 20060314
	594	! CALL WRFU_AlarmRingerOff( nest%alarms( COMPUTE_VORTEX_CENTER_ALARM ), rc=rc )
	595
	596	# else
	597	! this code, currently commented out, is an attempt to have the
	598	! vortex centering interval be set according to simulation start
	599	! time (rather than run start time) in case of a restart. But
	600	! there are other problems (the WRF clock is currently using
	601	! run-start as it's start time) so the alarm still would not fire
	602	! right if the model were started off-interval. Leave it here and
	603	! enable when the clock is changed to use sim-start for start time.
	604	! JM 20060314
	605	CALL nl_get_vortex_interval ( nest%id , vortex_interval )
	606	CALL WRFU_TimeIntervalSet( interval, M=vortex_interval, rc=rc )
	607
	608	CALL domain_clock_get( nest, timeSinceSimulationStart=TimeSinceStart )
	609
	610	CALL domain_alarm_create( nest, COMPUTE_VORTEX_CENTER_ALARM, interval )
	611	CALL WRFU_AlarmEnable( nest%alarms( COMPUTE_VORTEX_CENTER_ALARM ), rc=rc )
	612	n = WRFU_TimeIntervalDIVQuot( TimeSinceStart , interval )
	613	IF ( ( interval * n ) .EQ. TimeSinceStart ) THEN
	614	CALL WRFU_AlarmRingerOn( nest%alarms( COMPUTE_VORTEX_CENTER_ALARM ), rc=rc )
	615	ELSE
	616	CALL WRFU_AlarmRingerOff( nest%alarms( COMPUTE_VORTEX_CENTER_ALARM ), rc=rc )
	617	ENDIF
	618	# endif
	619	#endif
	620
	621	ENDIF
	622
	623	#endif
	624
	625	#if (NMM_CORE == 1 && NMM_NEST == 1)
	626	!===================================================================================
	627	! Added for the NMM core. This is gopal's doing.
	628	!===================================================================================
	629
	630	INTERFACE
	631
	632	SUBROUTINE med_nest_egrid_configure ( parent , nest )
	633	USE module_domain , ONLY : domain
	634	TYPE(domain) , POINTER :: parent , nest
	635	END SUBROUTINE med_nest_egrid_configure
	636
	637	SUBROUTINE med_construct_egrid_weights ( parent , nest )
	638	USE module_domain , ONLY : domain
	639	TYPE(domain) , POINTER :: parent , nest
	640	END SUBROUTINE med_construct_egrid_weights
	641
	642	SUBROUTINE BASE_STATE_PARENT ( Z3d,Q3d,T3d,PSTD, &
	643	PINT,T,Q,CWM, &
	644	FIS,QSH,PD,PDTOP,PTOP, &
	645	ETA1,ETA2, &
	646	DETA1,DETA2, &
	647	IDS,IDE,JDS,JDE,KDS,KDE, &
	648	IMS,IME,JMS,JME,KMS,KME, &
	649	IPS,IPE,JPS,JPE,KPS,KPE )
	650	!
	651
	652	USE MODULE_MODEL_CONSTANTS
	653	IMPLICIT NONE
	654	INTEGER, INTENT(IN ) :: IDS,IDE,JDS,JDE,KDS,KDE
	655	INTEGER, INTENT(IN ) :: IMS,IME,JMS,JME,KMS,KME
	656	INTEGER, INTENT(IN ) :: IPS,IPE,JPS,JPE,KPS,KPE
	657	REAL, INTENT(IN ) :: PDTOP,PTOP
	658	REAL, DIMENSION(KMS:KME), INTENT(IN) :: ETA1,ETA2,DETA1,DETA2
	659	REAL, DIMENSION(IMS:IME,JMS:JME), INTENT(IN) :: FIS,PD,QSH
	660	REAL, DIMENSION(IMS:IME,JMS:JME,KMS:KME), INTENT(IN) :: PINT,T,Q,CWM
	661	REAL, DIMENSION(KMS:KME) , INTENT(OUT):: PSTD
	662	REAL, DIMENSION(IMS:IME,JMS:JME,KMS:KME), INTENT(OUT):: Z3d,Q3d,T3d
	663
	664	END SUBROUTINE BASE_STATE_PARENT
	665
	666	SUBROUTINE NEST_TERRAIN ( nest, config_flags )
	667	USE module_domain , ONLY : domain
	668	USE module_configure , ONLY : grid_config_rec_type
	669	TYPE(domain) , POINTER :: nest
	670	TYPE(grid_config_rec_type) , INTENT(IN) :: config_flags
	671	END SUBROUTINE NEST_TERRAIN
	672
	673	SUBROUTINE med_interp_domain ( parent , nest )
	674	USE module_domain , ONLY : domain
	675	TYPE(domain) , POINTER :: parent , nest
	676	END SUBROUTINE med_interp_domain
	677
	678	SUBROUTINE med_init_domain_constants_nmm ( parent, nest )
	679	USE module_domain , ONLY : domain
	680	TYPE(domain) , POINTER :: parent , nest
	681	END SUBROUTINE med_init_domain_constants_nmm
	682
	683	SUBROUTINE start_domain ( grid , allowed_to_move )
	684	USE module_domain , ONLY : domain
	685	TYPE(domain) :: grid
	686	LOGICAL, INTENT(IN) :: allowed_to_move
	687	END SUBROUTINE start_domain
	688
	689	END INTERFACE
	690
	691	#ifdef HWRF
	692	!zhang's doing test
	693	if (config_flags%restart .or. nest%analysis) then
	694	nest%first_force = .true.
	695	else
	696	nest%first_force = .false.
	697	endif
	698	!end of zhang's doing
	699
	700	!zhang's doing for analysis option
	701	IF(.not. nest%analysis .and. .not. config_flags%restart)THEN ! initialize for cold-start
	702	#endif
	703
	704	!----------------------------------------------------------------------------
	705	! initialize nested domain configurations including setting up wbd,sbd, etc
	706	!----------------------------------------------------------------------------
	707
	708	CALL med_nest_egrid_configure ( parent , nest )
	709
	710	!-------------------------------------------------------------------------
	711	! initialize lat-lons and determine weights
	712	!-------------------------------------------------------------------------
	713
	714	CALL med_construct_egrid_weights ( parent, nest )
	715	!
	716	!
	717	! De-reference dimension information stored in the grid data structure.
	718	!
	719	! From the hybrid, construct the GPMs on isobaric surfaces and then interpolate those
	720	! values on to the nested domain. 23 standard prssure levels are assumed here. For
	721	! levels below ground, lapse rate atmosphere is assumed before the use of vertical
	722	! spline interpolation
	723	!
	724
	725
	726	IDS = parent%sd31
	727	IDE = parent%ed31
	728	JDS = parent%sd32
	729	JDE = parent%ed32
	730	KDS = parent%sd33
	731	KDE = parent%ed33
	732
	733	IMS = parent%sm31
	734	IME = parent%em31
	735	JMS = parent%sm32
	736	JME = parent%em32
	737	KMS = parent%sm33
	738	KME = parent%em33
	739
	740	IPS = parent%sp31
	741	IPE = parent%ep31
	742	JPS = parent%sp32
	743	JPE = parent%ep32
	744	KPS = parent%sp33
	745	KPE = parent%ep33
	746
	747	CALL BASE_STATE_PARENT ( parent%Z3d,parent%Q3d,parent%T3d,parent%PSTD, &
	748	parent%PINT,parent%T,parent%Q,parent%CWM, &
	749	parent%FIS,parent%QSH,parent%PD,parent%pdtop,parent%pt, &
	750	parent%ETA1,parent%ETA2, &
	751	parent%DETA1,parent%DETA2, &
	752	IDS,IDE,JDS,JDE,KDS,KDE, &
	753	IMS,IME,JMS,JME,KMS,KME, &
	754	IPS,IPE,JPS,JPE,KPS,KPE )
	755
	756	!
	757	! Set new terrain. Since some terrain adjustment is done within the interpolation calls
	758	! at the next step, the new terrain over the nested domain has to be called here.
	759	!
	760	IDS = nest%sd31
	761	IDE = nest%ed31
	762	JDS = nest%sd32
	763	JDE = nest%ed32
	764	KDS = nest%sd33
	765	KDE = nest%ed33
	766
	767	IMS = nest%sm31
	768	IME = nest%em31
	769	JMS = nest%sm32
	770	JME = nest%em32
	771	KMS = nest%sm33
	772	KME = nest%em33
	773
	774	IPS = nest%sp31
	775	IPE = nest%ep31
	776	JPS = nest%sp32
	777	JPE = nest%ep32
	778	KPS = nest%sp33
	779	KPE = nest%ep33
	780
	781
	782	CALL NEST_TERRAIN ( nest, config_flags )
	783
	784	! Initialize some more constants required especially for terrain adjustment processes
	785
	786	nest%PSTD=parent%PSTD
	787	nest%KZMAX=KME
	788	parent%KZMAX=KME ! just for safety
	789
	790	DO J = JPS, MIN(JPE,JDE-1)
	791	DO I = IPS, MIN(IPE,IDE-1)
	792	nest%fis(I,J)=nest%hres_fis(I,J)
	793	ENDDO
	794	ENDDO
	795
	796	!--------------------------------------------------------------------------
	797	! interpolation call
	798	!--------------------------------------------------------------------------
	799
	800	! initialize nest with interpolated data from the parent
	801
	802	nest%imask_nostag = 0
	803	nest%imask_xstag = 0
	804	nest%imask_ystag = 0
	805	nest%imask_xystag = 0
	806
	807	#ifdef HWRF
	808	CALL med_interp_domain( parent, nest )
	809	#else
	810	CALL domain_clock_get( parent, start_time=strt_time, current_time=cur_time )
	811
	812	IF ( .not. ( config_flags%restart .AND. strt_time .EQ. cur_time ) ) THEN
	813
	814	CALL med_interp_domain( parent, nest )
	815
	816	ELSE
	817
	818	CALL domain_clock_get( nest, current_timestr=timestr )
	819	CALL construct_filename2a ( rstname , config_flags%rst_inname , nest%id , 2 , timestr )
	820
	821	WRITE(message,*)'RESTART: nest, opening ',TRIM(rstname),' for reading'
	822	CALL wrf_message ( message )
	823	CALL model_to_grid_config_rec ( nest%id , model_config_rec , nest_config_flags )
	824	CALL open_r_dataset ( fid , TRIM(rstname) , nest , nest_config_flags , "DATASET=RESTART", ierr )
	825	IF ( ierr .NE. 0 ) THEN
	826	WRITE( message , '("program wrf: error opening ",A32," for reading")') TRIM(rstname)
	827	CALL WRF_ERROR_FATAL ( message )
	828	ENDIF
	829	CALL input_restart ( fid, nest , nest_config_flags , ierr )
	830	CALL close_dataset ( fid , nest_config_flags , "DATASET=RESTART" )
	831
	832	END IF
	833
	834	#endif
	835	!------------------------------------------------------------------------------
	836	! set up constants (module_initialize_real.F for nested nmm domain)
	837	!-----------------------------------------------------------------------------
	838
	839	CALL med_init_domain_constants_nmm ( parent, nest )
	840
	841	!--------------------------------------------------------------------------------------
	842	! set some other initial fields, fill out halos, etc.
	843	!--------------------------------------------------------------------------------------
	844
	845	CALL start_domain ( nest, .TRUE.)
	846
	847	#ifdef HWRF
	848	!zhang's doing: else for analysis or restart option
	849
	850	!zhang test
	851	CALL nl_set_isice ( nest%id , config_flags%isice )
	852	CALL nl_set_isoilwater ( nest%id , config_flags%isoilwater )
	853	CALL nl_set_isurban ( nest%id , config_flags%isurban )
	854	CALL nl_set_gmt ( nest%id , config_flags%gmt )
	855	CALL nl_set_julyr (nest%id, config_flags%julyr)
	856	CALL nl_set_julday ( nest%id , config_flags%julday )
	857	!zhang test ends
	858	CALL med_analysis_out ( nest, config_flags )
	859
	860	ELSE
	861
	862	!------------------------------------------------------------------------------------
	863	! read in analysis (equivalent of restart for the nested domains)
	864	!------------------------------------------------------------------------------------
	865
	866	!zhang's doing
	867	IF( nest%analysis .and. .not. config_flags%restart)THEN
	868	CALL med_analysis_in ( nest, config_flags )
	869	ELSE IF (config_flags%restart)THEN
	870	CALL med_restart_in ( nest, config_flags )
	871	ENDIF
	872	!end of zhang's doing
	873
	874	!----------------------------------------------------------------------------
	875	! initialize nested domain configurations including setting up wbd,sbd, etc
	876	!----------------------------------------------------------------------------
	877
	878	CALL med_nest_egrid_configure ( parent , nest )
	879
	880	!-------------------------------------------------------------------------
	881	! initialize lat-lons and determine weights (overwrite for safety)
	882	!-------------------------------------------------------------------------
	883
	884	CALL med_construct_egrid_weights ( parent, nest )
	885
	886	nest%imask_nostag = 0
	887	nest%imask_xstag = 0
	888	nest%imask_ystag = 0
	889	nest%imask_xystag = 0
	890
	891	!------------------------------------------------------------------------------
	892	! set up constants (module_initialize_real.F for nested nmm domain)
	893	!-----------------------------------------------------------------------------
	894
	895	CALL med_init_domain_constants_nmm ( parent, nest )
	896
	897	!--------------------------------------------------------------------------------------
	898	! set some other initial fields, fill out halos, etc. (again, safety sake only)
	899	! Also, in order to accomodate some physics initialization after nest move, set
	900	! analysis back to false for future use
	901	!--------------------------------------------------------------------------------------
	902
	903	CALL start_domain ( nest, .TRUE.)
	904
	905	nest%analysis=.FALSE.
	906	CALL nl_set_analysis( nest%id, nest%analysis)
	907
	908	ENDIF
	909
	910	#endif
	911
	912	!===================================================================================
	913	! Added for the NMM core. End of gopal's doing.
	914	!===================================================================================
	915	#endif
	916	RETURN
	917	END SUBROUTINE med_nest_initial
	918
	919	SUBROUTINE init_domain_constants ( parent , nest )
	920	USE module_domain , ONLY : domain
	921	IMPLICIT NONE
	922	TYPE(domain) :: parent , nest
	923	#if (EM_CORE == 1)
	924	CALL init_domain_constants_em ( parent, nest )
	925	#endif
	926	END SUBROUTINE init_domain_constants
	927
	928
	929	SUBROUTINE med_nest_force ( parent , nest )
	930	! Driver layer
	931	USE module_domain , ONLY : domain
	932	USE module_timing
	933	USE module_configure , ONLY : grid_config_rec_type
	934	! Model layer
	935	! External
	936	USE module_utility
	937
	938	IMPLICIT NONE
	939
	940	! Arguments
	941	TYPE(domain) , POINTER :: parent, nest
	942	! Local
	943	INTEGER :: idum1 , idum2 , fid, rc
	944
	945	#if (NMM_CORE == 1 && NMM_NEST == 1)
	946	INTEGER :: IDS,IDE,JDS,JDE,KDS,KDE ! gopal
	947	INTEGER :: IMS,IME,JMS,JME,KMS,KME
	948	INTEGER :: ITS,ITE,JTS,JTE,KTS,KTE
	949	#endif
	950
	951	INTERFACE
	952	SUBROUTINE med_force_domain ( parent , nest )
	953	USE module_domain , ONLY : domain
	954	TYPE(domain) , POINTER :: parent , nest
	955	END SUBROUTINE med_force_domain
	956	SUBROUTINE med_interp_domain ( parent , nest )
	957	USE module_domain , ONLY : domain
	958	TYPE(domain) , POINTER :: parent , nest
	959	END SUBROUTINE med_interp_domain
	960	#if (NMM_CORE == 1 && NMM_NEST == 1)
	961	!===================================================================================
	962	! Added for the NMM core. This is gopal's doing.
	963	!===================================================================================
	964
	965	SUBROUTINE BASE_STATE_PARENT ( Z3d,Q3d,T3d,PSTD, &
	966	PINT,T,Q,CWM, &
	967	FIS,QSH,PD,PDTOP,PTOP, &
	968	ETA1,ETA2, &
	969	DETA1,DETA2, &
	970	IDS,IDE,JDS,JDE,KDS,KDE, &
	971	IMS,IME,JMS,JME,KMS,KME, &
	972	ITS,ITE,JTS,JTE,KTS,KTE )
	973	!
	974
	975	USE MODULE_MODEL_CONSTANTS
	976	IMPLICIT NONE
	977	INTEGER, INTENT(IN ) :: IDS,IDE,JDS,JDE,KDS,KDE
	978	INTEGER, INTENT(IN ) :: IMS,IME,JMS,JME,KMS,KME
	979	INTEGER, INTENT(IN ) :: ITS,ITE,JTS,JTE,KTS,KTE
	980	REAL, INTENT(IN ) :: PDTOP,PTOP
	981	REAL, DIMENSION(KMS:KME), INTENT(IN) :: ETA1,ETA2,DETA1,DETA2
	982	REAL, DIMENSION(IMS:IME,JMS:JME), INTENT(IN) :: FIS,PD,QSH
	983	REAL, DIMENSION(IMS:IME,JMS:JME,KMS:KME), INTENT(IN) :: PINT,T,Q,CWM
	984	REAL, DIMENSION(KMS:KME) , INTENT(OUT):: PSTD
	985	REAL, DIMENSION(IMS:IME,JMS:JME,KMS:KME), INTENT(OUT):: Z3d,Q3d,T3d
	986
	987	END SUBROUTINE BASE_STATE_PARENT
	988
	989	#endif
	990	END INTERFACE
	991
	992	#if (NMM_CORE == 1 && NMM_NEST == 1)
	993
	994	! De-reference dimension information stored in the grid data structure.
	995
	996	IDS = parent%sd31
	997	IDE = parent%ed31
	998	JDS = parent%sd32
	999	JDE = parent%ed32
	1000	KDS = parent%sd33
	1001	KDE = parent%ed33
	1002
	1003	IMS = parent%sm31
	1004	IME = parent%em31
	1005	JMS = parent%sm32
	1006	JME = parent%em32
	1007	KMS = parent%sm33
	1008	KME = parent%em33
	1009
	1010	ITS = parent%sp31
	1011	ITE = parent%ep31
	1012	JTS = parent%sp32
	1013	JTE = parent%ep32
	1014	KTS = parent%sp33
	1015	KTE = parent%ep33
	1016
	1017
	1018	CALL BASE_STATE_PARENT ( parent%Z3d,parent%Q3d,parent%T3d,parent%PSTD, &
	1019	parent%PINT,parent%T,parent%Q,parent%CWM, &
	1020	parent%FIS,parent%QSH,parent%PD,parent%pdtop,parent%pt, &
	1021	parent%ETA1,parent%ETA2, &
	1022	parent%DETA1,parent%DETA2, &
	1023	IDS,IDE,JDS,JDE,KDS,KDE, &
	1024	IMS,IME,JMS,JME,KMS,KME, &
	1025	ITS,ITE,JTS,JTE,KTS,KTE )
	1026
	1027	#endif
	1028
	1029	IF ( .NOT. WRFU_ClockIsStopTime(nest%domain_clock ,rc=rc) ) THEN
	1030	! initialize nest with interpolated data from the parent
	1031	nest%imask_nostag = 1
	1032	nest%imask_xstag = 1
	1033	nest%imask_ystag = 1
	1034	nest%imask_xystag = 1
	1035	CALL med_force_domain( parent, nest )
	1036	ENDIF
	1037
	1038	! might also have calls here to do input from a file into the nest
	1039
	1040	RETURN
	1041	END SUBROUTINE med_nest_force
	1042
	1043	SUBROUTINE med_nest_feedback ( parent , nest , config_flags )
	1044	! Driver layer
	1045	USE module_domain , ONLY : domain , get_ijk_from_grid
	1046	USE module_timing
	1047	USE module_configure , ONLY : grid_config_rec_type
	1048	! Model layer
	1049	! External
	1050	USE module_utility
	1051	IMPLICIT NONE
	1052
	1053
	1054	! Arguments
	1055	TYPE(domain) , POINTER :: parent, nest
	1056	TYPE (grid_config_rec_type) , INTENT(IN) :: config_flags
	1057	! Local
	1058	INTEGER :: idum1 , idum2 , fid, rc
	1059	INTEGER :: ids , ide , jds , jde , kds , kde , &
	1060	ims , ime , jms , jme , kms , kme , &
	1061	ips , ipe , jps , jpe , kps , kpe
	1062	INTEGER i,j
	1063
	1064	INTERFACE
	1065	SUBROUTINE med_feedback_domain ( parent , nest )
	1066	USE module_domain , ONLY : domain
	1067	TYPE(domain) , POINTER :: parent , nest
	1068	END SUBROUTINE med_feedback_domain
	1069	END INTERFACE
	1070
	1071	! feedback nest to the parent
	1072	IF ( .NOT. WRFU_ClockIsStopTime(nest%domain_clock ,rc=rc) .AND. &
	1073	config_flags%feedback .NE. 0 ) THEN
	1074	CALL med_feedback_domain( parent, nest )
	1075	#ifdef MOVE_NESTS
	1076	CALL get_ijk_from_grid ( parent , &
	1077	ids, ide, jds, jde, kds, kde, &
	1078	ims, ime, jms, jme, kms, kme, &
	1079	ips, ipe, jps, jpe, kps, kpe )
	1080	! gopal's change- added ifdef
	1081	#if ( EM_CORE == 1 )
	1082	DO j = jps, MIN(jpe,jde-1)
	1083	DO i = ips, MIN(ipe,ide-1)
	1084	IF ( parent%nest_pos(i,j) .EQ. 9021000. ) THEN
	1085	parent%nest_pos(i,j) = parent%ht(i,j)*1.5 + 1000.
	1086	ELSE IF ( parent%ht(i,j) .NE. 0. ) THEN
	1087	parent%nest_pos(i,j) = parent%ht(i,j) + 500.
	1088	ELSE
	1089	parent%nest_pos(i,j) = 0.
	1090	ENDIF
	1091	ENDDO
	1092	ENDDO
	1093	#endif
	1094	#endif
	1095	END IF
	1096
	1097	RETURN
	1098	END SUBROUTINE med_nest_feedback
	1099
	1100	SUBROUTINE med_last_solve_io ( grid , config_flags )
	1101	! Driver layer
	1102	USE module_state_description
	1103	USE module_domain , ONLY : domain, domain_clock_get
	1104	USE module_configure , ONLY : grid_config_rec_type
	1105	USE module_utility
	1106	USE module_streams
	1107	! Model layer
	1108
	1109	IMPLICIT NONE
	1110
	1111	! Arguments
	1112	TYPE(domain) :: grid
	1113	TYPE (grid_config_rec_type) , INTENT(IN) :: config_flags
	1114	! Local
	1115	INTEGER :: rc
	1116	#ifdef HWRF
	1117	!zhang's doing
	1118	TYPE(WRFU_Time) :: CurrTime !zhang new
	1119	INTEGER :: hr, min, sec, ms,julyr,julday
	1120	REAL :: GMT
	1121	!end of zhang's doing
	1122	#endif
	1123
	1124	! #if (EM_CORE == 1)
	1125	IF( WRFU_AlarmIsRinging( grid%alarms( HISTORY_ALARM ), rc=rc ) .AND. &
	1126	(grid%dfi_write_dfi_history .OR. grid%dfi_stage == DFI_FST .OR. grid%dfi_opt == DFI_NODFI) ) THEN
	1127	! #else
	1128	! IF( WRFU_AlarmIsRinging( grid%alarms( HISTORY_ALARM ), rc=rc )) THEN
	1129	! #endif
	1130	CALL med_hist_out ( grid , HISTORY_ALARM , config_flags )
	1131	ENDIF
	1132
	1133	IF( WRFU_AlarmIsRinging( grid%alarms( INPUTOUT_ALARM ), rc=rc ) ) THEN
	1134	CALL med_filter_out ( grid , config_flags )
	1135	ENDIF
	1136
	1137	! registry-generated file of the following
	1138	! IF( WRFU_AlarmIsRinging( grid%alarms( AUXHIST1_ALARM ), rc=rc ) ) THEN
	1139	! CALL med_hist_out ( grid , AUXHIST1_ALARM , config_flags )
	1140	! ENDIF
	1141	#include "med_last_solve_io.inc"
	1142
	1143	! - RESTART OUTPUT
	1144	IF( WRFU_AlarmIsRinging( grid%alarms( RESTART_ALARM ), rc=rc ) ) THEN
	1145	#ifdef HWRF
	1146	!zhang's doing
	1147	!zhang new CALL ESMF_TimeGet( grid%current_time, YY=julyr, dayOfYear=julday, H=hr, M=min, S=sec, MS=ms, rc=rc)
	1148	CALL domain_clock_get( grid, current_time=CurrTime )
	1149	CALL WRFU_TimeGet( CurrTime, YY=julyr, dayOfYear=julday, H=hr, M=min, S=sec, MS=ms, rc=rc)
	1150	gmt=hr+real(min)/60.+real(sec)/3600.+real(ms)/(1000*3600)
	1151	if (grid%id .eq. 2) call med_namelist_out ( grid , config_flags )
	1152	!end of zhang's doing
	1153	#endif
	1154	IF ( grid%id .EQ. 1 ) THEN
	1155	CALL med_restart_out ( grid , config_flags )
	1156	ENDIF
	1157	ENDIF
	1158
	1159	! Write out time series
	1160	CALL write_ts( grid )
	1161
	1162	RETURN
	1163	END SUBROUTINE med_last_solve_io
	1164
	1165	#endif
	1166
	1167	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
	1168
	1169	#ifdef HWRF
	1170	!==================================================================================
	1171	! Added for the NMM 3d var. This is simply an extension of med_restart_out.
	1172	! The file is simply called wrfanal***. This is gopal's doing
	1173	!===================================================================================
	1174	!
	1175	SUBROUTINE med_analysis_in ( grid , config_flags )
	1176	! Driver layer
	1177	USE module_domain , ONLY : domain, domain_clock_get
	1178	USE module_io_domain
	1179	USE module_timing
	1180	! Model layer
	1181	USE module_configure , ONLY : grid_config_rec_type
	1182	USE module_bc_time_utilities
	1183	!zhang USE WRF_ESMF_MOD
	1184
	1185	IMPLICIT NONE
	1186
	1187	! Arguments
	1188	TYPE(domain) :: grid
	1189	TYPE (grid_config_rec_type) , INTENT(IN) :: config_flags
	1190
	1191	! Local
	1192	LOGICAL, EXTERNAL :: wrf_dm_on_monitor
	1193	CHARACTER*80 :: rstname , outname
	1194	INTEGER :: fid , rid
	1195	CHARACTER (LEN=256) :: message
	1196	INTEGER :: ierr
	1197	INTEGER :: myproc
	1198	!zhang old TYPE(ESMF_Time) :: CurrTime
	1199	TYPE(WRFU_Time) :: CurrTime
	1200	CHARACTER*80 :: timestr
	1201
	1202	IF ( wrf_dm_on_monitor() ) THEN
	1203	CALL start_timing
	1204	END IF
	1205
	1206	rid=grid%id
	1207
	1208	!zhang's doing CALL ESMF_ClockGet( grid%domain_clock, CurrTime=CurrTime, rc=ierr )
	1209	!zhang's doing CALL wrf_timetoa ( CurrTime, timestr )
	1210	CALL domain_clock_get( grid, current_timestr=timestr )
	1211	CALL construct_filename2a ( rstname ,config_flags%anl_outname, grid%id , 2 , timestr )
	1212
	1213	WRITE( message , '("med_analysis_in: opening ",A," for reading")' ) TRIM ( rstname )
	1214	CALL wrf_debug( 1 , message )
	1215	CALL open_r_dataset ( rid, TRIM(rstname), grid , &
	1216	config_flags , "DATASET=RESTART", ierr )
	1217
	1218	IF ( ierr .NE. 0 ) THEN
	1219	! CALL WRF_message( message )
	1220	CALL WRF_ERROR_FATAL('NESTED DOMAIN ERROR: FOR ANALYSIS SET TO TRUE, YOU NEED wrfanal FILE')
	1221	ENDIF
	1222	CALL input_restart ( rid, grid , config_flags , ierr )
	1223	IF ( wrf_dm_on_monitor() ) THEN
	1224	WRITE ( message , FMT = '("Reading restart for domain ",I8)' ) grid%id
	1225	CALL end_timing ( TRIM(message) )
	1226	END IF
	1227	CALL close_dataset ( rid , config_flags , "DATASET=RESTART" )
	1228	RETURN
	1229
	1230	END SUBROUTINE med_analysis_in
	1231	!=========================================================================================================
	1232	!=========================================================================================================
	1233	SUBROUTINE med_analysis_out ( grid , config_flags )
	1234	! Driver layer
	1235	USE module_domain , ONLY : domain, domain_clock_get
	1236	USE module_io_domain
	1237	USE module_timing
	1238	! Model layer
	1239	USE module_configure , ONLY : grid_config_rec_type
	1240	USE module_bc_time_utilities
	1241	!zhang USE WRF_ESMF_MOD
	1242
	1243	IMPLICIT NONE
	1244
	1245	! Arguments
	1246	TYPE(domain) :: grid
	1247	TYPE (grid_config_rec_type) , INTENT(IN) :: config_flags
	1248
	1249	! Local
	1250	LOGICAL, EXTERNAL :: wrf_dm_on_monitor
	1251	CHARACTER*80 :: rstname , outname
	1252	INTEGER :: fid , rid
	1253	CHARACTER (LEN=256) :: message
	1254	INTEGER :: ierr
	1255	INTEGER :: myproc
	1256	!zhang TYPE(ESMF_Time) :: CurrTime
	1257	TYPE(WRFU_Time) :: CurrTime
	1258	CHARACTER*80 :: timestr
	1259
	1260	IF ( wrf_dm_on_monitor() ) THEN
	1261	CALL start_timing
	1262	END IF
	1263
	1264	rid=grid%id
	1265
	1266	!zhang's doing CALL ESMF_ClockGet( grid%domain_clock, CurrTime=CurrTime, rc=ierr )
	1267	!zhang's doing CALL wrf_timetoa ( CurrTime, timestr )
	1268	CALL domain_clock_get( grid, current_timestr=timestr )
	1269	CALL construct_filename2a ( rstname ,config_flags%anl_outname, grid%id , 2 , timestr )
	1270
	1271	WRITE( message , '("med_analysis_out: opening ",A," for writing")' ) TRIM ( rstname )
	1272	CALL wrf_debug( 1 , message )
	1273	CALL open_w_dataset ( rid, TRIM(rstname), grid , &
	1274	config_flags , output_restart , "DATASET=RESTART", ierr )
	1275
	1276	IF ( ierr .NE. 0 ) THEN
	1277	CALL WRF_message( message )
	1278	ENDIF
	1279	CALL output_restart ( rid, grid , config_flags , ierr )
	1280	IF ( wrf_dm_on_monitor() ) THEN
	1281	WRITE ( message , FMT = '("Writing restart for domain ",I8)' ) grid%id
	1282	CALL end_timing ( TRIM(message) )
	1283	END IF
	1284	CALL close_dataset ( rid , config_flags , "DATASET=RESTART" )
	1285	RETURN
	1286	END SUBROUTINE med_analysis_out
	1287
	1288	#endif
	1289
	1290	RECURSIVE SUBROUTINE med_restart_out ( grid , config_flags )
	1291	! Driver layer
	1292	USE module_domain , ONLY : domain , domain_clock_get
	1293	USE module_io_domain
	1294	USE module_timing
	1295	USE module_configure , ONLY : grid_config_rec_type
	1296	! Model layer
	1297	USE module_bc_time_utilities
	1298	USE module_utility
	1299
	1300	IMPLICIT NONE
	1301
	1302	! Arguments
	1303	TYPE(domain) :: grid
	1304	TYPE (grid_config_rec_type) , INTENT(IN) :: config_flags
	1305
	1306	! Local
	1307	LOGICAL, EXTERNAL :: wrf_dm_on_monitor
	1308	CHARACTER*80 :: rstname , outname
	1309	INTEGER :: fid , rid, kid
	1310	CHARACTER (LEN=256) :: message
	1311	INTEGER :: ierr
	1312	INTEGER :: myproc
	1313	CHARACTER*80 :: timestr
	1314	TYPE (grid_config_rec_type) :: kid_config_flags
	1315
	1316	IF ( wrf_dm_on_monitor() ) THEN
	1317	CALL start_timing
	1318	END IF
	1319
	1320	! take this out - no effect - LPC
	1321	! rid=grid%id !zhang's doing
	1322
	1323	! write out this domains restart file first
	1324
	1325	CALL domain_clock_get( grid, current_timestr=timestr )
	1326	CALL construct_filename2a ( rstname , config_flags%rst_outname , grid%id , 2 , timestr )
	1327
	1328	WRITE( message , '("med_restart_out: opening ",A," for writing")' ) TRIM ( rstname )
	1329	CALL wrf_debug( 1 , message )
	1330	CALL open_w_dataset ( rid, TRIM(rstname), grid , &
	1331	config_flags , output_restart , "DATASET=RESTART", ierr )
	1332
	1333	IF ( ierr .NE. 0 ) THEN
	1334	CALL WRF_message( message )
	1335	ENDIF
	1336	CALL output_restart ( rid, grid , config_flags , ierr )
	1337	IF ( wrf_dm_on_monitor() ) THEN
	1338	WRITE ( message , FMT = '("Writing restart for domain ",I8)' ) grid%id
	1339	CALL end_timing ( TRIM(message) )
	1340	END IF
	1341	CALL close_dataset ( rid , config_flags , "DATASET=RESTART" )
	1342
	1343	! call recursively for children, (if any)
	1344	DO kid = 1, max_nests
	1345	IF ( ASSOCIATED( grid%nests(kid)%ptr ) ) THEN
	1346	CALL model_to_grid_config_rec ( grid%nests(kid)%ptr%id , model_config_rec , kid_config_flags )
	1347	CALL med_restart_out ( grid%nests(kid)%ptr , kid_config_flags )
	1348	ENDIF
	1349	ENDDO
	1350
	1351	RETURN
	1352	END SUBROUTINE med_restart_out
	1353
	1354	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
	1355
	1356	#ifdef HWRF
	1357	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
	1358	!zhang's doing
	1359	SUBROUTINE med_restart_in ( grid , config_flags )
	1360	! Driver layer
	1361	USE module_domain , ONLY : domain, domain_clock_get
	1362	USE module_io_domain
	1363	USE module_timing
	1364	! Model layer
	1365	USE module_configure , ONLY : grid_config_rec_type
	1366	USE module_bc_time_utilities
	1367
	1368	IMPLICIT NONE
	1369
	1370	! Arguments
	1371	TYPE(domain) :: grid
	1372	TYPE (grid_config_rec_type) , INTENT(IN) :: config_flags
	1373
	1374	! Local
	1375	LOGICAL, EXTERNAL :: wrf_dm_on_monitor
	1376	CHARACTER*80 :: rstname , outname
	1377	INTEGER :: fid , rid
	1378	CHARACTER (LEN=256) :: message
	1379	INTEGER :: ierr
	1380	INTEGER :: myproc
	1381	!zhang old TYPE(ESMF_Time) :: CurrTime
	1382	TYPE(WRFU_Time) :: CurrTime
	1383	CHARACTER*80 :: timestr
	1384
	1385	IF ( wrf_dm_on_monitor() ) THEN
	1386	CALL start_timing
	1387	END IF
	1388
	1389	rid=grid%id
	1390
	1391	!zhang's doing CALL ESMF_ClockGet( grid%domain_clock, CurrTime=CurrTime, rc=ierr )
	1392	!zhang's doing CALL wrf_timetoa ( CurrTime, timestr )
	1393	CALL domain_clock_get( grid, current_timestr=timestr )
	1394	CALL construct_filename2a ( rstname ,config_flags%rst_outname, grid%id , 2 , timestr )
	1395
	1396	WRITE( message , '("med_restart_in: opening ",A," for reading")' ) TRIM ( rstname )
	1397	CALL wrf_debug( 1 , message )
	1398	CALL open_r_dataset ( rid, TRIM(rstname), grid , &
	1399	config_flags , "DATASET=RESTART", ierr )
	1400
	1401	IF ( ierr .NE. 0 ) THEN
	1402	! CALL WRF_message( message )
	1403	CALL WRF_ERROR_FATAL('NESTED DOMAIN ERROR: FOR ANALYSIS SET TO TRUE, YOU NEED wrfanal FILE')
	1404	ENDIF
	1405	CALL input_restart ( rid, grid , config_flags , ierr )
	1406	IF ( wrf_dm_on_monitor() ) THEN
	1407	WRITE ( message , FMT = '("Reading restart for domain ",I8)' ) grid%id
	1408	CALL end_timing ( TRIM(message) )
	1409	END IF
	1410	CALL close_dataset ( rid , config_flags , "DATASET=RESTART" )
	1411	RETURN
	1412
	1413	END SUBROUTINE med_restart_in
	1414	!end of zhang's doing
	1415	#endif
	1416
	1417	SUBROUTINE med_hist_out ( grid , stream, config_flags )
	1418	! Driver layer
	1419	USE module_domain , ONLY : domain
	1420	USE module_timing
	1421	USE module_io_domain
	1422	USE module_configure , ONLY : grid_config_rec_type
	1423	USE module_bc_time_utilities
	1424	USE module_utility
	1425
	1426	IMPLICIT NONE
	1427	! Arguments
	1428	TYPE(domain) :: grid
	1429	TYPE (grid_config_rec_type) , INTENT(IN) :: config_flags
	1430	INTEGER , INTENT(IN) :: stream
	1431	! Local
	1432	LOGICAL, EXTERNAL :: wrf_dm_on_monitor
	1433	CHARACTER*80 :: fname, n2
	1434	CHARACTER (LEN=256) :: message
	1435	INTEGER :: ierr
	1436
	1437	IF ( wrf_dm_on_monitor() ) THEN
	1438	CALL start_timing
	1439	END IF
	1440
	1441	IF ( stream .LT. first_history .OR. stream .GT. last_auxhist ) THEN
	1442	WRITE(message,*)'med_hist_out: invalid history stream ',stream
	1443	CALL wrf_error_fatal( message )
	1444	ENDIF
	1445
	1446	SELECT CASE( stream )
	1447	CASE ( HISTORY_ALARM )
	1448	CALL open_hist_w( grid, config_flags, stream, HISTORY_ALARM, &
	1449	config_flags%history_outname, grid%oid, &
	1450	output_history, fname, n2, ierr )
	1451	CALL output_history ( grid%oid, grid , config_flags , ierr )
	1452
	1453	! registry-generated selections and calls top open_hist_w for aux streams
	1454	#include "med_hist_out_opens.inc"
	1455
	1456	END SELECT
	1457
	1458	WRITE(message,*)'med_hist_out: opened ',TRIM(fname),' as ',TRIM(n2)
	1459	CALL wrf_debug( 1, message )
	1460
	1461	grid%nframes(stream) = grid%nframes(stream) + 1
	1462
	1463	SELECT CASE( stream )
	1464	CASE ( HISTORY_ALARM )
	1465	IF ( grid%nframes(stream) >= config_flags%frames_per_outfile ) THEN
	1466	CALL close_dataset ( grid%oid , config_flags , n2 )
	1467	grid%oid = 0
	1468	grid%nframes(stream) = 0
	1469	ENDIF
	1470	! registry-generated selections and calls top close_dataset for aux streams
	1471	#include "med_hist_out_closes.inc"
	1472
	1473	END SELECT
	1474	IF ( wrf_dm_on_monitor() ) THEN
	1475	WRITE ( message , FMT = '("Writing ",A30," for domain ",I8)' )TRIM(fname),grid%id
	1476	CALL end_timing ( TRIM(message) )
	1477	END IF
	1478
	1479	RETURN
	1480	END SUBROUTINE med_hist_out
	1481
	1482	#if (DA_CORE != 1)
	1483	SUBROUTINE med_fddaobs_in ( grid , config_flags )
	1484	USE module_domain , ONLY : domain
	1485	USE module_configure , ONLY : grid_config_rec_type
	1486	IMPLICIT NONE
	1487	TYPE(domain) :: grid
	1488	TYPE (grid_config_rec_type) , INTENT(IN) :: config_flags
	1489	CALL wrf_fddaobs_in( grid, config_flags )
	1490	RETURN
	1491	END SUBROUTINE med_fddaobs_in
	1492	#endif
	1493
	1494	SUBROUTINE med_auxinput_in ( grid , stream, config_flags )
	1495	! Driver layer
	1496	USE module_domain , ONLY : domain
	1497	USE module_io_domain
	1498	! Model layer
	1499	USE module_configure , ONLY : grid_config_rec_type
	1500	USE module_bc_time_utilities
	1501	USE module_utility
	1502
	1503	IMPLICIT NONE
	1504	! Arguments
	1505	TYPE(domain) :: grid
	1506	TYPE (grid_config_rec_type) , INTENT(IN) :: config_flags
	1507	INTEGER , INTENT(IN) :: stream
	1508	! Local
	1509	CHARACTER (LEN=256) :: message
	1510	INTEGER :: ierr
	1511
	1512	IF ( stream .LT. first_auxinput .OR. stream .GT. last_auxinput ) THEN
	1513	WRITE(message,*)'med_auxinput_in: invalid input stream ',stream
	1514	CALL wrf_error_fatal( message )
	1515	ENDIF
	1516
	1517	grid%nframes(stream) = grid%nframes(stream) + 1
	1518
	1519	SELECT CASE( stream )
	1520	! registry-generated file of calls to open filename
	1521	! CASE ( AUXINPUT1_ALARM )
	1522	! CALL open_aux_u( grid, config_flags, stream, AUXINPUT1_ALARM, &
	1523	! config_flags%auxinput1_inname, grid%auxinput1_oid, &
	1524	! input_auxinput1, ierr )
	1525	! CALL input_auxinput1 ( grid%auxinput1_oid, grid , config_flags , ierr )
	1526	#include "med_auxinput_in.inc"
	1527	END SELECT
	1528
	1529	SELECT CASE( stream )
	1530	! registry-generated selections and calls top close_dataset for aux streams
	1531	#include "med_auxinput_in_closes.inc"
	1532	END SELECT
	1533
	1534	RETURN
	1535	END SUBROUTINE med_auxinput_in
	1536
	1537	SUBROUTINE med_filter_out ( grid , config_flags )
	1538	! Driver layer
	1539	USE module_domain , ONLY : domain , domain_clock_get
	1540	USE module_io_domain
	1541	USE module_timing
	1542	USE module_configure , ONLY : grid_config_rec_type
	1543	! Model layer
	1544	USE module_bc_time_utilities
	1545
	1546	IMPLICIT NONE
	1547
	1548	! Arguments
	1549	TYPE(domain) :: grid
	1550	TYPE (grid_config_rec_type) , INTENT(IN) :: config_flags
	1551
	1552	LOGICAL, EXTERNAL :: wrf_dm_on_monitor
	1553	CHARACTER*80 :: rstname , outname
	1554	INTEGER :: fid , rid
	1555	CHARACTER (LEN=256) :: message
	1556	INTEGER :: ierr
	1557	INTEGER :: myproc
	1558	CHARACTER*80 :: timestr
	1559
	1560	IF ( config_flags%write_input ) THEN
	1561
	1562	IF ( wrf_dm_on_monitor() ) THEN
	1563	CALL start_timing
	1564	END IF
	1565
	1566	CALL domain_clock_get( grid, current_timestr=timestr )
	1567	CALL construct_filename2a ( outname , config_flags%input_outname , grid%id , 2 , timestr )
	1568
	1569	WRITE ( message , '("med_filter_out 1: opening ",A," for writing. ")') TRIM ( outname )
	1570	CALL wrf_debug( 1, message )
	1571
	1572	CALL open_w_dataset ( fid, TRIM(outname), grid , &
	1573	config_flags , output_input , "DATASET=INPUT", ierr )
	1574	IF ( ierr .NE. 0 ) THEN
	1575	CALL wrf_error_fatal( message )
	1576	ENDIF
	1577
	1578	IF ( ierr .NE. 0 ) THEN
	1579	CALL wrf_error_fatal( message )
	1580	ENDIF
	1581
	1582	CALL output_input ( fid, grid , config_flags , ierr )
	1583	CALL close_dataset ( fid , config_flags , "DATASET=INPUT" )
	1584
	1585	IF ( wrf_dm_on_monitor() ) THEN
	1586	WRITE ( message , FMT = '("Writing filter output for domain ",I8)' ) grid%id
	1587	CALL end_timing ( TRIM(message) )
	1588	END IF
	1589	ENDIF
	1590
	1591	RETURN
	1592	END SUBROUTINE med_filter_out
	1593
	1594	SUBROUTINE med_latbound_in ( grid , config_flags )
	1595	! Driver layer
	1596	USE module_domain , ONLY : domain , domain_clock_get, head_grid
	1597	USE module_io_domain
	1598	USE module_timing
	1599	USE module_configure , ONLY : grid_config_rec_type
	1600	! Model layer
	1601	USE module_bc_time_utilities
	1602	USE module_utility
	1603
	1604	IMPLICIT NONE
	1605
	1606	#include <wrf_status_codes.h>
	1607
	1608	! Arguments
	1609	TYPE(domain) :: grid
	1610	TYPE (grid_config_rec_type) , INTENT(IN) :: config_flags
	1611
	1612	! Local data
	1613	LOGICAL, EXTERNAL :: wrf_dm_on_monitor
	1614	LOGICAL :: lbc_opened
	1615	INTEGER :: idum1 , idum2 , ierr , open_status , fid, rc
	1616	REAL :: bfrq
	1617	CHARACTER (LEN=256) :: message
	1618	CHARACTER (LEN=80) :: bdyname
	1619	Type (WRFU_Time ) :: startTime, stopTime, currentTime
	1620	Type (WRFU_TimeInterval ) :: stepTime
	1621	integer myproc,i,j,k
	1622
	1623	#include <wrf_io_flags.h>
	1624
	1625	CALL wrf_debug ( 200 , 'in med_latbound_in' )
	1626
	1627	! #if (EM_CORE == 1)
	1628	! Avoid trying to re-read the boundary conditions if we are doing DFI integration
	1629	! and do not expect to find boundary conditions for the current time
	1630	IF ( (grid%dfi_opt .EQ. DFI_DDFI .OR. grid%dfi_opt .EQ. DFI_TDFI) .AND. grid%dfi_stage .EQ. DFI_FWD ) RETURN
	1631	! #endif
	1632
	1633	IF ( grid%id .EQ. 1 .AND. config_flags%specified .AND. config_flags%io_form_boundary .GT. 0 ) THEN
	1634
	1635	CALL domain_clock_get( grid, current_time=currentTime, &
	1636	start_time=startTime, &
	1637	stop_time=stopTime, &
	1638	time_step=stepTime )
	1639
	1640	IF ( ( lbc_read_time( currentTime ) ) .AND. &
	1641	( currentTime + stepTime .GE. stopTime ) .AND. &
	1642	( currentTime .NE. startTime ) ) THEN
	1643	CALL wrf_debug( 100 , 'med_latbound_in: Skipping attempt to read lateral boundary file during last time step ' )
	1644
	1645	ELSE IF ( WRFU_AlarmIsRinging( grid%alarms( BOUNDARY_ALARM ), rc=rc ) ) THEN
	1646	CALL wrf_debug ( 100 , 'in med_latbound_in preparing to read' )
	1647	CALL WRFU_AlarmRingerOff( grid%alarms( BOUNDARY_ALARM ), rc=rc )
	1648	IF ( wrf_dm_on_monitor() ) CALL start_timing
	1649
	1650	! typically a <date> wouldn't be part of the bdy_inname, so just pass a dummy
	1651	CALL construct_filename2a ( bdyname , config_flags%bdy_inname , grid%id , 2 , 'dummydate' )
	1652
	1653	CALL wrf_inquire_opened(grid%lbc_fid , TRIM(bdyname) , open_status , ierr )
	1654	IF ( open_status .EQ. WRF_FILE_OPENED_FOR_READ ) THEN
	1655	lbc_opened = .TRUE.
	1656	ELSE
	1657	lbc_opened = .FALSE.
	1658	ENDIF
	1659	CALL wrf_dm_bcast_bytes ( lbc_opened , LWORDSIZE )
	1660	IF ( .NOT. lbc_opened ) THEN
	1661	CALL construct_filename1 ( bdyname , 'wrfbdy' , grid%id , 2 )
	1662	WRITE(message,*)'Opening: ',TRIM(bdyname)
	1663	CALL wrf_debug(100,TRIM(message))
	1664	CALL open_r_dataset ( grid%lbc_fid, TRIM(bdyname) , grid , config_flags , "DATASET=BOUNDARY", ierr )
	1665	IF ( ierr .NE. 0 ) THEN
	1666	WRITE( message, * ) 'med_latbound_in: error opening ',TRIM(bdyname), ' for reading. IERR = ',ierr
	1667	CALL WRF_ERROR_FATAL( message )
	1668	ENDIF
	1669	ELSE
	1670	CALL wrf_debug( 100 , bdyname // 'already opened' )
	1671	ENDIF
	1672	CALL wrf_debug( 100 , 'med_latbound_in: calling input_boundary ' )
	1673	CALL input_boundary ( grid%lbc_fid, grid , config_flags , ierr )
	1674
	1675	! #if (EM_CORE == 1)
	1676	IF ( (config_flags%dfi_opt .NE. DFI_NODFI) .AND. (head_grid%dfi_stage .NE. DFI_FST) ) THEN
	1677	CALL wrf_debug( 100 , 'med_latbound_in: closing boundary file ' )
	1678	CALL close_dataset ( grid%lbc_fid , config_flags , "DATASET=BOUNDARY" )
	1679	END IF
	1680	! #endif
	1681
	1682	CALL domain_clock_get( grid, current_time=currentTime )
	1683	DO WHILE (currentTime .GE. grid%next_bdy_time ) ! next_bdy_time is set by input_boundary from bdy file
	1684	CALL wrf_debug( 100 , 'med_latbound_in: calling input_boundary ' )
	1685	CALL input_boundary ( grid%lbc_fid, grid , config_flags , ierr )
	1686	ENDDO
	1687	CALL WRFU_AlarmSet( grid%alarms( BOUNDARY_ALARM ), RingTime=grid%next_bdy_time, rc=rc )
	1688
	1689	IF ( ierr .NE. 0 .and. ierr .NE. WRF_WARN_NETCDF ) THEN
	1690	WRITE( message, * ) 'med_latbound_in: error reading ',TRIM(bdyname), ' IERR = ',ierr
	1691	CALL WRF_ERROR_FATAL( message )
	1692	ENDIF
	1693	IF ( currentTime .EQ. grid%this_bdy_time ) grid%dtbc = 0.
	1694
	1695	IF ( wrf_dm_on_monitor() ) THEN
	1696	WRITE ( message , FMT = '("processing lateral boundary for domain ",I8)' ) grid%id
	1697	CALL end_timing ( TRIM(message) )
	1698	ENDIF
	1699	ENDIF
	1700	ENDIF
	1701	RETURN
	1702	END SUBROUTINE med_latbound_in
	1703
	1704	SUBROUTINE med_setup_step ( grid , config_flags )
	1705	! Driver layer
	1706	USE module_domain , ONLY : domain
	1707	USE module_configure , ONLY : grid_config_rec_type
	1708	! Model layer
	1709
	1710	IMPLICIT NONE
	1711	!<DESCRIPTION>
	1712	!
	1713	!The driver layer routine integrate() calls this mediation layer routine
	1714	!prior to initiating a time step on the domain specified by the argument
	1715	!grid. This provides the model-layer contributor an opportunity to make
	1716	!any pre-time-step initializations that pertain to a particular model
	1717	!domain. In WRF, this routine is used to call
	1718	!set_scalar_indices_from_config for the specified domain.
	1719	!
	1720	!</DESCRIPTION>
	1721
	1722	! Arguments
	1723	TYPE(domain) :: grid
	1724	TYPE (grid_config_rec_type) , INTENT(IN) :: config_flags
	1725	! Local
	1726	INTEGER :: idum1 , idum2
	1727
	1728	CALL set_scalar_indices_from_config ( grid%id , idum1 , idum2 )
	1729
	1730	RETURN
	1731
	1732	END SUBROUTINE med_setup_step
	1733
	1734	SUBROUTINE med_endup_step ( grid , config_flags )
	1735	! Driver layer
	1736	USE module_domain , ONLY : domain
	1737	USE module_configure , ONLY : grid_config_rec_type, model_config_rec
	1738	! Model layer
	1739
	1740	IMPLICIT NONE
	1741	!<DESCRIPTION>
	1742	!
	1743	!The driver layer routine integrate() calls this mediation layer routine
	1744	!prior to initiating a time step on the domain specified by the argument
	1745	!grid. This provides the model-layer contributor an opportunity to make
	1746	!any pre-time-step initializations that pertain to a particular model
	1747	!domain. In WRF, this routine is used to call
	1748	!set_scalar_indices_from_config for the specified domain.
	1749	!
	1750	!</DESCRIPTION>
	1751
	1752	! Arguments
	1753	TYPE(domain) :: grid
	1754	TYPE (grid_config_rec_type) , INTENT(OUT) :: config_flags
	1755	! Local
	1756	INTEGER :: idum1 , idum2
	1757
	1758	IF ( grid%id .EQ. 1 ) THEN
	1759	! turn off the restart flag after the first mother-domain step is finished
	1760	model_config_rec%restart = .FALSE.
	1761	config_flags%restart = .FALSE.
	1762	CALL nl_set_restart(1, .FALSE.)
	1763
	1764	ENDIF
	1765
	1766	RETURN
	1767
	1768	END SUBROUTINE med_endup_step
	1769
	1770	SUBROUTINE open_aux_u ( grid , config_flags, stream, alarm_id, &
	1771	auxinput_inname, oid, insub, ierr )
	1772	! Driver layer
	1773	USE module_domain , ONLY : domain , domain_clock_get
	1774	USE module_io_domain
	1775	! Model layer
	1776	USE module_configure , ONLY : grid_config_rec_type
	1777	USE module_bc_time_utilities
	1778	USE module_utility
	1779
	1780	IMPLICIT NONE
	1781	! Arguments
	1782	TYPE(domain) :: grid
	1783	TYPE (grid_config_rec_type) , INTENT(IN) :: config_flags
	1784	INTEGER , INTENT(IN) :: stream
	1785	INTEGER , INTENT(IN) :: alarm_id
	1786	CHARACTER() , INTENT(IN) :: auxinput_inname
	1787	INTEGER , INTENT(INOUT) :: oid
	1788	EXTERNAL insub
	1789	INTEGER , INTENT(OUT) :: ierr
	1790	! Local
	1791	CHARACTER*80 :: fname, n2
	1792	CHARACTER (LEN=256) :: message
	1793	CHARACTER*80 :: timestr
	1794	TYPE(WRFU_Time) :: ST,CT
	1795	LOGICAL :: adjust
	1796
	1797	IF ( stream .LT. first_stream .OR. stream .GT. last_stream ) THEN
	1798	WRITE(message,*)'open_aux_u: invalid input stream ',stream
	1799	CALL wrf_error_fatal( message )
	1800	ENDIF
	1801
	1802	ierr = 0
	1803
	1804	IF ( oid .eq. 0 ) THEN
	1805	CALL domain_clock_get( grid, current_time=CT, start_time=ST, &
	1806	current_timestr=timestr )
	1807	CALL nl_get_adjust_input_times( grid%id, adjust )
	1808	IF ( adjust ) THEN
	1809	CALL adjust_io_timestr( grid%io_intervals( alarm_id ), CT, ST, timestr )
	1810	ENDIF
	1811	CALL construct_filename2a ( fname , auxinput_inname, &
	1812	grid%id , 2 , timestr )
	1813	IF ( stream .EQ. 10 ) THEN
	1814	WRITE(n2,'("DATASET=AUXINPUT10")')
	1815	ELSE IF ( stream .EQ. 11 ) THEN
	1816	WRITE(n2,'("DATASET=AUXINPUT11")')
	1817	ELSE IF ( stream .GE. 10 ) THEN
	1818	WRITE(n2,'("DATASET=AUXINPUT",I2)')stream-first_input
	1819	ELSE
	1820	WRITE(n2,'("DATASET=AUXINPUT",I1)')stream-first_input
	1821	ENDIF
	1822	WRITE ( message , '("open_aux_u : opening ",A," for reading. ")') TRIM ( fname )
	1823	CALL wrf_debug( 1, message )
	1824	!<DESCRIPTION>
	1825	!
	1826	!Open_u_dataset is called rather than open_r_dataset to allow interfaces
	1827	!that can do blending or masking to update an existing field. (MCEL IO does this).
	1828	!No effect for other interfaces; open_u_dataset is equivalent to open_r_dataset
	1829	!in those cases.
	1830	!
	1831	!</DESCRIPTION>
	1832	CALL open_u_dataset ( oid, TRIM(fname), grid , &
	1833	config_flags , insub , n2, ierr )
	1834	ENDIF
	1835	IF ( ierr .NE. 0 ) THEN
	1836	WRITE ( message , '("open_aux_u : error opening ",A," for reading. ",I3)') &
	1837	TRIM ( fname ), ierr
	1838	CALL wrf_message( message )
	1839	ENDIF
	1840	RETURN
	1841	END SUBROUTINE open_aux_u
	1842
	1843	SUBROUTINE open_hist_w ( grid , config_flags, stream, alarm_id, &
	1844	hist_outname, oid, outsub, fname, n2, ierr )
	1845	! Driver layer
	1846	USE module_domain , ONLY : domain , domain_clock_get
	1847	USE module_io_domain
	1848	! Model layer
	1849	USE module_configure , ONLY : grid_config_rec_type
	1850	USE module_bc_time_utilities
	1851	USE module_utility
	1852
	1853	IMPLICIT NONE
	1854	! Arguments
	1855	TYPE(domain) :: grid
	1856	TYPE (grid_config_rec_type) , INTENT(IN) :: config_flags
	1857	INTEGER , INTENT(IN) :: stream
	1858	INTEGER , INTENT(IN) :: alarm_id
	1859	CHARACTER() , INTENT(IN) :: hist_outname
	1860	INTEGER , INTENT(INOUT) :: oid
	1861	EXTERNAL outsub
	1862	CHARACTER() , INTENT(OUT) :: fname, n2
	1863	INTEGER , INTENT(OUT) :: ierr
	1864	! Local
	1865	INTEGER :: len_n2
	1866	CHARACTER (LEN=256) :: message
	1867	CHARACTER*80 :: timestr
	1868	TYPE(WRFU_Time) :: ST,CT
	1869	LOGICAL :: adjust
	1870
	1871	IF ( stream .LT. first_history .OR. stream .GT. last_history ) THEN
	1872	WRITE(message,*)'open_hist_w: invalid history stream ',stream
	1873	CALL wrf_error_fatal( message )
	1874	ENDIF
	1875
	1876	ierr = 0
	1877
	1878	! Note that computation of fname and n2 are outside of the oid IF statement
	1879	! since they are OUT args and may be used by callers even if oid/=0.
	1880	CALL domain_clock_get( grid, current_time=CT, start_time=ST, &
	1881	current_timestr=timestr )
	1882	CALL nl_get_adjust_output_times( grid%id, adjust )
	1883	IF ( adjust ) THEN
	1884	CALL adjust_io_timestr( grid%io_intervals( alarm_id ), CT, ST, timestr )
	1885	ENDIF
	1886	CALL construct_filename2a ( fname , hist_outname, &
	1887	grid%id , 2 , timestr )
	1888	IF ( stream .EQ. history_only ) THEN
	1889	WRITE(n2,'("DATASET=HISTORY")')
	1890	ELSE IF ( stream .GE. 10 ) THEN
	1891	WRITE(n2,'("DATASET=AUXHIST",I2)')stream-first_history
	1892	ELSE
	1893	WRITE(n2,'("DATASET=AUXHIST",I1)')stream-first_history
	1894	ENDIF
	1895	#if (DA_CORE == 1)
	1896	len_n2 = LEN_TRIM(n2)
	1897	WRITE(n2(len_n2+1:len_n2+19),'(",REAL_OUTPUT_SIZE=4")')
	1898	#endif
	1899	IF ( oid .eq. 0 ) THEN
	1900	WRITE ( message , '("open_hist_w : opening ",A," for writing. ")') TRIM ( fname )
	1901	CALL wrf_debug( 1, message )
	1902	!<DESCRIPTION>
	1903	!
	1904	!Open_u_dataset is called rather than open_r_dataset to allow interfaces
	1905	!that can do blending or masking to update an existing field. (MCEL IO does this).
	1906	!No effect for other interfaces; open_u_dataset is equivalent to open_r_dataset
	1907	!in those cases.
	1908	!
	1909	!</DESCRIPTION>
	1910	CALL open_w_dataset ( oid, TRIM(fname), grid , &
	1911	config_flags , outsub , n2, ierr )
	1912	ENDIF
	1913	IF ( ierr .NE. 0 ) THEN
	1914	WRITE ( message , '("open_hist_w : error opening ",A," for writing. ",I3)') &
	1915	TRIM ( fname ), ierr
	1916	CALL wrf_message( message )
	1917	ENDIF
	1918	RETURN
	1919	END SUBROUTINE open_hist_w
	1920
	1921
	1922	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
	1923
	1924	#ifdef WRF_CHEM
	1925
	1926	SUBROUTINE med_read_wrf_chem_input ( grid , config_flags )
	1927	! Driver layer
	1928	USE module_domain , ONLY : domain , domain_clock_get
	1929	USE module_io_domain
	1930	USE module_timing
	1931	USE module_configure , ONLY : grid_config_rec_type
	1932	! Model layer
	1933	USE module_bc_time_utilities
	1934	#ifdef DM_PARALLEL
	1935	USE module_dm
	1936	#endif
	1937	USE module_date_time
	1938	USE module_utility
	1939
	1940	IMPLICIT NONE
	1941
	1942	! Arguments
	1943	TYPE(domain) :: grid
	1944
	1945	TYPE (grid_config_rec_type) , INTENT(IN) :: config_flags
	1946
	1947	! Local data
	1948	LOGICAL, EXTERNAL :: wrf_dm_on_monitor
	1949
	1950	INTEGER :: ierr, efid
	1951	REAL :: time, tupdate
	1952	real, allocatable :: dumc0(:,:,:)
	1953	CHARACTER (LEN=256) :: message, current_date_char, date_string
	1954	CHARACTER (LEN=80) :: inpname
	1955
	1956	#include <wrf_io_flags.h>
	1957	! IF ( grid%id .EQ. 1 ) THEN
	1958
	1959	CALL domain_clock_get( grid, current_timestr=current_date_char )
	1960
	1961	CALL construct_filename1 ( inpname , config_flags%auxinput12_inname , grid%id , 2 )
	1962	WRITE(message,*)'mediation_integrate: med_read_wrf_chem_input: Open file ',TRIM(inpname)
	1963	CALL wrf_message( TRIM(message) )
	1964
	1965	if( grid%auxinput12_oid .NE. 0 ) then
	1966	CALL close_dataset ( grid%auxinput12_oid , config_flags , "DATASET=AUXINPUT12" )
	1967	endif
	1968
	1969	CALL open_r_dataset ( grid%auxinput12_oid, TRIM(inpname) , grid , config_flags, &
	1970	"DATASET=AUXINPUT12", ierr )
	1971	IF ( ierr .NE. 0 ) THEN
	1972	WRITE( message , * ) 'med_read_wrf_chem_input error opening ', TRIM( inpname )
	1973	CALL wrf_error_fatal( TRIM( message ) )
	1974	ENDIF
	1975
	1976	WRITE(message,*)'mediation_integrate: med_read_wrf_chem_input: Read chemistry from wrfout at time ',&
	1977	TRIM(current_date_char)
	1978	CALL wrf_message( TRIM(message) )
	1979
	1980	CALL wrf_debug (100 , 'mediation_integrate: calling input_auxinput12' )
	1981	CALL input_auxinput12 ( grid%auxinput12_oid, grid , config_flags , ierr )
	1982
	1983	CALL close_dataset ( grid%auxinput12_oid , config_flags , "DATASET=AUXINPUT12" )
	1984
	1985	! ENDIF
	1986	CALL wrf_debug (100 , 'mediation_integrate: med_read_wrf_chem_input: exit' )
	1987
	1988	END SUBROUTINE med_read_wrf_chem_input
	1989	!------------------------------------------------------------------------
	1990	! Chemistry emissions input control. Three options are available and are
	1991	! set via the namelist variable io_style_emissions:
	1992	!
	1993	! 0 = Emissions are not read in from a file. They will contain their
	1994	! default values, which can be set in the Registry.
	1995	! (Intended for debugging of chem code)
	1996	!
	1997	! 1 = Emissions are read in from two 12 hour files that are cycled.
	1998	! With this choice, auxinput5_inname should be set to
	1999	! the value "wrfchemi_hhZ_d<domain>".
	2000	!
	2001	! 2 = Emissions are read in from files identified by date and that have
	2002	! a length defined by frames_per_auxinput5. Both
	2003	! auxinput5_inname should be set to
	2004	! "wrfchemi_d<domain>_<date>".
	2005	!------------------------------------------------------------------------
	2006	SUBROUTINE med_read_wrf_chem_emiss ( grid , config_flags )
	2007	! Driver layer
	2008	USE module_domain , ONLY : domain , domain_clock_get
	2009	USE module_io_domain
	2010	USE module_timing
	2011	USE module_configure , ONLY : grid_config_rec_type
	2012	! Model layer
	2013	USE module_bc_time_utilities
	2014	#ifdef DM_PARALLEL
	2015	USE module_dm
	2016	#endif
	2017	USE module_date_time
	2018	USE module_utility
	2019
	2020	IMPLICIT NONE
	2021
	2022	! Arguments
	2023	TYPE(domain) :: grid
	2024
	2025	! TYPE (grid_config_rec_type) , INTENT(IN) :: config_flags
	2026	TYPE (grid_config_rec_type) :: config_flags
	2027	Type (WRFU_Time ) :: stopTime, currentTime
	2028	Type (WRFU_TimeInterval ) :: stepTime
	2029
	2030	! Local data
	2031	LOGICAL, EXTERNAL :: wrf_dm_on_monitor
	2032
	2033	INTEGER :: ierr, efid
	2034	INTEGER :: ihr, ihrdiff, i
	2035	REAL :: time, tupdate
	2036	real, allocatable :: dumc0(:,:,:)
	2037	CHARACTER (LEN=256) :: message, current_date_char, date_string
	2038	CHARACTER (LEN=80) :: inpname
	2039
	2040	#include <wrf_io_flags.h>
	2041
	2042	CALL model_to_grid_config_rec ( grid%id , model_config_rec , config_flags )
	2043
	2044	! This "if" should be commented out when using emission files for nested
	2045	! domains. Also comment out the "ENDIF" line noted below.
	2046	! IF ( grid%id .EQ. 1 ) THEN
	2047
	2048	CALL domain_clock_get( grid, current_time=currentTime, &
	2049	current_timestr=current_date_char, &
	2050	stop_time=stopTime, &
	2051	time_step=stepTime )
	2052
	2053	time = float(grid%itimestep) * grid%dt
	2054
	2055	!---
	2056	! io_style_emissions option 0: no emissions read in...
	2057	!---
	2058	if( config_flags%io_style_emissions == 0 ) then
	2059	! Do nothing.
	2060	!---
	2061	! io_style_emissions option 1: cycle through two 12 hour input files...
	2062	!---
	2063	else if( config_flags%io_style_emissions == 1 ) then
	2064
	2065	tupdate = mod( time, (12. * 3600.) )
	2066	read(current_date_char(12:13),'(I2)') ihr
	2067	ihr = MOD(ihr,24)
	2068	ihrdiff = 0
	2069
	2070	IF( tupdate .LT. grid%dt ) THEN
	2071	tupdate = 0.
	2072	ENDIF
	2073	IF( ihr .EQ. 00 .OR. ihr .EQ. 12 ) THEN
	2074	tupdate = 0.
	2075	ENDIF
	2076
	2077	IF( currentTime + stepTime .GE. stopTime .AND. &
	2078	grid%auxinput5_oid .NE. 0 ) THEN
	2079	CALL close_dataset ( grid%auxinput5_oid , config_flags , "DATASET=AUXINPUT5" )
	2080	tupdate = 1.
	2081	ENDIF
	2082
	2083	! write(message,FMT='(A,F10.1,A)') ' EMISSIONS UPDATE TIME ',time,TRIM(current_date_char(12:13))
	2084	! CALL wrf_message( TRIM(message) )
	2085
	2086	IF ( tupdate .EQ. 0. .AND. ihr .LT. 12 ) THEN
	2087	ihrdiff = ihr
	2088	CALL construct_filename1 ( inpname , 'wrfchemi_00z' , grid%id , 2 )
	2089	WRITE(message,*)'mediation_integrate: med_read_wrf_chem_emissions: Open file ',TRIM(inpname)
	2090	CALL wrf_message( TRIM(message) )
	2091
	2092	if( grid%auxinput5_oid .NE. 0 ) then
	2093	CALL close_dataset ( grid%auxinput5_oid , config_flags , "DATASET=AUXINPUT5" )
	2094	endif
	2095
	2096	CALL open_r_dataset ( grid%auxinput5_oid, TRIM(inpname) , grid , config_flags, &
	2097	"DATASET=AUXINPUT5", ierr )
	2098	IF ( ierr .NE. 0 ) THEN
	2099	WRITE( message , * ) 'med_read_wrf_chem_emissions: error opening ', TRIM( inpname )
	2100	CALL wrf_error_fatal( TRIM( message ) )
	2101	ENDIF
	2102
	2103	ELSE IF ( tupdate .EQ. 0. .AND. ihr .GE. 12 ) THEN
	2104	ihrdiff = ihr - 12
	2105
	2106	CALL construct_filename1 ( inpname , 'wrfchemi_12z' , grid%id , 2 )
	2107	WRITE(message,*)'mediation_integrate: med_read_wrf_chem_emissions: Open file ',TRIM(inpname)
	2108	CALL wrf_message( TRIM(message) )
	2109
	2110	if( grid%auxinput5_oid .NE. 0 ) then
	2111	CALL close_dataset ( grid%auxinput5_oid , config_flags , "DATASET=AUXINPUT5" )
	2112	endif
	2113
	2114	CALL open_r_dataset ( grid%auxinput5_oid, TRIM(inpname) , grid , config_flags, &
	2115	"DATASET=AUXINPUT5", ierr )
	2116	IF ( ierr .NE. 0 ) THEN
	2117	WRITE( message , * ) 'med_read_wrf_chem_emissions: error opening ', TRIM( inpname )
	2118	CALL wrf_error_fatal( TRIM( message ) )
	2119	ENDIF
	2120	ENDIF
	2121
	2122	WRITE( message, '(A,2F10.1)' ) ' HOURLY EMISSIONS UPDATE TIME ',time,mod(time,3600.)
	2123	CALL wrf_message( TRIM(message) )
	2124	!
	2125	! hourly updates to emissions
	2126	IF ( ( mod( time, 3600. ) .LT. grid%dt ) .AND. &
	2127	( currentTime + stepTime .LT. stopTime ) ) THEN
	2128	! IF ( wrf_dm_on_monitor() ) CALL start_timing
	2129
	2130	WRITE(message,'(A,A)')'mediation_integrate: med_read_wrf_chem_emissions: Read emissions for time ',TRIM(current_date_char)
	2131	CALL wrf_message( TRIM(message) )
	2132
	2133	IF ( tupdate .EQ. 0. .AND. ihrdiff .GT. 0) THEN
	2134	IF( ihrdiff .GT. 12) THEN
	2135	WRITE(message,'(A)')'mediation_integrate: med_read_wrf_chem_emissions: Error in emissions time, skipping all times in file '
	2136	CALL wrf_message( TRIM(message) )
	2137	ENDIF
	2138	DO i=1,ihrdiff
	2139	WRITE(message,'(A,I4)')'mediation_integrate: med_read_wrf_chem_emissions: Skip emissions ',i
	2140	CALL wrf_message( TRIM(message) )
	2141	CALL input_auxinput5 ( grid%auxinput5_oid, grid , config_flags , ierr )
	2142	ENDDO
	2143	ENDIF
	2144
	2145	CALL wrf_debug (100 , 'mediation_integrate: calling input_auxinput5' )
	2146	CALL input_auxinput5 ( grid%auxinput5_oid, grid , config_flags , ierr )
	2147	ELSE
	2148	CALL wrf_debug (100 , 'mediation_integrate: med_read_wrf_chem_emissions: Do not read emissions' )
	2149	ENDIF
	2150
	2151	!---
	2152	! io_style_emissions option 2: use dated emission files whose length is
	2153	! set via frames_per_auxinput5...
	2154	!---
	2155	else if( config_flags%io_style_emissions == 2 ) then
	2156	WRITE(message,*)'mediation_integrate: med_read_wrf_chem_emissions: Read emissions for time ',TRIM(current_date_char)
	2157	CALL wrf_message( TRIM(message) )
	2158	!
	2159	! Code to read hourly emission files...
	2160	!
	2161	if( grid%auxinput5_oid == 0 ) then
	2162	CALL construct_filename2a(inpname , grid%emi_inname, grid%id , 2, current_date_char)
	2163	WRITE(message,*)'mediation_integrate: med_read_wrf_chem_emissions: Open file ',TRIM(inpname)
	2164	CALL wrf_message( TRIM(message) )
	2165	CALL open_r_dataset ( grid%auxinput5_oid, TRIM(inpname) , grid , config_flags, &
	2166	"DATASET=AUXINPUT5", ierr )
	2167	IF ( ierr .NE. 0 ) THEN
	2168	WRITE( message , * ) 'med_read_wrf_chem_emissions: error opening ', TRIM( inpname )
	2169	CALL wrf_error_fatal( TRIM( message ) )
	2170	ENDIF
	2171	end if
	2172	!
	2173	! Read the emissions data.
	2174	!
	2175	CALL wrf_debug (100 , 'mediation_integrate: calling input_auxinput5' )
	2176	CALL input_auxinput5 ( grid%auxinput5_oid, grid , config_flags , ierr )
	2177	!
	2178	! If reached the indicated number of frames in the emissions file, close it.
	2179	!
	2180	grid%emissframes = grid%emissframes + 1
	2181	IF ( grid%emissframes >= config_flags%frames_per_auxinput5 ) THEN
	2182	CALL close_dataset ( grid%auxinput5_oid , config_flags , "DATASET=AUXINPUT5" )
	2183	grid%emissframes = 0
	2184	grid%auxinput5_oid = 0
	2185	ENDIF
	2186
	2187	!---
	2188	! unknown io_style_emissions option...
	2189	!---
	2190	else
	2191	call wrf_error_fatal("Unknown emission style selected via io_style_emissions.")
	2192	end if
	2193
	2194	! The following line should be commented out when using emission files
	2195	! for nested domains. Also comment out the "if" noted above.
	2196	! ENDIF
	2197
	2198	CALL wrf_debug (100 , 'mediation_integrate: med_read_wrf_chem_emissions: exit' )
	2199
	2200	END SUBROUTINE med_read_wrf_chem_emiss
	2201
	2202	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
	2203	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
	2204
	2205	SUBROUTINE med_read_wrf_chem_bioemiss ( grid , config_flags )
	2206	! Driver layer
	2207	USE module_domain , ONLY : domain , domain_clock_get
	2208	USE module_io_domain
	2209	USE module_timing
	2210	USE module_configure , ONLY : grid_config_rec_type
	2211	! Model layer
	2212	USE module_bc_time_utilities
	2213	#ifdef DM_PARALLEL
	2214	USE module_dm
	2215	#endif
	2216	USE module_date_time
	2217	USE module_utility
	2218
	2219	IMPLICIT NONE
	2220
	2221	! Arguments
	2222	TYPE(domain) :: grid
	2223
	2224	TYPE (grid_config_rec_type) , INTENT(IN) :: config_flags
	2225
	2226	! Local data
	2227	LOGICAL, EXTERNAL :: wrf_dm_on_monitor
	2228
	2229	INTEGER :: ierr, efid
	2230	REAL :: time, tupdate
	2231	real, allocatable :: dumc0(:,:,:)
	2232	CHARACTER (LEN=256) :: message, current_date_char, date_string
	2233	CHARACTER (LEN=80) :: inpname
	2234
	2235	#include <wrf_io_flags.h>
	2236	! IF ( grid%id .EQ. 1 ) THEN
	2237
	2238	CALL domain_clock_get( grid, current_timestr=current_date_char )
	2239
	2240	CALL construct_filename1 ( inpname , 'wrfbiochemi' , grid%id , 2 )
	2241	WRITE(message,*)'mediation_integrate: med_read_wrf_chem_bioemissions: Open file ',TRIM(inpname)
	2242	CALL wrf_message( TRIM(message) )
	2243
	2244	if( grid%auxinput6_oid .NE. 0 ) then
	2245	CALL close_dataset ( grid%auxinput6_oid , config_flags , "DATASET=AUXINPUT6" )
	2246	endif
	2247
	2248	CALL open_r_dataset ( grid%auxinput6_oid, TRIM(inpname) , grid , config_flags, &
	2249	"DATASET=AUXINPUT6", ierr )
	2250	IF ( ierr .NE. 0 ) THEN
	2251	WRITE( message , * ) 'med_read_wrf_chem_bioemissions: error opening ', TRIM( inpname )
	2252	CALL wrf_error_fatal( TRIM( message ) )
	2253	ENDIF
	2254
	2255	WRITE(message,*)'mediation_integrate: med_read_wrf_chem_bioemissions: Read biogenic emissions at time ',&
	2256	TRIM(current_date_char)
	2257	CALL wrf_message( TRIM(message) )
	2258
	2259	CALL wrf_debug (100 , 'mediation_integrate: calling input_auxinput6' )
	2260	CALL input_auxinput6 ( grid%auxinput6_oid, grid , config_flags , ierr )
	2261
	2262	CALL close_dataset ( grid%auxinput6_oid , config_flags , "DATASET=AUXINPUT6" )
	2263
	2264	! ENDIF
	2265	CALL wrf_debug (100 , 'mediation_integrate: med_read_wrf_chem_bioemissions: exit' )
	2266
	2267	END SUBROUTINE med_read_wrf_chem_bioemiss
	2268	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
	2269	SUBROUTINE med_read_wrf_chem_emissopt4 ( grid , config_flags )
	2270	! Driver layer
	2271	USE module_domain , ONLY : domain , domain_clock_get
	2272	USE module_io_domain
	2273	USE module_timing
	2274	USE module_configure , ONLY : grid_config_rec_type
	2275	! Model layer
	2276	USE module_bc_time_utilities
	2277	#ifdef DM_PARALLEL
	2278	USE module_dm
	2279	#endif
	2280	USE module_date_time
	2281	USE module_utility
	2282
	2283	IMPLICIT NONE
	2284
	2285	! Arguments
	2286	TYPE(domain) :: grid
	2287
	2288	TYPE (grid_config_rec_type) , INTENT(IN) :: config_flags
	2289
	2290	! Local data
	2291	LOGICAL, EXTERNAL :: wrf_dm_on_monitor
	2292
	2293	INTEGER :: ierr, efid
	2294	REAL :: time, tupdate
	2295	real, allocatable :: dumc0(:,:,:)
	2296	CHARACTER (LEN=256) :: message, current_date_char, date_string
	2297	CHARACTER (LEN=80) :: inpname
	2298
	2299	#include <wrf_io_flags.h>
	2300	! IF ( grid%id .EQ. 1 ) THEN
	2301
	2302	CALL domain_clock_get( grid, current_timestr=current_date_char )
	2303
	2304	CALL construct_filename1 ( inpname , 'wrfchemi' , grid%id , 2 )
	2305	WRITE(message,*)'mediation_integrate: med_read_wrf_chem_emissions: Open file ',TRIM(inpname)
	2306	CALL wrf_message( TRIM(message) )
	2307
	2308	if( grid%auxinput5_oid .NE. 0 ) then
	2309	CALL close_dataset ( grid%auxinput5_oid , config_flags , "DATASET=AUXINPUT5" )
	2310	endif
	2311
	2312	CALL open_r_dataset ( grid%auxinput5_oid, TRIM(inpname) , grid , config_flags, &
	2313	"DATASET=AUXINPUT5", ierr )
	2314	IF ( ierr .NE. 0 ) THEN
	2315	WRITE( message , * ) 'med_read_wrf_chem_emissions: error opening ', TRIM( inpname )
	2316	CALL wrf_error_fatal( TRIM( message ) )
	2317	ENDIF
	2318
	2319	WRITE(message,*)'mediation_integrate: med_read_wrf_chem_emissions: Read biogenic emissions at time ',&
	2320	TRIM(current_date_char)
	2321	CALL wrf_message( TRIM(message) )
	2322
	2323	CALL wrf_debug (100 , 'mediation_integrate: calling input_auxinput5' )
	2324	CALL input_auxinput5 ( grid%auxinput5_oid, grid , config_flags , ierr )
	2325
	2326	CALL close_dataset ( grid%auxinput5_oid , config_flags , "DATASET=AUXINPUT5" )
	2327
	2328	! ENDIF
	2329	CALL wrf_debug (100 , 'mediation_integrate: med_read_wrf_chem_emissions: exit' )
	2330
	2331	END SUBROUTINE med_read_wrf_chem_emissopt4
	2332
	2333	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
	2334	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
	2335
	2336	SUBROUTINE med_read_wrf_chem_dms_emiss ( grid , config_flags )
	2337	! Driver layer
	2338	USE module_domain , ONLY : domain , domain_clock_get
	2339	USE module_io_domain
	2340	USE module_timing
	2341	USE module_configure , ONLY : grid_config_rec_type
	2342	! Model layer
	2343	USE module_bc_time_utilities
	2344	#ifdef DM_PARALLEL
	2345	USE module_dm
	2346	#endif
	2347	USE module_date_time
	2348	USE module_utility
	2349
	2350	IMPLICIT NONE
	2351
	2352	! Arguments
	2353	TYPE(domain) :: grid
	2354
	2355	TYPE (grid_config_rec_type) , INTENT(IN) :: config_flags
	2356
	2357	! Local data
	2358	LOGICAL, EXTERNAL :: wrf_dm_on_monitor
	2359
	2360	INTEGER :: ierr, efid
	2361	REAL :: time, tupdate
	2362	real, allocatable :: dumc0(:,:,:)
	2363	CHARACTER (LEN=256) :: message, current_date_char, date_string
	2364	CHARACTER (LEN=80) :: inpname
	2365
	2366	#include <wrf_io_flags.h>
	2367	! IF ( grid%id .EQ. 1 ) THEN
	2368
	2369	CALL domain_clock_get( grid, current_timestr=current_date_char )
	2370
	2371	CALL construct_filename1 ( inpname , 'wrfchemi_dms' , grid%id , 2 )
	2372	WRITE(message,*)'mediation_integrate: med_read_wrf_chem_dms_emiss: Open file ',TRIM(inpname)
	2373	CALL wrf_message( TRIM(message) )
	2374
	2375	if( grid%auxinput7_oid .NE. 0 ) then
	2376	CALL close_dataset ( grid%auxinput7_oid , config_flags , "DATASET=AUXINPUT7" )
	2377	endif
	2378
	2379	CALL open_r_dataset ( grid%auxinput7_oid, TRIM(inpname) , grid , config_flags, &
	2380	"DATASET=AUXINPUT7", ierr )
	2381	IF ( ierr .NE. 0 ) THEN
	2382	WRITE( message , * ) 'med_read_wrf_chem_dms_emiss: error opening ', TRIM( inpname )
	2383	CALL wrf_error_fatal( TRIM( message ) )
	2384	ENDIF
	2385
	2386	WRITE(message,*)'mediation_integrate: med_read_wrf_chem_dms_emiss: Read dms reference fields',&
	2387	TRIM(current_date_char)
	2388	CALL wrf_message( TRIM(message) )
	2389
	2390	CALL wrf_debug (100 , 'mediation_integrate: calling input_auxinput7' )
	2391	CALL input_auxinput7 ( grid%auxinput7_oid, grid , config_flags , ierr )
	2392
	2393	CALL close_dataset ( grid%auxinput7_oid , config_flags , "DATASET=AUXINPUT7" )
	2394
	2395	! ENDIF
	2396	CALL wrf_debug (100 , 'mediation_integrate: med_read_wrf_chem_dms_emiss: exit' )
	2397
	2398	END SUBROUTINE med_read_wrf_chem_dms_emiss
	2399
	2400	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
	2401	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
	2402
	2403	SUBROUTINE med_read_wrf_chem_gocart_bg ( grid , config_flags )
	2404	! Driver layer
	2405	USE module_domain , ONLY : domain , domain_clock_get
	2406	USE module_io_domain
	2407	USE module_timing
	2408	USE module_configure , ONLY : grid_config_rec_type
	2409	! Model layer
	2410	USE module_bc_time_utilities
	2411	#ifdef DM_PARALLEL
	2412	USE module_dm
	2413	#endif
	2414	USE module_date_time
	2415	USE module_utility
	2416
	2417	IMPLICIT NONE
	2418
	2419	! Arguments
	2420	TYPE(domain) :: grid
	2421
	2422	TYPE (grid_config_rec_type) , INTENT(IN) :: config_flags
	2423
	2424	! Local data
	2425	LOGICAL, EXTERNAL :: wrf_dm_on_monitor
	2426
	2427	INTEGER :: ierr, efid
	2428	REAL :: time, tupdate
	2429	real, allocatable :: dumc0(:,:,:)
	2430	CHARACTER (LEN=256) :: message, current_date_char, date_string
	2431	CHARACTER (LEN=80) :: inpname
	2432
	2433	#include <wrf_io_flags.h>
	2434	! IF ( grid%id .EQ. 1 ) THEN
	2435
	2436	CALL domain_clock_get( grid, current_timestr=current_date_char )
	2437
	2438	CALL construct_filename1 ( inpname , 'wrfchemi_gocart_bg' , grid%id , 2 )
	2439	WRITE(message,*)'mediation_integrate: med_read_wrf_chem_gocart_bg: Open file ',TRIM(inpname)
	2440	CALL wrf_message( TRIM(message) )
	2441
	2442	if( grid%auxinput8_oid .NE. 0 ) then
	2443	CALL close_dataset ( grid%auxinput8_oid , config_flags , "DATASET=AUXINPUT8" )
	2444	endif
	2445
	2446	CALL open_r_dataset ( grid%auxinput8_oid, TRIM(inpname) , grid , config_flags, &
	2447	"DATASET=AUXINPUT8", ierr )
	2448	IF ( ierr .NE. 0 ) THEN
	2449	WRITE( message , * ) 'med_read_wrf_chem_gocart_bg: error opening ', TRIM( inpname )
	2450	CALL wrf_error_fatal( TRIM( message ) )
	2451	ENDIF
	2452
	2453	WRITE(message,*)'mediation_integrate: med_read_wrf_chem_gocart_bg: Read gocart_bg at time ',&
	2454	TRIM(current_date_char)
	2455	CALL wrf_message( TRIM(message) )
	2456
	2457	CALL wrf_debug (100 , 'mediation_integrate: calling input_auxinput8' )
	2458	CALL input_auxinput8 ( grid%auxinput8_oid, grid , config_flags , ierr )
	2459
	2460	CALL close_dataset ( grid%auxinput8_oid , config_flags , "DATASET=AUXINPUT8" )
	2461
	2462	!
	2463	! CALL wrf_global_to_patch_real ( backg_no3_io , grid%backg_no3 , grid%domdesc, ' ' , 'xyz' , &
	2464	! ids, ide-1 , jds , jde-1 , kds , kde-1, &
	2465	! ims, ime , jms , jme , kms , kme , &
	2466	! ips, ipe , jps , jpe , kps , kpe )
	2467	!
	2468	! ENDIF
	2469	CALL wrf_debug (100 , 'mediation_integrate: med_read_wrf_chem_gocart_bg: exit' )
	2470
	2471	END SUBROUTINE med_read_wrf_chem_gocart_bg
	2472	SUBROUTINE med_read_wrf_chem_emissopt3 ( grid , config_flags )
	2473	! Driver layer
	2474	USE module_domain , ONLY : domain , domain_clock_get
	2475	USE module_io_domain
	2476	USE module_timing
	2477	USE module_configure , ONLY : grid_config_rec_type
	2478	! Model layer
	2479	USE module_bc_time_utilities
	2480	#ifdef DM_PARALLEL
	2481	USE module_dm
	2482	#endif
	2483	USE module_date_time
	2484	USE module_utility
	2485
	2486	IMPLICIT NONE
	2487
	2488	! Arguments
	2489	TYPE(domain) :: grid
	2490
	2491	TYPE (grid_config_rec_type) , INTENT(IN) :: config_flags
	2492
	2493	! Local data
	2494	LOGICAL, EXTERNAL :: wrf_dm_on_monitor
	2495
	2496	INTEGER :: ierr, efid
	2497	REAL :: time, tupdate
	2498	real, allocatable :: dumc0(:,:,:)
	2499	CHARACTER (LEN=256) :: message, current_date_char, date_string
	2500	CHARACTER (LEN=80) :: inpname
	2501
	2502	#include <wrf_io_flags.h>
	2503	! IF ( grid%id .EQ. 1 ) THEN
	2504
	2505	CALL domain_clock_get( grid, current_timestr=current_date_char )
	2506
	2507	CALL construct_filename1 ( inpname , 'wrffirechemi' , grid%id , 2 )
	2508	WRITE(message,*)'mediation_integrate: med_read_wrf_chem_fireemissions: Open file ',TRIM(inpname)
	2509	CALL wrf_message( TRIM(message) )
	2510
	2511	if( grid%auxinput7_oid .NE. 0 ) then
	2512	CALL close_dataset ( grid%auxinput7_oid , config_flags , "DATASET=AUXINPUT7" )
	2513	endif
	2514
	2515	CALL open_r_dataset ( grid%auxinput7_oid, TRIM(inpname) , grid , config_flags, &
	2516	"DATASET=AUXINPUT7", ierr )
	2517	IF ( ierr .NE. 0 ) THEN
	2518	WRITE( message , * ) 'med_read_wrf_chem_fireemissions: error opening ', TRIM( inpname )
	2519	CALL wrf_error_fatal( TRIM( message ) )
	2520	ENDIF
	2521
	2522	WRITE(message,*)'mediation_integrate: med_read_wrf_chem_fireemissions: Read fire emissions at time ',&
	2523	TRIM(current_date_char)
	2524	CALL wrf_message( TRIM(message) )
	2525
	2526	CALL wrf_debug (00 , 'mediation_integrate: calling input_auxinput7' )
	2527	CALL input_auxinput7 ( grid%auxinput7_oid, grid , config_flags , ierr )
	2528
	2529	CALL close_dataset ( grid%auxinput7_oid , config_flags , "DATASET=AUXINPUT7" )
	2530
	2531	! ENDIF
	2532	CALL wrf_debug (00 , 'mediation_integrate: med_read_wrf_chem_fireemissions: exit' )
	2533
	2534	END SUBROUTINE med_read_wrf_chem_emissopt3
	2535	#endif
	2536
	2537	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
	2538
	2539	#ifdef HWRF
	2540	!zhang's doing for outputing restart namelist parameters
	2541	RECURSIVE SUBROUTINE med_namelist_out ( grid , config_flags )
	2542	! Driver layer
	2543	USE module_domain , ONLY : domain, domain_clock_get
	2544	USE module_io_domain
	2545	USE module_timing
	2546	! Model layer
	2547	USE module_configure , ONLY : grid_config_rec_type
	2548	USE module_bc_time_utilities
	2549	!zhang new USE WRF_ESMF_MOD
	2550	USE module_utility
	2551	!zhang new ends
	2552
	2553	IMPLICIT NONE
	2554
	2555	! Arguments
	2556	TYPE(domain), INTENT(IN) :: grid
	2557	TYPE (grid_config_rec_type) , INTENT(IN) :: config_flags
	2558
	2559	! Local
	2560	!zhang new TYPE(ESMF_Time) :: CurrTime
	2561	TYPE(WRFU_Time) :: CurrTime
	2562	INTEGER :: nout,rc,kid
	2563	INTEGER :: hr, min, sec, ms,julyr,julday
	2564	REAL :: GMT
	2565	CHARACTER*80 :: prefix, outname
	2566	CHARACTER*80 :: timestr
	2567	LOGICAL :: exist
	2568	LOGICAL,EXTERNAL :: wrf_dm_on_monitor
	2569
	2570	TYPE (grid_config_rec_type) :: kid_config_flags
	2571
	2572	!zhang new
	2573	IF ( wrf_dm_on_monitor() ) THEN
	2574	CALL start_timing
	2575	END IF
	2576
	2577	prefix = "wrfnamelist_d<domain>_<date>"
	2578	nout = 99
	2579
	2580	!zhang new CALL ESMF_ClockGet( grid%domain_clock, CurrTime=CurrTime, rc=rc )
	2581	!zhang new CALL wrf_timetoa ( CurrTime, timestr )
	2582	CALL domain_clock_get( grid, current_timestr=timestr )
	2583	!zhang new ends
	2584	CALL construct_filename2a ( outname , prefix, grid%id , 2 , timestr )
	2585
	2586	IF ( wrf_dm_on_monitor() ) THEN
	2587
	2588	CLOSE (NOUT)
	2589	OPEN ( FILE = trim(outname) , UNIT = nout, STATUS = 'UNKNOWN', FORM = 'FORMATTED')
	2590	!zhang new CALL ESMF_TimeGet( grid%current_time, YY=julyr, dayOfYear=julday, H=hr, M=min, S=sec, MS=ms, rc=rc)
	2591	CALL domain_clock_get( grid, current_time=CurrTime )
	2592	CALL WRFU_TimeGet( CurrTime, YY=julyr, dayOfYear=julday, H=hr, M=min, S=sec, MS=ms, rc=rc)
	2593	!zhang new ends
	2594	gmt=hr+real(min)/60.+real(sec)/3600.+real(ms)/(1000*3600)
	2595	WRITE(NOUT,*) grid%i_parent_start
	2596	WRITE(NOUT,*) grid%j_parent_start
	2597	WRITE(NOUT,*) julyr
	2598	WRITE(NOUT,*) julday
	2599	WRITE(NOUT,*) gmt
	2600
	2601	CLOSE (NOUT)
	2602	ENDIF
	2603
	2604	! call recursively for children, (if any)
	2605	DO kid = 1, max_nests
	2606	IF ( ASSOCIATED( grid%nests(kid)%ptr ) ) THEN
	2607	CALL model_to_grid_config_rec ( grid%nests(kid)%ptr%id , model_config_rec , kid_config_flags )
	2608	CALL med_namelist_out ( grid%nests(kid)%ptr , kid_config_flags )
	2609	ENDIF
	2610	ENDDO
	2611
	2612	RETURN
	2613	END SUBROUTINE med_namelist_out
	2614	!end of zhang's doing
	2615	#endif

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

Download in other formats:

Original Format