[MITgcm-devel] problem with my last checkin

Jean-Michel Campin jmc at ocean.mit.edu
Tue Feb 5 09:48:20 EST 2008


Hi Martin,

It's a good idea to have this option to rotate the grid.
Regarding the fizhi problem,
will try to have a look at this today (so don't remove
your changes now).
Jean-Michel

On Tue, Feb 05, 2008 at 03:06:41PM +0100, Martin Losch wrote:
> 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
> 
> _______________________________________________
> MITgcm-devel mailing list
> MITgcm-devel at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-devel



More information about the MITgcm-devel mailing list