def METHOD_NAME(psy_trans, fortran_reader):
    ''' Check that the transformation correctly generates NEMO PSyIR for
    a simple case. '''
    code = '''subroutine basic_loop()
  integer, parameter :: jpi=16, jpj=16
  integer :: ji, jj
  real :: a(jpi, jpj), fconst
  do jj = 1, jpj
    do ji = 1, jpi
      a(ji) = fconst
    end do
  end do