id	summary	reporter	owner	description	type	status	priority	milestone	component	version	resolution	keywords	cc
15	Indices in ropp_fm_hybrid_phi_flv.f90	frcm	Huw Lewis	"In the routine for calculating geopotential height on hybrid vertical levels (ropp_fm_hybrid_phi_flv.f90), the following construction
{{{
! Calculate the integral (eq. 2.21)
! ---------------------------------

  do lvl = 1, size(p_hlv)
     phi_hlv(lvl) = phi_sfc - sum(del_phi(lvl:size(Tv_flv)))
  enddo
}}}
exists. however, the half levels (currently) have 61 elements, but there are only 60 full levels. Thus, for the calculation of phi_hlv(lvl), the result is 
{{{
phi_hlv(61) = phi_sfc - sum(del_phi(61:60))
}}}
where del_phi(61) isn't defined. Is that a bug in the implementation, or does it use a particular feature of Fortran 90?

Chris."	defect	closed	normal	1.2	ropp_fm	0.8	fixed