Skip to content

Helpful printout info: 1) gross check for d01 time step, 2) current and sim start date in printout #784

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
davegill merged 8 commits into from
Mar 18, 2019
Merged

Helpful printout info: 1) gross check for d01 time step, 2) current and sim start date in printout #784

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
8 commits
Select commit Hold shift + click to select a range
851dc07
Gross check for d01 timestep
davegill Feb 12, 2019
52b722a
Move the test from check_a_mundo to start_em, allows map_proj comparison
davegill Mar 9, 2019
7a4cb3b
removing debugging prints, oops
davegill Mar 9, 2019
a70169e
replacing blank line that was unintentionally removed
davegill Mar 9, 2019
9bb8a8a
sheesh with the spaces
davegill Mar 9, 2019
cc295f3
Add in prints towards top of printou: dt, dx, and ratio
davegill Mar 11, 2019
bf18bb9
Add date info to the time info at the top of the model print out
davegill Mar 12, 2019
5f7e72f
reasonable time step ratio in registry is now 6, not 10
davegill Mar 12, 2019
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
2 changes: 2 additions & 0 deletions Registry/Registry.EM_COMMON
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -1980,6 +1980,8 @@ rconfig integer time_step_fract_num namelist,domains 1 0
rconfig integer time_step_fract_den namelist,domains 1 1 ih "time_step_fract_den"
rconfig integer time_step_dfi namelist,domains 1 -1 ih "time_step_dfi"

rconfig real reasonable_time_step_ratio namelist,domains 1 6. ih "reasonable_time_step_ratio" "Any d01, real-data case with a time step ratio larger than this is stopped"

