Canyon test case: unable to open restart NetCDF file

[breed]

Hi,

I cannot get the canyon 2D test case to work, but weirdly the 3D case works fine. For 2D it's failing at the very start when it tries to create/open the restart ncfile. I checked the directory and it doesn't get created (whereas it does for the 3D case).

Code: Select all

 TIME-STEP YYYY-MM-DD hh:mm:ss.ss  KINETIC_ENRG   POTEN_ENRG    TOTAL_ENRG    NET_VOLUME
                       C => (i,j)       Cu            Cv           C Max       Max Speed

         1 0001-01-01 00:00:05.00  0.000000E+00  0.000000E+00  0.000000E+00  3.414983E+13
                          (00,00)  0.000000E+00  0.000000E+00  0.000000E+00  0.000000E+00
  DEF_RST_NF90     - inquiring restart file,           Grid 01: roms_rst.nc
 Found Error: 02   Line: 8960     Source: ROMS/Modules/mod_netcdf.F, netcdf_open

 NETCDF_OPEN - unable to open existing NetCDF file:
               roms_rst.nc
               call from:  ROMS/Utility/def_rst.F, def_rst_nf90
               No such file or directory                                                       
 Found Error: 03   Line: 1597     Source: ROMS/Utility/def_rst.F, def_rst_nf90

 DEF_RST_NF90 - unable to open restart NetCDF file: roms_rst.nc
 Found Error: 03   Line: 82       Source: ROMS/Utility/def_rst.F
 Found Error: 03   Line: 585      Source: ROMS/Nonlinear/output.F
 Found Error: 03   Line: 407      Source: ROMS/Nonlinear/main2d.F
 Found Error: 03   Line: 298      Source: ROMS/Drivers/nl_roms.h, ROMS_run

Does anyone have any ideas? Do I need to provide more information?

Thanks,
Brad

[MIDHILAVARNAV]

The reason for this is probably the location of the file directory you provided, ROMS could not locate the restart file location. Pl check if the location you provided in .in file is correct or not.

[wilkin]

That looks like a write permissions error.

ROMS is not reading the restart file ... because you have ANA_INITIAL and it creates the initial conditions internally.

The log shows the error occurs in DEF_RST_NF90, i.e., when ROMS is trying to define/create the restart file.

You may not have permissions to write in the output directory, or you have set a path to the restart file that does not exist.

[breed]

Thank you both for your suggestions.

That looks like a write permissions error.

ROMS is not reading the restart file ... because you have ANA_INITIAL and it creates the initial conditions internally.

The log shows the error occurs in DEF_RST_NF90, i.e., when ROMS is trying to define/create the restart file.

You may not have permissions to write in the output directory, or you have set a path to the restart file that does not exist.

What confuses me about it being a permissions issue is why the 3D works fine but the 2D doesn't. All model runs get written to the same parent directory.

It seems like the 2D run fails in def_rst_nf90 (line 1597) when trying to call netcdf_open - but this should be after it has already created the nc file earlier in the script. I don't quite understand the following section of code in def_rst.F (lines 144-148), but this appears to be the switch to ensure the netcdf file gets written (on line 163):

Code: Select all

IF (((nrrec(ng).eq.0).and.(iic(ng).eq.ntstart(ng))).or.           &
     &    ((nrrec(ng).ne.0).and.                                        &
     &     (TRIM(ncname).ne.TRIM(INI(ng)%name)))) THEN
        Ldefine=.TRUE.
      END IF

In the process of trying to diagnose this I noticed some other things that differ between the 2D and 3D cases that I don't understand. Please can you help clarify them if possible?

• In the 3D log file (attached) there is a 0 timestep, whilst there isn't for the 2D case

• In the 3D log file (attached) it says that it creates the history, average and restart files (in that order at timestep 0), whereas for 2D it just says 'inquiring restart file' (at timestep 1), without the his, avg first

[wilkin]

That does look like a bug.

There aren't a lot of of users running 2D ROMS, so it's possible some bug crept in with file handling.

[arango]

Yes, I fixed main2d.F some time ago in my research repositories, but I forgot to do it in the distributed svn and git repositories. Please update for code. Here is that trac ticket.

[breed]

Thank you for the update, works perfectly!