[MITgcm-devel] problem with my last checkin

[Martin Losch]

Hi there,

I have just checked in some functionality that is handy for rotated  
spherical grids (basically a few new scalar variables in PARAMS.h and  
a new subroutine that recomputes XC/YC/XG/YG in one special case. The  
new code does not change the verification experiments on my  
linux_ia32_g77-machine. But now I am rerunning testreport on  
hugo.csail.mit.edu with the same build_options_file and I get  
segmentation faults (linux_ia32_g77) for fizhi-cs-aqualev and fizhi- 
cs-32x32x40. Everything else is OK.

I have tried to find the problem, but the segmentation fault happens  
somewhere within fizhi, here' s the debugger output
>
> Program received signal SIGSEGV, Segmentation fault.
> 0x0805aca8 in solir_ (m=0xbf897574, n=0x8581468, ndim=0x8581468,  
> np=0x85810e0, wh=0xbf114b50, taucld=0xbf24cde0, tauclb=0xbf1c8b80,  
> tauclf=0xbf18cb70, reff=0xbf2c4df0, ict=0xbf897598, icb=0xbf89759c,
>     fcld=0xbf574f40, cc=0xbf204b90, taual=0xbf210dd0,  
> csm=0xbefdcad0, rsirbm=0xbf47aec0, rsirdf=0xbf478eb0,  
> flx=0xbf42ee40, flc=0xbf3f0e30, fdirir=0xbf472e80,  
> fdifir=0xbf470e70) at fizhi_swrad.f:1732
> 1732               ssaclt(i,j)=1.0
> Current language:  auto; currently fortran
>
I have no idea what's going on, and I can't even run g77 -fbounds- 
check, because in fizhi, there are so many assignments where this  
array bound check chockes, e.g. variable mndy(12,4) is access via DO  
I=1,48; mnc(I,1)=...; ENDDO, which is technically correct by makes it  
impossible to debug these files. What am I to do? Remove my changes  
again?

Martin