﻿id	summary	reporter	owner	description	type	status	priority	milestone	component	version	resolution	keywords	cc
54	ROPP FM Coding Problem	frae	frcm	"The g95 bails out in ropp_fm/common/ropp_fm_obs2roprof.f90:186

{{{
  do i = 1, int((int(sqrt(8.0_wp * size(y%cov%d)) + 1.0_wp) - 1.0_wp) / 2)
     ro_data%Lev2a%refrac_sigma = sqrt(y%cov%d(i + i*(i-1)/2)) ! pp, uplo = 'U'
  enddo
}}}

Not sure if I understand this construct correctly, the first int(...) in the for loop was introduced by me since the g95 compiler was complaining. But it seems that the size(y%cov%d) is actually 1, and the loop runs to 1 before it core dumps. Mabye the netcdf file should have a covariace in there, the ones I use at the moment have not (also the ones in ropp_fm/data don't seem to have one). So there should be an error message somewhere. Also, should errors be part of the FM, or should they be just in the 1DVar part?

In addition, the ro_data%Lev2a%refrac_sigma = puzzles me, why is there no index in there. Maybe my f90 knowledge is failing me here.

The same problem will likely also occur a few lines earlier when the same is done for bending angles.

Funny enough the Intel Fortran compiler does not seem to mind at all.
"	defect	closed	normal		ropp_fm	0.8	fixed		dave.offiler@… axel.vonengeln@…