rconfig integer min_time_step namelist,domains max_domains -1 h "min_time_step"
rconfig integer min_time_step_den namelist,domains max_domains 0 h "min_time_step denominator"
rconfig integer max_time_step namelist,domains max_domains -1 h "max_time_step"
Expand Down
62 changes: 50 additions & 12 deletions dyn_em/start_em.F
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -126,6 +126,8 @@ SUBROUTINE start_domain_em ( grid, allowed_to_read &
! CCN for MP=18 initializatio
REAL :: ccn_max_val

REAL :: dt_s, dx_km

REAL :: max_mf, max_rot_angle
CALL get_ijk_from_grid ( grid , &
ids, ide, jds, jde, kds, kde, &
Expand Down Expand Up @@ -233,23 +235,59 @@ SUBROUTINE start_domain_em ( grid, allowed_to_read &
first_trip_for_this_domain = .TRUE.
ENDIF

! Print out the maximum map scale factor on the coarse domain
! Print out the maximum map scale factor on the coarse domain

IF ( ( first_trip_for_this_domain ) .AND. ( grid%id .EQ. 1 ) .AND. &
( .NOT. config_flags%polar ) ) THEN
max_mf = grid%msft(its,jts)
DO j=jts,MIN(jde-1,jte)
DO i=its,MIN(ide-1,ite)
max_mf = MAX ( max_mf , grid%msft(i,j) )
END DO
IF ( ( first_trip_for_this_domain ) .AND. ( grid%id .EQ. 1 ) .AND. &
( .NOT. config_flags%polar ) ) THEN
max_mf = grid%msft(its,jts)
DO j=jts,MIN(jde-1,jte)
DO i=its,MIN(ide-1,ite)
max_mf = MAX ( max_mf , grid%msft(i,j) )
END DO
END DO
#if ( defined(DM_PARALLEL) && ! defined(STUBMPI) )
max_mf = wrf_dm_max_real ( max_mf )
max_mf = wrf_dm_max_real ( max_mf )
#endif
WRITE ( a_message , FMT='(A,F5.2,A)' ) 'Max map factor in domain 1 = ',max_mf, &
'. Scale the dt in the model accordingly.'
CALL wrf_message ( a_message )
WRITE ( a_message , FMT='(A,F5.2,A)' ) 'Max map factor in domain 1 = ',max_mf, &
'. Scale the dt in the model accordingly.'
CALL wrf_message ( a_message )
END IF

!-----------------------------------------------------------------------
! Do a gross check on the time step for the most coarse grid (not for
! adaptive time steps. If the value is significantly larger than reasonable,
! then stop immediately. Likely the user made a mistake.
! If the rule of thumb is 6 DX = DT, then anything > 10 is flagged.
! Only look at domain 1, only do this for the firt time into this routine,
! only look at real-data cases (where the map projection is not Cartesian),
! and only look at user-specified time steps (not the adaptive dt option).
!-----------------------------------------------------------------------
IF ( ( grid%id .EQ. 1 ) .AND. &
( first_trip_for_this_domain ) .AND. &
( config_flags%map_proj .NE. 0 ) .AND. &
( .NOT. config_flags%use_adaptive_time_step ) ) THEN
dt_s = REAL(config_flags%time_step) + &
REAL(config_flags%time_step_fract_num) / &
REAL(config_flags%time_step_fract_den)
dx_km = MIN ( config_flags%dx , config_flags%dy ) / 1000.
WRITE (message,*) 'D01: Time step = ',dt_s ,' (s)'
CALL wrf_message ( TRIM(message) )
WRITE (message,*) 'D01: Grid Distance = ',dx_km ,' (km)'
CALL wrf_message ( TRIM(message) )
WRITE (message,*) 'D01: Grid Distance Ratio dt/dx = ', dt_s / dx_km, ' (s/km)'
CALL wrf_message ( TRIM(message) )
WRITE (message,*) 'D01: Ratio Including Maximum Map Factor = ', dt_s / (dx_km / max_mf) , ' (s/km)'
CALL wrf_message ( TRIM(message) )
WRITE (message,*) 'D01: NML defined reasonable_time_step_ratio = ', config_flags%reasonable_time_step_ratio
CALL wrf_message ( TRIM(message) )
IF ( dt_s / dx_km > config_flags%reasonable_time_step_ratio ) THEN
CALL wrf_message ( 'The time step is probably too large for this grid distance, reduce it.' )
WRITE (message, * ) 'If you are sure of your settings, set reasonable_time_step_ratio in namelist.input > ' &
,dt_s / (dx_km / max_mf)
CALL wrf_message ( TRIM(message) )
CALL wrf_error_fatal ( '--- ERROR: Time step too large')
END IF
END IF

if(config_flags%cycling) then
! Clear the buckets for diagnostics at initial time
Expand Down
10 changes: 10 additions & 0 deletions share/input_wrf.F
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -430,6 +430,10 @@ SUBROUTINE input_wrf ( fid , grid , config_flags , switch , ierr )

#if (EM_CORE == 1)

IF ( switch .EQ. input_only ) then
CALL wrf_get_dom_ti_char ( fid , 'SIMULATION_START_DATE' , simulation_start_date , ierr )
END IF

! Test to make sure that the grid distances are the right size.

CALL wrf_get_dom_ti_real ( fid , 'DX' , dx_compare , 1 , icnt , ierr )
Expand Down Expand Up @@ -1140,6 +1144,12 @@ SUBROUTINE input_wrf ( fid , grid , config_flags , switch , ierr )
CALL domain_clock_get( grid, current_time=currtime, &
current_timestr=currtimestr )
#if (DA_CORE != 1)
IF ( (switch .EQ. input_only) .and. (grid%id .EQ. 1) .and. (grid%itimestep .EQ. 0) ) THEN
WRITE( wrf_err_message , * ) 'CURRENT DATE = ',TRIM( currtimestr )
CALL wrf_message ( trim(wrf_err_message) )
WRITE( wrf_err_message , * ) 'SIMULATION START DATE = ', simulation_start_date(1:19)
CALL wrf_message ( TRIM(wrf_err_message ) )
END IF
! Don't perform the check for WRFVAR, as we're not passing the right dates
! around
CALL domain_clockprint(150, grid, &
Expand Down