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.