New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 5282 for branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/NEMO/SAS_SRC/daymod.F90 – NEMO

Ignore:
Timestamp:
2015-05-18T17:19:50+02:00 (9 years ago)
Author:
diovino
Message:

Dev. branch CMCC4_simplification ticket #1456

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/NEMO/SAS_SRC/daymod.F90

    r4162 r5282  
    2020   !! 
    2121   !!   we suppose that the time step is deviding the number of second of in a day 
    22    !!             ---> MOD( rday, rdttra(1) ) == 0 
     22   !!             ---> MOD( rday, rdt ) == 0 
    2323   !! 
    2424   !!           ----------- WARNING ----------- 
     
    7171      !!---------------------------------------------------------------------- 
    7272      ! 
    73       ! all calendar staff is based on the fact that MOD( rday, rdttra(1) ) == 0 
    74       IF( MOD( rday     , rdttra(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' ) 
    7575      IF( MOD( rday     , 2.        ) /= 0. )   CALL ctl_stop( 'the number of second of in a day must be an even number'    ) 
    76       IF( MOD( rdttra(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'           ) 
    7777      nsecd   = NINT(rday           ) 
    7878      nsecd05 = NINT(0.5 * rday     ) 
    79       ndt     = NINT(      rdttra(1)) 
    80       ndt05   = NINT(0.5 * rdttra(1)) 
     79      ndt     = NINT(      rdt) 
     80      ndt05   = NINT(0.5 * rdt) 
    8181 
    8282      ! parameters corresponding to nit000 - 1 (as we start the step loop with a call to day) 
    8383      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 ) * rdttra(1) ) / rday 
     84      adatrj = ( REAL( nit000-1, wp ) * rdt ) / rday 
    8585      IF( ABS(adatrj  - REAL(NINT(adatrj),wp)) < 0.1 / rday )   adatrj = REAL(NINT(adatrj),wp)   ! avoid truncation error 
    8686      ! 
     
    226226      nsec_week  = nsec_week  + ndt 
    227227      nsec_day   = nsec_day   + ndt                 
    228       adatrj  = adatrj  + rdttra(1) / rday 
    229       fjulday = fjulday + rdttra(1) / rday 
     228      adatrj  = adatrj  + rdt / rday 
     229      fjulday = fjulday + rdt / rday 
    230230      IF( ABS(fjulday - REAL(NINT(fjulday),wp)) < zprec )   fjulday = REAL(NINT(fjulday),wp)   ! avoid truncation error 
    231231      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.