Opened 10 years ago
Last modified 10 years ago
#678 closed bug
IMPORTANT: Problem with multiple telescoping nested grids — at Initial Version
| Reported by: | arango | Owned by: | arango |
|---|---|---|---|
| Priority: | major | Milestone: | Release ROMS/TOMS 3.7 |
| Component: | Nonlinear | Version: | 3.7 |
| Keywords: | Cc: |
Description
There seems to be a problem when setting the time-stepping in multiple (more than one) telescoping grids. Recall that technically inside ROMS, a telescoping grid is a refined grid with RefineScale > 0 including another finer refined grid inside. For more information about telescoping grids check the following information in the forum.
The variable RefineStepsCounter in ntimestep.F need to be set to zero every time that the parent coarser grid steps for all of telescoped finer resolution children grids:
ELSE IF (RefinedGrid(ng).and.(RefineScale(ng).gt.0)) THEN
IF (step_counter(ng).le.(WindowSteps(ig)+1)) THEN
IF (Telescoping(ng)) THEN
my_Nsteps(ig)=1
step_counter(ng)=step_counter(ng)+1
RefineStepsCounter(ng)=RefineStepsCounter(ng)+1
DO il=nl+1,NestLayers ! When a parent steps,
gn=GridNumber(ig,il) ! set all its telescoped
IF (Telescoping(gn)) THEN ! children counters to
RefineStepsCounter(gn)=0 ! zero
END IF
END DO
ELSE
my_Nsteps(ig)=RefineSteps(ng)
step_counter(ng)=step_counter(ng)+RefineSteps(ng)
RefineStepsCounter(ng)=RefineStepsCounter(ng)+ &
& RefineSteps(ng)
END IF
ELSE
Many thanks to John Warner for correcting this problem.
