Nesting Run Problem
Nesting Run Problem
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
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
- Attachments
-
nest_header.h
- (612 Bytes) Downloaded 388 times
Re: Nesting Run Problem
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.
LBC(isFsur) == Cha Cha Cha Cha \
Nes Nes Nes Nes \
Nes Nes Nes Nes
and see if that makes it happy.
Re: Nesting Run Problem
Same error, just minus the comment!
INP_PAR - illegal lateral boundary condition keyword: Cha
Nes Nes Nes Nes \
INP_PAR - illegal lateral boundary condition keyword: Cha
Nes Nes Nes Nes \
Re: Nesting Run Problem
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
-
- Posts: 6
- Joined: Wed May 09, 2012 8:41 pm
- Location: University of South Carolina
Re: Nesting Run Problem
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
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
- Attachments
-
ocean_crocker_reef.in
- Input File
- (118.2 KiB) Downloaded 409 times
Re: Nesting Run Problem
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
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
-
- Posts: 6
- Joined: Wed May 09, 2012 8:41 pm
- Location: University of South Carolina
Re: Nesting Run Problem
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
Legna
Re: Nesting Run Problem
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.
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.
Re: Nesting Run Problem
no. you need N= the same number for each nested grid.