Opened 16 years ago
Closed 9 years ago
#201 closed Enhancement (fixed)
simplified definition of the model time step
Reported by: | gm | Owned by: | nemo |
---|---|---|---|
Priority: | low | Milestone: | Unscheduled |
Component: | OCE | Version: | trunk |
Severity: | Keywords: | ||
Cc: |
Description
the model time step is 2*rdt on dynamics and 2*rdtra on tracers, except at the first time step when starting from rest or when starting from a short restart file (i.e. containing only 1 time step info). In this cas an euler time step is used and thus rdt on dyn and 2*rdtra on tracers. Therefore in many places (tranxt, dynnxt, trazdf, dynzdf, dynspg...) lines similar to the following lines can be found:
! set time step length (Euler/Leapfrog) IF( neuler == 0 .AND. kt == nit000 ) THEN ; r2dt_t(:) = rdttra(:) ! at nit000 (Euler) ELSEIF( kt <= nit000 + 1 ) THEN ; r2dt_t(:) = 2. * rdttra(:) ! at nit000 or nit000+1 (Leapfrog) ENDIF
A significant simplification can be achieved if rdt_d and rdt_t variables are introduced in dom_oce. These variables will be set to the model time step at the beginning of step and used in all routines that need the model time step on dynamics or tracers
Commit History (0)
(No commits)
Change History (2)
comment:1 Changed 13 years ago by gm
- Milestone changed from NEMO improvements to Miscellaneous Items
comment:2 Changed 9 years ago by clevy
- Resolution set to fixed
- Status changed from new to closed