- Timestamp:
- 2016-01-08T10:35:19+01:00 (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2014/dev_r4704_NOC5_MPP_BDY_UPDATE/NEMOGCM/NEMO/OPA_SRC/SBC/sbctide.F90
- Property svn:keywords set to Id
r4292 r6225 6 6 !! History : 9.0 ! 2007 (O. Le Galloudec) Original code 7 7 !!---------------------------------------------------------------------- 8 USE oce ! ocean dynamics and tracers variables 9 USE dom_oce ! ocean space and time domain 10 USE phycst 11 USE daymod 12 USE dynspg_oce 13 USE tideini 8 USE oce ! ocean dynamics and tracers variables 9 USE dom_oce ! ocean space and time domain 10 USE phycst ! physical constant 11 USE daymod ! calandar 12 USE tideini ! 14 13 ! 15 USE i om16 USE i n_out_manager ! I/O units17 USE ioipsl 18 USE lbclnk 14 USE in_out_manager ! I/O units 15 USE iom ! xIOs server 16 USE ioipsl ! NetCDF IPSL library 17 USE lbclnk ! ocean lateral boundary conditions (or mpp link) 19 18 20 19 IMPLICIT NONE … … 36 35 !!---------------------------------------------------------------------- 37 36 !! NEMO/OPA 3.5 , NEMO Consortium (2013) 38 !! $Id :$37 !! $Id$ 39 38 !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 40 39 !!---------------------------------------------------------------------- … … 47 46 INTEGER, INTENT( in ) :: kt ! ocean time-step 48 47 INTEGER :: jk ! dummy loop index 48 INTEGER :: nsec_day_orig ! Temporary variable 49 49 !!---------------------------------------------------------------------- 50 51 IF( nsec_day == NINT(0.5_wp * rdt tra(1))) THEN ! start a new day50 51 IF( nsec_day == NINT(0.5_wp * rdt) .OR. kt == nit000 ) THEN ! start a new day 52 52 ! 53 53 IF( kt == nit000 ) THEN … … 60 60 pot_astro(:,:) = 0._wp 61 61 ! 62 ! If the run does not start from midnight then need to initialise tides 63 ! at the start of the current day (only occurs when kt==nit000) 64 ! Temporarily set nsec_day to beginning of day. 65 nsec_day_orig = nsec_day 66 IF ( nsec_day /= NINT(0.5_wp * rdt) ) THEN 67 kt_tide = kt - (nsec_day - 0.5_wp * rdt)/rdt 68 nsec_day = NINT(0.5_wp * rdt) 69 ELSE 70 kt_tide = kt 71 ENDIF 62 72 CALL tide_harmo( omega_tide, v0tide, utide, ftide, ntide, nb_harmo ) 63 73 ! 64 kt_tide = kt65 74 ! 66 75 IF(lwp) THEN … … 75 84 IF( ln_tide_pot ) CALL tide_init_potential 76 85 ! 86 ! Reset nsec_day 87 nsec_day = nsec_day_orig 77 88 ENDIF 78 89 !
Note: See TracChangeset
for help on using the changeset viewer.