verify_compatible_dim_sizes Subroutine

  subroutine verify_compatible_dim_sizes(ncid, dim_name, expected_size)
    use mp, only: mp_abort, iproc
    integer, intent(in) :: ncid, expected_size
    character(len = *), intent(in) :: dim_name
    integer :: actual_size
    actual_size = get_dim_size(ncid, dim_name)
    if (actual_size /= expected_size) then
       write(*,'("Restart error (",A,"): got/expected = ",I0,"/",I0," : ",I0)') &
            dim_name, actual_size, expected_size, iproc
       call mp_abort('Incompatible '//dim_name//' in restart', .true.)
    end if
  end subroutine verify_compatible_dim_sizes