Ignore:
Timestamp:
2020-02-19T10:22:15+01:00 (10 months ago)
Author:
davestorkey
Message:
  1. Rename the namelist timestep parameter: rn_rdt → rn_Dt
  2. Use the namelist parameter instead of the non-DOCTOR parameter: rdt → rn_Dt

This version has same SETTE results as the last version. Passes all tests but some
tests don't bit compare with the trunk.

Location:
NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2
Files:
135 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/cfgs/AGRIF_DEMO/EXPREF/1_namelist_cfg

    r12377 r12406  
    3030   ln_linssh   = .false.   !  =T  linear free surface  ==>>  model level are fixed in time 
    3131   ! 
    32    rn_rdt      = 5400.     !  time step for the dynamics and tracer 
     32   rn_Dt      = 5400.     !  time step for the dynamics and tracer 
    3333/ 
    3434!----------------------------------------------------------------------- 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/cfgs/AGRIF_DEMO/EXPREF/2_namelist_cfg

    r12377 r12406  
    3131   ln_linssh   = .false.   !  =T  linear free surface  ==>>  model level are fixed in time 
    3232   ! 
    33    rn_rdt      = 1350.     !  time step for the dynamics (and tracer if nn_acc=0) 
     33   rn_Dt      = 1350.     !  time step for the dynamics (and tracer if nn_acc=0) 
    3434/ 
    3535!----------------------------------------------------------------------- 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/cfgs/AGRIF_DEMO/EXPREF/3_namelist_cfg

    r12377 r12406  
    3131   ln_linssh   = .false.   !  =T  linear free surface  ==>>  model level are fixed in time 
    3232   ! 
    33    rn_rdt      = 450.     !  time step for the dynamics (and tracer if nn_acc=0) 
     33   rn_Dt      = 450.     !  time step for the dynamics (and tracer if nn_acc=0) 
    3434/ 
    3535!----------------------------------------------------------------------- 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/cfgs/AGRIF_DEMO/EXPREF/namelist_cfg

    r12377 r12406  
    3030   ln_linssh   = .false.   !  =T  linear free surface  ==>>  model level are fixed in time 
    3131   ! 
    32    rn_rdt      = 5400.     !  time step for the dynamics and tracer 
     32   rn_Dt      = 5400.     !  time step for the dynamics and tracer 
    3333/ 
    3434!----------------------------------------------------------------------- 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/cfgs/AMM12/EXPREF/namelist_cfg

    r12377 r12406  
    3333&namdom        !   time and space domain 
    3434!----------------------------------------------------------------------- 
    35    rn_rdt      =   600.    !  time step for the dynamics (and tracer if nn_acc=0) 
     35   rn_Dt      =   600.    !  time step for the dynamics (and tracer if nn_acc=0) 
    3636/ 
    3737!----------------------------------------------------------------------- 
     
    303303   ln_dynspg_ts = .true.   ! split-explicit free surface 
    304304   ln_bt_auto   = .false.  ! Number of sub-step defined from: 
    305    nn_baro      = 30       ! =F : the number of sub-step in rn_rdt seconds 
     305   nn_baro      = 30       ! =F : the number of sub-step in rn_Dt seconds 
    306306/ 
    307307!----------------------------------------------------------------------- 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/cfgs/C1D_PAPA/EXPREF/namelist_cfg

    r12377 r12406  
    4949&namdom        !   time and space domain 
    5050!----------------------------------------------------------------------- 
    51    rn_rdt      =  360.     !  time step for the dynamics and tracer 
     51   rn_Dt      =  360.     !  time step for the dynamics and tracer 
    5252/ 
    5353!----------------------------------------------------------------------- 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/cfgs/GYRE_BFM/EXPREF/namelist_cfg

    r12377 r12406  
    4545   ln_linssh   = .true.    !  =T  linear free surface  ==>>  model level are fixed in time 
    4646   ! 
    47    rn_rdt      = 7200.     !  time step for the dynamics 
     47   rn_Dt      = 7200.     !  time step for the dynamics 
    4848/ 
    4949 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/cfgs/GYRE_PISCES/EXPREF/namelist_cfg

    r12377 r12406  
    4545   ln_linssh   = .true.    !  =T  linear free surface  ==>>  model level are fixed in time 
    4646   ! 
    47    rn_rdt      = 7200.     !  time step for the dynamics 
     47   rn_Dt      = 7200.     !  time step for the dynamics 
    4848/ 
    4949!!====================================================================== 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/cfgs/ORCA2_ICE_ABL/EXPREF/namelist_cfg

    r12021 r12406  
    3131&namdom        !   time and space domain 
    3232!----------------------------------------------------------------------- 
    33    rn_rdt      = 5400.     !  time step for the dynamics and tracer 
     33   rn_Dt      = 5400.     !  time step for the dynamics and tracer 
    3434/ 
    3535!----------------------------------------------------------------------- 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/cfgs/ORCA2_ICE_PISCES/EXPREF/namelist_cfg

    r12377 r12406  
    2828&namdom        !   time and space domain 
    2929!----------------------------------------------------------------------- 
    30    rn_rdt      = 5400.     !  time step for the dynamics and tracer 
     30   rn_Dt      = 5400.     !  time step for the dynamics and tracer 
    3131/ 
    3232!----------------------------------------------------------------------- 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/cfgs/ORCA2_OFF_PISCES/EXPREF/namelist_cfg

    r12377 r12406  
    3434   ln_linssh   = .true.   !  =T  linear free surface  ==>>  model level are fixed in time 
    3535   ! 
    36    rn_rdt      = 21600.     !  time step for the dynamics and tracer 
     36   rn_Dt      = 21600.     !  time step for the dynamics and tracer 
    3737/ 
    3838!----------------------------------------------------------------------- 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/cfgs/ORCA2_OFF_TRC/EXPREF/namelist_cfg

    r12377 r12406  
    3434   ln_linssh   = .true.   !  =T  linear free surface  ==>>  model level are fixed in time 
    3535   ! 
    36    rn_rdt      = 21600.     !  time step for the dynamics and tracer 
     36   rn_Dt      = 21600.     !  time step for the dynamics and tracer 
    3737/ 
    3838!----------------------------------------------------------------------- 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/cfgs/SHARED/namelist_ref

    r12397 r12406  
    7272   ln_linssh   = .false.   !  =T  linear free surface  ==>>  model level are fixed in time 
    7373   ! 
    74    rn_rdt      = 5400.     !  time step for the dynamics and tracer 
     74   rn_Dt      = 5400.     !  time step for the dynamics and tracer 
    7575   rn_atfp     =    0.1    !  asselin time filter parameter 
    7676   ! 
     
    993993      ln_bt_auto    = .true.     ! Number of sub-step defined from: 
    994994         rn_bt_cmax   =  0.8        ! =T : the Maximum Courant Number allowed 
    995          nn_baro      = 30          ! =F : the number of sub-step in rn_rdt seconds 
     995         nn_baro      = 30          ! =F : the number of sub-step in rn_Dt seconds 
    996996      rn_bt_alpha   = 0.         ! Temporal diffusion parameter (if ln_bt_av=F) 
    997997/ 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/cfgs/SPITZ12/EXPREF/namelist_cfg

    r12377 r12406  
    2727&namdom        !   time and space domain 
    2828!----------------------------------------------------------------------- 
    29    rn_rdt      =  720.     !  time step for the dynamics and tracer 
     29   rn_Dt      =  720.     !  time step for the dynamics and tracer 
    3030   ! 
    3131/ 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/cfgs/WED025/EXPREF/namelist_cfg

    r12105 r12406  
    4545&namdom        !   time and space domain 
    4646!----------------------------------------------------------------------- 
    47    rn_rdt      = 1200.     !  time step for the dynamics and tracer 
     47   rn_Dt      = 1200.     !  time step for the dynamics and tracer 
    4848   ln_meshmask = .true.    !  =T create a mesh file 
    4949/ 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/ABL/abl.F90

    r11305 r12406  
    1111   USE dom_oce, ONLY: e1t, e1u, e1v, e1f         ! scale factors for horizontal grid 
    1212   USE dom_oce, ONLY: e2t, e2u, e2v, e2f         !  
    13    USE dom_oce, ONLY: rdt                        ! oceanic time-step 
     13   USE dom_oce, ONLY: rn_Dt                      ! oceanic time-step 
    1414   USE sbc_oce, ONLY: ght_abl, ghw_abl, e3t_abl, e3w_abl, jpka   ! scale factors and altitudes of ABL grid points in the vertical  
    1515   
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/ABL/ablmod.F90

    r12353 r12406  
    436436               zcff2 = jp_alp3_dyn * zsig**3 + jp_alp2_dyn * zsig**2   & 
    437437                  &  + jp_alp1_dyn * zsig    + jp_alp0_dyn 
    438                zcff  = (1._wp-zmsk) + zmsk * zcff2 * rdt   ! zcff = 1 for masked points 
    439                                                            ! rdt = rdt_abl / nn_fsbc                           
     438               zcff  = (1._wp-zmsk) + zmsk * zcff2 * rn_Dt   ! zcff = 1 for masked points 
     439                                                             ! rn_Dt = rdt_abl / nn_fsbc                           
    440440               zcff  = zcff * rest_eq(ji,jj) 
    441441               z_cft( ji, jj, jk ) = zcff 
     
    460460            zcff2 = jp_alp3_tra * zsig**3 + jp_alp2_tra * zsig**2   & 
    461461               &  + jp_alp1_tra * zsig    + jp_alp0_tra 
    462             zcff  = (1._wp-zmsk) + zmsk * zcff2 * rdt   ! zcff = 1 for masked points 
    463                                                         ! rdt = rdt_abl / nn_fsbc                           
     462            zcff  = (1._wp-zmsk) + zmsk * zcff2 * rn_Dt   ! zcff = 1 for masked points 
     463                                                          ! rn_Dt = rdt_abl / nn_fsbc                           
    464464            !z_cft( ji, jj, jk ) = zcff 
    465465            tq_abl( ji, jj, jk, nt_a, jp_ta ) = (1._wp - zcff ) * tq_abl( ji, jj, jk, nt_a, jp_ta )   & 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/ABL/sbcabl.F90

    r12353 r12406  
    202202 
    203203      ! ABL timestep 
    204       rdt_abl = nn_fsbc * rdt 
     204      rdt_abl = nn_fsbc * rn_Dt 
    205205 
    206206      ! Check parameters for dynamics 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/ICE/icedia.F90

    r12377 r12406  
    128128 
    129129      IF(  iom_use('ibgfrchfxtop') .OR. iom_use('ibgfrchfxbot') ) THEN 
    130          CALL iom_put( 'ibgfrchfxtop' , frc_temtop * z1_e1e2 * 1.e-20 * kt*rdt ) ! heat on top of ice/snw/ocean      (W/m2) 
    131          CALL iom_put( 'ibgfrchfxbot' , frc_tembot * z1_e1e2 * 1.e-20 * kt*rdt ) ! heat on top of ocean(below ice)   (W/m2)  
     130         CALL iom_put( 'ibgfrchfxtop' , frc_temtop * z1_e1e2 * 1.e-20 * kt*rn_Dt ) ! heat on top of ice/snw/ocean      (W/m2) 
     131         CALL iom_put( 'ibgfrchfxbot' , frc_tembot * z1_e1e2 * 1.e-20 * kt*rn_Dt ) ! heat on top of ocean(below ice)   (W/m2)  
    132132      ENDIF 
    133133       
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/ICE/icestp.F90

    r12377 r12406  
    338338      IF( ln_bdy .AND. ln_icediachk )   CALL ctl_warn('par_init: online conservation check does not work with BDY') 
    339339      ! 
    340       rdt_ice   = REAL(nn_fsbc) * rdt          !--- sea-ice timestep and its inverse 
     340      rdt_ice   = REAL(nn_fsbc) * rn_Dt          !--- sea-ice timestep and its inverse 
    341341      r1_rdtice = 1._wp / rdt_ice 
    342342      IF(lwp) WRITE(numout,*) 
    343       IF(lwp) WRITE(numout,*) '      ice timestep rdt_ice = nn_fsbc*rdt = ', rdt_ice 
     343      IF(lwp) WRITE(numout,*) '      ice timestep rdt_ice = nn_fsbc*rn_Dt = ', rdt_ice 
    344344      ! 
    345345      r1_nlay_i = 1._wp / REAL( nlay_i, wp )   !--- inverse of nlay_i and nlay_s 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/NST/agrif_oce_sponge.F90

    r12397 r12406  
    439439 
    440440         !* set relaxation time scale 
    441          IF( l_1st_euler .AND. lk_agrif_fstep ) THEN   ;   ztrelax =   rn_trelax_tra  / (        rdt ) 
    442          ELSE                                          ;   ztrelax =   rn_trelax_tra  / (2._wp * rdt ) 
     441         IF( l_1st_euler .AND. lk_agrif_fstep ) THEN   ;   ztrelax =   rn_trelax_tra  / (        rn_Dt ) 
     442         ELSE                                          ;   ztrelax =   rn_trelax_tra  / (2._wp * rn_Dt ) 
    443443         ENDIF 
    444444 
     
    596596#endif 
    597597         !* set relaxation time scale 
    598          IF( l_1st_euler .AND. lk_agrif_fstep ) THEN   ;   ztrelax =   rn_trelax_dyn  / (        rdt ) 
    599          ELSE                                          ;   ztrelax =   rn_trelax_dyn  / (2._wp * rdt ) 
     598         IF( l_1st_euler .AND. lk_agrif_fstep ) THEN   ;   ztrelax =   rn_trelax_dyn  / (        rn_Dt ) 
     599         ELSE                                          ;   ztrelax =   rn_trelax_dyn  / (2._wp * rn_Dt ) 
    600600         ENDIF 
    601601         ! 
     
    772772# endif 
    773773         !* set relaxation time scale 
    774          IF( l_1st_euler .AND. lk_agrif_fstep ) THEN   ;   ztrelax =   rn_trelax_dyn  / (        rdt ) 
    775          ELSE                                          ;   ztrelax =   rn_trelax_dyn  / (2._wp * rdt ) 
     774         IF( l_1st_euler .AND. lk_agrif_fstep ) THEN   ;   ztrelax =   rn_trelax_dyn  / (        rn_Dt ) 
     775         ELSE                                          ;   ztrelax =   rn_trelax_dyn  / (2._wp * rn_Dt ) 
    776776         ENDIF 
    777777         ! 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/NST/agrif_oce_update.F90

    r12397 r12406  
    11171117         IF (western_side) THEN 
    11181118            DO jj=j1,j2 
    1119                zcor = rdt * r1_e1e2t(i1  ,jj) * e2u(i1,jj) * (ub2_b(i1,jj)-tabres(i1,jj))  
     1119               zcor = rn_Dt * r1_e1e2t(i1  ,jj) * e2u(i1,jj) * (ub2_b(i1,jj)-tabres(i1,jj))  
    11201120               ssh(i1  ,jj,Kmm_a) = ssh(i1  ,jj,Kmm_a) + zcor 
    11211121               IF (.NOT.(lk_agrif_fstep.AND.(l_1st_euler))) ssh(i1  ,jj,Kbb_a) = ssh(i1  ,jj,Kbb_a) + atfp * zcor 
     
    11241124         IF (eastern_side) THEN 
    11251125            DO jj=j1,j2 
    1126                zcor = - rdt * r1_e1e2t(i2+1,jj) * e2u(i2,jj) * (ub2_b(i2,jj)-tabres(i2,jj)) 
     1126               zcor = - rn_Dt * r1_e1e2t(i2+1,jj) * e2u(i2,jj) * (ub2_b(i2,jj)-tabres(i2,jj)) 
    11271127               ssh(i2+1,jj,Kmm_a) = ssh(i2+1,jj,Kmm_a) + zcor 
    11281128               IF (.NOT.(lk_agrif_fstep.AND.(l_1st_euler))) ssh(i2+1,jj,Kbb_a) = ssh(i2+1,jj,Kbb_a) + atfp * zcor 
     
    12051205         IF (southern_side) THEN 
    12061206            DO ji=i1,i2 
    1207                zcor = rdt * r1_e1e2t(ji,j1  ) * e1v(ji,j1  ) * (vb2_b(ji,j1)-tabres(ji,j1)) 
     1207               zcor = rn_Dt * r1_e1e2t(ji,j1  ) * e1v(ji,j1  ) * (vb2_b(ji,j1)-tabres(ji,j1)) 
    12081208               ssh(ji,j1  ,Kmm_a) = ssh(ji,j1  ,Kmm_a) + zcor 
    12091209               IF (.NOT.(lk_agrif_fstep.AND.(l_1st_euler))) ssh(ji,j1  ,Kbb_a) = ssh(ji,j1,Kbb_a) + atfp * zcor 
     
    12121212         IF (northern_side) THEN                
    12131213            DO ji=i1,i2 
    1214                zcor = - rdt * r1_e1e2t(ji,j2+1) * e1v(ji,j2  ) * (vb2_b(ji,j2)-tabres(ji,j2)) 
     1214               zcor = - rn_Dt * r1_e1e2t(ji,j2+1) * e1v(ji,j2  ) * (vb2_b(ji,j2)-tabres(ji,j2)) 
    12151215               ssh(ji,j2+1,Kmm_a) = ssh(ji,j2+1,Kmm_a) + zcor 
    12161216               IF (.NOT.(lk_agrif_fstep.AND.(l_1st_euler))) ssh(ji,j2+1,Kbb_a) = ssh(ji,j2+1,Kbb_a) + atfp * zcor 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/NST/agrif_top_sponge.F90

    r12397 r12406  
    137137 
    138138         !* set relaxation time scale 
    139          IF( l_1st_euler .AND. lk_agrif_fstep ) THEN   ;   ztrelax =   rn_trelax_tra  / (        rdt ) 
    140          ELSE                                          ;   ztrelax =   rn_trelax_tra  / (2._wp * rdt ) 
     139         IF( l_1st_euler .AND. lk_agrif_fstep ) THEN   ;   ztrelax =   rn_trelax_tra  / (        rn_Dt ) 
     140         ELSE                                          ;   ztrelax =   rn_trelax_tra  / (2._wp * rn_Dt ) 
    141141         ENDIF 
    142142 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/NST/agrif_user.F90

    r12377 r12406  
    202202 
    203203         ! Check time steps            
    204          IF( NINT(Agrif_Rhot()) * NINT(rdt) .NE. Agrif_Parent(rdt) ) THEN 
    205             WRITE(cl_check1,*)  NINT(Agrif_Parent(rdt)) 
    206             WRITE(cl_check2,*)  NINT(rdt) 
    207             WRITE(cl_check3,*)  NINT(Agrif_Parent(rdt)/Agrif_Rhot()) 
     204         IF( NINT(Agrif_Rhot()) * NINT(rn_Dt) .NE. Agrif_Parent(rn_Dt) ) THEN 
     205            WRITE(cl_check1,*)  NINT(Agrif_Parent(rn_Dt)) 
     206            WRITE(cl_check2,*)  NINT(rn_Dt) 
     207            WRITE(cl_check3,*)  NINT(Agrif_Parent(rn_Dt)/Agrif_Rhot()) 
    208208            CALL ctl_stop( 'Incompatible time step between ocean grids',   & 
    209209                  &               'parent grid value : '//cl_check1    ,   &  
     
    613613      IF( check_namelist ) THEN 
    614614         ! Check time steps 
    615       IF( NINT(Agrif_Rhot()) * NINT(rdt) .NE. Agrif_Parent(rdt) ) THEN 
    616          WRITE(cl_check1,*)  Agrif_Parent(rdt) 
    617          WRITE(cl_check2,*)  rdt 
    618          WRITE(cl_check3,*)  rdt*Agrif_Rhot() 
     615      IF( NINT(Agrif_Rhot()) * NINT(rn_Dt) .NE. Agrif_Parent(rn_Dt) ) THEN 
     616         WRITE(cl_check1,*)  Agrif_Parent(rn_Dt) 
     617         WRITE(cl_check2,*)  rn_Dt 
     618         WRITE(cl_check3,*)  rn_Dt*Agrif_Rhot() 
    619619         CALL ctl_stop( 'incompatible time step between grids',   & 
    620620               &               'parent grid value : '//cl_check1    ,   &  
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OCE/ASM/asminc.F90

    r12397 r12406  
    534534            ! 
    535535            it = kt - nit000 + 1 
    536             zincwgt = wgtiau(it) / rdt   ! IAU weight for the current time step 
     536            zincwgt = wgtiau(it) / rn_Dt   ! IAU weight for the current time step 
    537537            ! 
    538538            IF(lwp) THEN 
     
    651651            ! 
    652652            it = kt - nit000 + 1 
    653             zincwgt = wgtiau(it) / rdt   ! IAU weight for the current time step 
     653            zincwgt = wgtiau(it) / rn_Dt   ! IAU weight for the current time step 
    654654            ! 
    655655            IF(lwp) THEN 
     
    722722            ! 
    723723            it = kt - nit000 + 1 
    724             zincwgt = wgtiau(it) / rdt   ! IAU weight for the current time step 
     724            zincwgt = wgtiau(it) / rn_Dt   ! IAU weight for the current time step 
    725725            ! 
    726726            IF(lwp) THEN 
     
    841841            it = kt - nit000 + 1 
    842842            zincwgt = wgtiau(it)      ! IAU weight for the current time step  
    843             ! note this is not a tendency so should not be divided by rdt (as with the tracer and other increments) 
     843            ! note this is not a tendency so should not be divided by rn_Dt (as with the tracer and other increments) 
    844844            ! 
    845845            IF(lwp) THEN 
     
    876876#if defined key_cice && defined key_asminc 
    877877            ! Sea-ice : CICE case. Pass ice increment tendency into CICE 
    878             ndaice_da(:,:) = seaice_bkginc(:,:) * zincwgt / rdt 
     878            ndaice_da(:,:) = seaice_bkginc(:,:) * zincwgt / rn_Dt 
    879879#endif 
    880880            ! 
     
    926926#if defined key_cice && defined key_asminc 
    927927            ! Sea-ice : CICE case. Pass ice increment tendency into CICE 
    928            ndaice_da(:,:) = seaice_bkginc(:,:) / rdt 
     928           ndaice_da(:,:) = seaice_bkginc(:,:) / rn_Dt 
    929929#endif 
    930930            IF ( .NOT. PRESENT(kindic) ) THEN 
     
    959959!           ! fwf : ice formation and melting 
    960960! 
    961 !                 zfons = ( -nfresh_da(ji,jj)*soce + nfsalt_da(ji,jj) )*rdt 
     961!                 zfons = ( -nfresh_da(ji,jj)*soce + nfsalt_da(ji,jj) )*rn_Dt 
    962962! 
    963963!           ! change salinity down to mixed layer depth 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OCE/BDY/bdylib.F90

    r12377 r12406  
    240240         ! Centred derivative is calculated as average of "left" and "right" derivatives for  
    241241         ! this reason.  
    242          ! Note no rdt factor in expression for zdt because it cancels in the expressions for  
     242         ! Note no rn_Dt factor in expression for zdt because it cancels in the expressions for  
    243243         ! zrx and zry. 
    244244         zdt   =     phia(iibm1   ,ijbm1   ) - phib(iibm1   ,ijbm1   ) 
     
    259259         zout = sign( 1., zrx ) 
    260260         zout = 0.5*( zout + abs(zout) ) 
    261          zwgt = 2.*rdt*( (1.-zout) * idx%nbd(jb,igrd) + zout * idx%nbdout(jb,igrd) ) 
     261         zwgt = 2.*rn_Dt*( (1.-zout) * idx%nbd(jb,igrd) + zout * idx%nbdout(jb,igrd) ) 
    262262         ! only apply radiation on outflow points  
    263263         if( ll_npo ) then     !! NPO version !! 
     
    425425            zout = sign( 1., zrx ) 
    426426            zout = 0.5*( zout + abs(zout) ) 
    427             zwgt = 2.*rdt*( (1.-zout) * idx%nbd(jb,igrd) + zout * idx%nbdout(jb,igrd) ) 
     427            zwgt = 2.*rn_Dt*( (1.-zout) * idx%nbd(jb,igrd) + zout * idx%nbdout(jb,igrd) ) 
    428428            ! only apply radiation on outflow points  
    429429            if( ll_npo ) then     !! NPO version !! 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OCE/BDY/bdytides.F90

    r12377 r12406  
    297297      ! Absolute time from model initialization:    
    298298      IF( PRESENT(kit) ) THEN   
    299          z_arg = ( REAL(kt, wp) + ( REAL(kit, wp) + zt_offset - 1. ) / REAL(nn_baro, wp) ) * rdt 
     299         z_arg = ( REAL(kt, wp) + ( REAL(kit, wp) + zt_offset - 1. ) / REAL(nn_baro, wp) ) * rn_Dt 
    300300      ELSE                               
    301          z_arg = ( REAL(kt, wp) + zt_offset ) * rdt 
     301         z_arg = ( REAL(kt, wp) + zt_offset ) * rn_Dt 
    302302      ENDIF 
    303303 
    304304      ! Linear ramp on tidal component at open boundaries  
    305305      zramp = 1. 
    306       IF (ln_tide_ramp) zramp = MIN(MAX( (z_arg - REAL(nit000,wp)*rdt)/(rn_tide_ramp_dt*rday),0.),1.) 
     306      IF (ln_tide_ramp) zramp = MIN(MAX( (z_arg - REAL(nit000,wp)*rn_Dt)/(rn_tide_ramp_dt*rday),0.),1.) 
    307307 
    308308      DO ib_bdy = 1,nb_bdy 
     
    319319            ! We refresh nodal factors every day below 
    320320            ! This should be done somewhere else 
    321             IF ( ( nsec_day == NINT(0.5_wp * rdt) .OR. kt==nit000 ) .AND. lk_first_btstp ) THEN 
    322                ! 
    323                kt_tide = kt - NINT((REAL(nsec_day,wp) - 0.5_wp * rdt)/rdt) 
     321            IF ( ( nsec_day == NINT(0.5_wp * rn_Dt) .OR. kt==nit000 ) .AND. lk_first_btstp ) THEN 
     322               ! 
     323               kt_tide = kt - NINT((REAL(nsec_day,wp) - 0.5_wp * rn_Dt)/rn_Dt) 
    324324               ! 
    325325               IF(lwp) THEN 
     
    333333               ! 
    334334            ENDIF 
    335             zoff = REAL(-kt_tide,wp) * rdt ! time offset relative to nodal factor computation time 
     335            zoff = REAL(-kt_tide,wp) * rn_Dt ! time offset relative to nodal factor computation time 
    336336            ! 
    337337            ! If time splitting, initialize arrays from slow varying open boundary data: 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OCE/DIA/dia25h.F90

    r12377 r12406  
    140140      ! ----------------- 
    141141      ! Define frequency of summing to create 25 h mean 
    142       IF( MOD( 3600,NINT(rdt) ) == 0 ) THEN 
    143          i_steps = 3600/NINT(rdt) 
     142      IF( MOD( 3600,NINT(rn_Dt) ) == 0 ) THEN 
     143         i_steps = 3600/NINT(rn_Dt) 
    144144      ELSE 
    145          CALL ctl_stop('STOP', 'dia_wri_tide: timestep must give MOD(3600,rdt) = 0 otherwise no hourly values are possible') 
     145         CALL ctl_stop('STOP', 'dia_wri_tide: timestep must give MOD(3600,rn_Dt) = 0 otherwise no hourly values are possible') 
    146146      ENDIF 
    147147 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OCE/DIA/diadetide.F90

    r12377 r12406  
    99   USE in_out_manager , ONLY :   lwp, numout 
    1010   USE iom            , ONLY :   iom_put 
    11    USE dom_oce        , ONLY :   rdt, nsec_day 
     11   USE dom_oce        , ONLY :   rn_Dt, nsec_day 
    1212   USE phycst         , ONLY :   rpi 
    1313   USE tide_mod 
     
    100100      zwght = 0.0_wp 
    101101      DO jn = 1, ndiadetide 
    102          ztmp = ( tdiadetide(jn) - REAL( nsec_day, KIND=wp ) ) / rdt 
     102         ztmp = ( tdiadetide(jn) - REAL( nsec_day, KIND=wp ) ) / rn_Dt 
    103103         IF ( ( ztmp < 0.5_wp ).AND.( ztmp >= -0.5_wp ) ) THEN 
    104104            zwght = zwght + 1.0_wp / REAL( ndiadetide, KIND=wp ) 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OCE/DIA/diahsb.F90

    r12377 r12406  
    121121      ENDIF 
    122122 
    123       frc_v = frc_v + z_frc_trd_v * rdt 
    124       frc_t = frc_t + z_frc_trd_t * rdt 
    125       frc_s = frc_s + z_frc_trd_s * rdt 
     123      frc_v = frc_v + z_frc_trd_v * rn_Dt 
     124      frc_t = frc_t + z_frc_trd_t * rn_Dt 
     125      frc_s = frc_s + z_frc_trd_s * rn_Dt 
    126126      !                                          ! Advection flux through fixed surface (z=0) 
    127127      IF( ln_linssh ) THEN 
    128          frc_wn_t = frc_wn_t + z_wn_trd_t * rdt 
    129          frc_wn_s = frc_wn_s + z_wn_trd_s * rdt 
     128         frc_wn_t = frc_wn_t + z_wn_trd_t * rn_Dt 
     129         frc_wn_s = frc_wn_s + z_wn_trd_s * rn_Dt 
    130130      ENDIF 
    131131 
     
    199199      CALL iom_put(   'bgfrctem' , frc_t    * rau0 * rcp * 1.e-20 )   ! hc  - surface forcing (1.e20 J)  
    200200      CALL iom_put(   'bgfrchfx' , frc_t    * rau0 * rcp /  &         ! hc  - surface forcing (W/m2)  
    201          &                       ( surf_tot * kt * rdt )        ) 
     201         &                       ( surf_tot * kt * rn_Dt )        ) 
    202202      CALL iom_put(   'bgfrcsal' , frc_s    * 1.e-9    )              ! sc  - surface forcing (psu*km3)  
    203203 
     
    207207         CALL iom_put( 'bgheatco' , zdiff_hc * 1.e-20 * rau0 * rcp )   ! Heat content drift    (1.e20 J)  
    208208         CALL iom_put( 'bgheatfx' , zdiff_hc * rau0 * rcp /  &         ! Heat flux drift       (W/m2)  
    209             &                       ( surf_tot * kt * rdt )        ) 
     209            &                       ( surf_tot * kt * rn_Dt )        ) 
    210210         CALL iom_put( 'bgsaltco' , zdiff_sc * 1.e-9    )              ! Salt content drift    (psu*km3) 
    211211         CALL iom_put( 'bgvolssh' , zdiff_v1 * 1.e-9    )              ! volume ssh drift      (km3)   
     
    227227         CALL iom_put( 'bgheatco' , zdiff_hc1 * 1.e-20 * rau0 * rcp )  ! Heat content drift    (1.e20 J)  
    228228         CALL iom_put( 'bgheatfx' , zdiff_hc1 * rau0 * rcp /  &        ! Heat flux drift       (W/m2)  
    229             &                       ( surf_tot * kt * rdt )         ) 
     229            &                       ( surf_tot * kt * rn_Dt )         ) 
    230230         CALL iom_put( 'bgsaltco' , zdiff_sc1 * 1.e-9    )             ! Salt content drift    (psu*km3) 
    231231         CALL iom_put( 'bgvolssh' , zdiff_v1 * 1.e-9    )              ! volume ssh drift      (km3)   
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OCE/DIA/dianam.F90

    r10068 r12406  
    7272 
    7373      IF( llfsec .OR. kfreq < 0 ) THEN   ;   inbsec = kfreq                       ! output frequency already in seconds 
    74       ELSE                               ;   inbsec = kfreq * NINT( rdt )   ! from time-step to seconds 
     74      ELSE                               ;   inbsec = kfreq * NINT( rn_Dt )   ! from time-step to seconds 
    7575      ENDIF 
    7676      iddss = NINT( rday          )                                         ! number of seconds in 1 day 
     
    116116      ! date of the beginning and the end of the run 
    117117 
    118       zdrun = rdt / rday * REAL( nitend - nit000, wp )                ! length of the run in days 
    119       zjul  = fjulday - rdt / rday 
     118      zdrun = rn_Dt / rday * REAL( nitend - nit000, wp )                ! length of the run in days 
     119      zjul  = fjulday - rn_Dt / rday 
    120120      CALL ju2ymds( zjul        , iyear1, imonth1, iday1, zsec1 )           ! year/month/day of the beginning of run 
    121121      CALL ju2ymds( zjul + zdrun, iyear2, imonth2, iday2, zsec2 )           ! year/month/day of the end       of run 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OCE/DIA/diawri.F90

    r12377 r12406  
    432432      clop = "x"         ! no use of the mask value (require less cpu time and otherwise the model crashes) 
    433433#if defined key_diainstant 
    434       zsto = nn_write * rdt 
     434      zsto = nn_write * rn_Dt 
    435435      clop = "inst("//TRIM(clop)//")" 
    436436#else 
    437       zsto=rdt 
     437      zsto=rn_Dt 
    438438      clop = "ave("//TRIM(clop)//")" 
    439439#endif 
    440       zout = nn_write * rdt 
    441       zmax = ( nitend - nit000 + 1 ) * rdt 
     440      zout = nn_write * rn_Dt 
     441      zmax = ( nitend - nit000 + 1 ) * rn_Dt 
    442442 
    443443      ! Define indices of the horizontal output zoom and vertical limit storage 
     
    460460 
    461461         ! Compute julian date from starting date of the run 
    462          CALL ymds2ju( nyear, nmonth, nday, rdt, zjulian ) 
     462         CALL ymds2ju( nyear, nmonth, nday, rn_Dt, zjulian ) 
    463463         zjulian = zjulian - adatrj   !   set calendar origin to the beginning of the experiment 
    464464         IF(lwp)WRITE(numout,*) 
     
    482482         CALL histbeg( clhstnam, jpi, glamt, jpj, gphit,           &  ! Horizontal grid: glamt and gphit 
    483483            &          iimi, iima-iimi+1, ijmi, ijma-ijmi+1,       & 
    484             &          nit000-1, zjulian, rdt, nh_T, nid_T, domain_id=nidom, snc4chunks=snc4set ) 
     484            &          nit000-1, zjulian, rn_Dt, nh_T, nid_T, domain_id=nidom, snc4chunks=snc4set ) 
    485485         CALL histvert( nid_T, "deptht", "Vertical T levels",      &  ! Vertical grid: gdept 
    486486            &           "m", ipk, gdept_1d, nz_T, "down" ) 
     
    518518         CALL histbeg( clhstnam, jpi, glamu, jpj, gphiu,           &  ! Horizontal grid: glamu and gphiu 
    519519            &          iimi, iima-iimi+1, ijmi, ijma-ijmi+1,       & 
    520             &          nit000-1, zjulian, rdt, nh_U, nid_U, domain_id=nidom, snc4chunks=snc4set ) 
     520            &          nit000-1, zjulian, rn_Dt, nh_U, nid_U, domain_id=nidom, snc4chunks=snc4set ) 
    521521         CALL histvert( nid_U, "depthu", "Vertical U levels",      &  ! Vertical grid: gdept 
    522522            &           "m", ipk, gdept_1d, nz_U, "down" ) 
     
    531531         CALL histbeg( clhstnam, jpi, glamv, jpj, gphiv,           &  ! Horizontal grid: glamv and gphiv 
    532532            &          iimi, iima-iimi+1, ijmi, ijma-ijmi+1,       & 
    533             &          nit000-1, zjulian, rdt, nh_V, nid_V, domain_id=nidom, snc4chunks=snc4set ) 
     533            &          nit000-1, zjulian, rn_Dt, nh_V, nid_V, domain_id=nidom, snc4chunks=snc4set ) 
    534534         CALL histvert( nid_V, "depthv", "Vertical V levels",      &  ! Vertical grid : gdept 
    535535            &          "m", ipk, gdept_1d, nz_V, "down" ) 
     
    544544         CALL histbeg( clhstnam, jpi, glamt, jpj, gphit,           &  ! Horizontal grid: glamt and gphit 
    545545            &          iimi, iima-iimi+1, ijmi, ijma-ijmi+1,       & 
    546             &          nit000-1, zjulian, rdt, nh_W, nid_W, domain_id=nidom, snc4chunks=snc4set ) 
     546            &          nit000-1, zjulian, rn_Dt, nh_W, nid_W, domain_id=nidom, snc4chunks=snc4set ) 
    547547         CALL histvert( nid_W, "depthw", "Vertical W levels",      &  ! Vertical grid: gdepw 
    548548            &          "m", ipk, gdepw_1d, nz_W, "down" ) 
     
    554554            CALL histbeg( clhstnam, jpi, glamt, jpj, gphit,           &  ! Horizontal grid: glamt and gphit 
    555555               &          iimi, iima-iimi+1, ijmi, ijma-ijmi+1,       & 
    556                &          nit000-1, zjulian, rdt, nh_A, nid_A, domain_id=nidom, snc4chunks=snc4set ) 
     556               &          nit000-1, zjulian, rn_Dt, nh_A, nid_A, domain_id=nidom, snc4chunks=snc4set ) 
    557557            CALL histvert( nid_A, "ght_abl", "Vertical T levels",      &  ! Vertical grid: gdept 
    558558               &           "m", ipka, ght_abl(2:jpka), nz_A, "up" ) 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OCE/DIU/diu_coolskin.F90

    r12377 r12406  
    6767 
    6868 
    69    SUBROUTINE diurnal_sst_coolskin_step(psqflux, pstauflux, psrho, rdt) 
     69   SUBROUTINE diurnal_sst_coolskin_step(psqflux, pstauflux, psrho, pDt) 
    7070      !!---------------------------------------------------------------------- 
    7171      !! *** ROUTINE diurnal_sst_takaya_step *** 
     
    8181      REAL(wp), INTENT(IN), DIMENSION(jpi,jpj) :: pstauflux   ! Wind stress (kg/ m s^2) 
    8282      REAL(wp), INTENT(IN), DIMENSION(jpi,jpj) :: psrho       ! Water density (kg/m^3) 
    83       REAL(wp), INTENT(IN) :: rdt                             ! Time-step 
     83      REAL(wp), INTENT(IN) :: pDt                             ! Time-step 
    8484      
    8585      ! Local variables 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OCE/DIU/diu_layers.F90

    r12377 r12406  
    3939      ! Cool skin 
    4040 
    41       CALL diurnal_sst_coolskin_step( qns, taum, rhop(:,:,1), rdt) 
     41      CALL diurnal_sst_coolskin_step( qns, taum, rhop(:,:,1), rn_Dt) 
    4242 
    4343      CALL iom_put( "sst_wl"   , x_dsst               )    ! warm layer (write out before update below). 
     
    4545 
    4646      ! Diurnal warm layer model        
    47       CALL diurnal_sst_takaya_step( kstp, qsr, qns, taum, rhop(:,:,1), rdt)  
     47      CALL diurnal_sst_takaya_step( kstp, qsr, qns, taum, rhop(:,:,1), rn_Dt)  
    4848 
    4949   END SUBROUTINE diurnal_layers   
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OCE/DOM/daymod.F90

    r12377 r12406  
    2020   !!                    ------------------------------- 
    2121   !!   sbcmod assume that the time step is dividing the number of second of  
    22    !!   in a day, i.e. ===> MOD( rday, rdt ) == 0  
     22   !!   in a day, i.e. ===> MOD( rday, rn_Dt ) == 0  
    2323   !!   except when user defined forcing is used (see sbcmod.F90) 
    2424   !!---------------------------------------------------------------------- 
     
    7373      ! 
    7474      ! max number of seconds between each restart 
    75       IF( REAL( nitend - nit000 + 1 ) * rdt > REAL( HUGE( nsec1jan000 ) ) ) THEN 
     75      IF( REAL( nitend - nit000 + 1 ) * rn_Dt > REAL( HUGE( nsec1jan000 ) ) ) THEN 
    7676         CALL ctl_stop( 'The number of seconds between each restart exceeds the integer 4 max value: 2^31-1. ',   & 
    7777            &           'You must do a restart at higher frequency (or remove this stop and recompile the code in I8)' ) 
     
    7979      nsecd   = NINT(       rday ) 
    8080      nsecd05 = NINT( 0.5 * rday ) 
    81       ndt     = NINT(       rdt  ) 
    82       ndt05   = NINT( 0.5 * rdt  ) 
     81      ndt     = NINT(       rn_Dt  ) 
     82      ndt05   = NINT( 0.5 * rn_Dt  ) 
    8383 
    8484      IF( .NOT. l_offline )   CALL day_rst( nit000, 'READ' ) 
     
    239239      nsec_monday  = nsec_monday  + ndt 
    240240      nsec_day   = nsec_day   + ndt 
    241       adatrj  = adatrj  + rdt / rday 
    242       fjulday = fjulday + rdt / rday 
     241      adatrj  = adatrj  + rn_Dt / rday 
     242      fjulday = fjulday + rn_Dt / rday 
    243243      IF( ABS(fjulday - REAL(NINT(fjulday),wp)) < zprec )   fjulday = REAL(NINT(fjulday),wp)   ! avoid truncation error 
    244244      IF( ABS(adatrj  - REAL(NINT(adatrj ),wp)) < zprec )   adatrj  = REAL(NINT(adatrj ),wp)   ! avoid truncation error 
     
    309309      !!       In both those options, the  exact duration of the experiment 
    310310      !!       since the beginning (cumulated duration of all previous restart runs) 
    311       !!       is not stored in the restart and is assumed to be (nit000-1)*rdt. 
     311      !!       is not stored in the restart and is assumed to be (nit000-1)*rn_Dt. 
    312312      !!       This is valid is the time step has remained constant. 
    313313      !! 
     
    379379               isecond = ( nhour * NINT(rhhmm) + nminute ) * NINT(rmmss) 
    380380               IF( isecond - ndt05 .lt. 0 )   ndastp = ndastp - 1      ! Start hour is specified in the namelist (default 0) 
    381                adatrj = ( REAL( nit000-1, wp ) * rdt ) / rday 
     381               adatrj = ( REAL( nit000-1, wp ) * rn_Dt ) / rday 
    382382               ! note this is wrong if time step has changed during run 
    383383            ENDIF 
     
    389389            isecond = ( nhour * NINT(rhhmm) + nminute ) * NINT(rmmss) 
    390390            IF( isecond - ndt05 .LT. 0 )   ndastp = ndastp - 1         ! Start hour is specified in the namelist (default 0) 
    391             adatrj = ( REAL( nit000-1, wp ) * rdt ) / rday 
     391            adatrj = ( REAL( nit000-1, wp ) * rn_Dt ) / rday 
    392392         ENDIF 
    393393         IF( ABS(adatrj  - REAL(NINT(adatrj),wp)) < 0.1 / rday )   adatrj = REAL(NINT(adatrj),wp)   ! avoid truncation error 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OCE/DOM/dom_oce.F90

    r12397 r12406  
    3333   LOGICAL , PUBLIC ::   ln_linssh      !: =T  linear free surface ==>> model level are fixed in time 
    3434   LOGICAL , PUBLIC ::   ln_meshmask    !: =T  create a mesh-mask file (mesh_mask.nc) 
    35    REAL(wp), PUBLIC ::   rn_rdt         !: time step for the dynamics and tracer 
     35   REAL(wp), PUBLIC ::   rn_Dt         !: time step for the dynamics and tracer 
    3636   REAL(wp), PUBLIC ::   rn_atfp        !: asselin time filter parameter 
    3737   LOGICAL , PUBLIC ::   ln_1st_euler   !: =T start with forward time step or not (=F) 
     
    4949   LOGICAL,  PUBLIC :: ln_bt_auto       !: Set number of barotropic iterations automatically 
    5050   INTEGER,  PUBLIC :: nn_bt_flt        !: Filter choice 
    51    INTEGER,  PUBLIC :: nn_baro          !: Number of barotropic iterations during one baroclinic step (rdt) 
     51   INTEGER,  PUBLIC :: nn_baro          !: Number of barotropic iterations during one baroclinic step (rn_Dt) 
    5252   REAL(wp), PUBLIC :: rn_bt_cmax       !: Maximum allowed courant number (used if ln_bt_auto=T) 
    5353   REAL(wp), PUBLIC :: rn_bt_alpha      !: Time stepping diffusion parameter 
     
    5656   !                                   !! old non-DOCTOR names still used in the model 
    5757   REAL(wp), PUBLIC ::   atfp           !: asselin time filter parameter 
    58    REAL(wp), PUBLIC ::   rdt            !: time step for the dynamics and tracer 
    5958 
    6059   !                                   !!! associated variables 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OCE/DOM/domain.F90

    r12397 r12406  
    289289         &             nn_stock, nn_write , ln_mskland  , ln_clobber   , nn_chunksz, ln_1st_euler  , & 
    290290         &             ln_cfmeta, ln_xios_read, nn_wxios 
    291       NAMELIST/namdom/ ln_linssh, rn_rdt, rn_atfp, ln_crs, ln_meshmask 
     291      NAMELIST/namdom/ ln_linssh, rn_Dt, rn_atfp, ln_crs, ln_meshmask 
    292292#if defined key_netcdf4 
    293293      NAMELIST/namnc4/ nn_nchunks_i, nn_nchunks_j, nn_nchunks_k, ln_nc4zip 
     
    408408         WRITE(numout,*) '      linear free surface (=T)                ln_linssh   = ', ln_linssh 
    409409         WRITE(numout,*) '      create mesh/mask file                   ln_meshmask = ', ln_meshmask 
    410          WRITE(numout,*) '      ocean time step                         rn_rdt      = ', rn_rdt 
     410         WRITE(numout,*) '      ocean time step                         rn_Dt      = ', rn_Dt 
    411411         WRITE(numout,*) '      asselin time filter parameter           rn_atfp     = ', rn_atfp 
    412412         WRITE(numout,*) '      online coarsening of dynamical fields   ln_crs      = ', ln_crs 
     
    415415      !          ! conversion DOCTOR names into model names (this should disappear soon) 
    416416      atfp = rn_atfp 
    417       rdt  = rn_rdt 
     417!!      Initialise current model timestep rDt = 2*rn_Dt if MLF or rDt = rn_Dt if RK3 
     418!!      rDt  = rn_Dt 
    418419 
    419420      IF( TRIM(Agrif_CFixed()) == '0' ) THEN 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OCE/DOM/domvvl.F90

    r12397 r12406  
    235235         IF( ln_vvl_ztilde_as_zstar ) THEN   ! z-star emulation using z-tile 
    236236            frq_rst_e3t(:,:) = 0._wp               !Ignore namelist settings 
    237             frq_rst_hdv(:,:) = 1._wp / rdt 
     237            frq_rst_hdv(:,:) = 1._wp / rn_Dt 
    238238         ENDIF 
    239239         IF ( ln_vvl_zstar_at_eqtor ) THEN   ! use z-star in vicinity of the Equator 
     
    247247                  ! values inside the equatorial band (ztilde as zstar) 
    248248                  frq_rst_e3t(ji,jj) =  0.0_wp 
    249                   frq_rst_hdv(ji,jj) =  1.0_wp / rdt 
     249                  frq_rst_hdv(ji,jj) =  1.0_wp / rn_Dt 
    250250               ELSE                                      ! transition band (2.5 to 6 degrees N/S) 
    251251                  !                                      ! (linearly transition from z-tilde to z-star) 
     
    253253                     &            * (  1.0_wp - COS( rad*(ABS(gphit(ji,jj))-2.5_wp)  & 
    254254                     &                                          * 180._wp / 3.5_wp ) ) 
    255                   frq_rst_hdv(ji,jj) = (1.0_wp / rdt)                                & 
    256                      &            + (  frq_rst_hdv(ji,jj)-(1.e0_wp / rdt) )*0.5_wp   & 
     255                  frq_rst_hdv(ji,jj) = (1.0_wp / rn_Dt)                                & 
     256                     &            + (  frq_rst_hdv(ji,jj)-(1.e0_wp / rn_Dt) )*0.5_wp   & 
    257257                     &            * (  1._wp  - COS( rad*(ABS(gphit(ji,jj))-2.5_wp)  & 
    258258                     &                                          * 180._wp / 3.5_wp ) ) 
     
    264264                  ij0 = 128   ;   ij1 = 135   ;    
    265265                  frq_rst_e3t( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) ) =  0.0_wp 
    266                   frq_rst_hdv( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) ) =  1.e0_wp / rdt 
     266                  frq_rst_hdv( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) ) =  1.e0_wp / rn_Dt 
    267267               ENDIF 
    268268            ENDIF 
     
    373373            IF( kt > nit000 ) THEN 
    374374               DO jk = 1, jpkm1 
    375                   hdiv_lf(:,:,jk) = hdiv_lf(:,:,jk) - rdt * frq_rst_hdv(:,:)   & 
     375                  hdiv_lf(:,:,jk) = hdiv_lf(:,:,jk) - rn_Dt * frq_rst_hdv(:,:)   & 
    376376                     &          * ( hdiv_lf(:,:,jk) - e3t(:,:,jk,Kmm) * ( hdiv(:,:,jk) - zhdiv(:,:) ) ) 
    377377               END DO 
     
    10011001            WRITE(numout,*) '                         rn_rst_e3t     = 0.e0' 
    10021002            WRITE(numout,*) '            hard-wired : z-tilde cutoff frequency of low-pass filter (days)' 
    1003             WRITE(numout,*) '                         rn_lf_cutoff   = 1.0/rdt' 
     1003            WRITE(numout,*) '                         rn_lf_cutoff   = 1.0/rn_Dt' 
    10041004         ELSE 
    10051005            WRITE(numout,*) '      z-tilde to zstar restoration timescale (days)        rn_rst_e3t   = ', rn_rst_e3t 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OCE/DYN/dynatf.F90

    r12397 r12406  
    195195            ! Add volume filter correction: compatibility with tracer advection scheme 
    196196            ! => time filter + conservation correction 
    197             zcoef = atfp * rdt * r1_rau0 
     197            zcoef = atfp * rn_Dt * r1_rau0 
    198198            zwfld(:,:) = emp_b(:,:) - emp(:,:) 
    199199            IF ( ln_rnf ) zwfld(:,:) =  zwfld(:,:) - ( rnf_b(:,:) - rnf(:,:) ) 
     
    207207            !     to manage rnf, isf and possibly in the futur icb, tide water glacier (...) 
    208208            !     ...(kt, coef, ktop, kbot, hz, fwf_b, fwf) 
    209             IF ( ln_isf ) CALL isf_dynatf( kt, Kmm, ze3t_f, atfp * rdt ) 
     209            IF ( ln_isf ) CALL isf_dynatf( kt, Kmm, ze3t_f, atfp * rn_Dt ) 
    210210            ! 
    211211            pe3t(:,:,1:jpkm1,Kmm) = ze3t_f(:,:,1:jpkm1)        ! filtered scale factor at T-points 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OCE/DYN/dynspg.F90

    r12377 r12406  
    6969      !!                apgu = - 1/rau0 di[apr] = 0.5*grav di[ssh_ib+ssh_ibb] 
    7070      !!                apgv = - 1/rau0 dj[apr] = 0.5*grav dj[ssh_ib+ssh_ibb] 
    71       !!             Note that as all external forcing a time averaging over a two rdt 
     71      !!             Note that as all external forcing a time averaging over a two rn_Dt 
    7272      !!             period is used to prevent the divergence of odd and even time step. 
    7373      !!---------------------------------------------------------------------- 
     
    114114            ! 
    115115            ! Update tide potential at the beginning of current time step 
    116             zt0step = REAL(nsec_day, wp)-0.5_wp*rdt 
     116            zt0step = REAL(nsec_day, wp)-0.5_wp*rn_Dt 
    117117            CALL upd_tide(zt0step, Kmm) 
    118118            ! 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OCE/DYN/dynspg_ts.F90

    r12397 r12406  
    429429         ! Update tide potential at the beginning of current time substep 
    430430         IF( ln_tide_pot .AND. ln_tide ) THEN 
    431             zt0substep = REAL(nsec_day, wp) - 0.5_wp*rdt + (jn + noffset - 1) * rdt / REAL(nn_baro, wp) 
     431            zt0substep = REAL(nsec_day, wp) - 0.5_wp*rn_Dt + (jn + noffset - 1) * rn_Dt / REAL(nn_baro, wp) 
    432432            CALL upd_tide(zt0substep, Kmm) 
    433433         END IF 
     
    976976 
    977977      ! Estimate number of iterations to satisfy a max courant number= rn_bt_cmax 
    978       IF( ln_bt_auto )   nn_baro = CEILING( rdt / rn_bt_cmax * zcmax) 
     978      IF( ln_bt_auto )   nn_baro = CEILING( rn_Dt / rn_bt_cmax * zcmax) 
    979979       
    980       rdtbt = rdt / REAL( nn_baro , wp ) 
     980      rdtbt = rn_Dt / REAL( nn_baro , wp ) 
    981981      zcmax = zcmax * rdtbt 
    982982      ! Print results 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OCE/DYN/sshwzv.F90

    r12397 r12406  
    223223      !!              from the filter, see Leclair and Madec 2010) and swap : 
    224224      !!                pssh(:,:,Kmm) = pssh(:,:,Kaa) + atfp * ( pssh(:,:,Kbb) -2 pssh(:,:,Kmm) + pssh(:,:,Kaa) ) 
    225       !!                            - atfp * rdt * ( emp_b - emp ) / rau0 
     225      !!                            - atfp * rn_Dt * ( emp_b - emp ) / rau0 
    226226      !! 
    227227      !! ** action  : - pssh(:,:,Kmm) time filtered 
     
    248248         pssh(:,:,Kmm) = pssh(:,:,Kmm) + atfp * ( pssh(:,:,Kbb) - 2 * pssh(:,:,Kmm) + pssh(:,:,Kaa) ) 
    249249         IF( .NOT.ln_linssh ) THEN                          ! "now" <-- with forcing removed 
    250             zcoef = atfp * rdt * r1_rau0 
     250            zcoef = atfp * rn_Dt * r1_rau0 
    251251            pssh(:,:,Kmm) = pssh(:,:,Kmm) - zcoef * (     emp_b(:,:) - emp   (:,:)   & 
    252252               &                             - rnf_b(:,:)        + rnf   (:,:)       & 
     
    255255 
    256256            ! ice sheet coupling 
    257             IF ( ln_isf .AND. ln_isfcpl .AND. kt == nit000+1) pssh(:,:,Kbb) = pssh(:,:,Kbb) - atfp * rdt * ( risfcpl_ssh(:,:) - 0.0 ) * ssmask(:,:) 
     257            IF ( ln_isf .AND. ln_isfcpl .AND. kt == nit000+1) pssh(:,:,Kbb) = pssh(:,:,Kbb) - atfp * rn_Dt * ( risfcpl_ssh(:,:) - 0.0 ) * ssmask(:,:) 
    258258 
    259259         ENDIF 
     
    306306         DO_3D_00_00( 1, jpkm1 ) 
    307307            z1_e3t = 1._wp / e3t(ji,jj,jk,Kmm) 
    308             ! 2*rdt and not r2dt (for restartability) 
    309             Cu_adv(ji,jj,jk) = 2._wp * rdt * ( ( MAX( ww(ji,jj,jk) , 0._wp ) - MIN( ww(ji,jj,jk+1) , 0._wp ) )                       &   
     308            ! 2*rn_Dt and not r2dt (for restartability) 
     309            Cu_adv(ji,jj,jk) = 2._wp * rn_Dt * ( ( MAX( ww(ji,jj,jk) , 0._wp ) - MIN( ww(ji,jj,jk+1) , 0._wp ) )                       &   
    310310               &                             + ( MAX( e2u(ji  ,jj)*e3u(ji  ,jj,jk,Kmm)*uu(ji  ,jj,jk,Kmm) + un_td(ji  ,jj,jk), 0._wp ) -   & 
    311311               &                                 MIN( e2u(ji-1,jj)*e3u(ji-1,jj,jk,Kmm)*uu(ji-1,jj,jk,Kmm) + un_td(ji-1,jj,jk), 0._wp ) )   & 
     
    319319         DO_3D_00_00( 1, jpkm1 ) 
    320320            z1_e3t = 1._wp / e3t(ji,jj,jk,Kmm) 
    321             ! 2*rdt and not r2dt (for restartability) 
    322             Cu_adv(ji,jj,jk) = 2._wp * rdt * ( ( MAX( ww(ji,jj,jk) , 0._wp ) - MIN( ww(ji,jj,jk+1) , 0._wp ) )   &  
     321            ! 2*rn_Dt and not r2dt (for restartability) 
     322            Cu_adv(ji,jj,jk) = 2._wp * rn_Dt * ( ( MAX( ww(ji,jj,jk) , 0._wp ) - MIN( ww(ji,jj,jk+1) , 0._wp ) )   &  
    323323               &                             + ( MAX( e2u(ji  ,jj)*e3u(ji  ,jj,jk,Kmm)*uu(ji  ,jj,jk,Kmm), 0._wp ) -   & 
    324324               &                                 MIN( e2u(ji-1,jj)*e3u(ji-1,jj,jk,Kmm)*uu(ji-1,jj,jk,Kmm), 0._wp ) )   & 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OCE/FLO/flo4rk.F90

    r12377 r12406  
    130130         ! computation of Runge-Kutta factor 
    131131         DO jfl = 1, jpnfl 
    132             zrkxfl(jfl,jind) = rdt*zufl(jfl) 
    133             zrkyfl(jfl,jind) = rdt*zvfl(jfl) 
    134             zrkzfl(jfl,jind) = rdt*zwfl(jfl) 
     132            zrkxfl(jfl,jind) = rn_Dt*zufl(jfl) 
     133            zrkyfl(jfl,jind) = rn_Dt*zvfl(jfl) 
     134            zrkzfl(jfl,jind) = rn_Dt*zwfl(jfl) 
    135135         END DO 
    136136         IF( jind /= 4 ) THEN 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OCE/FLO/floblk.F90

    r12377 r12406  
    233233            ! test to know if the "age" of the float is not bigger than the  
    234234            ! time step 
    235             IF( zagenewfl(jfl) > rdt ) THEN 
    236                zttfl(jfl) = (rdt-zagefl(jfl)) / zvol 
    237                zagenewfl(jfl) = rdt 
     235            IF( zagenewfl(jfl) > rn_Dt ) THEN 
     236               zttfl(jfl) = (rn_Dt-zagefl(jfl)) / zvol 
     237               zagenewfl(jfl) = rn_Dt 
    238238            ENDIF 
    239239             
     
    340340         ifin = 1 
    341341         DO jfl = 1, jpnfl 
    342             IF( zagefl(jfl) < rdt )   ifin = 0 
     342            IF( zagefl(jfl) < rn_Dt )   ifin = 0 
    343343            tpifl(jfl) = zgifl(jfl) + 0.5 
    344344            tpjfl(jfl) = zgjfl(jfl) + 0.5 
     
    347347         ifin = 1 
    348348         DO jfl = 1, jpnfl 
    349             IF( zagefl(jfl) < rdt )   ifin = 0 
     349            IF( zagefl(jfl) < rn_Dt )   ifin = 0 
    350350            tpifl(jfl) = zgifl(jfl) + 0.5 
    351351            tpjfl(jfl) = zgjfl(jfl) + 0.5 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OCE/FLO/flowri.F90

    r12377 r12406  
    245245            !------------------------------- 
    246246            irec =  INT( (kt-nn_it000+1)/nn_writefl ) +1 
    247             ztime = ( kt-nn_it000 + 1 ) * rdt 
     247            ztime = ( kt-nn_it000 + 1 ) * rn_Dt 
    248248 
    249249            CALL flioputv( numflo , 'time_counter', ztime , start=(/irec/) ) 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OCE/ICB/icbini.F90

    r12377 r12406  
    6060      !!              - setup either test icebergs or calving file 
    6161      !!---------------------------------------------------------------------- 
    62       REAL(wp), INTENT(in) ::   pdt   ! iceberg time-step (rdt*nn_fsbc) 
     62      REAL(wp), INTENT(in) ::   pdt   ! iceberg time-step (rn_Dt*nn_fsbc) 
    6363      INTEGER , INTENT(in) ::   kt    ! time step number 
    6464      ! 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OCE/ICB/icbtrj.F90

    r10068 r12406  
    7474 
    7575      ! compute end time step date 
    76       zfjulday = fjulday + rdt / rday * REAL( nitend - nit000 + 1 , wp) 
     76      zfjulday = fjulday + rn_Dt / rday * REAL( nitend - nit000 + 1 , wp) 
    7777      IF( ABS(zfjulday - REAL(NINT(zfjulday),wp)) < 0.1 / rday )   zfjulday = REAL(NINT(zfjulday),wp)   ! avoid truncation error 
    7878      CALL ju2ymds( zfjulday, iyear, imonth, iday, zsec ) 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OCE/IOM/iom.F90

    r12377 r12406  
    274274      ! 
    275275      ! set time step length 
    276       dtime%second = rdt 
     276      dtime%second = rn_Dt 
    277277      CALL xios_set_timestep( dtime ) 
    278278      ! 
     
    410410   IF(cdmdl == "OPA") THEN 
    411411!from restart.F90 
    412    CALL iom_set_rstw_var_active("rdt") 
     412   CALL iom_set_rstw_var_active("rn_Dt") 
    413413   IF ( .NOT. ln_diurnal_only ) THEN 
    414414        CALL iom_set_rstw_var_active('ub'  ) 
     
    448448 
    449449        i = 0 
    450         i = i + 1; fields(i)%vname="rdt";            fields(i)%grid="grid_scalar" 
     450        i = i + 1; fields(i)%vname="rn_Dt";            fields(i)%grid="grid_scalar" 
    451451        i = i + 1; fields(i)%vname="un";             fields(i)%grid="grid_N_3D" 
    452452        i = i + 1; fields(i)%vname="ub";             fields(i)%grid="grid_N_3D" 
     
    23582358            idx = INDEX(clname,'@startdate@') + INDEX(clname,'@STARTDATE@') 
    23592359            DO WHILE ( idx /= 0 )  
    2360                cldate = iom_sdate( fjulday - rdt / rday ) 
     2360               cldate = iom_sdate( fjulday - rn_Dt / rday ) 
    23612361               clname = clname(1:idx-1)//TRIM(cldate)//clname(idx+11:LEN_TRIM(clname)) 
    23622362               idx = INDEX(clname,'@startdate@') + INDEX(clname,'@STARTDATE@') 
     
    23652365            idx = INDEX(clname,'@startdatefull@') + INDEX(clname,'@STARTDATEFULL@') 
    23662366            DO WHILE ( idx /= 0 )  
    2367                cldate = iom_sdate( fjulday - rdt / rday, ldfull = .TRUE. ) 
     2367               cldate = iom_sdate( fjulday - rn_Dt / rday, ldfull = .TRUE. ) 
    23682368               clname = clname(1:idx-1)//TRIM(cldate)//clname(idx+15:LEN_TRIM(clname)) 
    23692369               idx = INDEX(clname,'@startdatefull@') + INDEX(clname,'@STARTDATEFULL@') 
     
    23722372            idx = INDEX(clname,'@enddate@') + INDEX(clname,'@ENDDATE@') 
    23732373            DO WHILE ( idx /= 0 )  
    2374                cldate = iom_sdate( fjulday + rdt / rday * REAL( nitend - nit000, wp ), ld24 = .TRUE. ) 
     2374               cldate = iom_sdate( fjulday + rn_Dt / rday * REAL( nitend - nit000, wp ), ld24 = .TRUE. ) 
    23752375               clname = clname(1:idx-1)//TRIM(cldate)//clname(idx+9:LEN_TRIM(clname)) 
    23762376               idx = INDEX(clname,'@enddate@') + INDEX(clname,'@ENDDATE@') 
     
    23792379            idx = INDEX(clname,'@enddatefull@') + INDEX(clname,'@ENDDATEFULL@') 
    23802380            DO WHILE ( idx /= 0 )  
    2381                cldate = iom_sdate( fjulday + rdt / rday * REAL( nitend - nit000, wp ), ld24 = .TRUE., ldfull = .TRUE. ) 
     2381               cldate = iom_sdate( fjulday + rn_Dt / rday * REAL( nitend - nit000, wp ), ld24 = .TRUE., ldfull = .TRUE. ) 
    23822382               clname = clname(1:idx-1)//TRIM(cldate)//clname(idx+13:LEN_TRIM(clname)) 
    23832383               idx = INDEX(clname,'@enddatefull@') + INDEX(clname,'@ENDDATEFULL@') 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OCE/IOM/restart.F90

    r12397 r12406  
    144144      !!---------------------------------------------------------------------- 
    145145                     IF(lwxios) CALL iom_swap(      cwxios_context          ) 
    146                      CALL iom_rstput( kt, nitrst, numrow, 'rdt'    , rdt       , ldxios = lwxios)   ! dynamics time step 
     146                     CALL iom_rstput( kt, nitrst, numrow, 'rdt'    , rn_Dt       , ldxios = lwxios)   ! dynamics time step 
    147147                     CALL iom_delay_rst( 'WRITE', 'OCE', numrow )   ! save only ocean delayed global communication variables 
    148148 
     
    247247      IF( iom_varid( numror, 'rdt', ldstop = .FALSE. ) > 0 )   THEN 
    248248         CALL iom_get( numror, 'rdt', zrdt, ldxios = lrxios ) 
    249          IF( zrdt /= rn_rdt ) THEN 
     249         IF( zrdt /= rn_Dt ) THEN 
    250250            IF(lwp) WRITE( numout,*) 
    251251            IF(lwp) WRITE( numout,*) 'rst_read:  rdt not equal to the read one' 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OCE/ISF/isfcpl.F90

    r12397 r12406  
    502502      ! compute run length 
    503503      nstp_iscpl  = nitend - nit000 + 1 
    504       rdt_iscpl   = nstp_iscpl * rn_rdt 
     504      rdt_iscpl   = nstp_iscpl * rn_Dt 
    505505      z1_rdtiscpl = 1._wp / rdt_iscpl  
    506506 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OCE/ISF/isfdynatf.F90

    r12372 r12406  
    3939      REAL(wp), DIMENSION(jpi,jpj,jpk), INTENT(inout) :: pe3t_f   ! time filtered scale factor to be corrected 
    4040      ! 
    41       REAL(wp)                        , INTENT(in   ) :: pcoef    ! atfp * rdt * r1_rau0 
     41      REAL(wp)                        , INTENT(in   ) :: pcoef    ! atfp * rn_Dt * r1_rau0 
    4242      !!-------------------------------------------------------------------- 
    4343      INTEGER :: jk  ! loop index 
     
    7070      REAL(wp), DIMENSION(jpi,jpj)    , INTENT(in   ) :: pfrac, phtbl    ! fraction of bottom cell included in tbl, tbl thickness 
    7171      REAL(wp), DIMENSION(jpi,jpj)    , INTENT(in   ) :: pfwf , pfwf_b   ! now/before fwf 
    72       REAL(wp),                         INTENT(in   ) :: pcoef           ! atfp * rdt * r1_rau0 
     72      REAL(wp),                         INTENT(in   ) :: pcoef           ! atfp * rn_Dt * r1_rau0 
    7373      !!---------------------------------------------------------------------- 
    7474      INTEGER :: ji,jj,jk 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OCE/LDF/ldfdyn.F90

    r12377 r12406  
    407407            zcmsmag   = (rn_csmc/rpi)**2                                            ! (C_smag/pi)^2 
    408408            zstabf_lo = rn_minfac * rn_minfac / ( 2._wp * 12._wp * 12._wp * zcmsmag ) ! lower limit stability factor scaling 
    409             zstabf_up = rn_maxfac / ( 4._wp * zcmsmag * 2._wp * rdt )               ! upper limit stability factor scaling 
     409            zstabf_up = rn_maxfac / ( 4._wp * zcmsmag * 2._wp * rn_Dt )               ! upper limit stability factor scaling 
    410410            IF( ln_dynldf_blp ) zstabf_lo = ( 16._wp / 9._wp ) * zstabf_lo          ! provide |U|L^3/12 lower limit instead  
    411411            !                                                                       ! of |U|L^3/16 in blp case 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OCE/OBS/diaobs.F90

    r12377 r12406  
    539539      ENDIF 
    540540 
    541       idaystp = NINT( rday / rdt ) 
     541      idaystp = NINT( rday / rn_Dt ) 
    542542 
    543543      !----------------------------------------------------------------------- 
     
    774774         & rday 
    775775      USE dom_oce, ONLY : &           ! Ocean space and time domain variables 
    776          & rdt 
     776         & rn_Dt 
    777777 
    778778      IMPLICIT NONE 
     
    805805      !! Compute number of days + number of hours + min since initial time 
    806806      !!---------------------------------------------------------------------- 
    807       zdayfrc = kstp * rdt / rday 
     807      zdayfrc = kstp * rn_Dt / rday 
    808808      zdayfrc = zdayfrc - aint(zdayfrc) 
    809809      imin = imin + int( zdayfrc * 24 * 60 )  
     
    816816        iday=iday+1 
    817817      END DO  
    818       iday = iday + kstp * rdt / rday  
     818      iday = iday + kstp * rn_Dt / rday  
    819819 
    820820      !----------------------------------------------------------------------- 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OCE/OBS/obs_prep.F90

    r12377 r12406  
    613613      !! * Modules used 
    614614      USE dom_oce, ONLY : &  ! Geographical information 
    615          & rdt 
     615         & rn_Dt 
    616616      USE phycst, ONLY : &   ! Physical constants 
    617617         & rday,  &              
     
    662662 
    663663      ! Intialize the number of time steps per day 
    664       idaystp = NINT( rday / rdt ) 
     664      idaystp = NINT( rday / rn_Dt ) 
    665665 
    666666      !--------------------------------------------------------------------- 
     
    732732 
    733733         ! Add in the number of time steps to the observation minute 
    734          zminstp = rmmss / rdt 
     734         zminstp = rmmss / rn_Dt 
    735735         zhoustp = rhhmm * zminstp 
    736736 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OCE/SBC/fldread.F90

    r12377 r12406  
    172172      ! Note that all varibles starting by nsec_* are shifted time by +1/2 time step to be centrered 
    173173      IF( PRESENT(kit) ) THEN   ! ignore kn_fsbc in this case 
    174          isecsbc = nsec_year + nsec1jan000 + NINT( (     REAL(      kit,wp) + zt_offset ) * rdt / REAL(nn_baro,wp) ) 
     174         isecsbc = nsec_year + nsec1jan000 + NINT( (     REAL(      kit,wp) + zt_offset ) * rn_Dt / REAL(nn_baro,wp) ) 
    175175      ELSE                      ! middle of sbc time step 
    176176         ! note: we use kn_fsbc-1 because nsec_year is defined at the middle of the current time step 
    177          isecsbc = nsec_year + nsec1jan000 + NINT( ( 0.5*REAL(kn_fsbc-1,wp) + zt_offset ) * rdt ) 
     177         isecsbc = nsec_year + nsec1jan000 + NINT( ( 0.5*REAL(kn_fsbc-1,wp) + zt_offset ) * rn_Dt ) 
    178178      ENDIF 
    179179      imf = SIZE( sd ) 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OCE/SBC/sbcblk.F90

    r12377 r12406  
    275275            ENDIF 
    276276            ! 
    277             IF( slf_i(jfpr)%freqh > 0. .AND. MOD( NINT(3600. * slf_i(jfpr)%freqh), nn_fsbc * NINT(rdt) ) /= 0 )   & 
    278                &  CALL ctl_warn( 'sbc_blk_init: sbcmod timestep rdt*nn_fsbc is NOT a submultiple of atmospheric forcing frequency.',   & 
    279                &                 '               This is not ideal. You should consider changing either rdt or nn_fsbc value...' ) 
     277            IF( slf_i(jfpr)%freqh > 0. .AND. MOD( NINT(3600. * slf_i(jfpr)%freqh), nn_fsbc * NINT(rn_Dt) ) /= 0 )   & 
     278               &  CALL ctl_warn( 'sbc_blk_init: sbcmod timestep rn_Dt*nn_fsbc is NOT a submultiple of atmospheric forcing frequency.',   & 
     279               &                 '               This is not ideal. You should consider changing either rn_Dt or nn_fsbc value...' ) 
    280280         ENDIF 
    281281      END DO 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OCE/SBC/sbcblk_skin_coare.F90

    r12377 r12406  
    199199         ! Okay test on updated absorbed flux: 
    200200         !#LB: remove??? has a strong influence !!! 
    201          IF( (.NOT. l_exit).AND.(Qnt_ac(ji,jj) + zQabs*rdt <= 0._wp) ) THEN 
     201         IF( (.NOT. l_exit).AND.(Qnt_ac(ji,jj) + zQabs*rn_Dt <= 0._wp) ) THEN 
    202202            l_exit       = .TRUE. 
    203203            l_destroy_wl = .TRUE. 
     
    211211            ! 2/ Regardless of WL formed (dT==0 or dT>0), we are in the process to initiate one or warm further it ! 
    212212 
    213             ztac = Tau_ac(ji,jj) + MAX(.002_wp , pTau(ji,jj))*rdt      ! updated momentum integral 
     213            ztac = Tau_ac(ji,jj) + MAX(.002_wp , pTau(ji,jj))*rn_Dt      ! updated momentum integral 
    214214            !PRINT *, '#LBD: updated value for Tac=',  REAL(ztac,4) 
    215215 
     
    218218            DO jl = 1, 5 
    219219               zQabs = frac_solar_abs(zHwl)*pQsw(ji,jj) + pQnsol(ji,jj) 
    220                zqac  = Qnt_ac(ji,jj) + zQabs*rdt ! updated heat absorbed 
     220               zqac  = Qnt_ac(ji,jj) + zQabs*rn_Dt ! updated heat absorbed 
    221221               IF( zqac <= 0._wp ) EXIT 
    222222               zHwl = MAX( MIN( Hwl_max , zcd1*ztac/SQRT(zqac)) , 0.1_wp ) ! Warm-layer depth 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OCE/SBC/sbcblk_skin_ecmwf.F90

    r12377 r12406  
    214214         zcst2 = zcst1 / ( 5._wp*zHwl*zusw2 )  !OR: zcst2 = zcst1*rNuwl0 / ( 5._wp*zHwl*zusw2 ) ??? 
    215215 
    216          zcst0 = rdt * (rNuwl0 + 1._wp) / zHwl 
     216         zcst0 = rn_Dt * (rNuwl0 + 1._wp) / zHwl 
    217217 
    218218         ZA = zcst0 * zQabs / ( rNuwl0 * zRhoCp_w ) 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OCE/SBC/sbccpl.F90

    r12377 r12406  
    11231123      !                                                      ! Receive all the atmos. fields (including ice information) 
    11241124      !                                                      ! ======================================================= ! 
    1125       isec = ( kt - nit000 ) * NINT( rdt )                      ! date of exchanges 
     1125      isec = ( kt - nit000 ) * NINT( rn_Dt )                      ! date of exchanges 
    11261126      DO jn = 1, jprcv                                          ! received fields sent by the atmosphere 
    11271127         IF( srcv(jn)%laction )   CALL cpl_rcv( jn, isec, frcv(jn)%z3, xcplmask(:,:,1:nn_cplmodel), nrcvinfo(jn) ) 
     
    21002100      !!---------------------------------------------------------------------- 
    21012101      ! 
    2102       isec = ( kt - nit000 ) * NINT( rdt )        ! date of exchanges 
     2102      isec = ( kt - nit000 ) * NINT( rn_Dt )        ! date of exchanges 
    21032103 
    21042104      zfr_l(:,:) = 1.- fr_i(:,:) 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OCE/SBC/sbcdcy.F90

    r12377 r12406  
    8686      ! -------------- 
    8787      ! When are we during the day (from 0 to 1) 
    88       zlo = ( REAL(nsec_day, wp) - 0.5_wp * rdt ) / rday 
    89       zup = zlo + ( REAL(nn_fsbc, wp)     * rdt ) / rday 
     88      zlo = ( REAL(nsec_day, wp) - 0.5_wp * rn_Dt ) / rday 
     89      zup = zlo + ( REAL(nn_fsbc, wp)     * rn_Dt ) / rday 
    9090      ! 
    9191      IF( nday_qsr == -1 ) THEN       ! first time step only 
     
    251251         END_2D 
    252252         ! 
    253          ztmp = rday / ( rdt * REAL(nn_fsbc, wp) ) 
     253         ztmp = rday / ( rn_Dt * REAL(nn_fsbc, wp) ) 
    254254         rscal(:,:) = rscal(:,:) * ztmp 
    255255         ! 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OCE/SBC/sbcfwb.F90

    r12377 r12406  
    126126         ENDIF    
    127127         !                                         ! Update fwfold if new year start 
    128          ikty = 365 * 86400 / rdt                  !!bug  use of 365 days leap year or 360d year !!!!!!! 
     128         ikty = 365 * 86400 / rn_Dt                  !!bug  use of 365 days leap year or 360d year !!!!!!! 
    129129         IF( MOD( kt, ikty ) == 0 ) THEN 
    130130            a_fwb_b = a_fwb                           ! mean sea level taking into account the ice+snow 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OCE/SBC/sbcmod.F90

    r12377 r12406  
    187187      ! 
    188188      IF( .NOT.ln_usr ) THEN     ! the model calendar needs some specificities (except in user defined case) 
    189          IF( MOD( rday , rdt ) /= 0. )   CALL ctl_stop( 'the time step must devide the number of second of in a day' ) 
     189         IF( MOD( rday , rn_Dt ) /= 0. )   CALL ctl_stop( 'the time step must devide the number of second of in a day' ) 
    190190         IF( MOD( rday , 2.  ) /= 0. )   CALL ctl_stop( 'the number of second of in a day must be an even number'    ) 
    191          IF( MOD( rdt  , 2.  ) /= 0. )   CALL ctl_stop( 'the time step (in second) must be an even number'           ) 
     191         IF( MOD( rn_Dt  , 2.  ) /= 0. )   CALL ctl_stop( 'the time step (in second) must be an even number'           ) 
    192192      ENDIF 
    193193      !                       !**  check option consistency 
     
    309309      !     SAS time-step has to be declared in OASIS (mandatory) -> nn_fsbc has to be modified accordingly 
    310310      IF( nn_components /= jp_iam_nemo ) THEN 
    311          IF( nn_components == jp_iam_opa )   nn_fsbc = cpl_freq('O_SFLX') / NINT(rdt) 
    312          IF( nn_components == jp_iam_sas )   nn_fsbc = cpl_freq('I_SFLX') / NINT(rdt) 
     311         IF( nn_components == jp_iam_opa )   nn_fsbc = cpl_freq('O_SFLX') / NINT(rn_Dt) 
     312         IF( nn_components == jp_iam_sas )   nn_fsbc = cpl_freq('I_SFLX') / NINT(rn_Dt) 
    313313         ! 
    314314         IF(lwp)THEN 
     
    331331      ENDIF 
    332332      ! 
    333       IF( MOD( rday, REAL(nn_fsbc, wp) * rdt ) /= 0 )   & 
     333      IF( MOD( rday, REAL(nn_fsbc, wp) * rn_Dt ) /= 0 )   & 
    334334         &  CALL ctl_warn( 'sbc_init : nn_fsbc is NOT a multiple of the number of time steps in a day' ) 
    335335      ! 
    336       IF( ln_dm2dc .AND. NINT(rday) / ( nn_fsbc * NINT(rdt) ) < 8  )   & 
     336      IF( ln_dm2dc .AND. NINT(rday) / ( nn_fsbc * NINT(rn_Dt) ) < 8  )   & 
    337337         &   CALL ctl_warn( 'sbc_init : diurnal cycle for qsr: the sampling of the diurnal cycle is too small...' ) 
    338338      ! 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OCE/TDE/tide_mod.F90

    r12343 r12406  
    171171      IF( ln_scal_load.AND.ln_read_load ) & 
    172172          &   CALL ctl_stop('Choose between ln_scal_load and ln_read_load') 
    173       IF( ln_tide_ramp.AND.((nitend-nit000+1)*rdt/rday < rn_tide_ramp_dt) )   & 
     173      IF( ln_tide_ramp.AND.((nitend-nit000+1)*rn_Dt/rday < rn_tide_ramp_dt) )   & 
    174174         &   CALL ctl_stop('rn_tide_ramp_dt must be lower than run duration') 
    175175      IF( ln_tide_ramp.AND.(rn_tide_ramp_dt<0.) ) & 
     
    424424      !!---------------------------------------------------------------------- 
    425425       
    426       IF( nsec_day == NINT(0.5_wp * rdt) .OR. kt == nit000 ) THEN      ! start a new day 
     426      IF( nsec_day == NINT(0.5_wp * rn_Dt) .OR. kt == nit000 ) THEN      ! start a new day 
    427427         ! 
    428428         CALL tide_harmo(tide_components, tide_harmonics, ndt05) ! Update oscillation parameters of tidal components for start of current day 
     
    441441         IF( ln_tide_pot )   CALL tide_init_potential 
    442442         ! 
    443          rn_tide_ramp_t = (kt - nit000)*rdt !   Elapsed time in seconds 
     443         rn_tide_ramp_t = (kt - nit000)*rn_Dt !   Elapsed time in seconds 
    444444      ENDIF 
    445445      ! 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OCE/TRA/traatf.F90

    r12397 r12406  
    123123         ENDIF 
    124124         ! total trend for the non-time-filtered variables.  
    125          zfact = 1.0 / rdt 
     125         zfact = 1.0 / rn_Dt 
    126126         ! G Nurser 23 Mar 2017. Recalculate trend as Delta(e3t*T)/e3tn; e3tn cancel from pts(Kmm) terms 
    127127         DO jk = 1, jpkm1 
     
    151151      ELSE                                            ! Leap-Frog + Asselin filter time stepping 
    152152         ! 
    153          IF( ln_linssh ) THEN   ;   CALL tra_atf_fix( kt, Kbb, Kmm, Kaa, nit000,      'TRA', pts, jpts )  ! linear free surface  
    154          ELSE                   ;   CALL tra_atf_vvl( kt, Kbb, Kmm, Kaa, nit000, rdt, 'TRA', pts, sbc_tsc, sbc_tsc_b, jpts )  ! non-linear free surface 
     153         IF( ln_linssh ) THEN   ;   CALL tra_atf_fix( kt, Kbb, Kmm, Kaa, nit000,        'TRA', pts, jpts )  ! linear free surface  
     154         ELSE                   ;   CALL tra_atf_vvl( kt, Kbb, Kmm, Kaa, nit000, rn_Dt, 'TRA', pts, sbc_tsc, sbc_tsc_b, jpts )  ! non-linear free surface 
    155155         ENDIF 
    156156         ! 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OCE/TRA/tranpc.F90

    r12377 r12406  
    301301         ! 
    302302         IF( l_trdtra ) THEN         ! send the Non penetrative mixing trends for diagnostic 
    303             z1_r2dt = 1._wp / (2._wp * rdt) 
     303            z1_r2dt = 1._wp / (2._wp * rn_Dt) 
    304304            ztrdt(:,:,:) = ( pts(:,:,:,jp_tem,Kaa) - ztrdt(:,:,:) ) * z1_r2dt 
    305305            ztrds(:,:,:) = ( pts(:,:,:,jp_sal,Kaa) - ztrds(:,:,:) ) * z1_r2dt 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OCE/TRD/trdvor.F90

    r12377 r12406  
    385385         ! III.1 compute total trend 
    386386         ! ------------------------ 
    387          zmean = 1._wp / (  REAL( nmoydpvor, wp ) * 2._wp * rdt  ) 
     387         zmean = 1._wp / (  REAL( nmoydpvor, wp ) * 2._wp * rn_Dt  ) 
    388388         vor_avrtot(:,:) = (  vor_avr(:,:) - vor_avrbn(:,:) + vor_avrb(:,:) - vor_avrbb(:,:) ) * zmean 
    389389 
     
    504504      ENDIF 
    505505#if defined key_diainstant 
    506       zsto = nn_write*rdt 
     506      zsto = nn_write*rn_Dt 
    507507      clop = "inst("//TRIM(clop)//")" 
    508508#else 
    509       zsto = rdt 
     509      zsto = rn_Dt 
    510510      clop = "ave("//TRIM(clop)//")" 
    511511#endif 
    512       zout = nn_trd*rdt 
     512      zout = nn_trd*rn_Dt 
    513513 
    514514      IF(lwp) WRITE(numout,*) '               netCDF initialization' 
     
    516516      ! II.2 Compute julian date from starting date of the run 
    517517      ! ------------------------ 
    518       CALL ymds2ju( nyear, nmonth, nday, rdt, zjulian ) 
     518      CALL ymds2ju( nyear, nmonth, nday, rn_Dt, zjulian ) 
    519519      zjulian = zjulian - adatrj   !   set calendar origin to the beginning of the experiment 
    520520      IF(lwp) WRITE(numout,*)' '   
     
    528528      IF(lwp) WRITE(numout,*) ' Name of NETCDF file ', clhstnam 
    529529      CALL histbeg( clhstnam, jpi, glamf, jpj, gphif,1, jpi,   &  ! Horizontal grid : glamt and gphit 
    530          &          1, jpj, nit000-1, zjulian, rdt, nh_t, nidvor, domain_id=nidom, snc4chunks=snc4set ) 
     530         &          1, jpj, nit000-1, zjulian, rn_Dt, nh_t, nidvor, domain_id=nidom, snc4chunks=snc4set ) 
    531531      CALL wheneq( jpi*jpj, fmask, 1, 1., ndexvor1, ndimvor1 )    ! surface 
    532532 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OCE/USR/usrdef_hgr.F90

    r12377 r12406  
    107107      IF( ln_bench ) THEN     ! benchmark: forced the resolution to be 106 km  
    108108         ze1 = 106000._wp     ! but keep (lat,lon) at the right nn_GYRE resolution 
    109          CALL ctl_warn( ' GYRE used as Benchmark: e1=e2=106km, no need to adjust rdt, ahm,aht ' ) 
     109         CALL ctl_warn( ' GYRE used as Benchmark: e1=e2=106km, no need to adjust rn_Dt, ahm,aht ' ) 
    110110      ENDIF 
    111111      IF( nprint==1 .AND. lwp )   THEN 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OCE/USR/usrdef_sbc.F90

    r12377 r12406  
    8989 
    9090      ! current day (in hours) since january the 1st of the current year 
    91       ztime = REAL( kt ) * rdt / (rmmss * rhhmm)   &       !  total incrementation (in hours) 
     91      ztime = REAL( kt ) * rn_Dt / (rmmss * rhhmm)   &       !  total incrementation (in hours) 
    9292         &      - (nyear  - 1) * rjjhh * zyydd             !  minus years since beginning of experiment (in hours) 
    9393 
     
    154154      !accumulates days of previous months of this year 
    155155      ! day (in hours) since january the 1st 
    156       ztime = FLOAT( kt ) * rdt / (rmmss * rhhmm)  &  ! incrementation in hour 
     156      ztime = FLOAT( kt ) * rn_Dt / (rmmss * rhhmm)  &  ! incrementation in hour 
    157157         &     - (nyear - 1) * rjjhh * zyydd          !  - nber of hours the precedent years 
    158158      ztimemax = ((5.*30.)+21.)* 24.               ! 21th june     in hours 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OCE/ZDF/zdfdrg.F90

    r12377 r12406  
    165165      !!--------------------------------------------------------------------- 
    166166      ! 
    167 !!gm bug : time step is only rdt (not 2 rdt if euler start !) 
    168       zm1_2dt = - 1._wp / ( 2._wp * rdt ) 
     167!!gm bug : time step is only rn_Dt (not 2 rn_Dt if euler start !) 
     168      zm1_2dt = - 1._wp / ( 2._wp * rn_Dt ) 
    169169 
    170170      IF( l_trddyn ) THEN      ! trends: store the input trends 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OCE/ZDF/zdfgls.F90

    r12377 r12406  
    267267         zd_up(ji,jj,jk) = zcof * ( p_avm(ji,jj,jk+1) + p_avm(ji,jj,jk  ) ) / ( e3t(ji,jj,jk  ,Kmm) * e3w(ji,jj,jk,Kmm) ) 
    268268         !                                        ! diagonal 
    269          zdiag(ji,jj,jk) = 1._wp - zd_lw(ji,jj,jk) - zd_up(ji,jj,jk)  + rdt * zdiss * wmask(ji,jj,jk)  
     269         zdiag(ji,jj,jk) = 1._wp - zd_lw(ji,jj,jk) - zd_up(ji,jj,jk)  + rn_Dt * zdiss * wmask(ji,jj,jk)  
    270270         !                                        ! right hand side in en 
    271          en(ji,jj,jk) = en(ji,jj,jk) + rdt * zesh2 * wmask(ji,jj,jk) 
     271         en(ji,jj,jk) = en(ji,jj,jk) + rn_Dt * zesh2 * wmask(ji,jj,jk) 
    272272      END_3D 
    273273      ! 
     
    477477         zd_up(ji,jj,jk) = zcof * ( p_avm(ji,jj,jk+1) + p_avm(ji,jj,jk  ) ) / ( e3t(ji,jj,jk  ,Kmm) * e3w(ji,jj,jk,Kmm) ) 
    478478         !                                               ! diagonal 
    479          zdiag(ji,jj,jk) = 1._wp - zd_lw(ji,jj,jk) - zd_up(ji,jj,jk) + rdt * zdiss * wmask(ji,jj,jk) 
     479         zdiag(ji,jj,jk) = 1._wp - zd_lw(ji,jj,jk) - zd_up(ji,jj,jk) + rn_Dt * zdiss * wmask(ji,jj,jk) 
    480480         !                                               ! right hand side in psi 
    481          psi(ji,jj,jk) = psi(ji,jj,jk) + rdt * zesh2 * wmask(ji,jj,jk) 
     481         psi(ji,jj,jk) = psi(ji,jj,jk) + rn_Dt * zesh2 * wmask(ji,jj,jk) 
    482482      END_3D 
    483483      ! 
     
    10071007      rc04  = rc03 * rc0 
    10081008      rsbc_tke1 = -3._wp/2._wp*rn_crban*ra_sf*rl_sf                      ! Dirichlet + Wave breaking 
    1009       rsbc_tke2 = rdt * rn_crban / rl_sf                                 ! Neumann + Wave breaking  
     1009      rsbc_tke2 = rn_Dt * rn_crban / rl_sf                                 ! Neumann + Wave breaking  
    10101010      zcr = MAX(rsmall, rsbc_tke1**(1./(-ra_sf*3._wp/2._wp))-1._wp ) 
    10111011      rtrans = 0.2_wp / zcr                                              ! Ad. inverse transition length between log and wave layer  
    10121012      rsbc_zs1  = rn_charn/grav                                          ! Charnock formula for surface roughness 
    10131013      rsbc_zs2  = rn_frac_hs / 0.85_wp / grav * 665._wp                  ! Rascle formula for surface roughness  
    1014       rsbc_psi1 = -0.5_wp * rdt * rc0**(rpp-2._wp*rmm) / rsc_psi 
    1015       rsbc_psi2 = -0.5_wp * rdt * rc0**rpp * rnn * vkarmn**rnn / rsc_psi ! Neumann + NO Wave breaking  
    1016       ! 
    1017       rfact_tke = -0.5_wp / rsc_tke * rdt                                ! Cst used for the Diffusion term of tke 
    1018       rfact_psi = -0.5_wp / rsc_psi * rdt                                ! Cst used for the Diffusion term of tke 
     1014      rsbc_psi1 = -0.5_wp * rn_Dt * rc0**(rpp-2._wp*rmm) / rsc_psi 
     1015      rsbc_psi2 = -0.5_wp * rn_Dt * rc0**rpp * rnn * vkarmn**rnn / rsc_psi ! Neumann + NO Wave breaking  
     1016      ! 
     1017      rfact_tke = -0.5_wp / rsc_tke * rn_Dt                                ! Cst used for the Diffusion term of tke 
     1018      rfact_psi = -0.5_wp / rsc_psi * rn_Dt                                ! Cst used for the Diffusion term of tke 
    10191019      ! 
    10201020      !                                !* Wall proximity function 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OCE/ZDF/zdfosm.F90

    r12377 r12406  
    441441                        &            + 0.135 * zla(ji,jj) * zwstrl(ji,jj)**3/hbl(ji,jj) ) 
    442442 
    443                    zvel_max =  - ( 1.0 + 1.0 * ( zwstrl(ji,jj)**3 + 0.5 * zwstrc(ji,jj)**3 )**pthird * rn_rdt / hbl(ji,jj) ) & 
     443                   zvel_max =  - ( 1.0 + 1.0 * ( zwstrl(ji,jj)**3 + 0.5 * zwstrc(ji,jj)**3 )**pthird * rn_Dt / hbl(ji,jj) ) & 
    444444                        &   * zwb_ent(ji,jj) / ( zwstrl(ji,jj)**3 + 0.5 * zwstrc(ji,jj)**3 )**pthird 
    445445! Entrainment including component due to shear turbulence. Modified Langmuir component, but gives same result for La=0.3 For testing uncomment. 
     
    447447!                           &            + ( 0.15 * ( 1.0 - EXP( -0.5 * zla(ji,jj) ) ) + 0.03 / zla(ji,jj)**2 ) * zustar(ji,jj)**3/hbl(ji,jj) ) 
    448448 
    449 !                      zvel_max = - ( 1.0 + 1.0 * ( zvstr(ji,jj)**3 + 0.5 * zwstrc(ji,jj)**3 )**pthird * rn_rdt / zhbl(ji,jj) ) * zwb_ent(ji,jj) / & 
     449!                      zvel_max = - ( 1.0 + 1.0 * ( zvstr(ji,jj)**3 + 0.5 * zwstrc(ji,jj)**3 )**pthird * rn_Dt / zhbl(ji,jj) ) * zwb_ent(ji,jj) / & 
    450450!                           &       ( zvstr(ji,jj)**3 + 0.5 * zwstrc(ji,jj)**3 )**pthird 
    451451                   zzdhdt = - zwb_ent(ji,jj) / ( zvel_max + MAX(zdb_bl(ji,jj),0.0) ) 
     
    458458                   IF ( zzdhdt < 0._wp ) THEN 
    459459                   ! For long timsteps factor in brackets slows the rapid collapse of the OSBL 
    460                       zpert   = 2.0 * ( 1.0 + 2.0 * zwstrl(ji,jj) * rn_rdt / hbl(ji,jj) ) * zwstrl(ji,jj)**2 / hbl(ji,jj) 
     460                      zpert   = 2.0 * ( 1.0 + 2.0 * zwstrl(ji,jj) * rn_Dt / hbl(ji,jj) ) * zwstrl(ji,jj)**2 / hbl(ji,jj) 
    461461                   ELSE 
    462                       zpert   = 2.0 * ( 1.0 + 2.0 * zwstrl(ji,jj) * rn_rdt / hbl(ji,jj) ) * zwstrl(ji,jj)**2 / hbl(ji,jj) & 
     462                      zpert   = 2.0 * ( 1.0 + 2.0 * zwstrl(ji,jj) * rn_Dt / hbl(ji,jj) ) * zwstrl(ji,jj)**2 / hbl(ji,jj) & 
    463463                           &  + MAX( zdb_bl(ji,jj), 0.0 ) 
    464464                   ENDIF 
     
    472472      ibld(:,:) = 3 
    473473 
    474       zhbl_t(:,:) = hbl(:,:) + (zdhdt(:,:) - ww(ji,jj,ibld(ji,jj)))* rn_rdt ! certainly need wb here, so subtract it 
     474      zhbl_t(:,:) = hbl(:,:) + (zdhdt(:,:) - ww(ji,jj,ibld(ji,jj)))* rn_Dt ! certainly need wb here, so subtract it 
    475475      zhbl_t(:,:) = MIN(zhbl_t(:,:), ht(:,:)) 
    476       zdhdt(:,:) = MIN(zdhdt(:,:), (zhbl_t(:,:) - hbl(:,:))/rn_rdt + ww(ji,jj,ibld(ji,jj))) ! adjustment to represent limiting by ocean bottom 
     476      zdhdt(:,:) = MIN(zdhdt(:,:), (zhbl_t(:,:) - hbl(:,:))/rn_Dt + ww(ji,jj,ibld(ji,jj))) ! adjustment to represent limiting by ocean bottom 
    477477 
    478478      DO_3D_00_00( 4, jpkm1 ) 
     
    496496            IF ( lconv(ji,jj) ) THEN 
    497497!unstable 
    498                zvel_max =  - ( 1.0 + 1.0 * ( zvstr(ji,jj)**3 + 0.5 * zwstrc(ji,jj)**3 )**pthird * rn_rdt / hbl(ji,jj) ) & 
     498               zvel_max =  - ( 1.0 + 1.0 * ( zvstr(ji,jj)**3 + 0.5 * zwstrc(ji,jj)**3 )**pthird * rn_Dt / hbl(ji,jj) ) & 
    499499                    &   * zwb_ent(ji,jj) / ( zvstr(ji,jj)**3 + 0.5 * zwstrc(ji,jj)**3 )**pthird 
    500500 
     
    503503                       & - zbeta * ( zs_bl(ji,jj) - ts(ji,jj,jm,jp_sal,Kmm) ) ), 0.0 ) + zvel_max 
    504504 
    505                   zhbl_s = zhbl_s + MIN( - zwb_ent(ji,jj) / zdb * rn_rdt / FLOAT(ibld(ji,jj)-imld(ji,jj) ), e3w(ji,jj,jk,Kmm) ) 
     505                  zhbl_s = zhbl_s + MIN( - zwb_ent(ji,jj) / zdb * rn_Dt / FLOAT(ibld(ji,jj)-imld(ji,jj) ), e3w(ji,jj,jk,Kmm) ) 
    506506                  zhbl_s = MIN(zhbl_s, ht(ji,jj)) 
    507507 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OCE/ZDF/zdftke.F90

    r12377 r12406  
    207207      ! 
    208208      zbbrau = rn_ebb / rau0       ! Local constant initialisation 
    209       zfact1 = -.5_wp * rdt  
    210       zfact2 = 1.5_wp * rdt * rn_ediss 
     209      zfact1 = -.5_wp * rn_Dt  
     210      zfact2 = 1.5_wp * rn_Dt * rn_ediss 
    211211      zfact3 = 0.5_wp       * rn_ediss 
    212212      ! 
     
    288288                  zwlc = rn_lc * SIN( rpi * gdepw(ji,jj,jk,Kmm) / zhlc(ji,jj) )   ! warning: optimization: zus^3 is in zfr_i 
    289289                  !                                           ! TKE Langmuir circulation source term 
    290                   en(ji,jj,jk) = en(ji,jj,jk) + rdt * zfr_i(ji,jj) * ( zwlc * zwlc * zwlc ) / zhlc(ji,jj) 
     290                  en(ji,jj,jk) = en(ji,jj,jk) + rn_Dt * zfr_i(ji,jj) * ( zwlc * zwlc * zwlc ) / zhlc(ji,jj) 
    291291               ENDIF 
    292292            ENDIF 
     
    325325         ! 
    326326         !                                   ! right hand side in en 
    327          en(ji,jj,jk) = en(ji,jj,jk) + rdt * (  p_sh2(ji,jj,jk)                          &   ! shear 
     327         en(ji,jj,jk) = en(ji,jj,jk) + rn_Dt * (  p_sh2(ji,jj,jk)                        &   ! shear 
    328328            &                                 - p_avt(ji,jj,jk) * rn2(ji,jj,jk)          &   ! stratification 
    329329            &                                 + zfact3 * dissl(ji,jj,jk) * en(ji,jj,jk)  &   ! dissipation 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OCE/nemogcm.F90

    r12397 r12406  
    482482 
    483483      !                                      ! Icebergs 
    484                            CALL icb_init( rdt, nit000)   ! initialise icebergs instance 
     484                           CALL icb_init( rn_Dt, nit000)   ! initialise icebergs instance 
    485485 
    486486                                                ! ice shelf 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OCE/step.F90

    r12397 r12406  
    105105      !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< 
    106106      ! 
    107       IF( l_1st_euler ) THEN   ;   r2dt =         rn_rdt   ! start or restart with Euler 1st time-step 
    108       ELSE                     ;   r2dt = 2._wp * rn_rdt   ! restart with leapfrog 
     107      IF( l_1st_euler ) THEN   ;   r2dt =         rn_Dt   ! start or restart with Euler 1st time-step 
     108      ELSE                     ;   r2dt = 2._wp * rn_Dt   ! restart with leapfrog 
    109109      ENDIF 
    110110      r1_2dt = 1._wp / r2dt 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OCE/timing.F90

    r11536 r12406  
    390390            WRITE(numtime,'(A28,F11.6,A2, F4.1,A3,A25,I8)') 'Waiting  global time : ',timing_glob(4*ji-1)   & 
    391391               &                                                         , ' (',      zperc,' %)',   ' on MPI rank : ', ji 
    392             zsypd = rn_rdt * REAL(nitend-nit000-1, wp) / (timing_glob(4*ji) * 365.) 
     392            zsypd = rn_Dt * REAL(nitend-nit000-1, wp) / (timing_glob(4*ji) * 365.) 
    393393            WRITE(numtime,'(A28,F11.6,A7,F10.3,A2,A15,I8)') 'Total           time : ',timing_glob(4*ji  )   & 
    394394               &                                                         , ' (SYPD: ', zsypd, ')',   ' on MPI rank : ', ji 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OFF/dtadyn.F90

    r12377 r12406  
    622622      !!---------------------------------------------------------------------- 
    623623      ! 
    624       z2dt = 2._wp * rdt 
     624      z2dt = 2._wp * rn_Dt 
    625625      ! 
    626626      zhdiv(:,:) = 0._wp 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/SAS/diawri.F90

    r12377 r12406  
    174174      ENDIF 
    175175#if defined key_diainstant 
    176       zsto = nn_write * rdt 
     176      zsto = nn_write * rn_Dt 
    177177      clop = "inst("//TRIM(clop)//")" 
    178178#else 
    179       zsto=rdt 
     179      zsto=rn_Dt 
    180180      clop = "ave("//TRIM(clop)//")" 
    181181#endif 
    182       zout = nn_write * rdt 
    183       zmax = ( nitend - nit000 + 1 ) * rdt 
     182      zout = nn_write * rn_Dt 
     183      zmax = ( nitend - nit000 + 1 ) * rn_Dt 
    184184 
    185185      ! Define indices of the horizontal output zoom and vertical limit storage 
     
    202202 
    203203         ! Compute julian date from starting date of the run 
    204          CALL ymds2ju( nyear, nmonth, nday, rdt, zjulian ) 
     204         CALL ymds2ju( nyear, nmonth, nday, rn_Dt, zjulian ) 
    205205         zjulian = zjulian - adatrj   !   set calendar origin to the beginning of the experiment 
    206206         IF(lwp)WRITE(numout,*) 
     
    224224         CALL histbeg( clhstnam, jpi, glamt, jpj, gphit,           &  ! Horizontal grid: glamt and gphit 
    225225            &          iimi, iima-iimi+1, ijmi, ijma-ijmi+1,       & 
    226             &          nit000-1, zjulian, rdt, nh_T, nid_T, domain_id=nidom, snc4chunks=snc4set ) 
     226            &          nit000-1, zjulian, rn_Dt, nh_T, nid_T, domain_id=nidom, snc4chunks=snc4set ) 
    227227         CALL histvert( nid_T, "deptht", "Vertical T levels",      &  ! Vertical grid: gdept 
    228228            &           "m", ipk, gdept_1d, nz_T, "down" ) 
     
    236236         CALL histbeg( clhstnam, jpi, glamu, jpj, gphiu,           &  ! Horizontal grid: glamu and gphiu 
    237237            &          iimi, iima-iimi+1, ijmi, ijma-ijmi+1,       & 
    238             &          nit000-1, zjulian, rdt, nh_U, nid_U, domain_id=nidom, snc4chunks=snc4set ) 
     238            &          nit000-1, zjulian, rn_Dt, nh_U, nid_U, domain_id=nidom, snc4chunks=snc4set ) 
    239239         CALL histvert( nid_U, "depthu", "Vertical U levels",      &  ! Vertical grid: gdept 
    240240            &           "m", ipk, gdept_1d, nz_U, "down" ) 
     
    248248         CALL histbeg( clhstnam, jpi, glamv, jpj, gphiv,           &  ! Horizontal grid: glamv and gphiv 
    249249            &          iimi, iima-iimi+1, ijmi, ijma-ijmi+1,       & 
    250             &          nit000-1, zjulian, rdt, nh_V, nid_V, domain_id=nidom, snc4chunks=snc4set ) 
     250            &          nit000-1, zjulian, rn_Dt, nh_V, nid_V, domain_id=nidom, snc4chunks=snc4set ) 
    251251         CALL histvert( nid_V, "depthv", "Vertical V levels",      &  ! Vertical grid : gdept 
    252252            &          "m", ipk, gdept_1d, nz_V, "down" ) 
     
    261261            CALL histbeg( clhstnam, jpi, glamt, jpj, gphit,           &  ! Horizontal grid: glamt and gphit 
    262262               &          iimi, iima-iimi+1, ijmi, ijma-ijmi+1,       & 
    263                &          nit000-1, zjulian, rdt, nh_A, nid_A, domain_id=nidom, snc4chunks=snc4set ) 
     263               &          nit000-1, zjulian, rn_Dt, nh_A, nid_A, domain_id=nidom, snc4chunks=snc4set ) 
    264264            CALL histvert( nid_A, "ght_abl", "Vertical T levels",      &  ! Vertical grid: gdept 
    265265               &           "m", ipka, ght_abl(2:jpka), nz_A, "up" ) 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/SAS/nemogcm.F90

    r12377 r12406  
    378378                           CALL bdy_init 
    379379      ! ==> 
    380                            CALL icb_init( rdt, nit000)   ! initialise icebergs instance 
     380                           CALL icb_init( rn_Dt, nit000)   ! initialise icebergs instance 
    381381      ! 
    382382      IF(lwp) WRITE(numout,cform_aaa)           ! Flag AAAAAAA 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/TOP/C14/trcatm_c14.F90

    r12377 r12406  
    223223      IF(kc14typ >= 1) THEN  ! Transient C14 & CO2 
    224224      ! 
    225          tyrc14_now = tyrc14_now + ( rdt / ( rday * nyear_len(1)) )    !  current time step in yr relative to tyrc14_beg 
     225         tyrc14_now = tyrc14_now + ( rn_Dt / ( rday * nyear_len(1)) )    !  current time step in yr relative to tyrc14_beg 
    226226      ! 
    227227      ! CO2 -------------------------------------------------------- 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/TOP/CFC/trcsms_cfc.F90

    r12377 r12406  
    166166 
    167167            ! cumulation of surface flux at each time step 
    168             qint_cfc(ji,jj,jl) = qint_cfc(ji,jj,jl) + qtr_cfc(ji,jj,jl) * rdt 
     168            qint_cfc(ji,jj,jl) = qint_cfc(ji,jj,jl) + qtr_cfc(ji,jj,jl) * rn_Dt 
    169169            !                                               !----------------! 
    170170         END_2D 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/TOP/PISCES/P2Z/p2zexp.F90

    r12397 r12406  
    9898         zwork = vsed * tr(ji,jj,ikt,jpdet,Kmm) 
    9999         zsedpoca(ji,jj) = ( zwork + dminl(ji,jj) * xksi(ji,jj)   & 
    100             &           - sedlam * sedpocn(ji,jj) - sedlostpoc * sedpocn(ji,jj) ) * rdt 
     100            &           - sedlam * sedpocn(ji,jj) - sedlostpoc * sedpocn(ji,jj) ) * rn_Dt 
    101101         zgeolpoc = zgeolpoc + sedlostpoc * sedpocn(ji,jj) * e1e2t(ji,jj) 
    102102      END_2D 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/TOP/PISCES/P4Z/p4zsms.F90

    r12397 r12406  
    106106         xfact = 1.e+3 * rfact2r 
    107107         IF(lwp) WRITE(numout,*)  
    108          IF(lwp) WRITE(numout,*) '    Passive Tracer  time step    rfact  = ', rfact, ' rdt = ', rdt 
     108         IF(lwp) WRITE(numout,*) '    Passive Tracer  time step    rfact  = ', rfact, ' rn_Dt = ', rn_Dt 
    109109         IF(lwp) write(numout,*) '    PISCES  Biology time step    rfact2 = ', rfact2 
    110110         IF(lwp) WRITE(numout,*) 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/TOP/PISCES/SED/oce_sed.F90

    r12377 r12406  
    1818   USE dom_oce , ONLY :   mbkt      =>   mbkt           !: vertical index of the bottom last T- ocean level 
    1919   USE dom_oce , ONLY :   tmask     =>   tmask          !: land/ocean mask at t-points 
    20    USE dom_oce , ONLY :   rdt       =>   rdt            !: time step for the dynamics 
     20   USE dom_oce , ONLY :   rn_Dt     =>   rn_Dt          !: time step for the dynamics 
    2121   USE dom_oce , ONLY :   nyear     =>   nyear          !: Current year 
    2222   USE dom_oce , ONLY :   ndastp    =>   ndastp         !: time step date in year/month/day aammjj 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/TOP/PISCES/SED/sedrst.F90

    r12377 r12406  
    330330      !!       In both those options, the  exact duration of the experiment 
    331331      !!       since the beginning (cumulated duration of all previous restart runs) 
    332       !!       is not stored in the restart and is assumed to be (nittrc000-1)*rdt. 
     332      !!       is not stored in the restart and is assumed to be (nittrc000-1)*rn_Dt. 
    333333      !!       This is valid is the time step has remained constant. 
    334334      !! 
     
    381381             ELSE 
    382382               ndastp = ndate0 - 1     ! ndate0 read in the namelist in dom_nam 
    383                adatrj = ( REAL( nittrc000-1, wp ) * rdt ) / rday 
     383               adatrj = ( REAL( nittrc000-1, wp ) * rn_Dt ) / rday 
    384384               ! note this is wrong if time step has changed during run 
    385385            ENDIF 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/TOP/TRP/trcrad.F90

    r12377 r12406  
    147147     ! 
    148148     IF( l_trdtrc )   ALLOCATE( ztrtrd(jpi,jpj,jpk) ) 
    149      zs2rdt = 1. / ( 2. * rdt ) 
     149     zs2rdt = 1. / ( 2. * rn_Dt ) 
    150150     ! 
    151151     DO jt = 1,2  ! Loop over time indices since exactly the same fix is applied to "now" and "after" fields 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/TOP/TRP/trdmxl_trc.F90

    r12377 r12406  
    401401         DO jn = 1, jptra 
    402402            IF( ln_trdtrc(jn) ) THEN 
    403                !-- Compute total trends    (use rdttrc instead of rdt ???) 
     403               !-- Compute total trends    (use rdttrc instead of rn_Dt ???) 
    404404               IF ( ln_trcadv_muscl .OR. ln_trcadv_muscl2 ) THEN  ! EULER-FORWARD schemes 
    405                   ztmltot(:,:,jn) =  ( tml_trc(:,:,jn) - tmlbn_trc(:,:,jn) )/rdt 
     405                  ztmltot(:,:,jn) =  ( tml_trc(:,:,jn) - tmlbn_trc(:,:,jn) )/rn_Dt 
    406406               ELSE                                                                     ! LEAP-FROG schemes 
    407                   ztmltot(:,:,jn) =  ( tml_trc(:,:,jn) - tmlbn_trc(:,:,jn) + tmlb_trc(:,:,jn) - tmlbb_trc(:,:,jn))/(2.*rdt) 
     407                  ztmltot(:,:,jn) =  ( tml_trc(:,:,jn) - tmlbn_trc(:,:,jn) + tmlb_trc(:,:,jn) - tmlbb_trc(:,:,jn))/(2.*rn_Dt) 
    408408               ENDIF 
    409409                
     
    439439            IF( ln_trdtrc(jn) ) THEN 
    440440               tml_sum_trc(:,:,jn) = tmlbn_trc(:,:,jn) + 2 * ( tml_sum_trc(:,:,jn) - tml_trc(:,:,jn) ) + tml_trc(:,:,jn) 
    441                ztmltot2   (:,:,jn) = ( tml_sum_trc(:,:,jn) - tml_sumb_trc(:,:,jn) ) /  ( 2.*rdt )    ! now tracer unit is /sec 
     441               ztmltot2   (:,:,jn) = ( tml_sum_trc(:,:,jn) - tml_sumb_trc(:,:,jn) ) /  ( 2.*rn_Dt )    ! now tracer unit is /sec 
    442442            ENDIF 
    443443         END DO 
     
    852852         CALL ctl_stop( 'STOP', 'trd_mxl_trc : this was never checked. Comment this line to proceed...' ) 
    853853      ENDIF 
    854       zsto = nn_trd_trc * rdt 
     854      zsto = nn_trd_trc * rn_Dt 
    855855      clop = "inst("//TRIM(clop)//")" 
    856856#  else 
    857857      IF( ln_trdmxl_trc_instant ) THEN 
    858          zsto = rdt                                               ! inst. diags : we use IOIPSL time averaging 
     858         zsto = rn_Dt                                               ! inst. diags : we use IOIPSL time averaging 
    859859      ELSE 
    860          zsto = nn_trd_trc * rdt                                    ! mean  diags : we DO NOT use any IOIPSL time averaging 
     860         zsto = nn_trd_trc * rn_Dt                                    ! mean  diags : we DO NOT use any IOIPSL time averaging 
    861861      ENDIF 
    862862      clop = "ave("//TRIM(clop)//")" 
    863863#  endif 
    864       zout = nn_trd_trc * rdt 
     864      zout = nn_trd_trc * rn_Dt 
    865865      iiter = nittrc000 - 1 
    866866 
     
    869869      ! II.2 Compute julian date from starting date of the run 
    870870      ! ------------------------------------------------------ 
    871       CALL ymds2ju( nyear, nmonth, nday, rdt, zjulian ) 
     871      CALL ymds2ju( nyear, nmonth, nday, rn_Dt, zjulian ) 
    872872      zjulian = zjulian - adatrj   !   set calendar origin to the beginning of the experiment 
    873873      IF(lwp) WRITE(numout,*)' '   
     
    901901            CALL dia_nam( clhstnam, nn_trd_trc, csuff ) 
    902902            CALL histbeg( clhstnam, jpi, glamt, jpj, gphit,                                            & 
    903                &        1, jpi, 1, jpj, iiter, zjulian, rdt, nh_t(jn), nidtrd(jn), domain_id=nidom, snc4chunks=snc4set ) 
     903               &        1, jpi, 1, jpj, iiter, zjulian, rn_Dt, nh_t(jn), nidtrd(jn), domain_id=nidom, snc4chunks=snc4set ) 
    904904       
    905905            !-- Define the ML depth variable 
     
    938938               CALL histdef(nidtrd(jn), trim(clvar)//trim(ctrd_trc(jl,2)), clmxl//" "//clvar//ctrd_trc(jl,1),                      &  
    939939                 &    cltrcu, jpi, jpj, nh_t(jn), 1  , 1, 1  , -99 , 32, clop, zsto, zout ) ! IOIPSL: time mean 
    940             END DO                                                                         ! if zsto=rdt above 
     940            END DO                                                                         ! if zsto=rn_Dt above 
    941941          
    942942            CALL histdef(nidtrd(jn), trim(clvar)//trim(ctrd_trc(jpmxl_trc_radb,2)), clmxl//" "//clvar//ctrd_trc(jpmxl_trc_radb,1), &  
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/TOP/trcnam.F90

    r12377 r12406  
    7676      ENDIF 
    7777      ! 
    78       rdttrc = rdt                              ! passive tracer time-step       
     78      rdttrc = rn_Dt                              ! passive tracer time-step       
    7979      !  
    8080      IF(lwp) THEN                              ! control print 
    8181        WRITE(numout,*)  
    82         WRITE(numout,*) '   ==>>>   Passive Tracer  time step    rdttrc = rdt = ', rdttrc 
     82        WRITE(numout,*) '   ==>>>   Passive Tracer  time step    rdttrc = rn_Dt = ', rdttrc 
    8383      ENDIF 
    8484      ! 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/TOP/trcrst.F90

    r12397 r12406  
    183183      !!       In both those options, the  exact duration of the experiment 
    184184      !!       since the beginning (cumulated duration of all previous restart runs) 
    185       !!       is not stored in the restart and is assumed to be (nittrc000-1)*rdt. 
     185      !!       is not stored in the restart and is assumed to be (nittrc000-1)*rn_Dt. 
    186186      !!       This is valid is the time step has remained constant. 
    187187      !! 
     
    263263               nminute = ( nn_time0 - nhour * 100 ) 
    264264               IF( nhour*3600+nminute*60-ndt05 .lt. 0 )  ndastp=ndastp-1      ! Start hour is specified in the namelist (default 0) 
    265                adatrj = ( REAL( nit000-1, wp ) * rdt ) / rday 
     265               adatrj = ( REAL( nit000-1, wp ) * rn_Dt ) / rday 
    266266               ! note this is wrong if time step has changed during run 
    267267            ENDIF 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/tests/BENCH/EXPREF/namelist_cfg_orca025_like

    r11536 r12406  
    3636&namdom        !   time and space domain 
    3737!----------------------------------------------------------------------- 
    38    rn_rdt      =    900.   !  time step for the dynamics (and tracer if nn_acc=0) 
     38   rn_Dt      =    900.   !  time step for the dynamics (and tracer if nn_acc=0) 
    3939   rn_atfp     =   0.05    !  asselin time filter parameter 
    4040/ 
     
    189189   ln_dynspg_ts   = .true.   ! split-explicit free surface 
    190190      ln_bt_auto    = .false.    ! Number of sub-step defined from: 
    191          nn_baro      =  30         ! =F : the number of sub-step in rn_rdt seconds 
     191         nn_baro      =  30         ! =F : the number of sub-step in rn_Dt seconds 
    192192/ 
    193193!----------------------------------------------------------------------- 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/tests/BENCH/EXPREF/namelist_cfg_orca12_like

    r11536 r12406  
    3636&namdom        !   time and space domain 
    3737!----------------------------------------------------------------------- 
    38    rn_rdt      =    300.   !  time step for the dynamics (and tracer if nn_acc=0) 
     38   rn_Dt      =    300.   !  time step for the dynamics (and tracer if nn_acc=0) 
    3939   rn_atfp     =   0.05    !  asselin time filter parameter 
    4040/ 
     
    188188   ln_dynspg_ts   = .true.   ! split-explicit free surface 
    189189      ln_bt_auto    = .false.    ! Number of sub-step defined from: 
    190          nn_baro      =  30         ! =F : the number of sub-step in rn_rdt seconds 
     190         nn_baro      =  30         ! =F : the number of sub-step in rn_Dt seconds 
    191191/ 
    192192!----------------------------------------------------------------------- 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/tests/BENCH/EXPREF/namelist_cfg_orca1_like

    r11536 r12406  
    3636&namdom        !   time and space domain 
    3737!----------------------------------------------------------------------- 
    38    rn_rdt      =   3600.   !  time step for the dynamics (and tracer if nn_acc=0) 
     38   rn_Dt      =   3600.   !  time step for the dynamics (and tracer if nn_acc=0) 
    3939   rn_atfp     =   0.05    !  asselin time filter parameter 
    4040/ 
     
    188188   ln_dynspg_ts   = .true.   ! split-explicit free surface 
    189189      ln_bt_auto    = .false.    ! Number of sub-step defined from: 
    190          nn_baro      =  30         ! =F : the number of sub-step in rn_rdt seconds 
     190         nn_baro      =  30         ! =F : the number of sub-step in rn_Dt seconds 
    191191/ 
    192192!----------------------------------------------------------------------- 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/tests/CANAL/EXPREF/namelist_cfg

    r12377 r12406  
    5959!----------------------------------------------------------------------- 
    6060   ln_linssh   =  .false.  !  =T  linear free surface  ==>>  model level are fixed in time 
    61    rn_rdt      =   1440.   !  time step for the dynamics (and tracer if nn_acc=0) 
     61   rn_Dt      =   1440.   !  time step for the dynamics (and tracer if nn_acc=0) 
    6262   rn_atfp     =   0.05    !  asselin time filter parameter 
    6363/ 
     
    211211         !                          !                     = 2 Boxcar over 2*nn_baro  "    " 
    212212      ln_bt_auto    = .false.    ! Number of sub-step defined from: 
    213          nn_baro      =  24         ! =F : the number of sub-step in rn_rdt seconds 
     213         nn_baro      =  24         ! =F : the number of sub-step in rn_Dt seconds 
    214214/ 
    215215!----------------------------------------------------------------------- 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/tests/CANAL/MY_SRC/diawri.F90

    r12377 r12406  
    563563      clop = "x"         ! no use of the mask value (require less cpu time and otherwise the model crashes) 
    564564#if defined key_diainstant 
    565       zsto = nn_write * rdt 
     565      zsto = nn_write * rn_Dt 
    566566      clop = "inst("//TRIM(clop)//")" 
    567567#else 
    568       zsto=rdt 
     568      zsto=rn_Dt 
    569569      clop = "ave("//TRIM(clop)//")" 
    570570#endif 
    571       zout = nn_write * rdt 
    572       zmax = ( nitend - nit000 + 1 ) * rdt 
     571      zout = nn_write * rn_Dt 
     572      zmax = ( nitend - nit000 + 1 ) * rn_Dt 
    573573 
    574574      ! Define indices of the horizontal output zoom and vertical limit storage 
     
    590590 
    591591         ! Compute julian date from starting date of the run 
    592          CALL ymds2ju( nyear, nmonth, nday, rdt, zjulian ) 
     592         CALL ymds2ju( nyear, nmonth, nday, rn_Dt, zjulian ) 
    593593         zjulian = zjulian - adatrj   !   set calendar origin to the beginning of the experiment 
    594594         IF(lwp)WRITE(numout,*) 
     
    612612         CALL histbeg( clhstnam, jpi, glamt, jpj, gphit,           &  ! Horizontal grid: glamt and gphit 
    613613            &          iimi, iima-iimi+1, ijmi, ijma-ijmi+1,       & 
    614             &          nit000-1, zjulian, rdt, nh_T, nid_T, domain_id=nidom, snc4chunks=snc4set ) 
     614            &          nit000-1, zjulian, rn_Dt, nh_T, nid_T, domain_id=nidom, snc4chunks=snc4set ) 
    615615         CALL histvert( nid_T, "deptht", "Vertical T levels",      &  ! Vertical grid: gdept 
    616616            &           "m", ipk, gdept_1d, nz_T, "down" ) 
     
    648648         CALL histbeg( clhstnam, jpi, glamu, jpj, gphiu,           &  ! Horizontal grid: glamu and gphiu 
    649649            &          iimi, iima-iimi+1, ijmi, ijma-ijmi+1,       & 
    650             &          nit000-1, zjulian, rdt, nh_U, nid_U, domain_id=nidom, snc4chunks=snc4set ) 
     650            &          nit000-1, zjulian, rn_Dt, nh_U, nid_U, domain_id=nidom, snc4chunks=snc4set ) 
    651651         CALL histvert( nid_U, "depthu", "Vertical U levels",      &  ! Vertical grid: gdept 
    652652            &           "m", ipk, gdept_1d, nz_U, "down" ) 
     
    661661         CALL histbeg( clhstnam, jpi, glamv, jpj, gphiv,           &  ! Horizontal grid: glamv and gphiv 
    662662            &          iimi, iima-iimi+1, ijmi, ijma-ijmi+1,       & 
    663             &          nit000-1, zjulian, rdt, nh_V, nid_V, domain_id=nidom, snc4chunks=snc4set ) 
     663            &          nit000-1, zjulian, rn_Dt, nh_V, nid_V, domain_id=nidom, snc4chunks=snc4set ) 
    664664         CALL histvert( nid_V, "depthv", "Vertical V levels",      &  ! Vertical grid : gdept 
    665665            &          "m", ipk, gdept_1d, nz_V, "down" ) 
     
    674674         CALL histbeg( clhstnam, jpi, glamt, jpj, gphit,           &  ! Horizontal grid: glamt and gphit 
    675675            &          iimi, iima-iimi+1, ijmi, ijma-ijmi+1,       & 
    676             &          nit000-1, zjulian, rdt, nh_W, nid_W, domain_id=nidom, snc4chunks=snc4set ) 
     676            &          nit000-1, zjulian, rn_Dt, nh_W, nid_W, domain_id=nidom, snc4chunks=snc4set ) 
    677677         CALL histvert( nid_W, "depthw", "Vertical W levels",      &  ! Vertical grid: gdepw 
    678678            &          "m", ipk, gdepw_1d, nz_W, "down" ) 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/tests/CANAL/MY_SRC/domvvl.F90

    r12397 r12406  
    209209         IF( ln_vvl_ztilde_as_zstar ) THEN   ! z-star emulation using z-tile 
    210210            frq_rst_e3t(:,:) = 0._wp               !Ignore namelist settings 
    211             frq_rst_hdv(:,:) = 1._wp / rdt 
     211            frq_rst_hdv(:,:) = 1._wp / rn_Dt 
    212212         ENDIF 
    213213         IF ( ln_vvl_zstar_at_eqtor ) THEN   ! use z-star in vicinity of the Equator 
     
    222222                     ! values inside the equatorial band (ztilde as zstar) 
    223223                     frq_rst_e3t(ji,jj) =  0.0_wp 
    224                      frq_rst_hdv(ji,jj) =  1.0_wp / rdt 
     224                     frq_rst_hdv(ji,jj) =  1.0_wp / rn_Dt 
    225225                  ELSE                                      ! transition band (2.5 to 6 degrees N/S) 
    226226                     !                                      ! (linearly transition from z-tilde to z-star) 
     
    228228                        &            * (  1.0_wp - COS( rad*(ABS(gphit(ji,jj))-2.5_wp)  & 
    229229                        &                                          * 180._wp / 3.5_wp ) ) 
    230                      frq_rst_hdv(ji,jj) = (1.0_wp / rdt)                                & 
    231                         &            + (  frq_rst_hdv(ji,jj)-(1.e0_wp / rdt) )*0.5_wp   & 
     230                     frq_rst_hdv(ji,jj) = (1.0_wp / rn_Dt)                                & 
     231                        &            + (  frq_rst_hdv(ji,jj)-(1.e0_wp / rn_Dt) )*0.5_wp   & 
    232232                        &            * (  1._wp  - COS( rad*(ABS(gphit(ji,jj))-2.5_wp)  & 
    233233                        &                                          * 180._wp / 3.5_wp ) ) 
     
    240240                  ij0 = 128   ;   ij1 = 135   ;    
    241241                  frq_rst_e3t( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) ) =  0.0_wp 
    242                   frq_rst_hdv( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) ) =  1.e0_wp / rdt 
     242                  frq_rst_hdv( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) ) =  1.e0_wp / rn_Dt 
    243243               ENDIF 
    244244            ENDIF 
     
    349349            IF( kt > nit000 ) THEN 
    350350               DO jk = 1, jpkm1 
    351                   hdiv_lf(:,:,jk) = hdiv_lf(:,:,jk) - rdt * frq_rst_hdv(:,:)   & 
     351                  hdiv_lf(:,:,jk) = hdiv_lf(:,:,jk) - rn_Dt * frq_rst_hdv(:,:)   & 
    352352                     &          * ( hdiv_lf(:,:,jk) - e3t(:,:,jk,Kmm) * ( hdiv(:,:,jk) - zhdiv(:,:) ) ) 
    353353               END DO 
     
    10101010            WRITE(numout,*) '                         rn_rst_e3t     = 0.e0' 
    10111011            WRITE(numout,*) '            hard-wired : z-tilde cutoff frequency of low-pass filter (days)' 
    1012             WRITE(numout,*) '                         rn_lf_cutoff   = 1.0/rdt' 
     1012            WRITE(numout,*) '                         rn_lf_cutoff   = 1.0/rn_Dt' 
    10131013         ELSE 
    10141014            WRITE(numout,*) '      z-tilde to zstar restoration timescale (days)        rn_rst_e3t   = ', rn_rst_e3t 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/tests/CANAL/MY_SRC/sbcmod.F90

    r12377 r12406  
    186186      ! 
    187187      IF( .NOT.ln_usr ) THEN     ! the model calendar needs some specificities (except in user defined case) 
    188          IF( MOD( rday , rdt ) /= 0. )   CALL ctl_stop( 'the time step must devide the number of second of in a day' ) 
    189          IF( MOD( rday , 2.  ) /= 0. )   CALL ctl_stop( 'the number of second of in a day must be an even number'    ) 
    190          IF( MOD( rdt  , 2.  ) /= 0. )   CALL ctl_stop( 'the time step (in second) must be an even number'           ) 
     188         IF( MOD( rday , rn_Dt ) /= 0. )   CALL ctl_stop( 'the time step must devide the number of second of in a day' ) 
     189         IF( MOD( rday , 2.    ) /= 0. )   CALL ctl_stop( 'the number of second of in a day must be an even number'    ) 
     190         IF( MOD( rn_Dt, 2.    ) /= 0. )   CALL ctl_stop( 'the time step (in second) must be an even number'           ) 
    191191      ENDIF 
    192192      !                       !**  check option consistency 
     
    298298      !     SAS time-step has to be declared in OASIS (mandatory) -> nn_fsbc has to be modified accordingly 
    299299      IF( nn_components /= jp_iam_nemo ) THEN 
    300          IF( nn_components == jp_iam_opa )   nn_fsbc = cpl_freq('O_SFLX') / NINT(rdt) 
    301          IF( nn_components == jp_iam_sas )   nn_fsbc = cpl_freq('I_SFLX') / NINT(rdt) 
     300         IF( nn_components == jp_iam_opa )   nn_fsbc = cpl_freq('O_SFLX') / NINT(rn_Dt) 
     301         IF( nn_components == jp_iam_sas )   nn_fsbc = cpl_freq('I_SFLX') / NINT(rn_Dt) 
    302302         ! 
    303303         IF(lwp)THEN 
     
    316316      ENDIF 
    317317      ! 
    318       IF( MOD( rday, REAL(nn_fsbc, wp) * rdt ) /= 0 )   & 
     318      IF( MOD( rday, REAL(nn_fsbc, wp) * rn_Dt ) /= 0 )   & 
    319319         &  CALL ctl_warn( 'sbc_init : nn_fsbc is NOT a multiple of the number of time steps in a day' ) 
    320320      ! 
    321       IF( ln_dm2dc .AND. NINT(rday) / ( nn_fsbc * NINT(rdt) ) < 8  )   & 
     321      IF( ln_dm2dc .AND. NINT(rday) / ( nn_fsbc * NINT(rn_Dt) ) < 8  )   & 
    322322         &   CALL ctl_warn( 'sbc_init : diurnal cycle for qsr: the sampling of the diurnal cycle is too small...' ) 
    323323      ! 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/tests/CANAL/MY_SRC/trazdf.F90

    r12377 r12406  
    6464      ENDIF 
    6565      ! 
    66       IF( neuler == 0 .AND. kt == nit000 ) THEN   ;   r2dt =      rdt   ! at nit000, =   rdt (restarting with Euler time stepping) 
    67       ELSEIF( kt <= nit000 + 1           ) THEN   ;   r2dt = 2. * rdt   ! otherwise, = 2 rdt (leapfrog) 
    68       ENDIF 
    69       ! 
    7066      IF( l_trdtra )   THEN                  !* Save ta and sa trends 
    7167         ALLOCATE( ztrdt(jpi,jpj,jpk) , ztrds(jpi,jpj,jpk) ) 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/tests/ICE_ADV1D/EXPREF/namelist_cfg

    r12377 r12406  
    3939   ln_linssh   = .true.    !  =T  linear free surface  ==>>  model level are fixed in time 
    4040   ! 
    41    rn_rdt      =   2.     !  time step for the dynamics (and tracer if nn_acc=0) 
     41   rn_Dt      =   2.     !  time step for the dynamics (and tracer if nn_acc=0) 
    4242   ! 
    4343   ln_meshmask = .true.   !  =T  create a mesh file 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/tests/ICE_ADV1D/EXPREF/namelist_cfg_120pts

    r12377 r12406  
    3939   ln_linssh   = .true.    !  =T  linear free surface  ==>>  model level are fixed in time 
    4040   ! 
    41    rn_rdt      =   2.     !  time step for the dynamics (and tracer if nn_acc=0) 
     41   rn_Dt      =   2.     !  time step for the dynamics (and tracer if nn_acc=0) 
    4242   ! 
    4343   ln_meshmask = .true.   !  =T  create a mesh file 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/tests/ICE_ADV1D/EXPREF/namelist_cfg_240pts

    r12377 r12406  
    3939   ln_linssh   = .true.    !  =T  linear free surface  ==>>  model level are fixed in time 
    4040   ! 
    41    rn_rdt      =   2.     !  time step for the dynamics (and tracer if nn_acc=0) 
     41   rn_Dt      =   2.     !  time step for the dynamics (and tracer if nn_acc=0) 
    4242   ! 
    4343   ln_meshmask = .true.   !  =T  create a mesh file 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/tests/ICE_ADV1D/EXPREF/namelist_cfg_60pts

    r12377 r12406  
    3939   ln_linssh   = .true.    !  =T  linear free surface  ==>>  model level are fixed in time 
    4040   ! 
    41    rn_rdt      =   2.     !  time step for the dynamics (and tracer if nn_acc=0) 
     41   rn_Dt      =   2.     !  time step for the dynamics (and tracer if nn_acc=0) 
    4242   ! 
    4343   ln_meshmask = .true.   !  =T  create a mesh file 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/tests/ICE_ADV2D/EXPREF/namelist_cfg

    r12377 r12406  
    3939   ln_linssh   = .true.    !  =T  linear free surface  ==>>  model level are fixed in time 
    4040   ! 
    41    rn_rdt      = 1200.     !  time step for the dynamics (and tracer if nn_acc=0) 
     41   rn_Dt      = 1200.     !  time step for the dynamics (and tracer if nn_acc=0) 
    4242   ! 
    4343   ln_meshmask = .true.   !  =T  create a mesh file 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/tests/ICE_AGRIF/EXPREF/1_namelist_cfg

    r12377 r12406  
    3939   ln_linssh   = .true.   !  =T  linear free surface  ==>>  model level are fixed in time 
    4040   ! 
    41    rn_rdt      = 400.     !  time step for the dynamics (and tracer if nn_acc=0) 
     41   rn_Dt      = 400.     !  time step for the dynamics (and tracer if nn_acc=0) 
    4242   ! 
    4343   ln_meshmask = .true.   !  =T  create a mesh file 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/tests/ICE_AGRIF/EXPREF/namelist_cfg

    r12377 r12406  
    3939   ln_linssh   = .true.    !  =T  linear free surface  ==>>  model level are fixed in time 
    4040   ! 
    41    rn_rdt      = 1200.     !  time step for the dynamics (and tracer if nn_acc=0) 
     41   rn_Dt      = 1200.     !  time step for the dynamics (and tracer if nn_acc=0) 
    4242   ! 
    4343   ln_meshmask = .true.   !  =T  create a mesh file 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/tests/ISOMIP+/EXPREF/namelist_cfg

    r12077 r12406  
    4444&namdom        !   time and space domain 
    4545!----------------------------------------------------------------------- 
    46 rn_rdt = 720. 
     46rn_Dt = 720. 
    4747/ 
    4848!----------------------------------------------------------------------- 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/tests/ISOMIP+/MY_SRC/sbcfwb.F90

    r12353 r12406  
    151151         ENDIF    
    152152         !                                         ! Update fwfold if new year start 
    153          ikty = 365 * 86400 / rdt                  !!bug  use of 365 days leap year or 360d year !!!!!!! 
     153         ikty = 365 * 86400 / rn_Dt               !!bug  use of 365 days leap year or 360d year !!!!!!! 
    154154         IF( MOD( kt, ikty ) == 0 ) THEN 
    155155            a_fwb_b = a_fwb                           ! mean sea level taking into account the ice+snow 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/tests/ISOMIP/EXPREF/namelist_cfg

    r12377 r12406  
    4444!----------------------------------------------------------------------- 
    4545   ln_linssh   = .false.   !  =T  linear free surface  ==>>  model level are fixed in time 
    46    rn_rdt      = 1800.     !  time step for the dynamics (and tracer if nn_acc=0) 
     46   rn_Dt      = 1800.     !  time step for the dynamics (and tracer if nn_acc=0) 
    4747/ 
    4848!----------------------------------------------------------------------- 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/tests/LOCK_EXCHANGE/EXPREF/namelist_FCT2_flux_cen2_cfg

    r12377 r12406  
    3333   ln_linssh   = .false.   !  =T  linear free surface  ==>>  model level are fixed in time 
    3434   ! 
    35    rn_rdt      =    1.     !  time step for the dynamics (and tracer if nn_acc=0) 
     35   rn_Dt      =    1.     !  time step for the dynamics (and tracer if nn_acc=0) 
    3636   rn_atfp     =    0.1    !  asselin time filter parameter 
    3737   ! 
     
    147147         !                          !                     = 2 Boxcar over 2*nn_baro  "    " 
    148148      ln_bt_auto    = .true.     ! Number of sub-step defined from: 
    149          nn_baro      =  1         ! =F : the number of sub-step in rn_rdt seconds 
     149         nn_baro      =  1         ! =F : the number of sub-step in rn_Dt seconds 
    150150/ 
    151151!----------------------------------------------------------------------- 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/tests/LOCK_EXCHANGE/EXPREF/namelist_FCT2_flux_ubs_cfg

    r12377 r12406  
    4747   ln_linssh   = .false.   !  =T  linear free surface  ==>>  model level are fixed in time 
    4848   ! 
    49    rn_rdt      =    1.     !  time step for the dynamics (and tracer if nn_acc=0) 
     49   rn_Dt      =    1.     !  time step for the dynamics (and tracer if nn_acc=0) 
    5050   rn_atfp     =    0.1    !  asselin time filter parameter 
    5151   ! 
     
    219219         !                          !                     = 2 Boxcar over 2*nn_baro  "    " 
    220220      ln_bt_auto    = .true.     ! Number of sub-step defined from: 
    221          nn_baro      =  1         ! =F : the number of sub-step in rn_rdt seconds 
     221         nn_baro      =  1         ! =F : the number of sub-step in rn_Dt seconds 
    222222/ 
    223223!----------------------------------------------------------------------- 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/tests/LOCK_EXCHANGE/EXPREF/namelist_FCT2_vect_eenH_cfg

    r12377 r12406  
    3333   ln_linssh   = .false.   !  =T  linear free surface  ==>>  model level are fixed in time 
    3434   ! 
    35    rn_rdt      =    1.     !  time step for the dynamics (and tracer if nn_acc=0) 
     35   rn_Dt      =    1.     !  time step for the dynamics (and tracer if nn_acc=0) 
    3636   rn_atfp     =    0.1    !  asselin time filter parameter 
    3737   ! 
     
    147147         !                          !                     = 2 Boxcar over 2*nn_baro  "    " 
    148148      ln_bt_auto    = .true.     ! Number of sub-step defined from: 
    149          nn_baro      =  1         ! =F : the number of sub-step in rn_rdt seconds 
     149         nn_baro      =  1         ! =F : the number of sub-step in rn_Dt seconds 
    150150/ 
    151151!----------------------------------------------------------------------- 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/tests/LOCK_EXCHANGE/EXPREF/namelist_FCT2_vect_een_cfg

    r12377 r12406  
    3333   ln_linssh   = .false.   !  =T  linear free surface  ==>>  model level are fixed in time 
    3434   ! 
    35    rn_rdt      =    1.     !  time step for the dynamics (and tracer if nn_acc=0) 
     35   rn_Dt      =    1.     !  time step for the dynamics (and tracer if nn_acc=0) 
    3636   rn_atfp     =    0.1    !  asselin time filter parameter 
    3737   ! 
     
    147147         !                          !                     = 2 Boxcar over 2*nn_baro  "    " 
    148148      ln_bt_auto    = .true.     ! Number of sub-step defined from: 
    149          nn_baro      =  1         ! =F : the number of sub-step in rn_rdt seconds 
     149         nn_baro      =  1         ! =F : the number of sub-step in rn_Dt seconds 
    150150/ 
    151151!----------------------------------------------------------------------- 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/tests/LOCK_EXCHANGE/EXPREF/namelist_FCT2_vect_ene_cfg

    r12377 r12406  
    3333   ln_linssh   = .false.   !  =T  linear free surface  ==>>  model level are fixed in time 
    3434   ! 
    35    rn_rdt      =    1.     !  time step for the dynamics (and tracer if nn_acc=0) 
     35   rn_Dt      =    1.     !  time step for the dynamics (and tracer if nn_acc=0) 
    3636   rn_atfp     =    0.1    !  asselin time filter parameter 
    3737   ! 
     
    147147         !                          !                     = 2 Boxcar over 2*nn_baro  "    " 
    148148      ln_bt_auto    = .true.     ! Number of sub-step defined from: 
    149          nn_baro      =  1         ! =F : the number of sub-step in rn_rdt seconds 
     149         nn_baro      =  1         ! =F : the number of sub-step in rn_Dt seconds 
    150150/ 
    151151!----------------------------------------------------------------------- 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/tests/LOCK_EXCHANGE/EXPREF/namelist_FCT2_vect_ens_cfg

    r12377 r12406  
    3333   ln_linssh   = .false.   !  =T  linear free surface  ==>>  model level are fixed in time 
    3434   ! 
    35    rn_rdt      =    1.     !  time step for the dynamics (and tracer if nn_acc=0) 
     35   rn_Dt      =    1.     !  time step for the dynamics (and tracer if nn_acc=0) 
    3636   rn_atfp     =    0.1    !  asselin time filter parameter 
    3737   ! 
     
    147147         !                          !                     = 2 Boxcar over 2*nn_baro  "    " 
    148148      ln_bt_auto    = .true.     ! Number of sub-step defined from: 
    149          nn_baro      =  1         ! =F : the number of sub-step in rn_rdt seconds 
     149         nn_baro      =  1         ! =F : the number of sub-step in rn_Dt seconds 
    150150/ 
    151151!----------------------------------------------------------------------- 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/tests/LOCK_EXCHANGE/EXPREF/namelist_FCT4_flux_cen2_cfg

    r12377 r12406  
    3333   ln_linssh   = .false.   !  =T  linear free surface  ==>>  model level are fixed in time 
    3434   ! 
    35    rn_rdt      =    1.     !  time step for the dynamics (and tracer if nn_acc=0) 
     35   rn_Dt      =    1.     !  time step for the dynamics (and tracer if nn_acc=0) 
    3636   rn_atfp     =    0.1    !  asselin time filter parameter 
    3737   ! 
     
    147147         !                          !                     = 2 Boxcar over 2*nn_baro  "    " 
    148148      ln_bt_auto    = .true.     ! Number of sub-step defined from: 
    149          nn_baro      =  1         ! =F : the number of sub-step in rn_rdt seconds 
     149         nn_baro      =  1         ! =F : the number of sub-step in rn_Dt seconds 
    150150/ 
    151151!----------------------------------------------------------------------- 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/tests/LOCK_EXCHANGE/EXPREF/namelist_FCT4_flux_ubs_cfg

    r12377 r12406  
    3333   ln_linssh   = .false.   !  =T  linear free surface  ==>>  model level are fixed in time 
    3434   ! 
    35    rn_rdt      =    1.     !  time step for the dynamics (and tracer if nn_acc=0) 
     35   rn_Dt      =    1.     !  time step for the dynamics (and tracer if nn_acc=0) 
    3636   rn_atfp     =    0.1    !  asselin time filter parameter 
    3737   ! 
     
    147147         !                          !                     = 2 Boxcar over 2*nn_baro  "    " 
    148148      ln_bt_auto    = .true.     ! Number of sub-step defined from: 
    149          nn_baro      =  1         ! =F : the number of sub-step in rn_rdt seconds 
     149         nn_baro      =  1         ! =F : the number of sub-step in rn_Dt seconds 
    150150/ 
    151151!----------------------------------------------------------------------- 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/tests/LOCK_EXCHANGE/EXPREF/namelist_FCT4_vect_eenH_cfg

    r12377 r12406  
    3333   ln_linssh   = .false.   !  =T  linear free surface  ==>>  model level are fixed in time 
    3434   ! 
    35    rn_rdt      =    1.     !  time step for the dynamics (and tracer if nn_acc=0) 
     35   rn_Dt      =    1.     !  time step for the dynamics (and tracer if nn_acc=0) 
    3636   rn_atfp     =    0.1    !  asselin time filter parameter 
    3737   ! 
     
    147147         !                          !                     = 2 Boxcar over 2*nn_baro  "    " 
    148148      ln_bt_auto    = .true.     ! Number of sub-step defined from: 
    149          nn_baro      =  1         ! =F : the number of sub-step in rn_rdt seconds 
     149         nn_baro      =  1         ! =F : the number of sub-step in rn_Dt seconds 
    150150/ 
    151151!----------------------------------------------------------------------- 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/tests/LOCK_EXCHANGE/EXPREF/namelist_FCT4_vect_een_cfg

    r12377 r12406  
    3333   ln_linssh   = .false.   !  =T  linear free surface  ==>>  model level are fixed in time 
    3434   ! 
    35    rn_rdt      =    1.     !  time step for the dynamics (and tracer if nn_acc=0) 
     35   rn_Dt      =    1.     !  time step for the dynamics (and tracer if nn_acc=0) 
    3636   rn_atfp     =    0.1    !  asselin time filter parameter 
    3737   ! 
     
    147147         !                          !                     = 2 Boxcar over 2*nn_baro  "    " 
    148148      ln_bt_auto    = .true.     ! Number of sub-step defined from: 
    149          nn_baro      =  1         ! =F : the number of sub-step in rn_rdt seconds 
     149         nn_baro      =  1         ! =F : the number of sub-step in rn_Dt seconds 
    150150/ 
    151151!----------------------------------------------------------------------- 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/tests/LOCK_EXCHANGE/EXPREF/namelist_FCT4_vect_ene_cfg

    r12377 r12406  
    3333   ln_linssh   = .false.   !  =T  linear free surface  ==>>  model level are fixed in time 
    3434   ! 
    35    rn_rdt      =    1.     !  time step for the dynamics (and tracer if nn_acc=0) 
     35   rn_Dt      =    1.     !  time step for the dynamics (and tracer if nn_acc=0) 
    3636   rn_atfp     =    0.1    !  asselin time filter parameter 
    3737   ! 
     
    147147         !                          !                     = 2 Boxcar over 2*nn_baro  "    " 
    148148      ln_bt_auto    = .true.     ! Number of sub-step defined from: 
    149          nn_baro      =  1         ! =F : the number of sub-step in rn_rdt seconds 
     149         nn_baro      =  1         ! =F : the number of sub-step in rn_Dt seconds 
    150150/ 
    151151!----------------------------------------------------------------------- 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/tests/LOCK_EXCHANGE/EXPREF/namelist_FCT4_vect_ens_cfg

    r12377 r12406  
    3333   ln_linssh   = .false.   !  =T  linear free surface  ==>>  model level are fixed in time 
    3434   ! 
    35    rn_rdt      =    1.     !  time step for the dynamics (and tracer if nn_acc=0) 
     35   rn_Dt      =    1.     !  time step for the dynamics (and tracer if nn_acc=0) 
    3636   rn_atfp     =    0.1    !  asselin time filter parameter 
    3737   ! 
     
    147147         !                          !                     = 2 Boxcar over 2*nn_baro  "    " 
    148148      ln_bt_auto    = .true.     ! Number of sub-step defined from: 
    149          nn_baro      =  1         ! =F : the number of sub-step in rn_rdt seconds 
     149         nn_baro      =  1         ! =F : the number of sub-step in rn_Dt seconds 
    150150/ 
    151151!----------------------------------------------------------------------- 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/tests/OVERFLOW/EXPREF/namelist_sco_FCT2_flux_cen-ahm1000_cfg

    r11848 r12406  
    3333&namdom        !   space and time domain (bathymetry, mesh, timestep) 
    3434!----------------------------------------------------------------------- 
    35    rn_rdt      =   10.     !  time step for the dynamics (and tracer if nn_acc=0) 
     35   rn_Dt      =   10.     !  time step for the dynamics (and tracer if nn_acc=0) 
    3636   rn_atfp     =    0.1    !  asselin time filter parameter 
    3737/ 
     
    158158         !                          !                     = 2 Boxcar over 2*nn_baro  "    " 
    159159      ln_bt_auto    = .true.     ! Number of sub-step defined from: 
    160          nn_baro      =  1         ! =F : the number of sub-step in rn_rdt seconds 
     160         nn_baro      =  1         ! =F : the number of sub-step in rn_Dt seconds 
    161161/ 
    162162!----------------------------------------------------------------------- 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/tests/OVERFLOW/EXPREF/namelist_sco_FCT2_flux_ubs_cfg

    r11848 r12406  
    3333&namdom        !   space and time domain (bathymetry, mesh, timestep) 
    3434!----------------------------------------------------------------------- 
    35    rn_rdt      =   10.     !  time step for the dynamics (and tracer if nn_acc=0) 
     35   rn_Dt      =   10.     !  time step for the dynamics (and tracer if nn_acc=0) 
    3636   rn_atfp     =    0.1    !  asselin time filter parameter 
    3737/ 
     
    158158         !                          !                     = 2 Boxcar over 2*nn_baro  "    " 
    159159      ln_bt_auto    = .true.     ! Number of sub-step defined from: 
    160          nn_baro      =  1         ! =F : the number of sub-step in rn_rdt seconds 
     160         nn_baro      =  1         ! =F : the number of sub-step in rn_Dt seconds 
    161161/ 
    162162!----------------------------------------------------------------------- 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/tests/OVERFLOW/EXPREF/namelist_sco_FCT4_flux_cen-ahm1000_cfg

    r10075 r12406  
    3333&namdom        !   space and time domain (bathymetry, mesh, timestep) 
    3434!----------------------------------------------------------------------- 
    35    rn_rdt      =   10.     !  time step for the dynamics (and tracer if nn_acc=0) 
     35   rn_Dt      =   10.     !  time step for the dynamics (and tracer if nn_acc=0) 
    3636   rn_atfp     =    0.1    !  asselin time filter parameter 
    3737/ 
     
    158158         !                          !                     = 2 Boxcar over 2*nn_baro  "    " 
    159159      ln_bt_auto    = .true.     ! Number of sub-step defined from: 
    160          nn_baro      =  1         ! =F : the number of sub-step in rn_rdt seconds 
     160         nn_baro      =  1         ! =F : the number of sub-step in rn_Dt seconds 
    161161/ 
    162162!----------------------------------------------------------------------- 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/tests/OVERFLOW/EXPREF/namelist_sco_FCT4_flux_ubs_cfg

    r11848 r12406  
    3333&namdom        !   space and time domain (bathymetry, mesh, timestep) 
    3434!----------------------------------------------------------------------- 
    35    rn_rdt      =   10.     !  time step for the dynamics (and tracer if nn_acc=0) 
     35   rn_Dt      =   10.     !  time step for the dynamics (and tracer if nn_acc=0) 
    3636   rn_atfp     =    0.1    !  asselin time filter parameter 
    3737/ 
     
    158158         !                          !                     = 2 Boxcar over 2*nn_baro  "    " 
    159159      ln_bt_auto    = .true.     ! Number of sub-step defined from: 
    160          nn_baro      =  1         ! =F : the number of sub-step in rn_rdt seconds 
     160         nn_baro      =  1         ! =F : the number of sub-step in rn_Dt seconds 
    161161/ 
    162162!----------------------------------------------------------------------- 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/tests/OVERFLOW/EXPREF/namelist_zps_FCT2_flux_ubs_cfg

    r11848 r12406  
    3333&namdom        !   space and time domain (bathymetry, mesh, timestep) 
    3434!----------------------------------------------------------------------- 
    35    rn_rdt      =   10.     !  time step for the dynamics (and tracer if nn_acc=0) 
     35   rn_Dt      =   10.     !  time step for the dynamics (and tracer if nn_acc=0) 
    3636   rn_atfp     =    0.1    !  asselin time filter parameter 
    3737/ 
     
    158158         !                          !                     = 2 Boxcar over 2*nn_baro  "    " 
    159159      ln_bt_auto    = .true.     ! Number of sub-step defined from: 
    160          nn_baro      =  1         ! =F : the number of sub-step in rn_rdt seconds 
     160         nn_baro      =  1         ! =F : the number of sub-step in rn_Dt seconds 
    161161/ 
    162162!----------------------------------------------------------------------- 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/tests/OVERFLOW/EXPREF/namelist_zps_FCT4_flux_ubs_cfg

    r12377 r12406  
    4646&namdom        !   time and space domain 
    4747!----------------------------------------------------------------------- 
    48    rn_rdt      =   10.     !  time step for the dynamics (and tracer if nn_acc=0) 
     48   rn_Dt      =   10.     !  time step for the dynamics (and tracer if nn_acc=0) 
    4949   rn_atfp     =    0.1    !  asselin time filter parameter 
    5050/ 
     
    219219         !                          !                     = 2 Boxcar over 2*nn_baro  "    " 
    220220      ln_bt_auto    = .true.     ! Number of sub-step defined from: 
    221          nn_baro      =  1         ! =F : the number of sub-step in rn_rdt seconds 
     221         nn_baro      =  1         ! =F : the number of sub-step in rn_Dt seconds 
    222222/ 
    223223!----------------------------------------------------------------------- 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/tests/OVERFLOW/EXPREF/namelist_zps_FCT4_vect_een_cfg

    r11848 r12406  
    3333&namdom        !   space and time domain (bathymetry, mesh, timestep) 
    3434!----------------------------------------------------------------------- 
    35    rn_rdt      =   10.     !  time step for the dynamics (and tracer if nn_acc=0) 
     35   rn_Dt      =   10.     !  time step for the dynamics (and tracer if nn_acc=0) 
    3636   rn_atfp     =    0.1    !  asselin time filter parameter 
    3737/ 
     
    158158         !                          !                     = 2 Boxcar over 2*nn_baro  "    " 
    159159      ln_bt_auto    = .true.     ! Number of sub-step defined from: 
    160          nn_baro      =  1         ! =F : the number of sub-step in rn_rdt seconds 
     160         nn_baro      =  1         ! =F : the number of sub-step in rn_Dt seconds 
    161161/ 
    162162!----------------------------------------------------------------------- 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/tests/STATION_ASF/EXPREF/namelist_coare3p6-noskin_cfg

    r12249 r12406  
    5757&namdom        !   time and space domain / STATION_ASF 
    5858!----------------------------------------------------------------------- 
    59    rn_rdt      =  3600.    !  time step for the dynamics and tracer => frequency of atmospheric forcing 
     59   rn_Dt      =  3600.    !  time step for the dynamics and tracer => frequency of atmospheric forcing 
    6060   ln_meshmask = .false. 
    6161   ln_linssh   = .true.    !  =T  linear free surface  ==>>  model level are fixed in time 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/tests/STATION_ASF/EXPREF/namelist_coare3p6_cfg

    r12249 r12406  
    5757&namdom        !   time and space domain / STATION_ASF 
    5858!----------------------------------------------------------------------- 
    59    rn_rdt      =  3600.    !  time step for the dynamics and tracer => frequency of atmospheric forcing 
     59   rn_Dt      =  3600.    !  time step for the dynamics and tracer => frequency of atmospheric forcing 
    6060   ln_meshmask = .false. 
    6161   ln_linssh   = .true.    !  =T  linear free surface  ==>>  model level are fixed in time 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/tests/STATION_ASF/EXPREF/namelist_ecmwf-noskin_cfg

    r12249 r12406  
    5757&namdom        !   time and space domain / STATION_ASF 
    5858!----------------------------------------------------------------------- 
    59    rn_rdt      =  3600.    !  time step for the dynamics and tracer => frequency of atmospheric forcing 
     59   rn_Dt      =  3600.    !  time step for the dynamics and tracer => frequency of atmospheric forcing 
    6060   ln_meshmask = .false. 
    6161   ln_linssh   = .true.    !  =T  linear free surface  ==>>  model level are fixed in time 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/tests/STATION_ASF/EXPREF/namelist_ecmwf_cfg

    r12249 r12406  
    5757&namdom        !   time and space domain / STATION_ASF 
    5858!----------------------------------------------------------------------- 
    59    rn_rdt      =  3600.    !  time step for the dynamics and tracer => frequency of atmospheric forcing 
     59   rn_Dt      =  3600.    !  time step for the dynamics and tracer => frequency of atmospheric forcing 
    6060   ln_meshmask = .false. 
    6161   ln_linssh   = .true.    !  =T  linear free surface  ==>>  model level are fixed in time 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/tests/STATION_ASF/EXPREF/namelist_ncar_cfg

    r12249 r12406  
    5757&namdom        !   time and space domain / STATION_ASF 
    5858!----------------------------------------------------------------------- 
    59    rn_rdt      =  3600.    !  time step for the dynamics and tracer => frequency of atmospheric forcing 
     59   rn_Dt      =  3600.    !  time step for the dynamics and tracer => frequency of atmospheric forcing 
    6060   ln_meshmask = .false. 
    6161   ln_linssh   = .true.    !  =T  linear free surface  ==>>  model level are fixed in time 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/tests/STATION_ASF/MY_SRC/diawri.F90

    r12353 r12406  
    168168      clop = "x"         ! no use of the mask value (require less cpu time and otherwise the model crashes) 
    169169#if defined key_diainstant 
    170       zsto = nn_write * rdt 
     170      zsto = nn_write * rn_Dt 
    171171      clop = "inst("//TRIM(clop)//")" 
    172172#else 
    173       zsto=rdt 
     173      zsto=rn_Dt 
    174174      clop = "ave("//TRIM(clop)//")" 
    175175#endif 
    176       zout = nn_write * rdt 
    177       zmax = ( nitend - nit000 + 1 ) * rdt 
     176      zout = nn_write * rn_Dt 
     177      zmax = ( nitend - nit000 + 1 ) * rn_Dt 
    178178 
    179179      ! Define indices of the horizontal output zoom and vertical limit storage 
     
    195195 
    196196         ! Compute julian date from starting date of the run 
    197          CALL ymds2ju( nyear, nmonth, nday, rdt, zjulian ) 
     197         CALL ymds2ju( nyear, nmonth, nday, rn_Dt, zjulian ) 
    198198         zjulian = zjulian - adatrj   !   set calendar origin to the beginning of the experiment 
    199199         IF(lwp)WRITE(numout,*) 
     
    217217         CALL histbeg( clhstnam, jpi, glamt, jpj, gphit,           &  ! Horizontal grid: glamt and gphit 
    218218            &          iimi, iima-iimi+1, ijmi, ijma-ijmi+1,       & 
    219             &          nit000-1, zjulian, rdt, nh_T, nid_T, domain_id=nidom, snc4chunks=snc4set ) 
     219            &          nit000-1, zjulian, rn_Dt, nh_T, nid_T, domain_id=nidom, snc4chunks=snc4set ) 
    220220         CALL histvert( nid_T, "deptht", "Vertical T levels",      &  ! Vertical grid: gdept 
    221221            &           "m", ipk, gdept_1d, nz_T, "down" ) 
     
    229229         CALL histbeg( clhstnam, jpi, glamu, jpj, gphiu,           &  ! Horizontal grid: glamu and gphiu 
    230230            &          iimi, iima-iimi+1, ijmi, ijma-ijmi+1,       & 
    231             &          nit000-1, zjulian, rdt, nh_U, nid_U, domain_id=nidom, snc4chunks=snc4set ) 
     231            &          nit000-1, zjulian, rn_Dt, nh_U, nid_U, domain_id=nidom, snc4chunks=snc4set ) 
    232232         CALL histvert( nid_U, "depthu", "Vertical U levels",      &  ! Vertical grid: gdept 
    233233            &           "m", ipk, gdept_1d, nz_U, "down" ) 
     
    241241         CALL histbeg( clhstnam, jpi, glamv, jpj, gphiv,           &  ! Horizontal grid: glamv and gphiv 
    242242            &          iimi, iima-iimi+1, ijmi, ijma-ijmi+1,       & 
    243             &          nit000-1, zjulian, rdt, nh_V, nid_V, domain_id=nidom, snc4chunks=snc4set ) 
     243            &          nit000-1, zjulian, rn_Dt, nh_V, nid_V, domain_id=nidom, snc4chunks=snc4set ) 
    244244         CALL histvert( nid_V, "depthv", "Vertical V levels",      &  ! Vertical grid : gdept 
    245245            &          "m", ipk, gdept_1d, nz_V, "down" ) 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/tests/VORTEX/EXPREF/1_namelist_cfg

    r12377 r12406  
    4040!----------------------------------------------------------------------- 
    4141   ln_linssh   =  .false.  !  =T  linear free surface  ==>>  model level are fixed in time 
    42    rn_rdt      =    480.   !  time step for the dynamics (and tracer if nn_acc=0) 
     42   rn_Dt      =    480.   !  time step for the dynamics (and tracer if nn_acc=0) 
    4343   rn_atfp     =   0.05    !  asselin time filter parameter 
    4444/ 
     
    216216         !                          !                     = 2 Boxcar over 2*nn_baro  "    " 
    217217      ln_bt_auto    = .false.    ! Number of sub-step defined from: 
    218          nn_baro      =  24         ! =F : the number of sub-step in rn_rdt seconds 
     218         nn_baro      =  24         ! =F : the number of sub-step in rn_Dt seconds 
    219219/ 
    220220!----------------------------------------------------------------------- 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/tests/VORTEX/EXPREF/namelist_cfg

    r12377 r12406  
    4040!----------------------------------------------------------------------- 
    4141   ln_linssh   =  .false.  !  =T  linear free surface  ==>>  model level are fixed in time 
    42    rn_rdt      =   1440.   !  time step for the dynamics (and tracer if nn_acc=0) 
     42   rn_Dt      =   1440.   !  time step for the dynamics (and tracer if nn_acc=0) 
    4343   rn_atfp     =   0.05    !  asselin time filter parameter 
    4444/ 
     
    207207         !                          !                     = 2 Boxcar over 2*nn_baro  "    " 
    208208      ln_bt_auto    = .false.    ! Number of sub-step defined from: 
    209          nn_baro      =  24         ! =F : the number of sub-step in rn_rdt seconds 
     209         nn_baro      =  24         ! =F : the number of sub-step in rn_Dt seconds 
    210210/ 
    211211!----------------------------------------------------------------------- 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/tests/VORTEX/MY_SRC/domvvl.F90

    r12397 r12406  
    237237         IF( ln_vvl_ztilde_as_zstar ) THEN   ! z-star emulation using z-tile 
    238238            frq_rst_e3t(:,:) = 0._wp               !Ignore namelist settings 
    239             frq_rst_hdv(:,:) = 1._wp / rdt 
     239            frq_rst_hdv(:,:) = 1._wp / rn_Dt 
    240240         ENDIF 
    241241         IF ( ln_vvl_zstar_at_eqtor ) THEN   ! use z-star in vicinity of the Equator 
     
    250250                     ! values inside the equatorial band (ztilde as zstar) 
    251251                     frq_rst_e3t(ji,jj) =  0.0_wp 
    252                      frq_rst_hdv(ji,jj) =  1.0_wp / rdt 
     252                     frq_rst_hdv(ji,jj) =  1.0_wp / rn_Dt 
    253253                  ELSE                                      ! transition band (2.5 to 6 degrees N/S) 
    254254                     !                                      ! (linearly transition from z-tilde to z-star) 
     
    256256                        &            * (  1.0_wp - COS( rad*(ABS(gphit(ji,jj))-2.5_wp)  & 
    257257                        &                                          * 180._wp / 3.5_wp ) ) 
    258                      frq_rst_hdv(ji,jj) = (1.0_wp / rdt)                                & 
    259                         &            + (  frq_rst_hdv(ji,jj)-(1.e0_wp / rdt) )*0.5_wp   & 
     258                     frq_rst_hdv(ji,jj) = (1.0_wp / rn_Dt)                                & 
     259                        &            + (  frq_rst_hdv(ji,jj)-(1.e0_wp / rn_Dt) )*0.5_wp   & 
    260260                        &            * (  1._wp  - COS( rad*(ABS(gphit(ji,jj))-2.5_wp)  & 
    261261                        &                                          * 180._wp / 3.5_wp ) ) 
     
    268268                  ij0 = 128   ;   ij1 = 135   ;    
    269269                  frq_rst_e3t( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) ) =  0.0_wp 
    270                   frq_rst_hdv( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) ) =  1.e0_wp / rdt 
     270                  frq_rst_hdv( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) ) =  1.e0_wp / rn_Dt 
    271271               ENDIF 
    272272            ENDIF 
     
    377377            IF( kt > nit000 ) THEN 
    378378               DO jk = 1, jpkm1 
    379                   hdiv_lf(:,:,jk) = hdiv_lf(:,:,jk) - rdt * frq_rst_hdv(:,:)   & 
     379                  hdiv_lf(:,:,jk) = hdiv_lf(:,:,jk) - rn_Dt * frq_rst_hdv(:,:)   & 
    380380                     &          * ( hdiv_lf(:,:,jk) - e3t(:,:,jk,Kmm) * ( hdiv(:,:,jk) - zhdiv(:,:) ) ) 
    381381               END DO 
     
    10341034            WRITE(numout,*) '                         rn_rst_e3t     = 0.e0' 
    10351035            WRITE(numout,*) '            hard-wired : z-tilde cutoff frequency of low-pass filter (days)' 
    1036             WRITE(numout,*) '                         rn_lf_cutoff   = 1.0/rdt' 
     1036            WRITE(numout,*) '                         rn_lf_cutoff   = 1.0/rn_Dt' 
    10371037         ELSE 
    10381038            WRITE(numout,*) '      z-tilde to zstar restoration timescale (days)        rn_rst_e3t   = ', rn_rst_e3t 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/tests/WAD/EXPREF/namelist_cfg

    r12377 r12406  
    4343   ! 
    4444   ln_meshmask = .true.    !  create (>0) a mesh file or not (=0) 
    45    rn_rdt      =    18.    !  time step for the dynamics  
     45   rn_Dt      =    18.    !  time step for the dynamics  
    4646/ 
    4747!----------------------------------------------------------------------- 
     
    344344   ln_dynspg_ts  = .true.  !  split-explicit free surface 
    345345   ln_bt_auto    = .false. ! Number of sub-step defined from: 
    346    nn_baro       =  12     ! =F : the number of sub-step in rn_rdt seconds 
     346   nn_baro       =  12     ! =F : the number of sub-step in rn_Dt seconds 
    347347/ 
    348348!----------------------------------------------------------------------- 
Note: See TracChangeset for help on using the changeset viewer.