Problem with netcdf file for Dogbone example

[stesp]

I am trying to run the Dogbone refinement grid example, but it is having problems reading dogbone_ngc_refined.nc . If I do ncdump on the file I get

ncdump: v1hpg.c:190: v1h_get_nc_type: Assertion `type == NC_BYTE || type == NC_CHAR || type == NC_SHORT || type == NC_INT || type == NC_FLOAT || type == NC_DOUBLE' failed.

I get a similar error running ROMS.

Maybe i'm using the wrong version of the netcdf library (4.1.3)?

Thanks,

Steve Spall

[kapitzah]

I'm having the same problem:

"ncdump -h dogbone_ngc_refined.nc" fails with the error message "ncdump: dogbone_ngc_refined.nc: Argument list too long".

The "file" command returns: "dogbone_ngc_refined.nc: data" while it should return something in the form of "...: NetCDF Data Format data"

So there is obviously something wrong with the file.

Cheers
Hartmut

[stesp]

Does anyone have any ideas on this? I can try to figure it out myself, but i'd rather be testing my own refinement grids rather than troubleshooting the standard example.

[turuncu]

Hi,

I am having same problem. Is this file corrupted somehow? Any idea to run the doggone test case? It could be related with NetCDF version but i am not sure.

Thanks,

--ufuk

[jcwarner]

can you try this file?

[turuncu]

Hi,

Thanks for your help. The file is okay now but dogbone_ngc_composite.nc has also same issue. If you don't mind could you regenerate that too. By the way, i could not see any difference in the definition of CPP flags in ticket 673

https://www.myroms.org/projects/src/ticket/673

Regards,

--ufuk

[turuncu]

Hi Again,

I am trying to generate 'dogbone_ngc_composite.nc' using Matlab toolbox using following Matlab commands

input = {'roms/Data/ROMS/Grid/dogbone_grd_left.nc'; ...
'roms/Data/ROMS/Grid/dogbone_grd_right.nc'};
output = 'dogbone_ngc_composite.nc'

[S,G] = contact(input, output);

but it gives errors as

output =

dogbone_ngc_composite.nc


Summary of Contact Regions Processed:

Grid 01: roms/Data/ROMS/Grid/dogbone_grd_left.nc
Grid 02: roms/Data/ROMS/Grid/dogbone_grd_right.nc

Contact Donor Receiver
Region Grid Grid

01 01 02
02 02 01

Error using sub2ind (line 52)
Out of range subscript.

Error in contact>Hweights (line 2453)
RindexB = sub2ind([Ir, Jr], ...

Error in contact (line 352)
S = Hweights(G, S, ImposeMask);

Error in dogbone (line 5)
[S,G] = contact(input, output);

Do you have any idea bout it?

--ufuk