program formatted_output
  use global_parameters, only : wp
  implicit none
  integer :: i
  real (kind=wp), dimension (2,3) :: a ! a is 2x3 matrix

! a is assingned by an array constructor
  a = RESHAPE((/1.11111e300_wp,-2.2222e12_wp,3.3333_wp,-4.0e7_wp,-5.5e300_wp, &
    6.0_wp/),(/2,3/))

! List-directed formatting
  print *, a
  print *
! Explicit formatting with es (scientific form) edit descriptor
  print '(3es12.2)', (a(i,:),i=1,2)
  print *
  print '(3en12.2)', TRANSPOSE(a)
end program formatted_output