Problem during nesting

[IBH123]

Dear All:

I'm trying to do a one-way nesting on a refinement grid, and I got the following error. Does anyone have an idea what might this error correspond to?

Code: Select all

At line 370 of file u2dbc_im.f90
Fortran runtime error: Index '-999' of dimension 1 of array 'ucontact%idg' below lower bound of 1

Error termination. Backtrace:
#0  0x10fd034 in __u2dbc_mod_MOD_u2dbc_tile
        at /gpfs/projects/WilsonGroup/ROMS/OIB_nwn/Build_romsG/u2dbc_im.f90:370
#1  0xc42c32 in ini_fields_tile
        at /gpfs/projects/WilsonGroup/ROMS/OIB_nwn/Build_romsG/ini_fields.f90:340
#2  0xc573a0 in __ini_fields_mod_MOD_ini_fields
        at /gpfs/projects/WilsonGroup/ROMS/OIB_nwn/Build_romsG/ini_fields.f90:82
#3  0x40c53e in main3d_
        at /gpfs/projects/WilsonGroup/ROMS/OIB_nwn/Build_romsG/main3d.f90:157
#4  0x407034 in __ocean_control_mod_MOD_roms_run
        at /gpfs/projects/WilsonGroup/ROMS/OIB_nwn/Build_romsG/ocean_control.f90:180
#5  0x404ee6 in ocean
        at /gpfs/projects/WilsonGroup/ROMS/OIB_nwn/Build_romsG/master.f90:108
#6  0x404f9e in main
        at /gpfs/projects/WilsonGroup/ROMS/OIB_nwn/Build_romsG/master.f90:50

The coarse grid can be run with no problem.
Does it means the connection file I use has problem? What is the best way to check that file? Or some other mistakes can also cause such error?

Thank you very much!

[IBH123]

Ah-ha, after playing with the examples for a while, I found this mistake is caused by masking out the boundary connect points on the land. If I turn off the Lmask option during the connectivity file generation, this mistake won't show up. But why do we have Lmask option for the connectivity file then?