Nesting Run Problem

[osua72]

Hello,

I am trying to create a nested ROMS model with two telescopic refined grids.

I have created the grids using the coarse2fine script with a Gfactor of 5 for both refined grids. I have also successfully created the contact.nc file using contact.m.

I created my frc.nc file for the outer (coarse) grid and set my boundaries as follows :

LBC(isFsur) == Cha Cha Cha Cha \ ! free-surface, Grid 1
Nes Nes Nes Nes \ ! free-surface, Grid 2
Nes Nes Nes Nes ! free-surface, Grid 3
LBC(isUbar) == Fla Fla Fla Fla \ ! 2D U-momentum, Grid 1
Nes Nes Nes Nes \ ! 2D U-momentum, Grid 2
Nes Nes Nes Nes ! 2D U-momentum, Grid 3
LBC(isVbar) == Fla Fla Fla Fla \ ! 2D V-momentum, Grid 1
Nes Nes Nes Nes \ ! 2D V-momentum, Grid 2
Nes Nes Nes Nes ! 2D V-momentum, Grid 3
LBC(isUvel) == Gra Gra Gra Gra \ ! 3D U-momentum, Grid 1
Nes Nes Nes Nes \ ! 3D U-momentum, Grid 2
Nes Nes Nes Nes ! 3D U-momentum, Grid 3
LBC(isVvel) == Gra Gra Gra Gra \ ! 3D V-momentum, Grid 1
Nes Nes Nes Nes \ ! 3D V-momentum, Grid 2
Nes Nes Nes Nes ! 3D V-momentum, Grid 3
LBC(isMtke) == Gra Gra Gra Gra \ ! mixing TKE, Grid 1
Nes Nes Nes Nes \ ! mixing TKE, Grid 2
Nes Nes Nes Nes ! mixing TKE, Grid 3

LBC(isTvar) == Gra Gra Gra Gra \ ! temperature, Grid 1
Gra Gra Gra Gra \ ! salinity, Grid 1
Nes Nes Nes Nes \ ! temperature, Grid 2
Nes Nes Nes Nes \ ! salinity, Grid 2
Nes Nes Nes Nes \ ! temperature, Grid 3
Nes Nes Nes Nes ! salinity, Grid 3

The model appears to compile fine but shortly after submitting the running script it breaks with the following error :

INP_PAR - illegal lateral boundary condition keyword: Cha
LBC(isFsur) == Cha Cha Cha Cha \ ! free-surface, Grid 1

Elapsed CPU time (seconds):


ROMS/TOMS - Output NetCDF summary for Grid 01:

ROMS/TOMS - Output NetCDF summary for Grid 02:

ROMS/TOMS - Output NetCDF summary for Grid 03:

ROMS/TOMS - Input error ............. exit_flag: 2


ERROR: Abnormal termination: NetCDF INPUT.
REASON: No error

Any clues to why this is happening would be much appreciated - see attached for header file if interested

[jcwarner]

just for fun, lets try

LBC(isFsur) == Cha Cha Cha Cha \
Nes Nes Nes Nes \
Nes Nes Nes Nes

and see if that makes it happy.

[osua72]

Same error, just minus the comment!

INP_PAR - illegal lateral boundary condition keyword: Cha
Nes Nes Nes Nes \

[osua72]

Not sure quite what I did - possibly just reduced the spacing between the inputs - but this problem has now been solved. Thank you John for your input

[TorresGarcia]

Hi,

I am trying to do a 4 Nested grid and I am having same problem when trying to run. I created all my parent grid and subsequently the childs unsing coarse2fine. I also create my contact file using contact.m

I was sucessfully able to run using two nested grid. However when adding the 3 grid or 4 grid I get the same error:

Input file:
! Number of nested grids.

Ngrids = 4

! Number of grid nesting layers. This parameter is used to allow refinement
! and composite grid combinations.

NestLayers = 4

! Number of grids in each nesting layer [1:NestLayers].

GridsInLayer = 1 1 1 1

! Grid dimension parameters. See notes below in the Glossary for how to set
! these parameters correctly.

Lm == 72 200 206 230 ! Number of I-direction INTERIOR RHO-points
Mm == 71 160 265 365 ! Number of J-direction INTERIOR RHO-points
N == 10 10 10 10 ! Number of vertical levels
ND == 0 0 0 0 ! Number of wave directional bins

Nbed = 1 1 1 1 ! Number of sediment bed layers

NAT = 2 2 2 2 ! Number of active tracers (usually, 2)
NPT = 0 0 0 0 ! Number of inactive passive tracers
NCS = 1 1 1 1 ! Number of cohesive (mud) sediment tracers
NNS = 0 0 0 0 ! Number of non-cohesive (sand) sediment tracers

