- Timestamp:
- 2015-05-18T17:19:50+02:00 (9 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/NEMO/SAS_SRC/daymod.F90
r4162 r5282 20 20 !! 21 21 !! we suppose that the time step is deviding the number of second of in a day 22 !! ---> MOD( rday, rdt tra(1)) == 022 !! ---> MOD( rday, rdt ) == 0 23 23 !! 24 24 !! ----------- WARNING ----------- … … 71 71 !!---------------------------------------------------------------------- 72 72 ! 73 ! all calendar staff is based on the fact that MOD( rday, rdt tra(1)) == 074 IF( MOD( rday , rdt tra(1)) /= 0. ) CALL ctl_stop( 'the time step must devide the number of second of in a day' )73 ! all calendar staff is based on the fact that MOD( rday, rdt ) == 0 74 IF( MOD( rday , rdt ) /= 0. ) CALL ctl_stop( 'the time step must devide the number of second of in a day' ) 75 75 IF( MOD( rday , 2. ) /= 0. ) CALL ctl_stop( 'the number of second of in a day must be an even number' ) 76 IF( MOD( rdt tra(1), 2. ) /= 0. ) CALL ctl_stop( 'the time step (in second) must be an even number' )76 IF( MOD( rdt, 2. ) /= 0. ) CALL ctl_stop( 'the time step (in second) must be an even number' ) 77 77 nsecd = NINT(rday ) 78 78 nsecd05 = NINT(0.5 * rday ) 79 ndt = NINT( rdt tra(1))80 ndt05 = NINT(0.5 * rdt tra(1))79 ndt = NINT( rdt) 80 ndt05 = NINT(0.5 * rdt) 81 81 82 82 ! parameters corresponding to nit000 - 1 (as we start the step loop with a call to day) 83 83 ndastp = ndate0 - 1 ! ndate0 read in the namelist in dom_nam, we assume that we start run at 00:00 84 adatrj = ( REAL( nit000-1, wp ) * rdt tra(1)) / rday84 adatrj = ( REAL( nit000-1, wp ) * rdt ) / rday 85 85 IF( ABS(adatrj - REAL(NINT(adatrj),wp)) < 0.1 / rday ) adatrj = REAL(NINT(adatrj),wp) ! avoid truncation error 86 86 ! … … 226 226 nsec_week = nsec_week + ndt 227 227 nsec_day = nsec_day + ndt 228 adatrj = adatrj + rdt tra(1)/ rday229 fjulday = fjulday + rdt tra(1)/ rday228 adatrj = adatrj + rdt / rday 229 fjulday = fjulday + rdt / rday 230 230 IF( ABS(fjulday - REAL(NINT(fjulday),wp)) < zprec ) fjulday = REAL(NINT(fjulday),wp) ! avoid truncation error 231 231 IF( ABS(adatrj - REAL(NINT(adatrj ),wp)) < zprec ) adatrj = REAL(NINT(adatrj ),wp) ! avoid truncation error
Note: See TracChangeset
for help on using the changeset viewer.