- Timestamp:
- 2018-07-13T09:28:50+02:00 (6 years ago)
- Location:
- NEMO/branches/2018/dev_r9838_ENHANCE04_RK3
- Files:
-
- 1 edited
- 1 copied
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2018/dev_r9838_ENHANCE04_RK3/src/OFF/dtadyn.F90
r9598 r9939 438 438 ENDIF 439 439 440 sshb(:,:) = sshn(:,:) + atfp * ( sshb(:,:) - 2 * sshn(:,:) + ssha(:,:))! before <-- now filtered440 sshb(:,:) = sshn(:,:) + rn_atfp * ( sshb(:,:) - 2 * sshn(:,:) + ssha(:,:) ) ! before <-- now filtered 441 441 sshn(:,:) = ssha(:,:) 442 442 … … 511 511 INTEGER :: jk 512 512 REAL(wp), DIMENSION(jpi,jpj) :: zhdiv 513 REAL(wp) :: z2dt 514 !!---------------------------------------------------------------------- 515 ! 516 z2dt = 2._wp * rdt 513 !!---------------------------------------------------------------------- 517 514 ! 518 515 zhdiv(:,:) = 0._wp … … 521 518 END DO 522 519 ! ! Sea surface elevation time-stepping 523 pssha(:,:) = ( psshb(:,:) - z2dt * ( r1_rau0 * pemp(:,:) + zhdiv(:,:) ) ) * ssmask(:,:)520 pssha(:,:) = ( psshb(:,:) - rDt * ( r1_rho0 * pemp(:,:) + zhdiv(:,:) ) ) * ssmask(:,:) 524 521 ! ! 525 522 ! ! After acale factors at t-points ( z_star coordinate )
Note: See TracChangeset
for help on using the changeset viewer.