! Domain decomposition parameters for serial, distributed-memory or
! shared-memory configurations used to determine tile horizontal range
! indices (Istr,Iend) and (Jstr,Jend), [1:Ngrids].

NtileI == 4 4 4 4 ! I-direction partition
NtileJ == 8 8 8 8 ! J-direction partition

! Set lateral boundary conditions keyword. Notice that a value is expected
! for each boundary segment per nested grid for each state variable.
!
! Each tracer variable requires [1:4,1:NAT+NPT,Ngrids] values. Otherwise,
! [1:4,1:Ngrids] values are expected for other variables. The boundary
! order is: 1=west, 2=south, 3=east, and 4=north. That is, anticlockwise
! starting at the western boundary.
!
! The keyword is case insensitive and usually has three characters. However,
! it is possible to have compound keywords, if applicable. For example, the
! keyword "RadNud" implies radiation boundary condition with nudging. This
! combination is usually used in active/passive radiation conditions.
!
! Keyword Lateral Boundary Condition Type
!
! Cha Chapman_implicit (free-surface)
! Che Chapman_explicit (free-surface)
! Cla Clamped
! Clo Closed
! Fla Flather (2D momentum) _____N_____ j=Mm
! Gra Gradient | 4 |
! Nes Nested (refinement) | |
! Nud Nudging 1 W E 3
! Per Periodic | |
! Rad Radiation |_____S_____|
! Red Reduced Physics (2D momentum) 2 j=1
! Shc Shchepetkin (2D momentum) i=1 i=Lm
!
! W S E N
! e o a o
! s u s r
! t t t t
! h h
!
! 1 2 3 4

LBC(isFsur) == Cha Cha Cha Clo \
Nes Nes Nes Nes \
Nes Nes Nes Nes \
Nes Nes Nes Nes ! free-surface
LBC(isUbar) == Fla Fla Fla Clo \
Nes Nes Nes Nes \
Nes Nes Nes Nes \
Nes Nes Nes Nes ! 2D U-momentum
LBC(isVbar) == Fla Fla Fla Clo \
Nes Nes Nes Nes \
Nes Nes Nes Nes \
Nes Nes Nes Nes ! 2D V-momentum
LBC(isUvel) == Rad Rad Rad Clo \
Nes Nes Nes Nes \
Nes Nes Nes Nes \
Nes Nes Nes Nes ! 3D U-momentum
LBC(isVvel) == Rad Rad Rad Clo \
Nes Nes Nes Nes \
Nes Nes Nes Nes \
Nes Nes Nes Nes ! 3D V-momentum
LBC(isMtke) == Rad Rad Rad Clo \
Nes Nes Nes Nes \
Nes Nes Nes Nes \
Nes Nes Nes Nes ! mixing TKE

LBC(isTvar) == Gra Gra Gra Clo \
Gra Gra Gra Clo ! temprature
Nes Nes Nes Nes \
Nes Nes Nes Nes ! salinity
Nes Nes Nes Nes \
Nes Nes Nes Nes
Nes Nes Nes Nes \
Nes Nes Nes Nes

! Wec boundary conditions

LBC(isU2Sd) == Gra Gra Gra Clo \
Nes Nes Nes Nes \
Nes Nes Nes Nes \
Nes Nes Nes Nes ! 2D U-stokes
LBC(isV2Sd) == Gra Gra Gra Clo \
Nes Nes Nes Nes \
Nes Nes Nes Nes \
Nes Nes Nes Nes ! 2D V-stokes
LBC(isU3Sd) == Gra Gra Gra Clo \
Nes Nes Nes Nes \
Nes Nes Nes Nes \
Nes Nes Nes Nes ! 3D U-stokes
LBC(isV3Sd) == Gra Gra Gra Clo \
Nes Nes Nes Nes \
Nes Nes Nes Nes \
Nes Nes Nes Nes ! 3D V-stokes




ERROR:
LOAD_LBC - illegal lateral boundary condition keyword: Cha
Nes Nes Nes Nes \

Elapsed CPU time (seconds):


ROMS/TOMS - Output NetCDF summary for Grid 01:

ROMS/TOMS - Output NetCDF summary for Grid 02:

ROMS/TOMS - Output NetCDF summary for Grid 03:

If anyone could help me if I am missing something in the text

[jcwarner]

your input files has a lot of 'tabs' in it.
No tabs allowed.
please remove all tabs and try again. only use the space bar. not sure if there is a way to show all the tab characters in emacs or someother software.

-j

[TorresGarcia]

Thank you. I used Notepad ++ and I could see clear all the tabs I had. That problem is solved, now it reads all my BC, input and output files.
Legna

[Vera]

hello~
I want to know that if I have two nested grids,can I set different N(number of vertical layers)for them?
Thanks in advance.

[jcwarner]

no. you need N= the same number for each nested grid.