- Timestamp:
- 2015-12-03T17:04:43+01:00 (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2015/dev_MetOffice_merge_2015/NEMOGCM/NEMO/OPA_SRC/SBC/sbctide.F90
r5930 r5997 46 46 INTEGER, INTENT( in ) :: kt ! ocean time-step 47 47 INTEGER :: jk ! dummy loop index 48 INTEGER :: nsec_day_orig ! Temporary variable 48 49 !!---------------------------------------------------------------------- 49 50 IF( nsec_day == NINT(0.5_wp * rdttra(1)) ) THEN ! start a new day50 51 IF( nsec_day == NINT(0.5_wp * rdttra(1)) .OR. kt == nit000 ) THEN ! start a new day 51 52 ! 52 53 IF( kt == nit000 ) THEN … … 59 60 pot_astro(:,:) = 0._wp 60 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 * rdttra(1)) ) THEN 67 kt_tide = kt - (nsec_day - 0.5_wp * rdttra(1))/rdttra(1) 68 nsec_day = NINT(0.5_wp * rdttra(1)) 69 ELSE 70 kt_tide = kt 71 ENDIF 61 72 CALL tide_harmo( omega_tide, v0tide, utide, ftide, ntide, nb_harmo ) 62 73 ! 63 kt_tide = kt64 74 ! 65 75 IF(lwp) THEN … … 74 84 IF( ln_tide_pot ) CALL tide_init_potential 75 85 ! 86 ! Reset nsec_day 87 nsec_day = nsec_day_orig 76 88 ENDIF 77 89 !
Note: See TracChangeset
for help on using the changeset viewer.