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

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 12443 for NEMO/branches/2020 – NEMO

Changeset 12443 for NEMO/branches/2020


Ignore:
Timestamp:
2020-02-24T14:00:21+01:00 (4 years ago)
Author:
davestorkey
Message:

2020/KERNEL-03_Storkey_Coward_RK3_stage2: More variable renaming:
atfp -> rn_atfp (use namelist parameter everywhere)
rdtbt -> rDt_e
nn_baro -> nn_e
rn_scal_load -> rn_load
rau0 -> rho0

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

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/cfgs/AMM12/EXPREF/namelist_cfg

    r12406 r12443  
    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_Dt seconds 
     305   nn_e      = 30       ! =F : the number of sub-step in rn_Dt seconds 
    306306/ 
    307307!----------------------------------------------------------------------- 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/cfgs/SHARED/namelist_ref

    r12406 r12443  
    651651         rn_tide_gamma = 0.7                   ! Tidal tilt factor 
    652652         ln_scal_load  = .false.               ! Use scalar approximation for 
    653             rn_scal_load = 0.094               !     load potential 
     653            rn_load = 0.094               !     load potential 
    654654         ln_read_load  = .false.               ! Or read load potential from file 
    655655            cn_tide_load = 'tide_LOAD_grid_T.nc'  ! filename for load potential 
     
    989989      ln_bt_av      = .true.     ! Time filtering of barotropic variables 
    990990         nn_bt_flt     = 1          ! Time filter choice  = 0 None 
    991          !                          !                     = 1 Boxcar over   nn_baro sub-steps 
    992          !                          !                     = 2 Boxcar over 2*nn_baro  "    " 
     991         !                          !                     = 1 Boxcar over   nn_e sub-steps 
     992         !                          !                     = 2 Boxcar over 2*nn_e  "    " 
    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_Dt seconds 
     995         nn_e      = 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/src/ICE/icedia.F90

    r12406 r12443  
    109109      ! ---------------------------! 
    110110      ! they must be kept outside an IF(iom_use) because of the call to dia_rst below 
    111       z_frc_volbot = r1_rau0 * glob_sum( 'icedia', -( wfx_ice(:,:) + wfx_snw(:,:) + wfx_err_sub(:,:) ) * e1e2t(:,:) ) * 1.e-9   ! freshwater flux ice/snow-ocean  
    112       z_frc_voltop = r1_rau0 * glob_sum( 'icedia', -( wfx_sub(:,:) + wfx_spr(:,:) )                    * e1e2t(:,:) ) * 1.e-9   ! freshwater flux ice/snow-atm 
    113       z_frc_sal    = r1_rau0 * glob_sum( 'icedia', -      sfx(:,:)                                     * e1e2t(:,:) ) * 1.e-9   ! salt fluxes ice/snow-ocean 
     111      z_frc_volbot = r1_rho0 * glob_sum( 'icedia', -( wfx_ice(:,:) + wfx_snw(:,:) + wfx_err_sub(:,:) ) * e1e2t(:,:) ) * 1.e-9   ! freshwater flux ice/snow-ocean  
     112      z_frc_voltop = r1_rho0 * glob_sum( 'icedia', -( wfx_sub(:,:) + wfx_spr(:,:) )                    * e1e2t(:,:) ) * 1.e-9   ! freshwater flux ice/snow-atm 
     113      z_frc_sal    = r1_rho0 * glob_sum( 'icedia', -      sfx(:,:)                                     * e1e2t(:,:) ) * 1.e-9   ! salt fluxes ice/snow-ocean 
    114114      z_frc_tembot =           glob_sum( 'icedia',  qt_oce_ai(:,:)                                     * e1e2t(:,:) ) * 1.e-20  ! heat on top of ocean (and below ice) 
    115115      z_frc_temtop =           glob_sum( 'icedia',  qt_atm_oi(:,:)                                     * e1e2t(:,:) ) * 1.e-20  ! heat on top of ice-coean 
     
    137137      IF(  iom_use('ibgvolume') .OR. iom_use('ibgsaltco') .OR. iom_use('ibgheatco') .OR. iom_use('ibgheatfx') ) THEN 
    138138             
    139          zdiff_vol = r1_rau0 * glob_sum( 'icedia', ( rhoi*vt_i(:,:) + rhos*vt_s(:,:) - vol_loc_ini(:,:) ) * e1e2t(:,:) ) * 1.e-9   ! freshwater trend (km3)  
    140          zdiff_sal = r1_rau0 * glob_sum( 'icedia', ( rhoi*st_i(:,:)                  - sal_loc_ini(:,:) ) * e1e2t(:,:) ) * 1.e-9   ! salt content trend (km3*pss) 
     139         zdiff_vol = r1_rho0 * glob_sum( 'icedia', ( rhoi*vt_i(:,:) + rhos*vt_s(:,:) - vol_loc_ini(:,:) ) * e1e2t(:,:) ) * 1.e-9   ! freshwater trend (km3)  
     140         zdiff_sal = r1_rho0 * glob_sum( 'icedia', ( rhoi*st_i(:,:)                  - sal_loc_ini(:,:) ) * e1e2t(:,:) ) * 1.e-9   ! salt content trend (km3*pss) 
    141141         zdiff_tem =           glob_sum( 'icedia', ( et_i(:,:) + et_s(:,:)           - tem_loc_ini(:,:) ) * e1e2t(:,:) ) * 1.e-20  ! heat content trend (1.e20 J) 
    142142         !                               + SUM( qevap_ice * a_i_b, dim=3 )       !! clem: I think this term should not be there (but needs a check) 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/ICE/icedyn_adv_pra.F90

    r12377 r12443  
    687687            IF ( pv_i(ji,jj,jl) > 0._wp ) THEN 
    688688               ! 
    689                zvs_excess = MAX( 0._wp, pv_s(ji,jj,jl) - pv_i(ji,jj,jl) * (rau0-rhoi) * r1_rhos ) 
     689               zvs_excess = MAX( 0._wp, pv_s(ji,jj,jl) - pv_i(ji,jj,jl) * (rho0-rhoi) * r1_rhos ) 
    690690               ! 
    691691               IF( zvs_excess > 0._wp ) THEN   ! snow-ice interface deplets below the ocean surface 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/ICE/icedyn_adv_umx.F90

    r12377 r12443  
    15051505            IF ( pv_i(ji,jj,jl) > 0._wp ) THEN 
    15061506               ! 
    1507                zvs_excess = MAX( 0._wp, pv_s(ji,jj,jl) - pv_i(ji,jj,jl) * (rau0-rhoi) * r1_rhos ) 
     1507               zvs_excess = MAX( 0._wp, pv_s(ji,jj,jl) - pv_i(ji,jj,jl) * (rho0-rhoi) * r1_rhos ) 
    15081508               ! 
    15091509               IF( zvs_excess > 0._wp ) THEN   ! snow-ice interface deplets below the ocean surface 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/ICE/icedyn_rhg_evp.F90

    r12377 r12443  
    116116      INTEGER ::   jter         ! local integers 
    117117      ! 
    118       REAL(wp) ::   zrhoco                                              ! rau0 * rn_cio 
     118      REAL(wp) ::   zrhoco                                              ! rho0 * rn_cio 
    119119      REAL(wp) ::   zdtevp, z1_dtevp                                    ! time step for subcycling 
    120120      REAL(wp) ::   ecc2, z1_ecc2                                       ! square of yield ellipse eccenticity 
     
    213213      ! 1) define some variables and initialize arrays 
    214214      !------------------------------------------------------------------------------! 
    215       zrhoco = rau0 * rn_cio  
     215      zrhoco = rho0 * rn_cio  
    216216 
    217217      ! ecc2: square of yield ellipse eccenticrity 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/ICE/iceistate.F90

    r12377 r12443  
    371371      IF( ln_ice_embd ) THEN            ! embedded sea-ice: deplete the initial ssh below sea-ice area 
    372372         ! 
    373          ssh(:,:,Kmm) = ssh(:,:,Kmm) - snwice_mass(:,:) * r1_rau0 
    374          ssh(:,:,Kbb) = ssh(:,:,Kbb) - snwice_mass(:,:) * r1_rau0 
     373         ssh(:,:,Kmm) = ssh(:,:,Kmm) - snwice_mass(:,:) * r1_rho0 
     374         ssh(:,:,Kbb) = ssh(:,:,Kbb) - snwice_mass(:,:) * r1_rho0 
    375375         ! 
    376376         IF( .NOT.ln_linssh ) THEN 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/ICE/icethd.F90

    r12377 r12443  
    116116      ELSE      !  if no ice dynamics => transmit directly the atmospheric stress to the ocean 
    117117         DO_2D_00_00 
    118             zfric(ji,jj) = r1_rau0 * SQRT( 0.5_wp *  & 
     118            zfric(ji,jj) = r1_rho0 * SQRT( 0.5_wp *  & 
    119119               &                         (  utau(ji,jj) * utau(ji,jj) + utau(ji-1,jj) * utau(ji-1,jj)   & 
    120120               &                          + vtau(ji,jj) * vtau(ji,jj) + vtau(ji,jj-1) * vtau(ji,jj-1) ) ) * tmask(ji,jj,1) 
     
    141141 
    142142         ! --- Energy needed to bring ocean surface layer until its freezing (mostly<0 but >0 if supercooling, J.m-2) --- ! 
    143          zqfr     = rau0 * rcp * e3t_m(ji,jj) * ( t_bo(ji,jj) - ( sst_m(ji,jj) + rt0 ) ) * tmask(ji,jj,1)  ! both < 0 (t_bo < sst) and > 0 (t_bo > sst) 
     143         zqfr     = rho0 * rcp * e3t_m(ji,jj) * ( t_bo(ji,jj) - ( sst_m(ji,jj) + rt0 ) ) * tmask(ji,jj,1)  ! both < 0 (t_bo < sst) and > 0 (t_bo > sst) 
    144144         zqfr_neg = MIN( zqfr , 0._wp )                                                                    ! only < 0 
    145145 
    146146         ! --- Sensible ocean-to-ice heat flux (mostly>0 but <0 if supercooling, W/m2) 
    147147         zfric_u            = MAX( SQRT( zfric(ji,jj) ), zfric_umin )  
    148          qsb_ice_bot(ji,jj) = rswitch * rau0 * rcp * zch * zfric_u * ( ( sst_m(ji,jj) + rt0 ) - t_bo(ji,jj) ) ! W.m-2 
     148         qsb_ice_bot(ji,jj) = rswitch * rho0 * rcp * zch * zfric_u * ( ( sst_m(ji,jj) + rt0 ) - t_bo(ji,jj) ) ! W.m-2 
    149149 
    150150         qsb_ice_bot(ji,jj) = rswitch * MIN( qsb_ice_bot(ji,jj), - zqfr_neg * r1_rdtice / MAX( at_i(ji,jj), epsi10 ) ) 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/ICE/icethd_dh.F90

    r10786 r12443  
    7676      REAL(wp) ::   zgrr         ! bottom growth rate 
    7777      REAL(wp) ::   zt_i_new     ! bottom formation temperature 
    78       REAL(wp) ::   z1_rho       ! 1/(rhos+rau0-rhoi) 
     78      REAL(wp) ::   z1_rho       ! 1/(rhos+rho0-rhoi) 
    7979 
    8080      REAL(wp) ::   zQm          ! enthalpy exchanged with the ocean (J/m2), >0 towards the ocean 
     
    571571      ! When snow load excesses Archimede's limit, snow-ice interface goes down under sea-level,  
    572572      ! flooding of seawater transforms snow into ice dh_snowice is positive for the ice 
    573       z1_rho = 1._wp / ( rhos+rau0-rhoi ) 
     573      z1_rho = 1._wp / ( rhos+rho0-rhoi ) 
    574574      DO ji = 1, npti 
    575575         ! 
    576          dh_snowice(ji) = MAX(  0._wp , ( rhos * h_s_1d(ji) + (rhoi-rau0) * h_i_1d(ji) ) * z1_rho ) 
     576         dh_snowice(ji) = MAX(  0._wp , ( rhos * h_s_1d(ji) + (rhoi-rho0) * h_i_1d(ji) ) * z1_rho ) 
    577577 
    578578         h_i_1d(ji)    = h_i_1d(ji) + dh_snowice(ji) 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/ICE/icethd_do.F90

    r12377 r12443  
    141141         ! Physical constants 
    142142         zhicrit = 0.04                                          ! frazil ice thickness 
    143          ztwogp  = 2. * rau0 / ( grav * 0.3 * ( rau0 - rhoi ) )  ! reduced grav 
     143         ztwogp  = 2. * rho0 / ( grav * 0.3 * ( rho0 - rhoi ) )  ! reduced grav 
    144144         zsqcd   = 1.0 / SQRT( 1.3 * zcai )                      ! 1/SQRT(airdensity*drag) 
    145145         zgamafr = 0.03 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/ICE/iceupdate.F90

    r12377 r12443  
    329329      ENDIF 
    330330 
    331       zrhoco = rau0 * rn_cio 
     331      zrhoco = rho0 * rn_cio 
    332332      ! 
    333333      IF( MOD( kt-1, nn_fsbc ) == 0 ) THEN     !==  Ice time-step only  ==!   (i.e. surface module time-step) 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/ICE/icevar.F90

    r12377 r12443  
    717717      !! ** Purpose :  compute the equivalent ssh in lead when sea ice is embedded 
    718718      !! 
    719       !! ** Method  :  ssh_lead = ssh + (Mice + Msnow) / rau0 
     719      !! ** Method  :  ssh_lead = ssh + (Mice + Msnow) / rho0 
    720720      !! 
    721721      !! ** Reference : Jean-Michel Campin, John Marshall, David Ferreira, 
     
    747747         zintb = REAL( nn_fsbc + 1 ) / REAL( nn_fsbc ) * 0.5_wp 
    748748         ! 
    749          zsnwiceload(:,:) = ( zintn * psnwice_mass(:,:) + zintb * psnwice_mass_b(:,:) ) * r1_rau0 
     749         zsnwiceload(:,:) = ( zintn * psnwice_mass(:,:) + zintb * psnwice_mass_b(:,:) ) * r1_rho0 
    750750         ! 
    751751      ELSE 
     
    937937               ! In case snow load is in excess that would lead to transformation from snow to ice 
    938938               ! Then, transfer the snow excess into the ice (different from icethd_dh) 
    939                zdh = MAX( 0._wp, ( rhos * ph_s(ji,jl) + ( rhoi - rau0 ) * ph_i(ji,jl) ) * r1_rau0 )  
     939               zdh = MAX( 0._wp, ( rhos * ph_s(ji,jl) + ( rhoi - rho0 ) * ph_i(ji,jl) ) * r1_rho0 )  
    940940               ! recompute h_i, h_s avoiding out of bounds values 
    941941               ph_i(ji,jl) = MIN( hi_max(jl), ph_i(ji,jl) + zdh ) 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/ICE/icewri.F90

    r12377 r12443  
    8787      ! Standard outputs 
    8888      !----------------- 
    89       zrho1 = ( rau0 - rhoi ) * r1_rau0 ; zrho2 = rhos * r1_rau0 
     89      zrho1 = ( rho0 - rhoi ) * r1_rho0 ; zrho2 = rhos * r1_rho0 
    9090      ! masks 
    9191      CALL iom_put( 'icemask'  , zmsk00 )   ! ice mask 0% 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/NST/agrif_oce_update.F90

    r12406 r12443  
    361361                           ztnu = tabres_child(ji,jj,jk,jn) * e3t(ji,jj,jk,Kmm_a) 
    362362                           ztno = ts(ji,jj,jk,jn,Kmm_a) * e3t(ji,jj,jk,Krhs_a) 
    363                            ts(ji,jj,jk,jn,Kbb_a) = ( ztb + atfp * ( ztnu - ztno) )  &  
     363                           ts(ji,jj,jk,jn,Kbb_a) = ( ztb + rn_atfp * ( ztnu - ztno) )  &  
    364364                                     &        * tmask(ji,jj,jk) / e3t(ji,jj,jk,Kbb_a) 
    365365                        ENDIF 
     
    432432                           ztnu = tabres(ji,jj,jk,jn) 
    433433                           ztno = ts(ji,jj,jk,jn,Kmm_a) * e3t(ji,jj,jk,Krhs_a) 
    434                            ts(ji,jj,jk,jn,Kbb_a) = ( ztb + atfp * ( ztnu - ztno) )  &  
     434                           ts(ji,jj,jk,jn,Kbb_a) = ( ztb + rn_atfp * ( ztnu - ztno) )  &  
    435435                                     &        * tmask(ji,jj,jk) / e3t(ji,jj,jk,Kbb_a) 
    436436                        ENDIF 
     
    555555                     zuno = uu(ji,jj,jk,Kmm_a) * e3u(ji,jj,jk,Krhs_a) 
    556556                     zunu = tabres_child(ji,jj,jk) * e3u(ji,jj,jk,Kmm_a) 
    557                      uu(ji,jj,jk,Kbb_a) = ( zub + atfp * ( zunu - zuno) ) &       
     557                     uu(ji,jj,jk,Kbb_a) = ( zub + rn_atfp * ( zunu - zuno) ) &       
    558558                                    & * umask(ji,jj,jk) / e3u(ji,jj,jk,Kbb_a) 
    559559                  ENDIF 
     
    601601                     zuno = uu(ji,jj,jk,Kmm_a) * e3u(ji,jj,jk,Krhs_a) 
    602602                     zunu = tabres(ji,jj,jk,1) 
    603                      uu(ji,jj,jk,Kbb_a) = ( zub + atfp * ( zunu - zuno) ) &       
     603                     uu(ji,jj,jk,Kbb_a) = ( zub + rn_atfp * ( zunu - zuno) ) &       
    604604                                    & * umask(ji,jj,jk) / e3u(ji,jj,jk,Kbb_a) 
    605605                  ENDIF 
     
    755755                     zvno = vv(ji,jj,jk,Kmm_a) * e3v(ji,jj,jk,Krhs_a) 
    756756                     zvnu = tabres_child(ji,jj,jk) * e3v(ji,jj,jk,Kmm_a) 
    757                      vv(ji,jj,jk,Kbb_a) = ( zvb + atfp * ( zvnu - zvno) ) &       
     757                     vv(ji,jj,jk,Kbb_a) = ( zvb + rn_atfp * ( zvnu - zvno) ) &       
    758758                                    & * vmask(ji,jj,jk) / e3v(ji,jj,jk,Kbb_a) 
    759759                  ENDIF 
     
    805805                     zvno = vv(ji,jj,jk,Kmm_a) * e3v(ji,jj,jk,Krhs_a) 
    806806                     zvnu = tabres(ji,jj,jk,1) 
    807                      vv(ji,jj,jk,Kbb_a) = ( zvb + atfp * ( zvnu - zvno) ) &       
     807                     vv(ji,jj,jk,Kbb_a) = ( zvb + rn_atfp * ( zvnu - zvno) ) &       
    808808                                    & * vmask(ji,jj,jk) / e3v(ji,jj,jk,Kbb_a) 
    809809                  ENDIF 
     
    909909                  IF (.NOT.(lk_agrif_fstep.AND.(l_1st_euler))) THEN ! Add asselin part 
    910910                     zcorr = (tabres(ji,jj) - uu_b(ji,jj,Kmm_a) * hu(ji,jj,Krhs_a)) * r1_hu(ji,jj,Kbb_a) 
    911                      uu_b(ji,jj,Kbb_a) = uu_b(ji,jj,Kbb_a) + atfp * zcorr * umask(ji,jj,1) 
     911                     uu_b(ji,jj,Kbb_a) = uu_b(ji,jj,Kbb_a) + rn_atfp * zcorr * umask(ji,jj,1) 
    912912                  END IF 
    913913               ENDIF     
     
    975975                  IF (.NOT.(lk_agrif_fstep.AND.(l_1st_euler))) THEN ! Add asselin part 
    976976                     zcorr = (tabres(ji,jj) - vv_b(ji,jj,Kmm_a) * hv(ji,jj,Krhs_a)) * r1_hv(ji,jj,Kbb_a) 
    977                      vv_b(ji,jj,Kbb_a) = vv_b(ji,jj,Kbb_a) + atfp * zcorr * vmask(ji,jj,1) 
     977                     vv_b(ji,jj,Kbb_a) = vv_b(ji,jj,Kbb_a) + rn_atfp * zcorr * vmask(ji,jj,1) 
    978978                  END IF 
    979979               ENDIF               
     
    10251025               DO ji=i1,i2 
    10261026                  ssh(ji,jj,Kbb_a) =   ssh(ji,jj,Kbb_a) & 
    1027                         & + atfp * ( tabres(ji,jj) - ssh(ji,jj,Kmm_a) ) * tmask(ji,jj,1) 
     1027                        & + rn_atfp * ( tabres(ji,jj) - ssh(ji,jj,Kmm_a) ) * tmask(ji,jj,1) 
    10281028               END DO 
    10291029            END DO 
     
    11191119               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 
    1121                IF (.NOT.(lk_agrif_fstep.AND.(l_1st_euler))) ssh(i1  ,jj,Kbb_a) = ssh(i1  ,jj,Kbb_a) + atfp * zcor 
     1121               IF (.NOT.(lk_agrif_fstep.AND.(l_1st_euler))) ssh(i1  ,jj,Kbb_a) = ssh(i1  ,jj,Kbb_a) + rn_atfp * zcor 
    11221122            END DO 
    11231123         ENDIF 
     
    11261126               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 
    1128                IF (.NOT.(lk_agrif_fstep.AND.(l_1st_euler))) ssh(i2+1,jj,Kbb_a) = ssh(i2+1,jj,Kbb_a) + atfp * zcor 
     1128               IF (.NOT.(lk_agrif_fstep.AND.(l_1st_euler))) ssh(i2+1,jj,Kbb_a) = ssh(i2+1,jj,Kbb_a) + rn_atfp * zcor 
    11291129            END DO 
    11301130         ENDIF 
     
    12071207               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 
    1209                IF (.NOT.(lk_agrif_fstep.AND.(l_1st_euler))) ssh(ji,j1  ,Kbb_a) = ssh(ji,j1,Kbb_a) + atfp * zcor 
     1209               IF (.NOT.(lk_agrif_fstep.AND.(l_1st_euler))) ssh(ji,j1  ,Kbb_a) = ssh(ji,j1,Kbb_a) + rn_atfp * zcor 
    12101210            END DO 
    12111211         ENDIF 
     
    12141214               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 
    1216                IF (.NOT.(lk_agrif_fstep.AND.(l_1st_euler))) ssh(ji,j2+1,Kbb_a) = ssh(ji,j2+1,Kbb_a) + atfp * zcor 
     1216               IF (.NOT.(lk_agrif_fstep.AND.(l_1st_euler))) ssh(ji,j2+1,Kbb_a) = ssh(ji,j2+1,Kbb_a) + rn_atfp * zcor 
    12171217            END DO 
    12181218         ENDIF 
     
    13641364                  DO ji=i1,i2 
    13651365                     e3t(ji,jj,jk,Kbb_a) =  e3t(ji,jj,jk,Kbb_a) & 
    1366                            & + atfp * ( ptab(ji,jj,jk) - e3t(ji,jj,jk,Kmm_a) ) 
     1366                           & + rn_atfp * ( ptab(ji,jj,jk) - e3t(ji,jj,jk,Kmm_a) ) 
    13671367                  END DO 
    13681368               END DO 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/NST/agrif_top_update.F90

    r12397 r12443  
    135135                           ztnu = tabres_child(ji,jj,jk,jn) * e3t(ji,jj,jk,Kmm_a) 
    136136                           ztno = tr(ji,jj,jk,jn,Kmm_a) * e3t(ji,jj,jk,Krhs_a) 
    137                            tr(ji,jj,jk,jn,Kbb_a) = ( ztb + atfp * ( ztnu - ztno) )  &  
     137                           tr(ji,jj,jk,jn,Kbb_a) = ( ztb + rn_atfp * ( ztnu - ztno) )  &  
    138138                                     &        * tmask(ji,jj,jk) / e3t(ji,jj,jk,Kbb_a) 
    139139                        ENDIF 
     
    209209                           ztnu = tabres(ji,jj,jk,jn) 
    210210                           ztno = tr(ji,jj,jk,jn,Kmm_a) * e3t(ji,jj,jk,Krhs_a) 
    211                            tr(ji,jj,jk,jn,Kbb_a) = ( ztb + atfp * ( ztnu - ztno) )  &  
     211                           tr(ji,jj,jk,jn,Kbb_a) = ( ztb + rn_atfp * ( ztnu - ztno) )  &  
    212212                                     &        * tmask(ji,jj,jk) / e3t(ji,jj,jk,Kbb_a) 
    213213                        ENDIF 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OCE/BDY/bdyice.F90

    r11536 r12443  
    179179 
    180180            ! Then, a) transfer the snow excess into the ice (different from icethd_dh) 
    181             zdh = MAX( 0._wp, ( rhos * h_s(ji,jj,jl) + ( rhoi - rau0 ) * h_i(ji,jj,jl) ) * r1_rau0 ) 
     181            zdh = MAX( 0._wp, ( rhos * h_s(ji,jj,jl) + ( rhoi - rho0 ) * h_i(ji,jj,jl) ) * r1_rho0 ) 
    182182            ! Or, b) transfer all the snow into ice (if incoming ice is likely to melt as it comes into a warmer environment) 
    183183            !zdh = MAX( 0._wp, h_s(ji,jj,jl) * rhos / rhoi ) 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OCE/BDY/bdytides.F90

    r12406 r12443  
    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) ) * rn_Dt 
     299         z_arg = ( REAL(kt, wp) + ( REAL(kit, wp) + zt_offset - 1. ) / REAL(nn_e, wp) ) * rn_Dt 
    300300      ELSE                               
    301301         z_arg = ( REAL(kt, wp) + zt_offset ) * rn_Dt 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OCE/BDY/bdyvol.F90

    r12377 r12443  
    7777      ! Calculate the cumulate surface Flux z_cflxemp (m3/s) over all the domain 
    7878      ! ----------------------------------------------------------------------- 
    79       IF ( kc == 1 ) z_cflxemp = glob_sum( 'bdyvol', ( emp(:,:) - rnf(:,:) + fwfisf_cav(:,:) + fwfisf_par(:,:) ) * bdytmask(:,:) * e1e2t(:,:)  ) / rau0 
     79      IF ( kc == 1 ) z_cflxemp = glob_sum( 'bdyvol', ( emp(:,:) - rnf(:,:) + fwfisf_cav(:,:) + fwfisf_par(:,:) ) * bdytmask(:,:) * e1e2t(:,:)  ) / rho0 
    8080 
    8181      ! Compute bdy surface each cycle if non linear free surface 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OCE/DIA/diaar5.F90

    r12377 r12443  
    103103         END DO 
    104104         CALL iom_put( 'volcello'  , zrhd(:,:,:)  )  ! WARNING not consistent with CMIP DR where volcello is at ca. 2000 
    105          CALL iom_put( 'masscello' , rau0 * e3t(:,:,:,Kmm) * tmask(:,:,:) )  ! ocean mass 
     105         CALL iom_put( 'masscello' , rho0 * e3t(:,:,:,Kmm) * tmask(:,:,:) )  ! ocean mass 
    106106      ENDIF  
    107107      ! 
     
    181181         CALL iom_put( 'sshsteric', zssh_steric ) 
    182182         !                                         ! ocean bottom pressure 
    183          zztmp = rau0 * grav * 1.e-4_wp               ! recover pressure from pressure anomaly and cover to dbar = 1.e4 Pa 
     183         zztmp = rho0 * grav * 1.e-4_wp               ! recover pressure from pressure anomaly and cover to dbar = 1.e4 Pa 
    184184         zbotpres(:,:) = zztmp * ( zbotpres(:,:) + ssh(:,:,Kmm) + thick0(:,:) ) 
    185185         CALL iom_put( 'botpres', zbotpres ) 
     
    213213         ztemp = glob_sum( 'diaar5', ztsn(:,:,1,jp_tem) ) 
    214214         zsal  = glob_sum( 'diaar5', ztsn(:,:,1,jp_sal) ) 
    215          zmass = rau0 * ( zarho + zvol )       
     215         zmass = rho0 * ( zarho + zvol )       
    216216         ! 
    217217         CALL iom_put( 'masstot', zmass ) 
     
    251251               z2d(:,:) = 0._wp 
    252252               DO_3D_11_11( 1, jpkm1 ) 
    253                   z2d(ji,jj) = z2d(ji,jj) + rau0 * e3t(ji,jj,jk,Kmm) *  ztpot(ji,jj,jk) 
     253                  z2d(ji,jj) = z2d(ji,jj) + rho0 * e3t(ji,jj,jk,Kmm) *  ztpot(ji,jj,jk) 
    254254               END_3D 
    255255               CALL iom_put( 'tosmint_pot', z2d )  
     
    285285          ELSE 
    286286            DO_3D_11_11( 1, jpk ) 
    287                zpe(ji,jj) = zpe(ji,jj) + avt(ji,jj,jk) * MIN(0._wp,rn2(ji,jj,jk)) * rau0 * e3w(ji,jj,jk,Kmm) 
     287               zpe(ji,jj) = zpe(ji,jj) + avt(ji,jj,jk) * MIN(0._wp,rn2(ji,jj,jk)) * rho0 * e3w(ji,jj,jk,Kmm) 
    288288            END_3D 
    289289         ENDIF 
     
    325325       CALL lbc_lnk( 'diaar5', z2d, 'U', -1. ) 
    326326       IF( cptr == 'adv' ) THEN 
    327           IF( ktra == jp_tem ) CALL iom_put( 'uadv_heattr' , rau0_rcp * z2d )  ! advective heat transport in i-direction 
    328           IF( ktra == jp_sal ) CALL iom_put( 'uadv_salttr' , rau0     * z2d )  ! advective salt transport in i-direction 
     327          IF( ktra == jp_tem ) CALL iom_put( 'uadv_heattr' , rho0_rcp * z2d )  ! advective heat transport in i-direction 
     328          IF( ktra == jp_sal ) CALL iom_put( 'uadv_salttr' , rho0     * z2d )  ! advective salt transport in i-direction 
    329329       ENDIF 
    330330       IF( cptr == 'ldf' ) THEN 
    331           IF( ktra == jp_tem ) CALL iom_put( 'udiff_heattr' , rau0_rcp * z2d ) ! diffusive heat transport in i-direction 
    332           IF( ktra == jp_sal ) CALL iom_put( 'udiff_salttr' , rau0     * z2d ) ! diffusive salt transport in i-direction 
     331          IF( ktra == jp_tem ) CALL iom_put( 'udiff_heattr' , rho0_rcp * z2d ) ! diffusive heat transport in i-direction 
     332          IF( ktra == jp_sal ) CALL iom_put( 'udiff_salttr' , rho0     * z2d ) ! diffusive salt transport in i-direction 
    333333       ENDIF 
    334334       ! 
     
    339339       CALL lbc_lnk( 'diaar5', z2d, 'V', -1. ) 
    340340       IF( cptr == 'adv' ) THEN 
    341           IF( ktra == jp_tem ) CALL iom_put( 'vadv_heattr' , rau0_rcp * z2d )  ! advective heat transport in j-direction 
    342           IF( ktra == jp_sal ) CALL iom_put( 'vadv_salttr' , rau0     * z2d )  ! advective salt transport in j-direction 
     341          IF( ktra == jp_tem ) CALL iom_put( 'vadv_heattr' , rho0_rcp * z2d )  ! advective heat transport in j-direction 
     342          IF( ktra == jp_sal ) CALL iom_put( 'vadv_salttr' , rho0     * z2d )  ! advective salt transport in j-direction 
    343343       ENDIF 
    344344       IF( cptr == 'ldf' ) THEN 
    345           IF( ktra == jp_tem ) CALL iom_put( 'vdiff_heattr' , rau0_rcp * z2d ) ! diffusive heat transport in j-direction 
    346           IF( ktra == jp_sal ) CALL iom_put( 'vdiff_salttr' , rau0     * z2d ) ! diffusive salt transport in j-direction 
     345          IF( ktra == jp_tem ) CALL iom_put( 'vdiff_heattr' , rho0_rcp * z2d ) ! diffusive heat transport in j-direction 
     346          IF( ktra == jp_sal ) CALL iom_put( 'vdiff_salttr' , rho0     * z2d ) ! diffusive salt transport in j-direction 
    347347       ENDIF 
    348348           
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OCE/DIA/diadct.F90

    r12377 r12443  
    676676                  zsn   = interp(Kmm,k%I,k%J,jk,'V',ts(:,:,:,jp_sal,Kmm) )  
    677677                  zrhop = interp(Kmm,k%I,k%J,jk,'V',rhop)  
    678                   zrhoi = interp(Kmm,k%I,k%J,jk,'V',rhd*rau0+rau0)  
     678                  zrhoi = interp(Kmm,k%I,k%J,jk,'V',rhd*rho0+rho0)  
    679679                  zsshn =  0.5*( ssh(k%I,k%J,Kmm) + ssh(k%I,k%J+1,Kmm)    ) * vmask(k%I,k%J,1)  
    680680               CASE(2,3)  
     
    682682                  zsn   = interp(Kmm,k%I,k%J,jk,'U',ts(:,:,:,jp_sal,Kmm) )  
    683683                  zrhop = interp(Kmm,k%I,k%J,jk,'U',rhop)  
    684                   zrhoi = interp(Kmm,k%I,k%J,jk,'U',rhd*rau0+rau0)  
     684                  zrhoi = interp(Kmm,k%I,k%J,jk,'U',rhd*rho0+rho0)  
    685685                  zsshn =  0.5*( ssh(k%I,k%J,Kmm) + ssh(k%I+1,k%J,Kmm)    ) * umask(k%I,k%J,1)   
    686686               END SELECT  
     
    849849                 zsn   = interp(Kmm,k%I,k%J,jk,'V',ts(:,:,:,jp_sal,Kmm) )  
    850850                 zrhop = interp(Kmm,k%I,k%J,jk,'V',rhop)  
    851                  zrhoi = interp(Kmm,k%I,k%J,jk,'V',rhd*rau0+rau0)  
     851                 zrhoi = interp(Kmm,k%I,k%J,jk,'V',rhd*rho0+rho0)  
    852852 
    853853              CASE(2,3)  
     
    855855                 zsn   = interp(Kmm,k%I,k%J,jk,'U',ts(:,:,:,jp_sal,Kmm) )  
    856856                 zrhop = interp(Kmm,k%I,k%J,jk,'U',rhop)  
    857                  zrhoi = interp(Kmm,k%I,k%J,jk,'U',rhd*rau0+rau0)  
     857                 zrhoi = interp(Kmm,k%I,k%J,jk,'U',rhd*rho0+rho0)  
    858858                 zsshn =  0.5*( ssh(k%I,k%J,Kmm)    + ssh(k%I+1,k%J,Kmm)    ) * umask(k%I,k%J,1)   
    859859              END SELECT  
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OCE/DIA/diahsb.F90

    r12406 r12443  
    9191      ! 1 - Trends due to forcing ! 
    9292      ! ------------------------- ! 
    93       z_frc_trd_v = r1_rau0 * glob_sum( 'diahsb', - ( emp(:,:) - rnf(:,:) + fwfisf_cav(:,:) + fwfisf_par(:,:) ) * surf(:,:) )   ! volume fluxes 
     93      z_frc_trd_v = r1_rho0 * glob_sum( 'diahsb', - ( emp(:,:) - rnf(:,:) + fwfisf_cav(:,:) + fwfisf_par(:,:) ) * surf(:,:) )   ! volume fluxes 
    9494      z_frc_trd_t =           glob_sum( 'diahsb', sbc_tsc(:,:,jp_tem) * surf(:,:) )                       ! heat fluxes 
    9595      z_frc_trd_s =           glob_sum( 'diahsb', sbc_tsc(:,:,jp_sal) * surf(:,:) )                       ! salt fluxes 
     
    101101         &                          + glob_sum( 'diahsb', ( risf_cav_tsc(:,:,jp_tem) + risf_par_tsc(:,:,jp_tem) ) * surf(:,:) ) 
    102102      !                    ! Add penetrative solar radiation 
    103       IF( ln_traqsr )   z_frc_trd_t = z_frc_trd_t + r1_rau0_rcp * glob_sum( 'diahsb', qsr     (:,:) * surf(:,:) ) 
     103      IF( ln_traqsr )   z_frc_trd_t = z_frc_trd_t + r1_rho0_rcp * glob_sum( 'diahsb', qsr     (:,:) * surf(:,:) ) 
    104104      !                    ! Add geothermal heat flux 
    105105      IF( ln_trabbc )   z_frc_trd_t = z_frc_trd_t +               glob_sum( 'diahsb', qgh_trd0(:,:) * surf(:,:) ) 
     
    197197 
    198198      CALL iom_put(   'bgfrcvol' , frc_v    * 1.e-9    )              ! vol - surface forcing (km3)  
    199       CALL iom_put(   'bgfrctem' , frc_t    * rau0 * rcp * 1.e-20 )   ! hc  - surface forcing (1.e20 J)  
    200       CALL iom_put(   'bgfrchfx' , frc_t    * rau0 * rcp /  &         ! hc  - surface forcing (W/m2)  
     199      CALL iom_put(   'bgfrctem' , frc_t    * rho0 * rcp * 1.e-20 )   ! hc  - surface forcing (1.e20 J)  
     200      CALL iom_put(   'bgfrchfx' , frc_t    * rho0 * rcp /  &         ! hc  - surface forcing (W/m2)  
    201201         &                       ( surf_tot * kt * rn_Dt )        ) 
    202202      CALL iom_put(   'bgfrcsal' , frc_s    * 1.e-9    )              ! sc  - surface forcing (psu*km3)  
     
    205205         CALL iom_put( 'bgtemper' , zdiff_hc / zvol_tot )              ! Temperature drift     (C)  
    206206         CALL iom_put( 'bgsaline' , zdiff_sc / zvol_tot )              ! Salinity    drift     (PSU) 
    207          CALL iom_put( 'bgheatco' , zdiff_hc * 1.e-20 * rau0 * rcp )   ! Heat content drift    (1.e20 J)  
    208          CALL iom_put( 'bgheatfx' , zdiff_hc * rau0 * rcp /  &         ! Heat flux drift       (W/m2)  
     207         CALL iom_put( 'bgheatco' , zdiff_hc * 1.e-20 * rho0 * rcp )   ! Heat content drift    (1.e20 J)  
     208         CALL iom_put( 'bgheatfx' , zdiff_hc * rho0 * rcp /  &         ! Heat flux drift       (W/m2)  
    209209            &                       ( surf_tot * kt * rn_Dt )        ) 
    210210         CALL iom_put( 'bgsaltco' , zdiff_sc * 1.e-9    )              ! Salt content drift    (psu*km3) 
     
    225225         CALL iom_put( 'bgtemper' , zdiff_hc1 / zvol_tot)              ! Heat content drift    (C)  
    226226         CALL iom_put( 'bgsaline' , zdiff_sc1 / zvol_tot)              ! Salt content drift    (PSU) 
    227          CALL iom_put( 'bgheatco' , zdiff_hc1 * 1.e-20 * rau0 * rcp )  ! Heat content drift    (1.e20 J)  
    228          CALL iom_put( 'bgheatfx' , zdiff_hc1 * rau0 * rcp /  &        ! Heat flux drift       (W/m2)  
     227         CALL iom_put( 'bgheatco' , zdiff_hc1 * 1.e-20 * rho0 * rcp )  ! Heat content drift    (1.e20 J)  
     228         CALL iom_put( 'bgheatfx' , zdiff_hc1 * rho0 * rcp /  &        ! Heat flux drift       (W/m2)  
    229229            &                       ( surf_tot * kt * rn_Dt )         ) 
    230230         CALL iom_put( 'bgsaltco' , zdiff_sc1 * 1.e-9    )             ! Salt content drift    (psu*km3) 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OCE/DIA/diahth.F90

    r12377 r12443  
    261261            zzdep = 300. 
    262262            CALL  dia_hth_htc( Kmm, zzdep, ts(:,:,:,jp_tem,Kmm), htc3 ) 
    263             CALL iom_put( 'hc300', rau0_rcp * htc3 )  ! vertically integrated heat content (J/m2) 
     263            CALL iom_put( 'hc300', rho0_rcp * htc3 )  ! vertically integrated heat content (J/m2) 
    264264         ENDIF 
    265265         ! 
     
    270270            zzdep = 700. 
    271271            CALL  dia_hth_htc( Kmm, zzdep, ts(:,:,:,jp_tem,Kmm), htc7 ) 
    272             CALL iom_put( 'hc700', rau0_rcp * htc7 )  ! vertically integrated heat content (J/m2) 
     272            CALL iom_put( 'hc700', rho0_rcp * htc7 )  ! vertically integrated heat content (J/m2) 
    273273   
    274274         ENDIF 
     
    280280            zzdep = 2000. 
    281281            CALL  dia_hth_htc( Kmm, zzdep, ts(:,:,:,jp_tem,Kmm), htc20 ) 
    282             CALL iom_put( 'hc2000', rau0_rcp * htc20 )  ! vertically integrated heat content (J/m2)   
     282            CALL iom_put( 'hc2000', rho0_rcp * htc20 )  ! vertically integrated heat content (J/m2)   
    283283         ENDIF 
    284284         ! 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OCE/DIA/diaptr.F90

    r12377 r12443  
    5050 
    5151   REAL(wp) ::   rc_sv    = 1.e-6_wp   ! conversion from m3/s to Sverdrup 
    52    REAL(wp) ::   rc_pwatt = 1.e-15_wp  ! conversion from W    to PW (further x rau0 x Cp) 
    53    REAL(wp) ::   rc_ggram = 1.e-9_wp   ! conversion from g    to Gg  (further x rau0) 
     52   REAL(wp) ::   rc_pwatt = 1.e-15_wp  ! conversion from W    to PW (further x rho0 x Cp) 
     53   REAL(wp) ::   rc_ggram = 1.e-9_wp   ! conversion from g    to Gg  (further x rho0) 
    5454 
    5555   REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: btmsk   ! T-point basin interior masks 
     
    346346         IF( dia_ptr_alloc() /= 0 )   CALL ctl_stop( 'STOP', 'dia_ptr_init : unable to allocate arrays' ) 
    347347 
    348          rc_pwatt = rc_pwatt * rau0_rcp          ! conversion from K.s-1 to PetaWatt 
    349          rc_ggram = rc_ggram * rau0              ! conversion from m3/s to Gg/s 
     348         rc_pwatt = rc_pwatt * rho0_rcp          ! conversion from K.s-1 to PetaWatt 
     349         rc_ggram = rc_ggram * rho0              ! conversion from m3/s to Gg/s 
    350350 
    351351         IF( lk_mpp )   CALL mpp_ini_znl( numout )     ! Define MPI communicator for zonal sum 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OCE/DIA/diawri.F90

    r12406 r12443  
    173173 
    174174      IF ( iom_use("taubot") ) THEN                ! bottom stress 
    175          zztmp = rau0 * 0.25 
     175         zztmp = rho0 * 0.25 
    176176         z2d(:,:) = 0._wp 
    177177         DO_2D_00_00 
     
    212212      IF( iom_use('w_masstr') .OR. iom_use('w_masstr2') ) THEN   ! vertical mass transport & its square value 
    213213         ! Caution: in the VVL case, it only correponds to the baroclinic mass transport. 
    214          z2d(:,:) = rau0 * e1e2t(:,:) 
     214         z2d(:,:) = rho0 * e1e2t(:,:) 
    215215         DO jk = 1, jpk 
    216216            z3d(:,:,jk) = ww(:,:,jk) * z2d(:,:) 
     
    249249            z2d(ji,jj) = z2d(ji,jj) + e3t(ji,jj,jk,Kmm) * ts(ji,jj,jk,jp_tem,Kmm) * tmask(ji,jj,jk) 
    250250         END_3D 
    251          CALL iom_put( "heatc", rau0_rcp * z2d )   ! vertically integrated heat content (J/m2) 
     251         CALL iom_put( "heatc", rho0_rcp * z2d )   ! vertically integrated heat content (J/m2) 
    252252      ENDIF 
    253253 
     
    257257            z2d(ji,jj) = z2d(ji,jj) + e3t(ji,jj,jk,Kmm) * ts(ji,jj,jk,jp_sal,Kmm) * tmask(ji,jj,jk) 
    258258         END_3D 
    259          CALL iom_put( "saltc", rau0 * z2d )          ! vertically integrated salt content (PSU*kg/m2) 
     259         CALL iom_put( "saltc", rho0 * z2d )          ! vertically integrated salt content (PSU*kg/m2) 
    260260      ENDIF 
    261261      ! 
     
    279279         z2d(:,:) = 0.e0 
    280280         DO jk = 1, jpkm1 
    281             z3d(:,:,jk) = rau0 * uu(:,:,jk,Kmm) * e2u(:,:) * e3u(:,:,jk,Kmm) * umask(:,:,jk) 
     281            z3d(:,:,jk) = rho0 * uu(:,:,jk,Kmm) * e2u(:,:) * e3u(:,:,jk,Kmm) * umask(:,:,jk) 
    282282            z2d(:,:) = z2d(:,:) + z3d(:,:,jk) 
    283283         END DO 
     
    308308         z3d(:,:,jpk) = 0.e0 
    309309         DO jk = 1, jpkm1 
    310             z3d(:,:,jk) = rau0 * vv(:,:,jk,Kmm) * e1v(:,:) * e3v(:,:,jk,Kmm) * vmask(:,:,jk) 
     310            z3d(:,:,jk) = rho0 * vv(:,:,jk,Kmm) * e1v(:,:) * e3v(:,:,jk,Kmm) * vmask(:,:,jk) 
    311311         END DO 
    312312         CALL iom_put( "v_masstr", z3d )              ! mass transport in j-direction 
     
    337337         END_3D 
    338338         CALL lbc_lnk( 'diawri', z2d, 'T', -1. ) 
    339          CALL iom_put( "tosmint", rau0 * z2d )        ! Vertical integral of temperature 
     339         CALL iom_put( "tosmint", rho0 * z2d )        ! Vertical integral of temperature 
    340340      ENDIF 
    341341      IF( iom_use("somint") ) THEN 
     
    345345         END_3D 
    346346         CALL lbc_lnk( 'diawri', z2d, 'T', -1. ) 
    347          CALL iom_put( "somint", rau0 * z2d )         ! Vertical integral of salinity 
     347         CALL iom_put( "somint", rho0 * z2d )         ! Vertical integral of salinity 
    348348      ENDIF 
    349349 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OCE/DOM/dom_oce.F90

    r12424 r12443  
    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_Dt         !: 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 (rn_Dt) 
     51   INTEGER,  PUBLIC :: nn_e          !: 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 
    5454 
    55  
    56    !                                   !! old non-DOCTOR names still used in the model 
    57    REAL(wp), PUBLIC ::   atfp           !: asselin time filter parameter 
    5855 
    5956   !                                   !!! associated variables 
    6057   LOGICAL , PUBLIC ::   l_1st_euler    !: Euler 1st time-step flag (=T if ln_restart=F or ln_1st_euler=T) 
    6158   REAL(wp), PUBLIC ::   rDt, r1_Dt     !: Current model timestep and reciprocal 
    62                                         !: = 2 * rn_Dt if leapfrog and l_1st_euler = F 
    63                                         !:  =     rn_Dt if leapfrog and l_1st_euler = T 
    64                                         !:  =     rn_Dt if RK3 
     59                                        !: rDt = 2 * rn_Dt if leapfrog and l_1st_euler = F 
     60                                        !:     =     rn_Dt if leapfrog and l_1st_euler = T 
     61                                        !:     =     rn_Dt if RK3 
    6562 
    6663   !!---------------------------------------------------------------------- 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OCE/DOM/domain.F90

    r12424 r12443  
    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_Dt      = ', rn_Dt 
     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 
    413413      ENDIF 
    414414      ! 
    415       !          ! conversion DOCTOR names into model names (this should disappear soon) 
    416       atfp = rn_atfp 
    417415      !! Initialise current model timestep rDt = 2*rn_Dt if MLF or rDt = rn_Dt if RK3 
    418416      rDt  = 2._wp * rn_Dt 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OCE/DOM/domvvl.F90

    r12424 r12443  
    621621         ELSE 
    622622            tilde_e3t_b(:,:,:) = tilde_e3t_n(:,:,:) &  
    623             &         + atfp * ( tilde_e3t_b(:,:,:) - 2.0_wp * tilde_e3t_n(:,:,:) + tilde_e3t_a(:,:,:) ) 
     623            &         + rn_atfp * ( tilde_e3t_b(:,:,:) - 2.0_wp * tilde_e3t_n(:,:,:) + tilde_e3t_a(:,:,:) ) 
    624624         ENDIF 
    625625         tilde_e3t_n(:,:,:) = tilde_e3t_a(:,:,:) 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OCE/DOM/phycst.F90

    r10068 r12443  
    3939   REAL(wp), PUBLIC ::   rt0      = 273.15_wp        !: freezing point of fresh water [Kelvin] 
    4040 
    41    REAL(wp), PUBLIC ::   rau0                        !: volumic mass of reference     [kg/m3] 
    42    REAL(wp), PUBLIC ::   r1_rau0                     !: = 1. / rau0                   [m3/kg] 
     41   REAL(wp), PUBLIC ::   rho0                        !: volumic mass of reference     [kg/m3] 
     42   REAL(wp), PUBLIC ::   r1_rho0                     !: = 1. / rho0                   [m3/kg] 
    4343   REAL(wp), PUBLIC ::   rcp                         !: ocean specific heat           [J/Kelvin] 
    4444   REAL(wp), PUBLIC ::   r1_rcp                      !: = 1. / rcp                    [Kelvin/J] 
    45    REAL(wp), PUBLIC ::   rau0_rcp                    !: = rau0 * rcp  
    46    REAL(wp), PUBLIC ::   r1_rau0_rcp                 !: = 1. / ( rau0 * rcp ) 
     45   REAL(wp), PUBLIC ::   rho0_rcp                    !: = rho0 * rcp  
     46   REAL(wp), PUBLIC ::   r1_rho0_rcp                 !: = 1. / ( rho0 * rcp ) 
    4747 
    4848   REAL(wp), PUBLIC ::   emic     =    0.97_wp       !: emissivity of snow or ice (not used?) 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OCE/DYN/dynatf.F90

    r12424 r12443  
    8787      !!             arrays to start the next time step: 
    8888      !!                (puu(Kmm),pvv(Kmm)) = (puu(Kmm),pvv(Kmm))  
    89       !!                                    + atfp [ (puu(Kbb),pvv(Kbb)) + (puu(Kaa),pvv(Kaa)) - 2 (puu(Kmm),pvv(Kmm)) ] 
     89      !!                                    + rn_atfp [ (puu(Kbb),pvv(Kbb)) + (puu(Kaa),pvv(Kaa)) - 2 (puu(Kmm),pvv(Kmm)) ] 
    9090      !!             Note that with flux form advection and non linear free surface, 
    9191      !!             the time filter is applied on thickness weighted velocity. 
     
    181181            !                             ! =============! 
    182182            DO_3D_11_11( 1, jpkm1 ) 
    183                puu(ji,jj,jk,Kmm) = puu(ji,jj,jk,Kmm) + atfp * ( puu(ji,jj,jk,Kbb) - 2._wp * puu(ji,jj,jk,Kmm) + puu(ji,jj,jk,Kaa) ) 
    184                pvv(ji,jj,jk,Kmm) = pvv(ji,jj,jk,Kmm) + atfp * ( pvv(ji,jj,jk,Kbb) - 2._wp * pvv(ji,jj,jk,Kmm) + pvv(ji,jj,jk,Kaa) ) 
     183               puu(ji,jj,jk,Kmm) = puu(ji,jj,jk,Kmm) + rn_atfp * ( puu(ji,jj,jk,Kbb) - 2._wp * puu(ji,jj,jk,Kmm) + puu(ji,jj,jk,Kaa) ) 
     184               pvv(ji,jj,jk,Kmm) = pvv(ji,jj,jk,Kmm) + rn_atfp * ( pvv(ji,jj,jk,Kbb) - 2._wp * pvv(ji,jj,jk,Kmm) + pvv(ji,jj,jk,Kaa) ) 
    185185            END_3D 
    186186            !                             ! ================! 
     
    191191            ALLOCATE( ze3t_f(jpi,jpj,jpk), zwfld(jpi,jpj) ) 
    192192            DO jk = 1, jpkm1 
    193                ze3t_f(:,:,jk) = pe3t(:,:,jk,Kmm) + atfp * ( pe3t(:,:,jk,Kbb) - 2._wp * pe3t(:,:,jk,Kmm) + pe3t(:,:,jk,Kaa) ) 
     193               ze3t_f(:,:,jk) = pe3t(:,:,jk,Kmm) + rn_atfp * ( pe3t(:,:,jk,Kbb) - 2._wp * pe3t(:,:,jk,Kmm) + pe3t(:,:,jk,Kaa) ) 
    194194            END DO 
    195195            ! Add volume filter correction: compatibility with tracer advection scheme 
    196196            ! => time filter + conservation correction 
    197             zcoef = atfp * rn_Dt * r1_rau0 
     197            zcoef = rn_atfp * rn_Dt * r1_rho0 
    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 * rn_Dt ) 
     209            IF ( ln_isf ) CALL isf_dynatf( kt, Kmm, ze3t_f, rn_atfp * rn_Dt ) 
    210210            ! 
    211211            pe3t(:,:,1:jpkm1,Kmm) = ze3t_f(:,:,1:jpkm1)        ! filtered scale factor at T-points 
     
    216216               CALL dom_vvl_interpol( pe3t(:,:,:,Kmm), pe3v(:,:,:,Kmm), 'V' ) 
    217217               DO_3D_11_11( 1, jpkm1 ) 
    218                   puu(ji,jj,jk,Kmm) = puu(ji,jj,jk,Kmm) + atfp * ( puu(ji,jj,jk,Kbb) - 2._wp * puu(ji,jj,jk,Kmm) + puu(ji,jj,jk,Kaa) ) 
    219                   pvv(ji,jj,jk,Kmm) = pvv(ji,jj,jk,Kmm) + atfp * ( pvv(ji,jj,jk,Kbb) - 2._wp * pvv(ji,jj,jk,Kmm) + pvv(ji,jj,jk,Kaa) ) 
     218                  puu(ji,jj,jk,Kmm) = puu(ji,jj,jk,Kmm) + rn_atfp * ( puu(ji,jj,jk,Kbb) - 2._wp * puu(ji,jj,jk,Kmm) + puu(ji,jj,jk,Kaa) ) 
     219                  pvv(ji,jj,jk,Kmm) = pvv(ji,jj,jk,Kmm) + rn_atfp * ( pvv(ji,jj,jk,Kbb) - 2._wp * pvv(ji,jj,jk,Kmm) + pvv(ji,jj,jk,Kaa) ) 
    220220               END_3D 
    221221               ! 
     
    234234                  zve3b = pe3v(ji,jj,jk,Kbb) * pvv(ji,jj,jk,Kbb) 
    235235                  ! 
    236                   puu(ji,jj,jk,Kmm) = ( zue3n + atfp * ( zue3b - 2._wp * zue3n  + zue3a ) ) / ze3u_f(ji,jj,jk) 
    237                   pvv(ji,jj,jk,Kmm) = ( zve3n + atfp * ( zve3b - 2._wp * zve3n  + zve3a ) ) / ze3v_f(ji,jj,jk) 
     236                  puu(ji,jj,jk,Kmm) = ( zue3n + rn_atfp * ( zue3b - 2._wp * zue3n  + zue3a ) ) / ze3u_f(ji,jj,jk) 
     237                  pvv(ji,jj,jk,Kmm) = ( zve3n + rn_atfp * ( zve3b - 2._wp * zve3n  + zve3a ) ) / ze3v_f(ji,jj,jk) 
    238238               END_3D 
    239239               pe3u(:,:,1:jpkm1,Kmm) = ze3u_f(:,:,1:jpkm1)   
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OCE/DYN/dynspg.F90

    r12406 r12443  
    6767      !!              ln_apr_dyn=T : the atmospheric pressure forcing is applied  
    6868      !!             as the gradient of the inverse barometer ssh: 
    69       !!                apgu = - 1/rau0 di[apr] = 0.5*grav di[ssh_ib+ssh_ibb] 
    70       !!                apgv = - 1/rau0 dj[apr] = 0.5*grav dj[ssh_ib+ssh_ibb] 
     69      !!                apgu = - 1/rho0 di[apr] = 0.5*grav di[ssh_ib+ssh_ibb] 
     70      !!                apgv = - 1/rho0 dj[apr] = 0.5*grav dj[ssh_ib+ssh_ibb] 
    7171      !!             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. 
     
    7878      ! 
    7979      INTEGER  ::   ji, jj, jk                   ! dummy loop indices 
    80       REAL(wp) ::   z2dt, zg_2, zintp, zgrau0r, zld   ! local scalars 
     80      REAL(wp) ::   z2dt, zg_2, zintp, zgrho0r, zld   ! local scalars 
    8181      REAL(wp), ALLOCATABLE, DIMENSION(:,:)   ::   zpice 
    8282      REAL(wp), ALLOCATABLE, DIMENSION(:,:,:) ::   ztrdu, ztrdv 
     
    123123            ! 
    124124            IF (ln_scal_load) THEN 
    125                zld = rn_scal_load * grav 
     125               zld = rn_load * grav 
    126126               DO_2D_00_00 
    127127                  spgu(ji,jj) = spgu(ji,jj) + zld * ( pssh(ji+1,jj,Kmm) - pssh(ji,jj,Kmm) ) * r1_e1u(ji,jj) 
     
    134134            ALLOCATE( zpice(jpi,jpj) ) 
    135135            zintp = REAL( MOD( kt-1, nn_fsbc ) ) / REAL( nn_fsbc ) 
    136             zgrau0r     = - grav * r1_rau0 
    137             zpice(:,:) = (  zintp * snwice_mass(:,:) + ( 1.- zintp ) * snwice_mass_b(:,:)  ) * zgrau0r 
     136            zgrho0r     = - grav * r1_rho0 
     137            zpice(:,:) = (  zintp * snwice_mass(:,:) + ( 1.- zintp ) * snwice_mass_b(:,:)  ) * zgrho0r 
    138138            DO_2D_00_00 
    139139               spgu(ji,jj) = spgu(ji,jj) + ( zpice(ji+1,jj) - zpice(ji,jj) ) * r1_e1u(ji,jj) 
     
    183183      NAMELIST/namdyn_spg/ ln_dynspg_exp       , ln_dynspg_ts,   & 
    184184      &                    ln_bt_fw, ln_bt_av  , ln_bt_auto  ,   & 
    185       &                    nn_baro , rn_bt_cmax, nn_bt_flt, rn_bt_alpha 
     185      &                    nn_e , rn_bt_cmax, nn_bt_flt, rn_bt_alpha 
    186186      !!---------------------------------------------------------------------- 
    187187      ! 
     
    222222      ! 
    223223      IF( nspg == np_TS ) THEN   ! split-explicit scheme initialisation 
    224          CALL dyn_spg_ts_init          ! do it first: set nn_baro used to allocate some arrays later on 
     224         CALL dyn_spg_ts_init          ! do it first: set nn_e used to allocate some arrays later on 
    225225      ENDIF 
    226226      ! 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OCE/DYN/dynspg_exp.F90

    r12377 r12443  
    4949      !!              momentum trend the surface pressure gradient : 
    5050      !!                      (uu(rhs),vv(rhs)) = (uu(rhs),vv(rhs)) + (spgu,spgv) 
    51       !!              where spgu = -1/rau0 d/dx(ps) = -g/e1u di( ssh(now) ) 
    52       !!                    spgv = -1/rau0 d/dy(ps) = -g/e2v dj( ssh(now) ) 
     51      !!              where spgu = -1/rho0 d/dx(ps) = -g/e1u di( ssh(now) ) 
     52      !!                    spgv = -1/rho0 d/dy(ps) = -g/e2v dj( ssh(now) ) 
    5353      !! 
    5454      !! ** Action :   (puu(:,:,:,Krhs),pvv(:,:,:,Krhs))   trend of horizontal velocity increased by  
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OCE/DYN/dynspg_ts.F90

    r12424 r12443  
    7272   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) ::   un_adv , vn_adv   !: Advection vel. at "now" barocl. step 
    7373   ! 
    74    INTEGER, SAVE :: icycle      ! Number of barotropic sub-steps for each internal step nn_baro <= 2.5 nn_baro 
    75    REAL(wp),SAVE :: rdtbt       ! Barotropic time step 
     74   INTEGER, SAVE :: icycle      ! Number of barotropic sub-steps for each internal step nn_e <= 2.5 nn_e 
     75   REAL(wp),SAVE :: rDt_e       ! Barotropic time step 
    7676   ! 
    7777   REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:)   ::   wgtbtp1, wgtbtp2   ! 1st & 2nd weights used in time filtering of barotropic fields 
     
    102102      ierr(:) = 0 
    103103      ! 
    104       ALLOCATE( wgtbtp1(3*nn_baro), wgtbtp2(3*nn_baro), zwz(jpi,jpj), STAT=ierr(1) ) 
     104      ALLOCATE( wgtbtp1(3*nn_e), wgtbtp2(3*nn_e), zwz(jpi,jpj), STAT=ierr(1) ) 
    105105      IF( ln_dynvor_een .OR. ln_dynvor_eeT )   & 
    106106         &     ALLOCATE( ftnw(jpi,jpj) , ftne(jpi,jpj) , ftsw(jpi,jpj) , ftse(jpi,jpj), STAT=ierr(2)   ) 
     
    185185      ll_fw_start = .FALSE. 
    186186      !                                         ! time offset in steps for bdy data update 
    187       IF( .NOT.ln_bt_fw ) THEN   ;   noffset = - nn_baro 
     187      IF( .NOT.ln_bt_fw ) THEN   ;   noffset = - nn_e 
    188188      ELSE                       ;   noffset =   0  
    189189      ENDIF 
     
    303303      IF( ln_bt_fw ) THEN                        ! Add wind forcing 
    304304         DO_2D_00_00 
    305             zu_frc(ji,jj) =  zu_frc(ji,jj) + r1_rau0 * utau(ji,jj) * r1_hu(ji,jj,Kmm) 
    306             zv_frc(ji,jj) =  zv_frc(ji,jj) + r1_rau0 * vtau(ji,jj) * r1_hv(ji,jj,Kmm) 
     305            zu_frc(ji,jj) =  zu_frc(ji,jj) + r1_rho0 * utau(ji,jj) * r1_hu(ji,jj,Kmm) 
     306            zv_frc(ji,jj) =  zv_frc(ji,jj) + r1_rho0 * vtau(ji,jj) * r1_hv(ji,jj,Kmm) 
    307307         END_2D 
    308308      ELSE 
    309          zztmp = r1_rau0 * r1_2 
     309         zztmp = r1_rho0 * r1_2 
    310310         DO_2D_00_00 
    311311            zu_frc(ji,jj) =  zu_frc(ji,jj) + zztmp * ( utau_b(ji,jj) + utau(ji,jj) ) * r1_hu(ji,jj,Kmm) 
     
    320320      !                                   ! ---------------------------------------------------  ! 
    321321      IF (ln_bt_fw) THEN                          ! FORWARD integration: use kt+1/2 fluxes (NOW+1/2) 
    322          zssh_frc(:,:) = r1_rau0 * ( emp(:,:) - rnf(:,:) + fwfisf_cav(:,:) + fwfisf_par(:,:) ) 
     322         zssh_frc(:,:) = r1_rho0 * ( emp(:,:) - rnf(:,:) + fwfisf_cav(:,:) + fwfisf_par(:,:) ) 
    323323      ELSE                                        ! CENTRED integration: use kt-1/2 + kt+1/2 fluxes (NOW) 
    324          zztmp = r1_rau0 * r1_2 
     324         zztmp = r1_rho0 * r1_2 
    325325         zssh_frc(:,:) = zztmp * (  emp(:,:)        + emp_b(:,:)                    & 
    326326                &                 - rnf(:,:)        - rnf_b(:,:)                    & 
     
    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*rn_Dt + (jn + noffset - 1) * rn_Dt / REAL(nn_baro, wp) 
     431            zt0substep = REAL(nsec_day, wp) - 0.5_wp*rn_Dt + (jn + noffset - 1) * rn_Dt / REAL(nn_e, wp) 
    432432            CALL upd_tide(zt0substep, Kmm) 
    433433         END IF 
     
    495495         IF( .NOT.Agrif_Root() .AND. ln_bt_fw ) CALL agrif_dyn_ts_flux( jn, zhU, zhV ) 
    496496#endif 
    497          IF( ln_wd_il )   CALL wad_lmt_bt(zhU, zhV, sshn_e, zssh_frc, rdtbt)    !!gm wad_lmt_bt use of lbc_lnk on zhU, zhV 
     497         IF( ln_wd_il )   CALL wad_lmt_bt(zhU, zhV, sshn_e, zssh_frc, rDt_e)    !!gm wad_lmt_bt use of lbc_lnk on zhU, zhV 
    498498 
    499499         IF( ln_wd_dl ) THEN           ! un_e and vn_e are set to zero at faces where  
     
    510510         DO_2D_00_00 
    511511            zhdiv = (   zhU(ji,jj) - zhU(ji-1,jj) + zhV(ji,jj) - zhV(ji,jj-1)   ) * r1_e1e2t(ji,jj) 
    512             ssha_e(ji,jj) = (  sshn_e(ji,jj) - rdtbt * ( zssh_frc(ji,jj) + zhdiv )  ) * ssmask(ji,jj) 
     512            ssha_e(ji,jj) = (  sshn_e(ji,jj) - rDt_e * ( zssh_frc(ji,jj) + zhdiv )  ) * ssmask(ji,jj) 
    513513         END_2D 
    514514         ! 
     
    552552         ! 
    553553         !                             ! Surface pressure gradient 
    554          zldg = ( 1._wp - rn_scal_load ) * grav    ! local factor 
     554         zldg = ( 1._wp - rn_load ) * grav    ! local factor 
    555555         DO_2D_00_00 
    556556            zu_spg(ji,jj) = - zldg * ( zsshp2_e(ji+1,jj) - zsshp2_e(ji,jj) ) * r1_e1u(ji,jj) 
     
    600600            DO_2D_00_00 
    601601               ua_e(ji,jj) = (                                 un_e(ji,jj)   &  
    602                          &     + rdtbt * (                   zu_spg(ji,jj)   & 
     602                         &     + rDt_e * (                   zu_spg(ji,jj)   & 
    603603                         &                                 + zu_trd(ji,jj)   & 
    604604                         &                                 + zu_frc(ji,jj) ) &  
     
    606606 
    607607               va_e(ji,jj) = (                                 vn_e(ji,jj)   & 
    608                          &     + rdtbt * (                   zv_spg(ji,jj)   & 
     608                         &     + rDt_e * (                   zv_spg(ji,jj)   & 
    609609                         &                                 + zv_trd(ji,jj)   & 
    610610                         &                                 + zv_frc(ji,jj) ) & 
     
    625625               ! 
    626626               ua_e(ji,jj) = (               hu_e  (ji,jj) *   un_e (ji,jj)      &  
    627                     &            + rdtbt * (  zhu_bck        * zu_spg (ji,jj)  &   ! 
     627                    &            + rDt_e * (  zhu_bck        * zu_spg (ji,jj)  &   ! 
    628628                    &                       + zhup2_e(ji,jj) * zu_trd (ji,jj)  &   ! 
    629629                    &                       +  hu(ji,jj,Kmm) * zu_frc (ji,jj)  )   ) * z1_hu 
    630630               ! 
    631631               va_e(ji,jj) = (               hv_e  (ji,jj) *   vn_e (ji,jj)      & 
    632                     &            + rdtbt * (  zhv_bck        * zv_spg (ji,jj)  &   ! 
     632                    &            + rDt_e * (  zhv_bck        * zv_spg (ji,jj)  &   ! 
    633633                    &                       + zhvp2_e(ji,jj) * zv_trd (ji,jj)  &   ! 
    634634                    &                       +  hv(ji,jj,Kmm) * zv_frc (ji,jj)  )   ) * z1_hv 
     
    638638         IF ( ll_wd ) THEN ! revert to explicit for bit comparison tests in non wad runs 
    639639            DO_2D_00_00 
    640                   ua_e(ji,jj) =  ua_e(ji,jj) /(1.0 -   rdtbt * zCdU_u(ji,jj) * hur_e(ji,jj)) 
    641                   va_e(ji,jj) =  va_e(ji,jj) /(1.0 -   rdtbt * zCdU_v(ji,jj) * hvr_e(ji,jj)) 
     640                  ua_e(ji,jj) =  ua_e(ji,jj) /(1.0 -   rDt_e * zCdU_u(ji,jj) * hur_e(ji,jj)) 
     641                  va_e(ji,jj) =  va_e(ji,jj) /(1.0 -   rDt_e * zCdU_v(ji,jj) * hvr_e(ji,jj)) 
    642642            END_2D 
    643643         ENDIF 
     
    707707               zvn_save = vn_adv(ji,jj) 
    708708               !                          ! apply the previously computed correction  
    709                un_adv(ji,jj) = r1_2 * ( ub2_b(ji,jj) + zun_save - atfp * un_bf(ji,jj) ) 
    710                vn_adv(ji,jj) = r1_2 * ( vb2_b(ji,jj) + zvn_save - atfp * vn_bf(ji,jj) ) 
     709               un_adv(ji,jj) = r1_2 * ( ub2_b(ji,jj) + zun_save - rn_atfp * un_bf(ji,jj) ) 
     710               vn_adv(ji,jj) = r1_2 * ( vb2_b(ji,jj) + zvn_save - rn_atfp * vn_bf(ji,jj) ) 
    711711               !                          ! Update corrective fluxes for next time step 
    712                un_bf(ji,jj)  = atfp * un_bf(ji,jj) + ( zun_save - ub2_b(ji,jj) ) 
    713                vn_bf(ji,jj)  = atfp * vn_bf(ji,jj) + ( zvn_save - vb2_b(ji,jj) ) 
     712               un_bf(ji,jj)  = rn_atfp * un_bf(ji,jj) + ( zun_save - ub2_b(ji,jj) ) 
     713               vn_bf(ji,jj)  = rn_atfp * vn_bf(ji,jj) + ( zvn_save - vb2_b(ji,jj) ) 
    714714               !                          ! Save integrated transport for next computation 
    715715               ub2_b(ji,jj) = zun_save 
     
    809809      LOGICAL, INTENT(in) ::   ll_fw      ! forward time splitting =.true. 
    810810      INTEGER, INTENT(inout) :: jpit      ! cycle length     
    811       REAL(wp), DIMENSION(3*nn_baro), INTENT(inout) ::   zwgt1, & ! Primary weights 
     811      REAL(wp), DIMENSION(3*nn_e), INTENT(inout) ::   zwgt1, & ! Primary weights 
    812812                                                         zwgt2    ! Secondary weights 
    813813       
     
    821821      ! Set time index when averaged value is requested 
    822822      IF (ll_fw) THEN  
    823          jic = nn_baro 
     823         jic = nn_e 
    824824      ELSE 
    825          jic = 2 * nn_baro 
     825         jic = 2 * nn_e 
    826826      ENDIF 
    827827 
     
    829829      IF (ll_av) THEN 
    830830           ! Define simple boxcar window for primary weights  
    831            ! (width = nn_baro, centered around jic)      
     831           ! (width = nn_e, centered around jic)      
    832832         SELECT CASE ( nn_bt_flt ) 
    833833              CASE( 0 )  ! No averaging 
     
    835835                 jpit = jic 
    836836 
    837               CASE( 1 )  ! Boxcar, width = nn_baro 
    838                  DO jn = 1, 3*nn_baro 
    839                     za1 = ABS(float(jn-jic))/float(nn_baro)  
     837              CASE( 1 )  ! Boxcar, width = nn_e 
     838                 DO jn = 1, 3*nn_e 
     839                    za1 = ABS(float(jn-jic))/float(nn_e)  
    840840                    IF (za1 < 0.5_wp) THEN 
    841841                      zwgt1(jn) = 1._wp 
     
    844844                 ENDDO 
    845845 
    846               CASE( 2 )  ! Boxcar, width = 2 * nn_baro 
    847                  DO jn = 1, 3*nn_baro 
    848                     za1 = ABS(float(jn-jic))/float(nn_baro)  
     846              CASE( 2 )  ! Boxcar, width = 2 * nn_e 
     847                 DO jn = 1, 3*nn_e 
     848                    za1 = ABS(float(jn-jic))/float(nn_e)  
    849849                    IF (za1 < 1._wp) THEN 
    850850                      zwgt1(jn) = 1._wp 
     
    976976 
    977977      ! Estimate number of iterations to satisfy a max courant number= rn_bt_cmax 
    978       IF( ln_bt_auto )   nn_baro = CEILING( rn_Dt / rn_bt_cmax * zcmax) 
     978      IF( ln_bt_auto )   nn_e = CEILING( rn_Dt / rn_bt_cmax * zcmax) 
    979979       
    980       rdtbt = rn_Dt / REAL( nn_baro , wp ) 
    981       zcmax = zcmax * rdtbt 
     980      rDt_e = rn_Dt / REAL( nn_e , wp ) 
     981      zcmax = zcmax * rDt_e 
    982982      ! Print results 
    983983      IF(lwp) WRITE(numout,*) 
     
    985985      IF(lwp) WRITE(numout,*) '~~~~~~~~~~~~~~~' 
    986986      IF( ln_bt_auto ) THEN 
    987          IF(lwp) WRITE(numout,*) '     ln_ts_auto =.true. Automatically set nn_baro ' 
     987         IF(lwp) WRITE(numout,*) '     ln_ts_auto =.true. Automatically set nn_e ' 
    988988         IF(lwp) WRITE(numout,*) '     Max. courant number allowed: ', rn_bt_cmax 
    989989      ELSE 
    990          IF(lwp) WRITE(numout,*) '     ln_ts_auto=.false.: Use nn_baro in namelist   nn_baro = ', nn_baro 
     990         IF(lwp) WRITE(numout,*) '     ln_ts_auto=.false.: Use nn_e in namelist   nn_e = ', nn_e 
    991991      ENDIF 
    992992 
    993993      IF(ln_bt_av) THEN 
    994          IF(lwp) WRITE(numout,*) '     ln_bt_av =.true.  ==> Time averaging over nn_baro time steps is on ' 
     994         IF(lwp) WRITE(numout,*) '     ln_bt_av =.true.  ==> Time averaging over nn_e time steps is on ' 
    995995      ELSE 
    996996         IF(lwp) WRITE(numout,*) '     ln_bt_av =.false. => No time averaging of barotropic variables ' 
     
    10121012      SELECT CASE ( nn_bt_flt ) 
    10131013         CASE( 0 )      ;   IF(lwp) WRITE(numout,*) '           Dirac' 
    1014          CASE( 1 )      ;   IF(lwp) WRITE(numout,*) '           Boxcar: width = nn_baro' 
    1015          CASE( 2 )      ;   IF(lwp) WRITE(numout,*) '           Boxcar: width = 2*nn_baro'  
     1014         CASE( 1 )      ;   IF(lwp) WRITE(numout,*) '           Boxcar: width = nn_e' 
     1015         CASE( 2 )      ;   IF(lwp) WRITE(numout,*) '           Boxcar: width = 2*nn_e'  
    10161016         CASE DEFAULT   ;   CALL ctl_stop( 'unrecognised value for nn_bt_flt: should 0,1, or 2' ) 
    10171017      END SELECT 
    10181018      ! 
    10191019      IF(lwp) WRITE(numout,*) ' ' 
    1020       IF(lwp) WRITE(numout,*) '     nn_baro = ', nn_baro 
    1021       IF(lwp) WRITE(numout,*) '     Barotropic time step [s] is :', rdtbt 
     1020      IF(lwp) WRITE(numout,*) '     nn_e = ', nn_e 
     1021      IF(lwp) WRITE(numout,*) '     Barotropic time step [s] is :', rDt_e 
    10221022      IF(lwp) WRITE(numout,*) '     Maximum Courant number is   :', zcmax 
    10231023      ! 
     
    10311031      ENDIF 
    10321032      IF( zcmax>0.9_wp ) THEN 
    1033          CALL ctl_stop( 'dynspg_ts ERROR: Maximum Courant number is greater than 0.9: Inc. nn_baro !' )           
     1033         CALL ctl_stop( 'dynspg_ts ERROR: Maximum Courant number is greater than 0.9: Inc. nn_e !' )           
    10341034      ENDIF 
    10351035      ! 
     
    14301430      ! 
    14311431      IF( ln_wd_il ) THEN      ! W/D : use the "clipped" bottom friction   !!gm   explain WHY, please ! 
    1432          zztmp = -1._wp / rdtbt 
     1432         zztmp = -1._wp / rDt_e 
    14331433         DO_2D_00_00 
    14341434            pu_RHSi(ji,jj) = pu_RHSi(ji,jj) + zu_i(ji,jj) *  wdrampu(ji,jj) * MAX(                                 &  
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OCE/DYN/dynzdf.F90

    r12424 r12443  
    261261         ze3ua =  ( 1._wp - r_vvl ) * e3u(ji,jj,1,Kmm) + r_vvl * e3u(ji,jj,1,Kaa)  
    262262         puu(ji,jj,1,Kaa) = puu(ji,jj,1,Kaa) + rDt * 0.5_wp * ( utau_b(ji,jj) + utau(ji,jj) )   & 
    263             &                                      / ( ze3ua * rau0 ) * umask(ji,jj,1)  
     263            &                                      / ( ze3ua * rho0 ) * umask(ji,jj,1)  
    264264      END_2D 
    265265      DO_3D_00_00( 2, jpkm1 ) 
     
    385385         ze3va =  ( 1._wp - r_vvl ) * e3v(ji,jj,1,Kmm) + r_vvl * e3v(ji,jj,1,Kaa)  
    386386         pvv(ji,jj,1,Kaa) = pvv(ji,jj,1,Kaa) + rDt * 0.5_wp * ( vtau_b(ji,jj) + vtau(ji,jj) )   & 
    387             &                                      / ( ze3va * rau0 ) * vmask(ji,jj,1)  
     387            &                                      / ( ze3va * rho0 ) * vmask(ji,jj,1)  
    388388      END_2D 
    389389      DO_3D_00_00( 2, jpkm1 ) 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OCE/DYN/sshwzv.F90

    r12424 r12443  
    8888      ENDIF 
    8989      ! 
    90       zcoef = 0.5_wp * r1_rau0 
     90      zcoef = 0.5_wp * r1_rho0 
    9191 
    9292      !                                           !------------------------------! 
     
    222222      !! ** Method  : - apply Asselin time fiter to now ssh (excluding the forcing 
    223223      !!              from the filter, see Leclair and Madec 2010) and swap : 
    224       !!                pssh(:,:,Kmm) = pssh(:,:,Kaa) + atfp * ( pssh(:,:,Kbb) -2 pssh(:,:,Kmm) + pssh(:,:,Kaa) ) 
    225       !!                            - atfp * rn_Dt * ( emp_b - emp ) / rau0 
     224      !!                pssh(:,:,Kmm) = pssh(:,:,Kaa) + rn_atfp * ( pssh(:,:,Kbb) -2 pssh(:,:,Kmm) + pssh(:,:,Kaa) ) 
     225      !!                            - rn_atfp * rn_Dt * ( emp_b - emp ) / rho0 
    226226      !! 
    227227      !! ** action  : - pssh(:,:,Kmm) time filtered 
     
    246246      IF ( .NOT.( l_1st_euler ) ) THEN   ! Only do time filtering for leapfrog timesteps 
    247247         !                                                  ! filtered "now" field 
    248          pssh(:,:,Kmm) = pssh(:,:,Kmm) + atfp * ( pssh(:,:,Kbb) - 2 * pssh(:,:,Kmm) + pssh(:,:,Kaa) ) 
     248         pssh(:,:,Kmm) = pssh(:,:,Kmm) + rn_atfp * ( pssh(:,:,Kbb) - 2 * pssh(:,:,Kmm) + pssh(:,:,Kaa) ) 
    249249         IF( .NOT.ln_linssh ) THEN                          ! "now" <-- with forcing removed 
    250             zcoef = atfp * rn_Dt * r1_rau0 
     250            zcoef = rn_atfp * rn_Dt * r1_rho0 
    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 * rn_Dt * ( risfcpl_ssh(:,:) - 0.0 ) * ssmask(:,:) 
     257            IF ( ln_isf .AND. ln_isfcpl .AND. kt == nit000+1) pssh(:,:,Kbb) = pssh(:,:,Kbb) - rn_atfp * rn_Dt * ( risfcpl_ssh(:,:) - 0.0 ) * ssmask(:,:) 
    258258 
    259259         ENDIF 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OCE/DYN/wet_dry.F90

    r12377 r12443  
    270270 
    271271 
    272    SUBROUTINE wad_lmt_bt( zflxu, zflxv, sshn_e, zssh_frc, rdtbt ) 
     272   SUBROUTINE wad_lmt_bt( zflxu, zflxv, sshn_e, zssh_frc, rDt_e ) 
    273273      !!---------------------------------------------------------------------- 
    274274      !!                  ***  ROUTINE wad_lmt  *** 
     
    280280      !! ** Action  : - calculate flux limiter and W/D flag 
    281281      !!---------------------------------------------------------------------- 
    282       REAL(wp)                , INTENT(in   ) ::   rdtbt    ! ocean time-step index 
     282      REAL(wp)                , INTENT(in   ) ::   rDt_e    ! ocean time-step index 
    283283      REAL(wp), DIMENSION(:,:), INTENT(inout) ::   zflxu,  zflxv, sshn_e, zssh_frc   
    284284      ! 
     
    299299      zdepwd = 50._wp   ! maximum depth that ocean cells can have W/D processes 
    300300      ! 
    301       z2dt = rdtbt    
     301      z2dt = rDt_e    
    302302      ! 
    303303      zflxp(:,:)   = 0._wp 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OCE/FLO/flowri.F90

    r12406 r12443  
    122122               ztem(jfl) = ts(iafloc,ibfloc,icfl,jp_tem,Kmm) 
    123123               zsal (jfl) = ts(iafloc,ibfloc,icfl,jp_sal,Kmm) 
    124                zrho (jfl) = (rhd(iafloc,ibfloc,icfl)+1)*rau0 
     124               zrho (jfl) = (rhd(iafloc,ibfloc,icfl)+1)*rho0 
    125125 
    126126            ENDIF 
     
    142142            ztem(jfl) = ts(iafloc,ibfloc,icfl,jp_tem,Kmm) 
    143143            zsal(jfl) = ts(iafloc,ibfloc,icfl,jp_sal,Kmm) 
    144             zrho(jfl) = (rhd(iafloc,ibfloc,icfl)+1)*rau0 
     144            zrho(jfl) = (rhd(iafloc,ibfloc,icfl)+1)*rho0 
    145145           
    146146         ENDIF 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OCE/IOM/restart.F90

    r12406 r12443  
    262262      IF ( ln_diurnal_only ) THEN  
    263263         IF(lwp) WRITE( numout, * ) & 
    264          &   "rst_read:- ln_diurnal_only set, setting rhop=rau0"  
    265          rhop = rau0 
     264         &   "rst_read:- ln_diurnal_only set, setting rhop=rho0"  
     265         rhop = rho0 
    266266         CALL iom_get( numror, jpdom_autoglo, 'tn'     , w3d, ldxios = lrxios )  
    267267         ts(:,:,1,jp_tem,Kmm) = w3d(:,:,1) 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OCE/ISF/isfcav.F90

    r12343 r12443  
    2424   USE oce      , ONLY: ts                              ! ocean tracers 
    2525   USE par_oce  , ONLY: jpi,jpj                         ! ocean space and time domain 
    26    USE phycst   , ONLY: grav,rau0,rau0_rcp,r1_rau0_rcp  ! physical constants 
     26   USE phycst   , ONLY: grav,rho0,rho0_rcp,r1_rho0_rcp  ! physical constants 
    2727   USE eosbn2   , ONLY: ln_teos10                       ! use ln_teos10 or not 
    2828   ! 
     
    8585      ! 
    8686      ! initialisation 
    87       IF (TRIM(cn_gammablk) == 'vel_stab' ) zqoce_b (:,:) = ptsc(:,:,jp_tem) * rau0_rcp ! last time step total heat fluxes (to speed up convergence) 
     87      IF (TRIM(cn_gammablk) == 'vel_stab' ) zqoce_b (:,:) = ptsc(:,:,jp_tem) * rho0_rcp ! last time step total heat fluxes (to speed up convergence) 
    8888      ! 
    8989      ! compute ice shelf melting 
     
    142142      ! 
    143143      ! set temperature content 
    144       ptsc(:,:,jp_tem) = - zqh(:,:) * r1_rau0_rcp 
     144      ptsc(:,:,jp_tem) = - zqh(:,:) * r1_rho0_rcp 
    145145      ! 
    146146      ! write restart variables (qoceisf, qhcisf, fwfisf for now and before) 
     
    215215            risf_lamb1 =-0.0564_wp 
    216216            risf_lamb2 = 0.0773_wp 
    217             risf_lamb3 =-7.8633e-8 * grav * rau0 
     217            risf_lamb3 =-7.8633e-8 * grav * rho0 
    218218         ELSE                  ! linearisation from table 4 (Asay-Davis et al., 2015) 
    219219            risf_lamb1 =-0.0573_wp 
    220220            risf_lamb2 = 0.0832_wp 
    221             risf_lamb3 =-7.5300e-8 * grav * rau0 
     221            risf_lamb3 =-7.5300e-8 * grav * rho0 
    222222         ENDIF 
    223223 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OCE/ISF/isfcavmlt.F90

    r12340 r12443  
    1717 
    1818   USE dom_oce                            ! ocean space and time domain 
    19    USE phycst , ONLY: rcp, rau0, rau0_rcp ! physical constants 
     19   USE phycst , ONLY: rcp, rho0, rho0_rcp ! physical constants 
    2020   USE eosbn2 , ONLY: eos_fzp             ! equation of state 
    2121 
     
    161161      ! 
    162162      ! compute ocean-ice heat flux and then derive fwf assuming that ocean heat flux equal latent heat 
    163       pqfwf(:,:) = - pgt(:,:) * rau0_rcp * zthd(:,:) / rLfusisf    ! fresh water flux  ( > 0 out ) 
     163      pqfwf(:,:) = - pgt(:,:) * rho0_rcp * zthd(:,:) / rLfusisf    ! fresh water flux  ( > 0 out ) 
    164164      pqoce(:,:) = - pqfwf(:,:) * rLfusisf                         ! ocea-ice flux     ( > 0 out ) 
    165165      pqhc (:,:) =   pqfwf(:,:) * ztfrz(:,:) * rcp                 ! heat content flux ( > 0 out ) 
     
    213213         ! 
    214214         ! compute coeficient to solve the 2nd order equation 
    215          zeps1 = rau0_rcp * pgt(ji,jj) 
    216          zeps2 = rLfusisf * rau0 * pgs(ji,jj) 
     215         zeps1 = rho0_rcp * pgt(ji,jj) 
     216         zeps2 = rLfusisf * rho0 * pgs(ji,jj) 
    217217         zeps3 = rhoisf * rcpisf * rkappa / MAX(risfdep(ji,jj),zeps) 
    218218         zeps4 = risf_lamb2 + risf_lamb3 * risfdep(ji,jj) 
     
    238238         ! 
    239239         ! compute the upward water and heat flux (eq. 24 and eq. 26) 
    240          pqfwf(ji,jj) = rau0     * pgs(ji,jj) * ( zsfrz - pstbl(ji,jj) ) / MAX(zsfrz,zeps) ! fresh water flux    (> 0 out) 
    241          pqoce(ji,jj) = rau0_rcp * pgt(ji,jj) * zthd (ji,jj)                               ! ocean-ice heat flux (> 0 out) 
     240         pqfwf(ji,jj) = rho0     * pgs(ji,jj) * ( zsfrz - pstbl(ji,jj) ) / MAX(zsfrz,zeps) ! fresh water flux    (> 0 out) 
     241         pqoce(ji,jj) = rho0_rcp * pgt(ji,jj) * zthd (ji,jj)                               ! ocean-ice heat flux (> 0 out) 
    242242         pqhc (ji,jj) = rcp      * pqfwf(ji,jj) * ztfrz(ji,jj)                             ! heat content   flux (> 0 out) 
    243243         ! 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OCE/ISF/isfdynatf.F90

    r12406 r12443  
    1313   USE isf_oce 
    1414 
    15    USE phycst , ONLY: r1_rau0         ! physical constant 
     15   USE phycst , ONLY: r1_rho0         ! physical constant 
    1616   USE dom_oce, ONLY: tmask, ssmask, ht, e3t, r1_e1e2t   ! time and space domain 
    1717 
     
    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 * rn_Dt * r1_rau0 
     41      REAL(wp)                        , INTENT(in   ) :: pcoef    ! rn_atfp * rn_Dt * r1_rho0 
    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 * rn_Dt * r1_rau0 
     72      REAL(wp),                         INTENT(in   ) :: pcoef           ! rn_atfp * rn_Dt * r1_rho0 
    7373      !!---------------------------------------------------------------------- 
    7474      INTEGER :: ji,jj,jk 
     
    7777      ! 
    7878      ! compute fwf conservation correction 
    79       zfwfinc(:,:) = pcoef * ( pfwf_b(:,:) - pfwf(:,:) ) / ( ht(:,:) + 1._wp - ssmask(:,:) ) * r1_rau0 
     79      zfwfinc(:,:) = pcoef * ( pfwf_b(:,:) - pfwf(:,:) ) / ( ht(:,:) + 1._wp - ssmask(:,:) ) * r1_rho0 
    8080      ! 
    8181      ! add the increment 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OCE/ISF/isfhdiv.F90

    r12340 r12443  
    1616 
    1717   USE dom_oce                ! time and space domain 
    18    USE phycst , ONLY: r1_rau0 ! physical constant 
     18   USE phycst , ONLY: r1_rho0 ! physical constant 
    1919   USE in_out_manager         ! 
    2020 
     
    9696      ! 
    9797      ! compute integrated divergence correction 
    98       zhdiv(:,:) = 0.5_wp * ( pfwf(:,:) + pfwf_b(:,:) ) * r1_rau0 / phtbl(:,:) 
     98      zhdiv(:,:) = 0.5_wp * ( pfwf(:,:) + pfwf_b(:,:) ) * r1_rho0 / phtbl(:,:) 
    9999      ! 
    100100      ! update divergence at each level affected by ice shelf top boundary layer 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OCE/ISF/isfpar.F90

    r12077 r12443  
    2424   USE dom_oce  , ONLY: bathy          ! ocean space and time domain 
    2525   USE par_oce  , ONLY: jpi,jpj        ! ocean space and time domain 
    26    USE phycst   , ONLY: r1_rau0_rcp    ! physical constants 
     26   USE phycst   , ONLY: r1_rho0_rcp    ! physical constants 
    2727   ! 
    2828   USE in_out_manager ! I/O manager 
     
    8888      ! 
    8989      ! set temperature content 
    90       ptsc(:,:,jp_tem) = zqh(:,:) * r1_rau0_rcp 
     90      ptsc(:,:,jp_tem) = zqh(:,:) * r1_rho0_rcp 
    9191      ! 
    9292      ! write restart variables (qoceisf, qhcisf, fwfisf for now and before) 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OCE/ISF/isfparmlt.F90

    r12077 r12443  
    1313   USE dom_oce                  ! ocean space and time domain 
    1414   USE oce    , ONLY: ts        ! ocean dynamics and tracers 
    15    USE phycst , ONLY: rcp, rau0 ! physical constants 
     15   USE phycst , ONLY: rcp, rho0 ! physical constants 
    1616   USE eosbn2 , ONLY: eos_fzp   ! equation of state 
    1717 
     
    148148      ! 
    149149      ! 2. ------------Net heat flux and fresh water flux due to the ice shelf 
    150       pqoce(:,:) =   rau0 * rcp * rn_gammat0 * risfLeff(:,:) * e1t(:,:) * ( ztavg(:,:) - ztfrz(:,:) ) * r1_e1e2t(:,:) 
     150      pqoce(:,:) =   rho0 * rcp * rn_gammat0 * risfLeff(:,:) * e1t(:,:) * ( ztavg(:,:) - ztfrz(:,:) ) * r1_e1e2t(:,:) 
    151151      pqfwf(:,:) = - pqoce(:,:) / rLfusisf             ! derived from the latent heat flux 
    152152      pqhc (:,:) =   pqfwf(:,:) * ztfrz(:,:) * rcp     ! heat content flux  
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OCE/LDF/ldftra.F90

    r12377 r12443  
    820820      ! 
    821821      IF( iom_use('weiv_masstr') ) THEN   ! vertical mass transport & its square value 
    822          zw2d(:,:) = rau0 * e1e2t(:,:) 
     822         zw2d(:,:) = rho0 * e1e2t(:,:) 
    823823         DO jk = 1, jpk 
    824824            zw3d(:,:,jk) = zw3d(:,:,jk) * zw2d(:,:) 
     
    830830         zw3d(:,:,:) = 0.e0 
    831831         DO jk = 1, jpkm1 
    832             zw3d(:,:,jk) = rau0 * ( psi_uw(:,:,jk+1) - psi_uw(:,:,jk) )  
     832            zw3d(:,:,jk) = rho0 * ( psi_uw(:,:,jk+1) - psi_uw(:,:,jk) )  
    833833         END DO 
    834834         CALL iom_put( "ueiv_masstr", zw3d )                  ! mass transport in i-direction 
    835835      ENDIF 
    836836      ! 
    837       zztmp = 0.5_wp * rau0 * rcp  
     837      zztmp = 0.5_wp * rho0 * rcp  
    838838      IF( iom_use('ueiv_heattr') .OR. iom_use('ueiv_heattr3d') ) THEN 
    839839        zw2d(:,:)   = 0._wp  
     
    853853         zw3d(:,:,:) = 0.e0 
    854854         DO jk = 1, jpkm1 
    855             zw3d(:,:,jk) = rau0 * ( psi_vw(:,:,jk+1) - psi_vw(:,:,jk) )  
     855            zw3d(:,:,jk) = rho0 * ( psi_vw(:,:,jk+1) - psi_vw(:,:,jk) )  
    856856         END DO 
    857857         CALL iom_put( "veiv_masstr", zw3d )                  ! mass transport in i-direction 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OCE/SBC/fldread.F90

    r12406 r12443  
    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 ) * rn_Dt / REAL(nn_baro,wp) ) 
     174         isecsbc = nsec_year + nsec1jan000 + NINT( (     REAL(      kit,wp) + zt_offset ) * rn_Dt / REAL(nn_e,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 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OCE/SBC/sbcapr.F90

    r12377 r12443  
    3636    
    3737   REAL(wp) ::   tarea                ! whole domain mean masked ocean surface 
    38    REAL(wp) ::   r1_grau              ! = 1.e0 / (grav * rau0) 
     38   REAL(wp) ::   r1_grau              ! = 1.e0 / (grav * rho0) 
    3939    
    4040   TYPE(FLD), ALLOCATABLE, DIMENSION(:) ::   sf_apr   ! structure of input fields (file informations, fields read) 
     
    9898      ENDIF 
    9999      ! 
    100       r1_grau = 1.e0 / (grav * rau0)               !* constant for optimization 
     100      r1_grau = 1.e0 / (grav * rho0)               !* constant for optimization 
    101101      ! 
    102102      !                                            !* control check 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OCE/SBC/sbccpl.F90

    r12406 r12443  
    193193 
    194194   REAL(wp) ::   rpref = 101000._wp   ! reference atmospheric pressure[N/m2]  
    195    REAL(wp) ::   r1_grau              ! = 1.e0 / (grav * rau0)  
     195   REAL(wp) ::   r1_grau              ! = 1.e0 / (grav * rho0)  
    196196 
    197197   INTEGER , ALLOCATABLE, SAVE, DIMENSION(:) ::   nrcvinfo           ! OASIS info argument 
     
    12501250          IF( kt /= nit000 )   ssh_ibb(:,:) = ssh_ib(:,:)    !* Swap of ssh_ib fields  
    12511251 
    1252           r1_grau = 1.e0 / (grav * rau0)               !* constant for optimization  
     1252          r1_grau = 1.e0 / (grav * rho0)               !* constant for optimization  
    12531253          ssh_ib(:,:) = - ( frcv(jpr_mslp)%z3(:,:,1) - rpref ) * r1_grau    ! equivalent ssh (inverse barometer)  
    12541254          apr   (:,:) =     frcv(jpr_mslp)%z3(:,:,1)                         !atmospheric pressure  
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OCE/SBC/sbcfwb.F90

    r12406 r12443  
    130130            a_fwb_b = a_fwb                           ! mean sea level taking into account the ice+snow 
    131131                                                      ! sum over the global domain 
    132             a_fwb   = glob_sum( 'sbcfwb', e1e2t(:,:) * ( ssh(:,:,Kmm) + snwice_mass(:,:) * r1_rau0 ) ) 
     132            a_fwb   = glob_sum( 'sbcfwb', e1e2t(:,:) * ( ssh(:,:,Kmm) + snwice_mass(:,:) * r1_rho0 ) ) 
    133133            a_fwb   = a_fwb * 1.e+3 / ( area * rday * 365. )     ! convert in Kg/m3/s = mm/s 
    134134!!gm        !                                                      !!bug 365d year  
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OCE/SBC/sbcice_cice.F90

    r12377 r12443  
    1313   USE dom_oce         ! ocean space and time domain 
    1414   USE domvvl 
    15    USE phycst, only : rcp, rau0, r1_rau0, rhos, rhoi 
     15   USE phycst, only : rcp, rho0, r1_rho0, rhos, rhoi 
    1616   USE in_out_manager  ! I/O manager 
    1717   USE iom, ONLY : iom_put,iom_use              ! I/O manager library !!Joakim edit 
     
    228228      IF( .NOT.ln_rstart ) THEN 
    229229         IF( ln_ice_embd ) THEN            ! embedded sea-ice: deplete the initial ssh below sea-ice area 
    230             ssh(:,:,Kmm) = ssh(:,:,Kmm) - snwice_mass(:,:) * r1_rau0 
    231             ssh(:,:,Kbb) = ssh(:,:,Kbb) - snwice_mass(:,:) * r1_rau0 
     230            ssh(:,:,Kmm) = ssh(:,:,Kmm) - snwice_mass(:,:) * r1_rho0 
     231            ssh(:,:,Kbb) = ssh(:,:,Kbb) - snwice_mass(:,:) * r1_rho0 
    232232 
    233233!!gm This should be put elsewhere....   (same remark for limsbc) 
     
    417417! Freezing/melting potential 
    418418! Calculated over NEMO leapfrog timestep (hence 2*dt) 
    419       nfrzmlt(:,:) = rau0 * rcp * e3t_m(:,:) * ( Tocnfrz-sst_m(:,:) ) / ( 2.0*dt ) 
     419      nfrzmlt(:,:) = rho0 * rcp * e3t_m(:,:) * ( Tocnfrz-sst_m(:,:) ) / ( 2.0*dt ) 
    420420 
    421421      ztmp(:,:) = nfrzmlt(:,:) 
     
    450450         zintb = REAL( nn_fsbc + 1 ) / REAL( nn_fsbc ) * 0.5_wp 
    451451          ! 
    452          zpice(:,:) = ssh_m(:,:) + (  zintn * snwice_mass(:,:) +  zintb * snwice_mass_b(:,:)  ) * r1_rau0 
     452         zpice(:,:) = ssh_m(:,:) + (  zintn * snwice_mass(:,:) +  zintb * snwice_mass_b(:,:)  ) * r1_rho0 
    453453          ! 
    454454         ! 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OCE/SBC/sbcrnf.F90

    r12377 r12443  
    137137         !                                                           ! set temperature & salinity content of runoffs 
    138138         IF( ln_rnf_tem ) THEN                                       ! use runoffs temperature data 
    139             rnf_tsc(:,:,jp_tem) = ( sf_t_rnf(1)%fnow(:,:,1) ) * rnf(:,:) * r1_rau0 
     139            rnf_tsc(:,:,jp_tem) = ( sf_t_rnf(1)%fnow(:,:,1) ) * rnf(:,:) * r1_rho0 
    140140            CALL eos_fzp( sss_m(:,:), ztfrz(:,:) ) 
    141141            WHERE( sf_t_rnf(1)%fnow(:,:,1) == -999._wp )             ! if missing data value use SST as runoffs temperature 
    142                rnf_tsc(:,:,jp_tem) = sst_m(:,:) * rnf(:,:) * r1_rau0 
     142               rnf_tsc(:,:,jp_tem) = sst_m(:,:) * rnf(:,:) * r1_rho0 
    143143            END WHERE 
    144144         ELSE                                                        ! use SST as runoffs temperature 
    145145            !CEOD River is fresh water so must at least be 0 unless we consider ice 
    146             rnf_tsc(:,:,jp_tem) = MAX( sst_m(:,:), 0.0_wp ) * rnf(:,:) * r1_rau0 
     146            rnf_tsc(:,:,jp_tem) = MAX( sst_m(:,:), 0.0_wp ) * rnf(:,:) * r1_rho0 
    147147         ENDIF 
    148148         !                                                           ! use runoffs salinity data 
    149          IF( ln_rnf_sal )   rnf_tsc(:,:,jp_sal) = ( sf_s_rnf(1)%fnow(:,:,1) ) * rnf(:,:) * r1_rau0 
     149         IF( ln_rnf_sal )   rnf_tsc(:,:,jp_sal) = ( sf_s_rnf(1)%fnow(:,:,1) ) * rnf(:,:) * r1_rho0 
    150150         !                                                           ! else use S=0 for runoffs (done one for all in the init) 
    151151                                         CALL iom_put( 'runoffs'     , rnf(:,:)                         )   ! output runoff mass flux 
    152          IF( iom_use('hflx_rnf_cea') )   CALL iom_put( 'hflx_rnf_cea', rnf_tsc(:,:,jp_tem) * rau0 * rcp )   ! output runoff sensible heat (W/m2) 
     152         IF( iom_use('hflx_rnf_cea') )   CALL iom_put( 'hflx_rnf_cea', rnf_tsc(:,:,jp_tem) * rho0 * rcp )   ! output runoff sensible heat (W/m2) 
    153153      ENDIF 
    154154      ! 
     
    210210            DO_2D_11_11 
    211211               DO jk = 1, nk_rnf(ji,jj) 
    212                   phdivn(ji,jj,jk) = phdivn(ji,jj,jk) - ( rnf(ji,jj) + rnf_b(ji,jj) ) * zfact * r1_rau0 / h_rnf(ji,jj) 
     212                  phdivn(ji,jj,jk) = phdivn(ji,jj,jk) - ( rnf(ji,jj) + rnf_b(ji,jj) ) * zfact * r1_rho0 / h_rnf(ji,jj) 
    213213               END DO 
    214214            END_2D 
     
    221221               !                          ! apply the runoff input flow 
    222222               DO jk = 1, nk_rnf(ji,jj) 
    223                   phdivn(ji,jj,jk) = phdivn(ji,jj,jk) - ( rnf(ji,jj) + rnf_b(ji,jj) ) * zfact * r1_rau0 / h_rnf(ji,jj) 
     223                  phdivn(ji,jj,jk) = phdivn(ji,jj,jk) - ( rnf(ji,jj) + rnf_b(ji,jj) ) * zfact * r1_rho0 / h_rnf(ji,jj) 
    224224               END DO 
    225225            END_2D 
     
    227227      ELSE                       !==   runoff put only at the surface   ==! 
    228228         h_rnf (:,:)   = e3t (:,:,1,Kmm)        ! update h_rnf to be depth of top box 
    229          phdivn(:,:,1) = phdivn(:,:,1) - ( rnf(:,:) + rnf_b(:,:) ) * zfact * r1_rau0 / e3t(:,:,1,Kmm) 
     229         phdivn(:,:,1) = phdivn(:,:,1) - ( rnf(:,:) + rnf_b(:,:) ) * zfact * r1_rho0 / e3t(:,:,1,Kmm) 
    230230      ENDIF 
    231231      ! 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OCE/TDE/tide_mod.F90

    r12406 r12443  
    6464   INTEGER , PUBLIC ::   nb_harmo            !: Number of active tidal components 
    6565   REAL(wp), PUBLIC ::   rn_tide_ramp_dt     !: 
    66    REAL(wp), PUBLIC ::   rn_scal_load        !: 
     66   REAL(wp), PUBLIC ::   rn_load        !: 
    6767   CHARACTER(lc), PUBLIC ::   cn_tide_load   !:  
    6868   REAL(wp)         ::   rn_tide_gamma       ! Tidal tilt factor 
     
    124124      NAMELIST/nam_tide/ln_tide, nn_tide_var, ln_tide_dia, ln_tide_pot, rn_tide_gamma, & 
    125125         &              ln_scal_load, ln_read_load, cn_tide_load,         & 
    126          &              ln_tide_ramp, rn_scal_load, rn_tide_ramp_dt,      & 
     126         &              ln_tide_ramp, rn_load, rn_tide_ramp_dt,      & 
    127127         &              sn_tide_cnames 
    128128      !!---------------------------------------------------------------------- 
     
    153153            WRITE(numout,*) '         Read load potential from file           ln_read_load    = ', ln_read_load 
    154154            WRITE(numout,*) '         Apply ramp on tides at startup          ln_tide_ramp    = ', ln_tide_ramp 
    155             WRITE(numout,*) '         Fraction of SSH used in scal. approx.   rn_scal_load    = ', rn_scal_load 
     155            WRITE(numout,*) '         Fraction of SSH used in scal. approx.   rn_load    = ', rn_load 
    156156            WRITE(numout,*) '         Duration (days) of ramp                 rn_tide_ramp_dt = ', rn_tide_ramp_dt 
    157157         ENDIF 
    158158      ELSE 
    159          rn_scal_load = 0._wp  
     159         rn_load = 0._wp  
    160160 
    161161         IF(lwp) WRITE(numout,*) 
     
    193193      IF( nb_harmo == 0 )   CALL ctl_stop( 'tide_init : No tidal components set in nam_tide' ) 
    194194      ! 
    195       IF (.NOT.ln_scal_load ) rn_scal_load = 0._wp 
     195      IF (.NOT.ln_scal_load ) rn_load = 0._wp 
    196196      ! 
    197197      ALLOCATE( amp_pot(jpi,jpj,nb_harmo),                      & 
     
    754754      ! 
    755755      IF( ln_tide_dia ) THEN          ! Output total tidal potential (incl. load potential) 
    756          IF ( iom_use( "tide_pot" ) ) CALL iom_put( "tide_pot", pot_astro(:,:) + rn_scal_load * ssh(:,:,Kmm) ) 
     756         IF ( iom_use( "tide_pot" ) ) CALL iom_put( "tide_pot", pot_astro(:,:) + rn_load * ssh(:,:,Kmm) ) 
    757757      END IF 
    758758      ! 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OCE/TRA/eosbn2.F90

    r12377 r12443  
    191191      !!                   ***  ROUTINE eos_insitu  *** 
    192192      !! 
    193       !! ** Purpose :   Compute the in situ density (ratio rho/rau0) from 
     193      !! ** Purpose :   Compute the in situ density (ratio rho/rho0) from 
    194194      !!       potential temperature and salinity using an equation of state 
    195195      !!       selected in the nameos namelist 
    196196      !! 
    197       !! ** Method  :   prd(t,s,z) = ( rho(t,s,z) - rau0 ) / rau0 
     197      !! ** Method  :   prd(t,s,z) = ( rho(t,s,z) - rho0 ) / rho0 
    198198      !!         with   prd    in situ density anomaly      no units 
    199199      !!                t      TEOS10: CT or EOS80: PT      Celsius 
     
    201201      !!                z      depth                        meters 
    202202      !!                rho    in situ density              kg/m^3 
    203       !!                rau0   reference density            kg/m^3 
     203      !!                rho0   reference density            kg/m^3 
    204204      !! 
    205205      !!     ln_teos10 : polynomial TEOS-10 equation of state is used for rho(t,s,z). 
     
    210210      !! 
    211211      !!     ln_seos : simplified equation of state 
    212       !!              prd(t,s,z) = ( -a0*(1+lambda/2*(T-T0)+mu*z+nu*(S-S0))*(T-T0) + b0*(S-S0) ) / rau0 
     212      !!              prd(t,s,z) = ( -a0*(1+lambda/2*(T-T0)+mu*z+nu*(S-S0))*(T-T0) + b0*(S-S0) ) / rho0 
    213213      !!              linear case function of T only: rn_alpha<>0, other coefficients = 0 
    214214      !!              linear eos function of T and S: rn_alpha and rn_beta<>0, other coefficients=0 
     
    267267            zn  = ( ( zn3 * zh + zn2 ) * zh + zn1 ) * zh + zn0 
    268268            ! 
    269             prd(ji,jj,jk) = (  zn * r1_rau0 - 1._wp  ) * ztm  ! density anomaly (masked) 
     269            prd(ji,jj,jk) = (  zn * r1_rho0 - 1._wp  ) * ztm  ! density anomaly (masked) 
    270270            ! 
    271271         END_3D 
     
    283283               &  - rn_nu * zt * zs 
    284284               !                                  
    285             prd(ji,jj,jk) = zn * r1_rau0 * ztm                ! density anomaly (masked) 
     285            prd(ji,jj,jk) = zn * r1_rho0 * ztm                ! density anomaly (masked) 
    286286         END_3D 
    287287         ! 
     
    299299      !!                  ***  ROUTINE eos_insitu_pot  *** 
    300300      !! 
    301       !! ** Purpose :   Compute the in situ density (ratio rho/rau0) and the 
     301      !! ** Purpose :   Compute the in situ density (ratio rho/rho0) and the 
    302302      !!      potential volumic mass (Kg/m3) from potential temperature and 
    303303      !!      salinity fields using an equation of state selected in the 
     
    379379                  prhop(ji,jj,jk) = prhop(ji,jj,jk) + zn0_sto(jsmp)                      ! potential density referenced at the surface 
    380380                  ! 
    381                   prd(ji,jj,jk) = prd(ji,jj,jk) + (  zn_sto(jsmp) * r1_rau0 - 1._wp  )   ! density anomaly (masked) 
     381                  prd(ji,jj,jk) = prd(ji,jj,jk) + (  zn_sto(jsmp) * r1_rho0 - 1._wp  )   ! density anomaly (masked) 
    382382               END DO 
    383383               prhop(ji,jj,jk) = 0.5_wp * prhop(ji,jj,jk) * ztm / nn_sto_eos 
     
    419419               prhop(ji,jj,jk) = zn0 * ztm                           ! potential density referenced at the surface 
    420420               ! 
    421                prd(ji,jj,jk) = (  zn * r1_rau0 - 1._wp  ) * ztm      ! density anomaly (masked) 
     421               prd(ji,jj,jk) = (  zn * r1_rho0 - 1._wp  ) * ztm      ! density anomaly (masked) 
    422422            END_3D 
    423423         ENDIF 
     
    434434               &  + rn_b0 * ( 1._wp - 0.5_wp*rn_lambda2*zs ) * zs   & 
    435435               &  - rn_nu * zt * zs 
    436             prhop(ji,jj,jk) = ( rau0 + zn ) * ztm 
     436            prhop(ji,jj,jk) = ( rho0 + zn ) * ztm 
    437437            !                                                     ! density anomaly (masked) 
    438438            zn = zn - ( rn_a0 * rn_mu1 * zt + rn_b0 * rn_mu2 * zs ) * zh 
    439             prd(ji,jj,jk) = zn * r1_rau0 * ztm 
     439            prd(ji,jj,jk) = zn * r1_rho0 * ztm 
    440440            ! 
    441441         END_3D 
     
    454454      !!                  ***  ROUTINE eos_insitu_2d  *** 
    455455      !! 
    456       !! ** Purpose :   Compute the in situ density (ratio rho/rau0) from 
     456      !! ** Purpose :   Compute the in situ density (ratio rho/rho0) from 
    457457      !!      potential temperature and salinity using an equation of state 
    458458      !!      selected in the nameos namelist. * 2D field case 
     
    508508            zn  = ( ( zn3 * zh + zn2 ) * zh + zn1 ) * zh + zn0 
    509509            ! 
    510             prd(ji,jj) = zn * r1_rau0 - 1._wp               ! unmasked in situ density anomaly 
     510            prd(ji,jj) = zn * r1_rho0 - 1._wp               ! unmasked in situ density anomaly 
    511511            ! 
    512512         END_2D 
     
    524524               &  - rn_nu * zt * zs 
    525525               ! 
    526             prd(ji,jj) = zn * r1_rau0               ! unmasked in situ density anomaly 
     526            prd(ji,jj) = zn * r1_rho0               ! unmasked in situ density anomaly 
    527527            ! 
    528528         END_2D 
     
    588588            zn  = ( ( zn3 * zh + zn2 ) * zh + zn1 ) * zh + zn0 
    589589            ! 
    590             pab(ji,jj,jk,jp_tem) = zn * r1_rau0 * ztm 
     590            pab(ji,jj,jk,jp_tem) = zn * r1_rho0 * ztm 
    591591            ! 
    592592            ! beta 
     
    609609            zn  = ( ( zn3 * zh + zn2 ) * zh + zn1 ) * zh + zn0 
    610610            ! 
    611             pab(ji,jj,jk,jp_sal) = zn / zs * r1_rau0 * ztm 
     611            pab(ji,jj,jk,jp_sal) = zn / zs * r1_rho0 * ztm 
    612612            ! 
    613613         END_3D 
     
    622622            ! 
    623623            zn  = rn_a0 * ( 1._wp + rn_lambda1*zt + rn_mu1*zh ) + rn_nu*zs 
    624             pab(ji,jj,jk,jp_tem) = zn * r1_rau0 * ztm   ! alpha 
     624            pab(ji,jj,jk,jp_tem) = zn * r1_rho0 * ztm   ! alpha 
    625625            ! 
    626626            zn  = rn_b0 * ( 1._wp - rn_lambda2*zs - rn_mu2*zh ) - rn_nu*zt 
    627             pab(ji,jj,jk,jp_sal) = zn * r1_rau0 * ztm   ! beta 
     627            pab(ji,jj,jk,jp_sal) = zn * r1_rho0 * ztm   ! beta 
    628628            ! 
    629629         END_3D 
     
    694694            zn  = ( ( zn3 * zh + zn2 ) * zh + zn1 ) * zh + zn0 
    695695            ! 
    696             pab(ji,jj,jp_tem) = zn * r1_rau0 
     696            pab(ji,jj,jp_tem) = zn * r1_rho0 
    697697            ! 
    698698            ! beta 
     
    715715            zn  = ( ( zn3 * zh + zn2 ) * zh + zn1 ) * zh + zn0 
    716716            ! 
    717             pab(ji,jj,jp_sal) = zn / zs * r1_rau0 
     717            pab(ji,jj,jp_sal) = zn / zs * r1_rho0 
    718718            ! 
    719719            ! 
     
    729729            ! 
    730730            zn  = rn_a0 * ( 1._wp + rn_lambda1*zt + rn_mu1*zh ) + rn_nu*zs 
    731             pab(ji,jj,jp_tem) = zn * r1_rau0   ! alpha 
     731            pab(ji,jj,jp_tem) = zn * r1_rho0   ! alpha 
    732732            ! 
    733733            zn  = rn_b0 * ( 1._wp - rn_lambda2*zs - rn_mu2*zh ) - rn_nu*zt 
    734             pab(ji,jj,jp_sal) = zn * r1_rau0   ! beta 
     734            pab(ji,jj,jp_sal) = zn * r1_rho0   ! beta 
    735735            ! 
    736736         END_2D 
     
    799799         zn  = ( ( zn3 * zh + zn2 ) * zh + zn1 ) * zh + zn0 
    800800         ! 
    801          pab(jp_tem) = zn * r1_rau0 
     801         pab(jp_tem) = zn * r1_rho0 
    802802         ! 
    803803         ! beta 
     
    820820         zn  = ( ( zn3 * zh + zn2 ) * zh + zn1 ) * zh + zn0 
    821821         ! 
    822          pab(jp_sal) = zn / zs * r1_rau0 
     822         pab(jp_sal) = zn / zs * r1_rho0 
    823823         ! 
    824824         ! 
     
    831831         ! 
    832832         zn  = rn_a0 * ( 1._wp + rn_lambda1*zt + rn_mu1*zh ) + rn_nu*zs 
    833          pab(jp_tem) = zn * r1_rau0   ! alpha 
     833         pab(jp_tem) = zn * r1_rho0   ! alpha 
    834834         ! 
    835835         zn  = rn_b0 * ( 1._wp - rn_lambda2*zs - rn_mu2*zh ) - rn_nu*zt 
    836          pab(jp_sal) = zn * r1_rau0   ! beta 
     836         pab(jp_sal) = zn * r1_rho0   ! beta 
    837837         ! 
    838838      CASE DEFAULT 
     
    10521052      !! ** Method  :   PE is defined analytically as the vertical  
    10531053      !!                   primitive of EOS times -g integrated between 0 and z>0. 
    1054       !!                pen is the nonlinear bsq-PE anomaly: pen = ( PE - rau0 gz ) / rau0 gz - rd 
     1054      !!                pen is the nonlinear bsq-PE anomaly: pen = ( PE - rho0 gz ) / rho0 gz - rd 
    10551055      !!                                                      = 1/z * /int_0^z rd dz - rd  
    10561056      !!                                where rd is the density anomaly (see eos_rhd function) 
    10571057      !!                ab_pe are partial derivatives of PE anomaly with respect to T and S: 
    1058       !!                    ab_pe(1) = - 1/(rau0 gz) * dPE/dT + drd/dT = - d(pen)/dT 
    1059       !!                    ab_pe(2) =   1/(rau0 gz) * dPE/dS + drd/dS =   d(pen)/dS 
     1058      !!                    ab_pe(1) = - 1/(rho0 gz) * dPE/dT + drd/dT = - d(pen)/dT 
     1059      !!                    ab_pe(2) =   1/(rho0 gz) * dPE/dS + drd/dS =   d(pen)/dS 
    10601060      !! 
    10611061      !! ** Action  : - pen         : PE anomaly given at T-points 
     
    11031103            zn  = ( zn2 * zh + zn1 ) * zh + zn0 
    11041104            ! 
    1105             ppen(ji,jj,jk)  = zn * zh * r1_rau0 * ztm 
     1105            ppen(ji,jj,jk)  = zn * zh * r1_rho0 * ztm 
    11061106            ! 
    11071107            ! alphaPE non-linear anomaly 
     
    11181118            zn  = ( zn2 * zh + zn1 ) * zh + zn0 
    11191119            !                               
    1120             pab_pe(ji,jj,jk,jp_tem) = zn * zh * r1_rau0 * ztm 
     1120            pab_pe(ji,jj,jk,jp_tem) = zn * zh * r1_rho0 * ztm 
    11211121            ! 
    11221122            ! betaPE non-linear anomaly 
     
    11331133            zn  = ( zn2 * zh + zn1 ) * zh + zn0 
    11341134            !                               
    1135             pab_pe(ji,jj,jk,jp_sal) = zn / zs * zh * r1_rau0 * ztm 
     1135            pab_pe(ji,jj,jk,jp_sal) = zn / zs * zh * r1_rho0 * ztm 
    11361136            ! 
    11371137         END_3D 
     
    11441144            zh  = gdept(ji,jj,jk,Kmm)              ! depth in meters  at t-point 
    11451145            ztm = tmask(ji,jj,jk)                ! tmask 
    1146             zn  = 0.5_wp * zh * r1_rau0 * ztm 
     1146            zn  = 0.5_wp * zh * r1_rho0 * ztm 
    11471147            !                                    ! Potential Energy 
    11481148            ppen(ji,jj,jk) = ( rn_a0 * rn_mu1 * zt + rn_b0 * rn_mu2 * zs ) * zn 
     
    11861186      IF(lwm) WRITE( numond, nameos ) 
    11871187      ! 
    1188       rau0        = 1026._wp                 !: volumic mass of reference     [kg/m3] 
     1188      rho0        = 1026._wp                 !: volumic mass of reference     [kg/m3] 
    11891189      rcp         = 3991.86795711963_wp      !: heat capacity     [J/K] 
    11901190      ! 
     
    15981598            WRITE(numout,*) '   ==>>>   use of simplified eos:    ' 
    15991599            WRITE(numout,*) '              rhd(dT=T-10,dS=S-35,Z) = [-a0*(1+lambda1/2*dT+mu1*Z)*dT ' 
    1600             WRITE(numout,*) '                                       + b0*(1+lambda2/2*dT+mu2*Z)*dS - nu*dT*dS] / rau0' 
     1600            WRITE(numout,*) '                                       + b0*(1+lambda2/2*dT+mu2*Z)*dS - nu*dT*dS] / rho0' 
    16011601            WRITE(numout,*) '              with the following coefficients :' 
    16021602            WRITE(numout,*) '                 thermal exp. coef.    rn_a0      = ', rn_a0 
     
    16171617      END SELECT 
    16181618      ! 
    1619       rau0_rcp    = rau0 * rcp  
    1620       r1_rau0     = 1._wp / rau0 
     1619      rho0_rcp    = rho0 * rcp  
     1620      r1_rho0     = 1._wp / rho0 
    16211621      r1_rcp      = 1._wp / rcp 
    1622       r1_rau0_rcp = 1._wp / rau0_rcp  
     1622      r1_rho0_rcp = 1._wp / rho0_rcp  
    16231623      ! 
    16241624      IF(lwp) THEN 
     
    16351635      IF(lwp) WRITE(numout,*) 
    16361636      IF(lwp) WRITE(numout,*) '   Associated physical constant' 
    1637       IF(lwp) WRITE(numout,*) '      volumic mass of reference           rau0  = ', rau0   , ' kg/m^3' 
    1638       IF(lwp) WRITE(numout,*) '      1. / rau0                        r1_rau0  = ', r1_rau0, ' m^3/kg' 
     1637      IF(lwp) WRITE(numout,*) '      volumic mass of reference           rho0  = ', rho0   , ' kg/m^3' 
     1638      IF(lwp) WRITE(numout,*) '      1. / rho0                        r1_rho0  = ', r1_rho0, ' m^3/kg' 
    16391639      IF(lwp) WRITE(numout,*) '      ocean specific heat                 rcp   = ', rcp    , ' J/Kelvin' 
    1640       IF(lwp) WRITE(numout,*) '      rau0 * rcp                       rau0_rcp = ', rau0_rcp 
    1641       IF(lwp) WRITE(numout,*) '      1. / ( rau0 * rcp )           r1_rau0_rcp = ', r1_rau0_rcp 
     1640      IF(lwp) WRITE(numout,*) '      rho0 * rcp                       rho0_rcp = ', rho0_rcp 
     1641      IF(lwp) WRITE(numout,*) '      1. / ( rho0 * rcp )           r1_rho0_rcp = ', r1_rho0_rcp 
    16421642      ! 
    16431643   END SUBROUTINE eos_init 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OCE/TRA/traadv_fct.F90

    r12377 r12443  
    2020   USE diaptr         ! poleward transport diagnostics 
    2121   USE diaar5         ! AR5 diagnostics 
    22    USE phycst  , ONLY : rau0_rcp 
     22   USE phycst  , ONLY : rho0_rcp 
    2323   USE zdf_oce , ONLY : ln_zad_Aimp 
    2424   ! 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OCE/TRA/traatf.F90

    r12424 r12443  
    214214            ztd = pt(ji,jj,jk,jn,Kaa) - 2._wp * ztn + pt(ji,jj,jk,jn,Kbb)  ! time laplacian on tracers 
    215215            ! 
    216             pt(ji,jj,jk,jn,Kmm) = ztn + atfp * ztd                      ! pt <-- filtered pt 
     216            pt(ji,jj,jk,jn,Kmm) = ztn + rn_atfp * ztd                      ! pt <-- filtered pt 
    217217         END_3D 
    218218         ! 
     
    229229      !!  
    230230      !! ** Method  : - Apply a thickness weighted Asselin time filter on now fields. 
    231       !!             pt(Kmm)  = ( e3t(Kmm)*pt(Kmm) + atfp*[ e3t(Kbb)*pt(Kbb) - 2 e3t(Kmm)*pt(Kmm) + e3t_a*pt(Kaa) ] ) 
    232       !!                       /( e3t(Kmm)         + atfp*[ e3t(Kbb)         - 2 e3t(Kmm)         + e3t(Kaa)    ] ) 
     231      !!             pt(Kmm)  = ( e3t(Kmm)*pt(Kmm) + rn_atfp*[ e3t(Kbb)*pt(Kbb) - 2 e3t(Kmm)*pt(Kmm) + e3t_a*pt(Kaa) ] ) 
     232      !!                       /( e3t(Kmm)         + rn_atfp*[ e3t(Kbb)         - 2 e3t(Kmm)         + e3t(Kaa)    ] ) 
    233233      !! 
    234234      !! ** Action  : - pt(Kmm) ready for the next time step 
     
    272272      ENDIF 
    273273      zfact = 1._wp / p2dt 
    274       zfact1 = atfp * p2dt 
    275       zfact2 = zfact1 * r1_rau0 
     274      zfact1 = rn_atfp * p2dt 
     275      zfact2 = zfact1 * r1_rho0 
    276276      DO jn = 1, kjpt       
    277277         DO_3D_00_00( 1, jpkm1 ) 
     
    287287            ztc_d  = ztc_a  - 2. * ztc_n  + ztc_b 
    288288            ! 
    289             ze3t_f = ze3t_n + atfp * ze3t_d 
    290             ztc_f  = ztc_n  + atfp * ztc_d 
     289            ze3t_f = ze3t_n + rn_atfp * ze3t_d 
     290            ztc_f  = ztc_n  + rn_atfp * ztc_d 
    291291            ! 
    292292            ! Add asselin correction on scale factors: 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OCE/TRA/trabbc.F90

    r12377 r12443  
    6666      !!       ocean bottom can be computed once and is added to the temperature 
    6767      !!       trend juste above the bottom at each time step: 
    68       !!            ta = ta + Qsf / (rau0 rcp e3T) for k= mbkt 
     68      !!            ta = ta + Qsf / (rho0 rcp e3T) for k= mbkt 
    6969      !!       Where Qsf is the geothermal heat flux. 
    7070      !! 
     
    102102      ENDIF 
    103103      ! 
    104       CALL iom_put ( "hfgeou" , rau0_rcp * qgh_trd0(:,:) ) 
     104      CALL iom_put ( "hfgeou" , rho0_rcp * qgh_trd0(:,:) ) 
    105105      IF(sn_cfctl%l_prtctl)   CALL prt_ctl( tab3d_1=pts(:,:,:,jp_tem,Krhs), clinfo1=' bbc  - Ta: ', mask1=tmask, clinfo3='tra-ta' ) 
    106106      ! 
     
    162162         CASE ( 1 )                          !* constant flux 
    163163            IF(lwp) WRITE(numout,*) '   ==>>>   constant heat flux  =   ', rn_geoflx_cst 
    164             qgh_trd0(:,:) = r1_rau0_rcp * rn_geoflx_cst 
     164            qgh_trd0(:,:) = r1_rho0_rcp * rn_geoflx_cst 
    165165            ! 
    166166         CASE ( 2 )                          !* variable geothermal heat flux : read the geothermal fluxes in mW/m2 
     
    179179 
    180180            CALL fld_read( nit000, 1, sf_qgh )                         ! Read qgh data 
    181             qgh_trd0(:,:) = r1_rau0_rcp * sf_qgh(1)%fnow(:,:,1) * 1.e-3 ! conversion in W/m2 
     181            qgh_trd0(:,:) = r1_rho0_rcp * sf_qgh(1)%fnow(:,:,1) * 1.e-3 ! conversion in W/m2 
    182182            ! 
    183183         CASE DEFAULT 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OCE/TRA/tramle.F90

    r12377 r12443  
    4141 
    4242   REAL(wp) ::   r5_21 = 5.e0 / 21.e0   ! factor used in mle streamfunction computation 
    43    REAL(wp) ::   rb_c                   ! ML buoyancy criteria = g rho_c /rau0 where rho_c is defined in zdfmld 
     43   REAL(wp) ::   rb_c                   ! ML buoyancy criteria = g rho_c /rho0 where rho_c is defined in zdfmld 
    4444   REAL(wp) ::   rc_f                   ! MLE coefficient (= rn_ce / (5 km * fo) ) in nn_mle=1 case 
    4545 
     
    112112         zc = e3t(ji,jj,jk,Kmm) * REAL( MIN( MAX( 0, inml_mle(ji,jj)-jk ) , 1  )  )    ! zc being 0 outside the ML t-points 
    113113         zmld(ji,jj) = zmld(ji,jj) + zc 
    114          zbm (ji,jj) = zbm (ji,jj) + zc * (rau0 - rhop(ji,jj,jk) ) * r1_rau0 
     114         zbm (ji,jj) = zbm (ji,jj) + zc * (rho0 - rhop(ji,jj,jk) ) * r1_rho0 
    115115         zn2 (ji,jj) = zn2 (ji,jj) + zc * (rn2(ji,jj,jk)+rn2(ji,jj,jk+1))*0.5_wp 
    116116      END_3D 
     
    273273      IF( ln_mle ) THEN                ! MLE initialisation 
    274274         ! 
    275          rb_c = grav * rn_rho_c_mle /rau0        ! Mixed Layer buoyancy criteria 
     275         rb_c = grav * rn_rho_c_mle /rho0        ! Mixed Layer buoyancy criteria 
    276276         IF(lwp) WRITE(numout,*) 
    277277         IF(lwp) WRITE(numout,*) '      ML buoyancy criteria = ', rb_c, ' m/s2 ' 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OCE/TRA/traqsr.F90

    r12397 r12443  
    8787      !!         I(k) = Qsr*( rn_abs*EXP(z(k)/rn_si0) + (1.-rn_abs)*EXP(z(k)/rn_si1) ) 
    8888      !!         The temperature trend associated with the solar radiation penetration  
    89       !!         is given by : zta = 1/e3t dk[ I ] / (rau0*Cp) 
     89      !!         is given by : zta = 1/e3t dk[ I ] / (rho0*Cp) 
    9090      !!         At the bottom, boudary condition for the radiation is no flux : 
    9191      !!      all heat which has not been absorbed in the above levels is put 
     
    155155         ! 
    156156         DO jk = 1, nksr 
    157             qsr_hc(:,:,jk) = r1_rau0_rcp * ( etot3(:,:,jk) - etot3(:,:,jk+1) ) 
     157            qsr_hc(:,:,jk) = r1_rho0_rcp * ( etot3(:,:,jk) - etot3(:,:,jk+1) ) 
    158158         END DO 
    159159         ! 
     
    227227         ! 
    228228         DO_3D_00_00( 1, nksr ) 
    229             qsr_hc(ji,jj,jk) = r1_rau0_rcp * ( zea(ji,jj,jk) - zea(ji,jj,jk+1) ) 
     229            qsr_hc(ji,jj,jk) = r1_rho0_rcp * ( zea(ji,jj,jk) - zea(ji,jj,jk+1) ) 
    230230         END_3D 
    231231         ! 
     
    234234      CASE( np_2BD  )            !==  2-bands fluxes  ==! 
    235235         ! 
    236          zz0 =        rn_abs   * r1_rau0_rcp      ! surface equi-partition in 2-bands 
    237          zz1 = ( 1. - rn_abs ) * r1_rau0_rcp 
     236         zz0 =        rn_abs   * r1_rho0_rcp      ! surface equi-partition in 2-bands 
     237         zz1 = ( 1. - rn_abs ) * r1_rho0_rcp 
    238238         DO_3D_00_00( 1, nksr ) 
    239239            zc0 = zz0 * EXP( -gdepw(ji,jj,jk  ,Kmm)*xsi0r ) + zz1 * EXP( -gdepw(ji,jj,jk  ,Kmm)*xsi1r ) 
     
    252252      ! sea-ice: store the 1st ocean level attenuation coefficient 
    253253      DO_2D_00_00 
    254          IF( qsr(ji,jj) /= 0._wp ) THEN   ;   fraqsr_1lev(ji,jj) = qsr_hc(ji,jj,1) / ( r1_rau0_rcp * qsr(ji,jj) ) 
     254         IF( qsr(ji,jj) /= 0._wp ) THEN   ;   fraqsr_1lev(ji,jj) = qsr_hc(ji,jj,1) / ( r1_rho0_rcp * qsr(ji,jj) ) 
    255255         ELSE                             ;   fraqsr_1lev(ji,jj) = 1._wp 
    256256         ENDIF 
     
    262262         zetot(:,:,nksr+1:jpk) = 0._wp     ! below ~400m set to zero 
    263263         DO jk = nksr, 1, -1 
    264             zetot(:,:,jk) = zetot(:,:,jk+1) + qsr_hc(:,:,jk) * rau0_rcp 
     264            zetot(:,:,jk) = zetot(:,:,jk+1) + qsr_hc(:,:,jk) * rho0_rcp 
    265265         END DO          
    266266         CALL iom_put( 'qsr3d', zetot )   ! 3D distribution of shortwave Radiation 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OCE/TRA/trasbc.F90

    r12377 r12443  
    124124      !                             !==  Now sbc tracer content fields  ==! 
    125125      DO_2D_01_00 
    126          sbc_tsc(ji,jj,jp_tem) = r1_rau0_rcp * qns(ji,jj)   ! non solar heat flux 
    127          sbc_tsc(ji,jj,jp_sal) = r1_rau0     * sfx(ji,jj)   ! salt flux due to freezing/melting 
     126         sbc_tsc(ji,jj,jp_tem) = r1_rho0_rcp * qns(ji,jj)   ! non solar heat flux 
     127         sbc_tsc(ji,jj,jp_sal) = r1_rho0     * sfx(ji,jj)   ! salt flux due to freezing/melting 
    128128      END_2D 
    129129      IF( ln_linssh ) THEN                !* linear free surface   
    130130         DO_2D_01_00 
    131             sbc_tsc(ji,jj,jp_tem) = sbc_tsc(ji,jj,jp_tem) + r1_rau0 * emp(ji,jj) * pts(ji,jj,1,jp_tem,Kmm) 
    132             sbc_tsc(ji,jj,jp_sal) = sbc_tsc(ji,jj,jp_sal) + r1_rau0 * emp(ji,jj) * pts(ji,jj,1,jp_sal,Kmm) 
     131            sbc_tsc(ji,jj,jp_tem) = sbc_tsc(ji,jj,jp_tem) + r1_rho0 * emp(ji,jj) * pts(ji,jj,1,jp_tem,Kmm) 
     132            sbc_tsc(ji,jj,jp_sal) = sbc_tsc(ji,jj,jp_sal) + r1_rho0 * emp(ji,jj) * pts(ji,jj,1,jp_sal,Kmm) 
    133133         END_2D 
    134134         IF( iom_use('emp_x_sst') )   CALL iom_put( "emp_x_sst", emp (:,:) * pts(:,:,1,jp_tem,Kmm) ) 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OCE/TRD/trddyn.F90

    r12377 r12443  
    140140                              !                                    ! wind stress trends 
    141141                              ALLOCATE( z2dx(jpi,jpj) , z2dy(jpi,jpj) ) 
    142                               z2dx(:,:) = ( utau_b(:,:) + utau(:,:) ) / ( e3u(:,:,1,Kmm) * rau0 ) 
    143                               z2dy(:,:) = ( vtau_b(:,:) + vtau(:,:) ) / ( e3v(:,:,1,Kmm) * rau0 ) 
     142                              z2dx(:,:) = ( utau_b(:,:) + utau(:,:) ) / ( e3u(:,:,1,Kmm) * rho0 ) 
     143                              z2dy(:,:) = ( vtau_b(:,:) + vtau(:,:) ) / ( e3v(:,:,1,Kmm) * rho0 ) 
    144144                              CALL iom_put( "utrd_tau", z2dx ) 
    145145                              CALL iom_put( "vtrd_tau", z2dy ) 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OCE/TRD/trdglo.F90

    r12377 r12443  
    7676      INTEGER ::   ji, jj, jk      ! dummy loop indices 
    7777      INTEGER ::   ikbu, ikbv      ! local integers 
    78       REAL(wp)::   zvm, zvt, zvs, z1_2rau0   ! local scalars 
     78      REAL(wp)::   zvm, zvt, zvs, z1_2rho0   ! local scalars 
    7979      REAL(wp), DIMENSION(jpi,jpj)  :: ztswu, ztswv, z2dx, z2dy   ! 2D workspace  
    8080      !!---------------------------------------------------------------------- 
     
    125125            !                  
    126126            IF( ktrd == jpdyn_zdf ) THEN      ! zdf trend: compute separately the surface forcing trend 
    127                z1_2rau0 = 0.5_wp / rau0 
     127               z1_2rho0 = 0.5_wp / rho0 
    128128               DO_2D_10_10 
    129129                  zvt = ( utau_b(ji,jj) + utau(ji,jj) ) * tmask_i(ji+1,jj) * tmask_i(ji,jj) * umask(ji,jj,jk)   & 
    130                      &                                                     * z1_2rau0       * e1e2u(ji,jj) 
     130                     &                                                     * z1_2rho0       * e1e2u(ji,jj) 
    131131                  zvs = ( vtau_b(ji,jj) + vtau(ji,jj) ) * tmask_i(ji,jj+1) * tmask_i(ji,jj) * vmask(ji,jj,jk)   & 
    132                      &                                                     * z1_2rau0       * e1e2v(ji,jj) 
     132                     &                                                     * z1_2rho0       * e1e2v(ji,jj) 
    133133                  umo(jpdyn_tau) = umo(jpdyn_tau) + zvt 
    134134                  vmo(jpdyn_tau) = vmo(jpdyn_tau) + zvs 
     
    141141!               ! 
    142142!               IF( ln_drgimp ) THEN                   ! implicit drag case: compute separately the bottom friction  
    143 !                  z1_2rau0 = 0.5_wp / rau0 
     143!                  z1_2rho0 = 0.5_wp / rho0 
    144144!                  DO jj = 1, jpjm1 
    145145!                     DO ji = 1, jpim1 
     
    203203         CALL eos( ts(:,:,:,:,Kmm), rhd, rhop )       ! now potential density 
    204204 
    205          zcof = 0.5_wp / rau0             ! Density flux at w-point 
     205         zcof = 0.5_wp / rho0             ! Density flux at w-point 
    206206         zkz(:,:,1) = 0._wp 
    207207         DO jk = 2, jpk 
     
    209209         END DO 
    210210          
    211          zcof   = 0.5_wp / rau0           ! Density flux at u and v-points 
     211         zcof   = 0.5_wp / rho0           ! Density flux at u and v-points 
    212212         DO_3D_10_10( 1, jpkm1 ) 
    213213            zkx(ji,jj,jk) = zcof * e2u(ji,jj) * e3u(ji,jj,jk,Kmm) * uu(ji,jj,jk,Kmm) * ( rhop(ji,jj,jk) + rhop(ji+1,jj,jk) ) 
     
    347347 9546    FORMAT(' 0 < horizontal diffusion                                  : ', e20.13) 
    348348 9547    FORMAT(' 0 < vertical diffusion                                    : ', e20.13) 
    349  9548    FORMAT(' pressure gradient u2 = - 1/rau0 u.dz(rhop)                : ', e20.13, '  u.dz(rhop) =', e20.13) 
     349 9548    FORMAT(' pressure gradient u2 = - 1/rho0 u.dz(rhop)                : ', e20.13, '  u.dz(rhop) =', e20.13) 
    350350         ! 
    351351         ! Save potential to kinetic energy conversion for next time step 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OCE/TRD/trdken.F90

    r12377 r12443  
    102102      zke(:,1, : ) = 0._wp 
    103103      DO_3D_01_01( 1, jpkm1 ) 
    104          zke(ji,jj,jk) = 0.5_wp * rau0 *( uu(ji  ,jj,jk,Kmm) * putrd(ji  ,jj,jk) * bu(ji  ,jj,jk)  & 
     104         zke(ji,jj,jk) = 0.5_wp * rho0 *( uu(ji  ,jj,jk,Kmm) * putrd(ji  ,jj,jk) * bu(ji  ,jj,jk)  & 
    105105            &                           + uu(ji-1,jj,jk,Kmm) * putrd(ji-1,jj,jk) * bu(ji-1,jj,jk)  & 
    106106            &                           + vv(ji,jj  ,jk,Kmm) * pvtrd(ji,jj  ,jk) * bv(ji,jj  ,jk)  & 
     
    123123                           zke2d(1,:) = 0._wp   ;   zke2d(:,1) = 0._wp 
    124124                           DO_2D_01_01 
    125                               zke2d(ji,jj) = r1_rau0 * 0.5_wp * (   z2dx(ji,jj) + z2dx(ji-1,jj)   & 
     125                              zke2d(ji,jj) = r1_rho0 * 0.5_wp * (   z2dx(ji,jj) + z2dx(ji-1,jj)   & 
    126126                              &                                   + z2dy(ji,jj) + z2dy(ji,jj-1)   ) * r1_bt(ji,jj,1) 
    127127                           END_2D 
     
    207207      ! 
    208208      ! Local constant initialization  
    209       zcoef = - rau0 * grav * 0.5_wp       
     209      zcoef = - rho0 * grav * 0.5_wp       
    210210       
    211211      !  Surface value (also valid in partial step case) 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OCE/TRD/trdvor.F90

    r12406 r12443  
    105105         ztswu(:,:) = 0.e0   ;   ztswv(:,:) = 0.e0 
    106106         DO_2D_00_00 
    107             ztswu(ji,jj) = 0.5 * ( utau_b(ji,jj) + utau(ji,jj) ) / ( e3u(ji,jj,1,Kmm) * rau0 ) 
    108             ztswv(ji,jj) = 0.5 * ( vtau_b(ji,jj) + vtau(ji,jj) ) / ( e3v(ji,jj,1,Kmm) * rau0 ) 
     107            ztswu(ji,jj) = 0.5 * ( utau_b(ji,jj) + utau(ji,jj) ) / ( e3u(ji,jj,1,Kmm) * rho0 ) 
     108            ztswv(ji,jj) = 0.5 * ( vtau_b(ji,jj) + vtau(ji,jj) ) / ( e3v(ji,jj,1,Kmm) * rho0 ) 
    109109         END_2D 
    110110         ! 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OCE/ZDF/zdfgls.F90

    r12406 r12443  
    170170         ! 
    171171         ! surface friction 
    172          ustar2_surf(ji,jj) = r1_rau0 * taum(ji,jj) * tmask(ji,jj,1) 
     172         ustar2_surf(ji,jj) = r1_rho0 * taum(ji,jj) * tmask(ji,jj,1) 
    173173         !    
    174174!!gm Rq we may add here r_ke0(_top/_bot) ?  ==>> think about that... 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OCE/ZDF/zdfiwm.F90

    r12377 r12443  
    8787      !!              This is divided into three components: 
    8888      !!                 1. Bottom-intensified low-mode dissipation at critical slopes 
    89       !!                     zemx_iwm(z) = ( ecri_iwm / rau0 ) * EXP( -(H-z)/hcri_iwm ) 
     89      !!                     zemx_iwm(z) = ( ecri_iwm / rho0 ) * EXP( -(H-z)/hcri_iwm ) 
    9090      !!                                   / ( 1. - EXP( - H/hcri_iwm ) ) * hcri_iwm 
    9191      !!              where hcri_iwm is the characteristic length scale of the bottom  
    9292      !!              intensification, ecri_iwm a map of available power, and H the ocean depth. 
    9393      !!                 2. Pycnocline-intensified low-mode dissipation 
    94       !!                     zemx_iwm(z) = ( epyc_iwm / rau0 ) * ( sqrt(rn2(z))^nn_zpyc ) 
     94      !!                     zemx_iwm(z) = ( epyc_iwm / rho0 ) * ( sqrt(rn2(z))^nn_zpyc ) 
    9595      !!                                   / SUM( sqrt(rn2(z))^nn_zpyc * e3w(z) ) 
    9696      !!              where epyc_iwm is a map of available power, and nn_zpyc 
     
    9898      !!              energy dissipation. 
    9999      !!                 3. WKB-height dependent high mode dissipation 
    100       !!                     zemx_iwm(z) = ( ebot_iwm / rau0 ) * rn2(z) * EXP(-z_wkb(z)/hbot_iwm) 
     100      !!                     zemx_iwm(z) = ( ebot_iwm / rho0 ) * rn2(z) * EXP(-z_wkb(z)/hbot_iwm) 
    101101      !!                                   / SUM( rn2(z) * EXP(-z_wkb(z)/hbot_iwm) * e3w(z) ) 
    102102      !!              where hbot_iwm is the characteristic length scale of the WKB bottom  
     
    151151      DO_2D_11_11 
    152152         zhdep(ji,jj) = gdepw_0(ji,jj,mbkt(ji,jj)+1)       ! depth of the ocean 
    153          zfact(ji,jj) = rau0 * (  1._wp - EXP( -zhdep(ji,jj) / hcri_iwm(ji,jj) )  ) 
     153         zfact(ji,jj) = rho0 * (  1._wp - EXP( -zhdep(ji,jj) / hcri_iwm(ji,jj) )  ) 
    154154         IF( zfact(ji,jj) /= 0._wp )   zfact(ji,jj) = ecri_iwm(ji,jj) / zfact(ji,jj) 
    155155      END_2D 
     
    181181         ! 
    182182         DO_2D_11_11 
    183             IF( zfact(ji,jj) /= 0 )   zfact(ji,jj) = epyc_iwm(ji,jj) / ( rau0 * zfact(ji,jj) ) 
     183            IF( zfact(ji,jj) /= 0 )   zfact(ji,jj) = epyc_iwm(ji,jj) / ( rho0 * zfact(ji,jj) ) 
    184184         END_2D 
    185185         ! 
     
    196196         ! 
    197197         DO_2D_11_11 
    198             IF( zfact(ji,jj) /= 0 )   zfact(ji,jj) = epyc_iwm(ji,jj) / ( rau0 * zfact(ji,jj) ) 
     198            IF( zfact(ji,jj) /= 0 )   zfact(ji,jj) = epyc_iwm(ji,jj) / ( rho0 * zfact(ji,jj) ) 
    199199         END_2D 
    200200         ! 
     
    243243      ! 
    244244      DO_2D_11_11 
    245          IF( zfact(ji,jj) /= 0 )   zfact(ji,jj) = ebot_iwm(ji,jj) / ( rau0 * zfact(ji,jj) ) 
     245         IF( zfact(ji,jj) /= 0 )   zfact(ji,jj) = ebot_iwm(ji,jj) / ( rho0 * zfact(ji,jj) ) 
    246246      END_2D 
    247247      ! 
     
    255255      ! Calculate molecular kinematic viscosity 
    256256      znu_t(:,:,:) = 1.e-4_wp * (  17.91_wp - 0.53810_wp * ts(:,:,:,jp_tem,Kmm) + 0.00694_wp * ts(:,:,:,jp_tem,Kmm) * ts(:,:,:,jp_tem,Kmm)  & 
    257          &                                  + 0.02305_wp * ts(:,:,:,jp_sal,Kmm)  ) * tmask(:,:,:) * r1_rau0 
     257         &                                  + 0.02305_wp * ts(:,:,:,jp_sal,Kmm)  ) * tmask(:,:,:) * r1_rho0 
    258258      DO jk = 2, jpkm1 
    259259         znu_w(:,:,jk) = 0.5_wp * ( znu_t(:,:,jk-1) + znu_t(:,:,jk) ) * wmask(:,:,jk) 
     
    293293         END_3D 
    294294         CALL mpp_sum( 'zdfiwm', zztmp ) 
    295          zztmp = rau0 * zztmp ! Global integral of rauo * Kz * N^2 = power contributing to mixing  
     295         zztmp = rho0 * zztmp ! Global integral of rauo * Kz * N^2 = power contributing to mixing  
    296296         ! 
    297297         IF(lwp) THEN 
     
    337337                                    !* output useful diagnostics: Kz*N^2 ,  
    338338!!gm Kz*N2 should take into account the ratio avs/avt if it is used.... (see diaar5) 
    339                                     !  vertical integral of rau0 * Kz * N^2 , energy density (zemx_iwm) 
     339                                    !  vertical integral of rho0 * Kz * N^2 , energy density (zemx_iwm) 
    340340      IF( iom_use("bflx_iwm") .OR. iom_use("pcmap_iwm") ) THEN 
    341341         ALLOCATE( z2d(jpi,jpj) , z3d(jpi,jpj,jpk) ) 
     
    345345            z2d(:,:) = z2d(:,:) + e3w(:,:,jk,Kmm) * z3d(:,:,jk) * wmask(:,:,jk) 
    346346         END DO 
    347          z2d(:,:) = rau0 * z2d(:,:) 
     347         z2d(:,:) = rho0 * z2d(:,:) 
    348348         CALL iom_put( "bflx_iwm", z3d ) 
    349349         CALL iom_put( "pcmap_iwm", z2d ) 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OCE/ZDF/zdfmxl.F90

    r12377 r12443  
    9797      nmln(:,:)  = nlb10               ! Initialization to the number of w ocean point 
    9898      hmlp(:,:)  = 0._wp               ! here hmlp used as a dummy variable, integrating vertically N^2 
    99       zN2_c = grav * rho_c * r1_rau0   ! convert density criteria into N^2 criteria 
     99      zN2_c = grav * rho_c * r1_rho0   ! convert density criteria into N^2 criteria 
    100100      DO_3D_11_11( nlb10, jpkm1 ) 
    101101         ikt = mbkt(ji,jj) 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OCE/ZDF/zdfosm.F90

    r12406 r12443  
    300300     DO_2D_00_00 
    301301        ! Surface downward irradiance (so always +ve) 
    302         zrad0(ji,jj) = qsr(ji,jj) * r1_rau0_rcp 
     302        zrad0(ji,jj) = qsr(ji,jj) * r1_rho0_rcp 
    303303        ! Downwards irradiance at base of boundary layer 
    304304        zradh(ji,jj) = zrad0(ji,jj) * ( zz0 * EXP( -hbl(ji,jj)/rn_si0 ) + zz1 * EXP( -hbl(ji,jj)/rn_si1) ) 
     
    312312        zbeta    = rab_n(ji,jj,1,jp_sal) 
    313313        ! Upwards surface Temperature flux for non-local term 
    314         zwth0(ji,jj) = - qns(ji,jj) * r1_rau0_rcp * tmask(ji,jj,1) 
     314        zwth0(ji,jj) = - qns(ji,jj) * r1_rho0_rcp * tmask(ji,jj,1) 
    315315        ! Upwards surface salinity flux for non-local term 
    316         zws0(ji,jj) = - ( ( emp(ji,jj)-rnf(ji,jj) ) * ts(ji,jj,1,jp_sal,Kmm)  + sfx(ji,jj) ) * r1_rau0 * tmask(ji,jj,1) 
     316        zws0(ji,jj) = - ( ( emp(ji,jj)-rnf(ji,jj) ) * ts(ji,jj,1,jp_sal,Kmm)  + sfx(ji,jj) ) * r1_rho0 * tmask(ji,jj,1) 
    317317        ! Non radiative upwards surface buoyancy flux 
    318318        zwb0(ji,jj) = grav * zthermal * zwth0(ji,jj) -  grav * zbeta * zws0(ji,jj) 
     
    324324        zwbav(ji,jj) = grav  * zthermal * zwthav(ji,jj) - grav  * zbeta * zwsav(ji,jj) 
    325325        ! Surface upward velocity fluxes 
    326         zuw0(ji,jj) = -utau(ji,jj) * r1_rau0 * tmask(ji,jj,1) 
    327         zvw0(ji,jj) = -vtau(ji,jj) * r1_rau0 * tmask(ji,jj,1) 
     326        zuw0(ji,jj) = -utau(ji,jj) * r1_rho0 * tmask(ji,jj,1) 
     327        zvw0(ji,jj) = -vtau(ji,jj) * r1_rho0 * tmask(ji,jj,1) 
    328328        ! Friction velocity (zustar), at T-point : LMD94 eq. 2 
    329329        zustar(ji,jj) = MAX( SQRT( SQRT( zuw0(ji,jj) * zuw0(ji,jj) + zvw0(ji,jj) * zvw0(ji,jj) ) ), 1.0e-8 ) 
     
    12501250            IF ( iom_use("us_x") ) CALL iom_put( "us_x", tmask(:,:,1)*zustke*zcos_wind )   ! x surface Stokes drift 
    12511251            IF ( iom_use("us_y") ) CALL iom_put( "us_y", tmask(:,:,1)*zustke*zsin_wind )  ! y surface Stokes drift 
    1252             IF ( iom_use("wind_wave_abs_power") ) CALL iom_put( "wind_wave_abs_power", 1000.*rau0*tmask(:,:,1)*zustar**2*zustke ) 
     1252            IF ( iom_use("wind_wave_abs_power") ) CALL iom_put( "wind_wave_abs_power", 1000.*rho0*tmask(:,:,1)*zustar**2*zustke ) 
    12531253         ! Stokes drift read in from sbcwave  (=2). 
    12541254         CASE(2) 
    12551255            IF ( iom_use("us_x") ) CALL iom_put( "us_x", ut0sd )               ! x surface Stokes drift 
    12561256            IF ( iom_use("us_y") ) CALL iom_put( "us_y", vt0sd )               ! y surface Stokes drift 
    1257             IF ( iom_use("wind_wave_abs_power") ) CALL iom_put( "wind_wave_abs_power", 1000.*rau0*tmask(:,:,1)*zustar**2* & 
     1257            IF ( iom_use("wind_wave_abs_power") ) CALL iom_put( "wind_wave_abs_power", 1000.*rho0*tmask(:,:,1)*zustar**2* & 
    12581258                 & SQRT(ut0sd**2 + vt0sd**2 ) ) 
    12591259         END SELECT 
     
    12711271         IF ( iom_use("zwstrl") ) CALL iom_put( "zwstrl", tmask(:,:,1)*zwstrl )         ! Langmuir velocity scale 
    12721272         IF ( iom_use("zustar") ) CALL iom_put( "zustar", tmask(:,:,1)*zustar )         ! friction velocity scale 
    1273          IF ( iom_use("wind_power") ) CALL iom_put( "wind_power", 1000.*rau0*tmask(:,:,1)*zustar**3 ) ! BL depth internal to zdf_osm routine 
    1274          IF ( iom_use("wind_wave_power") ) CALL iom_put( "wind_wave_power", 1000.*rau0*tmask(:,:,1)*zustar**2*zustke ) 
     1273         IF ( iom_use("wind_power") ) CALL iom_put( "wind_power", 1000.*rho0*tmask(:,:,1)*zustar**3 ) ! BL depth internal to zdf_osm routine 
     1274         IF ( iom_use("wind_wave_power") ) CALL iom_put( "wind_wave_power", 1000.*rho0*tmask(:,:,1)*zustar**2*zustke ) 
    12751275         IF ( iom_use("zhbl") ) CALL iom_put( "zhbl", tmask(:,:,1)*zhbl )               ! BL depth internal to zdf_osm routine 
    12761276         IF ( iom_use("zhml") ) CALL iom_put( "zhml", tmask(:,:,1)*zhml )               ! ML depth internal to zdf_osm routine 
     
    15001500     imld_rst(:,:)  = nlb10         ! Initialization to the number of w ocean point 
    15011501     hbl(:,:)  = 0._wp              ! here hbl used as a dummy variable, integrating vertically N^2 
    1502      zN2_c = grav * rho_c * r1_rau0 ! convert density criteria into N^2 criteria 
     1502     zN2_c = grav * rho_c * r1_rho0 ! convert density criteria into N^2 criteria 
    15031503     ! 
    15041504     hbl(:,:)  = 0._wp              ! here hbl used as a dummy variable, integrating vertically N^2 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OCE/ZDF/zdfric.F90

    r12377 r12443  
    174174         ! 
    175175         DO_2D_00_00 
    176             zustar = SQRT( taum(ji,jj) * r1_rau0 ) 
     176            zustar = SQRT( taum(ji,jj) * r1_rho0 ) 
    177177            zhek   = rn_ekmfc * zustar / ( ABS( ff_t(ji,jj) ) + rsmall )   ! Ekman depth 
    178178            zh_ekm(ji,jj) = MAX(  rn_mldmin , MIN( zhek , rn_mldmax )  )   ! set allowed range 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OCE/ZDF/zdftke.F90

    r12406 r12443  
    206206      !!-------------------------------------------------------------------- 
    207207      ! 
    208       zbbrau = rn_ebb / rau0       ! Local constant initialisation 
     208      zbbrau = rn_ebb / rho0       ! Local constant initialisation 
    209209      zfact1 = -.5_wp * rn_Dt  
    210210      zfact2 = 1.5_wp * rn_Dt * rn_ediss 
     
    228228      !                     !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< 
    229229      ! 
    230       !   en(bot)   = (ebb0/rau0)*0.5*sqrt(u_botfr^2+v_botfr^2) (min value rn_emin) 
     230      !   en(bot)   = (ebb0/rho0)*0.5*sqrt(u_botfr^2+v_botfr^2) (min value rn_emin) 
    231231      ! where ebb0 does not includes surface wave enhancement (i.e. ebb0=3.75) 
    232232      ! Note that stress averaged is done using an wet-only calculation of u and v at t-point like in zdfsh2 
     
    237237            zmsku = ( 2. - umask(ji-1,jj,mbkt(ji,jj)) * umask(ji,jj,mbkt(ji,jj)) ) 
    238238            zmskv = ( 2. - vmask(ji,jj-1,mbkt(ji,jj)) * vmask(ji,jj,mbkt(ji,jj)) ) 
    239             !                       ! where 0.001875 = (rn_ebb0/rau0) * 0.5 = 3.75*0.5/1000. (CAUTION CdU<0) 
     239            !                       ! where 0.001875 = (rn_ebb0/rho0) * 0.5 = 3.75*0.5/1000. (CAUTION CdU<0) 
    240240            zebot = - 0.001875_wp * rCdU_bot(ji,jj) * SQRT(  ( zmsku*( uu(ji,jj,mbkt(ji,jj),Kbb)+uu(ji-1,jj,mbkt(ji,jj),Kbb) ) )**2  & 
    241241               &                                           + ( zmskv*( vv(ji,jj,mbkt(ji,jj),Kbb)+vv(ji,jj-1,mbkt(ji,jj),Kbb) ) )**2  ) 
     
    246246               zmsku = ( 2. - umask(ji-1,jj,mikt(ji,jj)) * umask(ji,jj,mikt(ji,jj)) ) 
    247247               zmskv = ( 2. - vmask(ji,jj-1,mikt(ji,jj)) * vmask(ji,jj,mikt(ji,jj)) ) 
    248                !                             ! where 0.001875 = (rn_ebb0/rau0) * 0.5 = 3.75*0.5/1000.  (CAUTION CdU<0) 
     248               !                             ! where 0.001875 = (rn_ebb0/rho0) * 0.5 = 3.75*0.5/1000.  (CAUTION CdU<0) 
    249249               zetop = - 0.001875_wp * rCdU_top(ji,jj) * SQRT(  ( zmsku*( uu(ji,jj,mikt(ji,jj),Kbb)+uu(ji-1,jj,mikt(ji,jj),Kbb) ) )**2  & 
    250250                  &                                           + ( zmskv*( vv(ji,jj,mikt(ji,jj),Kbb)+vv(ji,jj-1,mikt(ji,jj),Kbb) ) )**2  ) 
     
    439439      zmxld(:,:,:)  = rmxl_min 
    440440      ! 
    441       IF( ln_mxl0 ) THEN            ! surface mixing length = F(stress) : l=vkarmn*2.e5*taum/(rau0*g) 
    442          zraug = vkarmn * 2.e5_wp / ( rau0 * grav ) 
     441      IF( ln_mxl0 ) THEN            ! surface mixing length = F(stress) : l=vkarmn*2.e5*taum/(rho0*g) 
     442         zraug = vkarmn * 2.e5_wp / ( rho0 * grav ) 
    443443         DO_2D_00_00 
    444444            zmxlm(ji,jj,1) = MAX( rn_mxl0, zraug * taum(ji,jj) * tmask(ji,jj,1) ) 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OCE/oce.F90

    r12377 r12443  
    2828   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:)     ::   rn2b ,  rn2    !: brunt-vaisala frequency**2     [s-2] 
    2929   ! 
    30    REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) ::   rhd    !: in situ density anomalie rhd=(rho-rau0)/rau0  [no units] 
     30   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) ::   rhd    !: in situ density anomalie rhd=(rho-rho0)/rho0  [no units] 
    3131   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) ::   rhop   !: potential volumic mass                           [kg/m3] 
    3232   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) ::   Cu_adv                   !: vertical Courant number (adaptive-implicit) 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OFF/dtadyn.F90

    r12406 r12443  
    546546      ENDIF 
    547547 
    548       ssh(:,:,Kmm) = ssh(:,:,Kmm) + atfp * ( ssh(:,:,Kbb) - 2 * ssh(:,:,Kmm) + ssh(:,:,Kaa))   
     548      ssh(:,:,Kmm) = ssh(:,:,Kmm) + rn_atfp * ( ssh(:,:,Kbb) - 2 * ssh(:,:,Kmm) + ssh(:,:,Kaa))   
    549549 
    550550      !! Do we also need to time filter e3t?? 
     
    629629      END DO 
    630630      !                                                ! Sea surface  elevation time-stepping 
    631       pssha(:,:) = ( psshb(:,:) - z2dt * ( r1_rau0 * pemp(:,:)  + zhdiv(:,:) ) ) * ssmask(:,:) 
     631      pssha(:,:) = ( psshb(:,:) - z2dt * ( r1_rho0 * pemp(:,:)  + zhdiv(:,:) ) ) * ssmask(:,:) 
    632632      !                                                 !  
    633633      !                                                 ! After acale factors at t-points ( z_star coordinate ) 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/TOP/PISCES/P2Z/p2zexp.F90

    r12406 r12443  
    122122        DO_2D_11_11 
    123123           zsedpocd = zsedpoca(ji,jj) - 2. * sedpocn(ji,jj) + sedpocb(ji,jj)      ! time laplacian on tracers 
    124            sedpocb(ji,jj) = sedpocn(ji,jj) + atfp * zsedpocd                     ! sedpocb <-- filtered sedpocn 
     124           sedpocb(ji,jj) = sedpocn(ji,jj) + rn_atfp * zsedpocd                     ! sedpocb <-- filtered sedpocn 
    125125           sedpocn(ji,jj) = zsedpoca(ji,jj)                                       ! sedpocn <-- sedpoca 
    126126        END_2D 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/TOP/TRP/trcatf.F90

    r12397 r12443  
    7171      !!      the divergence of two consecutive time-steps and tr arrays 
    7272      !!      to prepare the next time_step: 
    73       !!         (tr(Kmm)) = (tr(Kmm)) + atfp [ (tr(Kbb)) + (tr(Kaa)) - 2 (tr(Kmm)) ] 
     73      !!         (tr(Kmm)) = (tr(Kmm)) + rn_atfp [ (tr(Kbb)) + (tr(Kaa)) - 2 (tr(Kmm)) ] 
    7474      !! 
    7575      !! 
     
    200200      !!                  /( e3t(:,:,:,Kmm)    + rbcp*[ e3t(:,:,:,Kbb)    - 2 e3t(:,:,:,Kmm)    + e3t(:,:,:,Kaa)    ] )    
    201201      !!             ztm = 0                                                       otherwise 
    202       !!             tb  = ( e3t_n*tn + atfp*[ e3t_b*tb - 2 e3t_n*tn + e3t_a*ta ] ) 
    203       !!                  /( e3t(:,:,:,Kmm)    + atfp*[ e3t(:,:,:,Kbb)    - 2 e3t(:,:,:,Kmm)    + e3t(:,:,:,Kaa)    ] ) 
     202      !!             tb  = ( e3t_n*tn + rn_atfp*[ e3t_b*tb - 2 e3t_n*tn + e3t_a*ta ] ) 
     203      !!                  /( e3t(:,:,:,Kmm)    + rn_atfp*[ e3t(:,:,:,Kbb)    - 2 e3t(:,:,:,Kmm)    + e3t(:,:,:,Kaa)    ] ) 
    204204      !!             tn  = ta  
    205205      !!             ta  = zt        (NB: reset to 0 after eos_bn2 call) 
     
    222222         IF(lwp) WRITE(numout,*) '~~~~~~~~~~~' 
    223223         IF( .NOT. ln_linssh ) THEN 
    224             rfact1 = atfp * rdttrc 
    225             rfact2 = rfact1 / rau0 
     224            rfact1 = rn_atfp * rdttrc 
     225            rfact2 = rfact1 / rho0 
    226226         ENDIF 
    227227        !   
     
    241241            ztc_d  = ztc_a  - 2. * ztc_n  + ztc_b 
    242242            ! 
    243             ze3t_f = ze3t_n + atfp * ze3t_d 
    244             ztc_f  = ztc_n  + atfp * ztc_d 
     243            ze3t_f = ze3t_n + rn_atfp * ze3t_d 
     244            ztc_f  = ztc_n  + rn_atfp * ztc_d 
    245245            ! 
    246246            IF( .NOT. ln_linssh .AND. jk == mikt(ji,jj) ) THEN           ! first level  
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/TOP/TRP/trcsbc.F90

    r12377 r12443  
    121121         DO jn = 1, jptra 
    122122            DO_2D_01_00 
    123                sbc_trc(ji,jj,jn) = zsfx(ji,jj) * r1_rau0 * ptr(ji,jj,1,jn,Kmm) 
     123               sbc_trc(ji,jj,jn) = zsfx(ji,jj) * r1_rho0 * ptr(ji,jj,1,jn,Kmm) 
    124124            END_2D 
    125125         END DO 
     
    129129         DO jn = 1, jptra 
    130130            DO_2D_01_00 
    131                sbc_trc(ji,jj,jn) = ( zsfx(ji,jj) + fmmflx(ji,jj) ) * r1_rau0 * ptr(ji,jj,1,jn,Kmm) 
     131               sbc_trc(ji,jj,jn) = ( zsfx(ji,jj) + fmmflx(ji,jj) ) * r1_rho0 * ptr(ji,jj,1,jn,Kmm) 
    132132            END_2D 
    133133         END DO 
     
    145145               ztfx  = zftra                        ! net tracer flux 
    146146               ! 
    147                zdtra = r1_rau0 * ( ztfx + ( zsfx(ji,jj) + fmmflx(ji,jj) ) * ptr(ji,jj,1,jn,Kmm) )  
     147               zdtra = r1_rho0 * ( ztfx + ( zsfx(ji,jj) + fmmflx(ji,jj) ) * ptr(ji,jj,1,jn,Kmm) )  
    148148               IF ( zdtra < 0. ) THEN 
    149149                  zdtra  = MAX(zdtra, -ptr(ji,jj,1,jn,Kmm) * e3t(ji,jj,1,Kmm) / r2dttrc )   ! avoid negative concentrations to arise 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/TOP/oce_trc.F90

    r12377 r12443  
    3939   USE oce , ONLY :   ts     =>    ts     !: 4D array contaning ( tn, sn ) 
    4040   USE oce , ONLY :   rhop   =>    rhop   !: potential volumic mass (kg m-3)  
    41    USE oce , ONLY :   rhd    =>    rhd    !: in situ density anomalie rhd=(rho-rau0)/rau0 (no units) 
     41   USE oce , ONLY :   rhd    =>    rhd    !: in situ density anomalie rhd=(rho-rho0)/rho0 (no units) 
    4242   USE oce , ONLY :   hdiv   =>    hdiv   !: horizontal divergence (1/s) 
    4343   USE oce , ONLY :   ssh    =>    ssh    !: sea surface height at t-point [m]    
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/TOP/trcbc.F90

    r12377 r12443  
    416416            DO_2D_01_00 
    417417               DO jk = 1, nk_rnf(ji,jj) 
    418                   zrnf = (rnf(ji,jj) + rnf_b(ji,jj)) * 0.5_wp * r1_rau0 / h_rnf(ji,jj) 
     418                  zrnf = (rnf(ji,jj) + rnf_b(ji,jj)) * 0.5_wp * r1_rho0 / h_rnf(ji,jj) 
    419419                  ptr(ji,jj,jk,jn,Krhs) = ptr(ji,jj,jk,jn,Krhs)  + (ptr(ji,jj,jk,jn,Kmm) * zrnf) 
    420420               END DO 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/tests/BENCH/EXPREF/namelist_cfg_orca025_like

    r12406 r12443  
    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_Dt seconds 
     191         nn_e      =  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

    r12406 r12443  
    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_Dt seconds 
     190         nn_e      =  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

    r12406 r12443  
    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_Dt seconds 
     190         nn_e      =  30         ! =F : the number of sub-step in rn_Dt seconds 
    191191/ 
    192192!----------------------------------------------------------------------- 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/tests/BENCH/MY_SRC/diawri.F90

    r12377 r12443  
    169169 
    170170      IF ( iom_use("taubot") ) THEN                ! bottom stress 
    171          zztmp = rau0 * 0.25 
     171         zztmp = rho0 * 0.25 
    172172         z2d(:,:) = 0._wp 
    173173         DO jj = 2, jpjm1 
     
    214214      IF( iom_use('w_masstr') .OR. iom_use('w_masstr2') ) THEN   ! vertical mass transport & its square value 
    215215         ! Caution: in the VVL case, it only correponds to the baroclinic mass transport. 
    216          z2d(:,:) = rau0 * e1e2t(:,:) 
     216         z2d(:,:) = rho0 * e1e2t(:,:) 
    217217         DO jk = 1, jpk 
    218218            z3d(:,:,jk) = ww(:,:,jk) * z2d(:,:) 
     
    257257            END DO 
    258258         END DO 
    259          CALL iom_put( "heatc", rau0_rcp * z2d )   ! vertically integrated heat content (J/m2) 
     259         CALL iom_put( "heatc", rho0_rcp * z2d )   ! vertically integrated heat content (J/m2) 
    260260      ENDIF 
    261261 
     
    269269            END DO 
    270270         END DO 
    271          CALL iom_put( "saltc", rau0 * z2d )          ! vertically integrated salt content (PSU*kg/m2) 
     271         CALL iom_put( "saltc", rho0 * z2d )          ! vertically integrated salt content (PSU*kg/m2) 
    272272      ENDIF 
    273273      ! 
     
    295295         z2d(:,:) = 0.e0 
    296296         DO jk = 1, jpkm1 
    297             z3d(:,:,jk) = rau0 * uu(:,:,jk,Kmm) * e2u(:,:) * e3u(:,:,jk,Kmm) * umask(:,:,jk) 
     297            z3d(:,:,jk) = rho0 * uu(:,:,jk,Kmm) * e2u(:,:) * e3u(:,:,jk,Kmm) * umask(:,:,jk) 
    298298            z2d(:,:) = z2d(:,:) + z3d(:,:,jk) 
    299299         END DO 
     
    332332         z3d(:,:,jpk) = 0.e0 
    333333         DO jk = 1, jpkm1 
    334             z3d(:,:,jk) = rau0 * vv(:,:,jk,Kmm) * e1v(:,:) * e3v(:,:,jk,Kmm) * vmask(:,:,jk) 
     334            z3d(:,:,jk) = rho0 * vv(:,:,jk,Kmm) * e1v(:,:) * e3v(:,:,jk,Kmm) * vmask(:,:,jk) 
    335335         END DO 
    336336         CALL iom_put( "v_masstr", z3d )              ! mass transport in j-direction 
     
    373373         END DO 
    374374         CALL lbc_lnk( 'diawri', z2d, 'T', -1. ) 
    375          CALL iom_put( "tosmint", rau0 * z2d )        ! Vertical integral of temperature 
     375         CALL iom_put( "tosmint", rho0 * z2d )        ! Vertical integral of temperature 
    376376      ENDIF 
    377377      IF( iom_use("somint") ) THEN 
     
    385385         END DO 
    386386         CALL lbc_lnk( 'diawri', z2d, 'T', -1. ) 
    387          CALL iom_put( "somint", rau0 * z2d )         ! Vertical integral of salinity 
     387         CALL iom_put( "somint", rho0 * z2d )         ! Vertical integral of salinity 
    388388      ENDIF 
    389389 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/tests/BENCH/MY_SRC/zdfiwm.F90

    r12377 r12443  
    8585      !!              This is divided into three components: 
    8686      !!                 1. Bottom-intensified low-mode dissipation at critical slopes 
    87       !!                     zemx_iwm(z) = ( ecri_iwm / rau0 ) * EXP( -(H-z)/hcri_iwm ) 
     87      !!                     zemx_iwm(z) = ( ecri_iwm / rho0 ) * EXP( -(H-z)/hcri_iwm ) 
    8888      !!                                   / ( 1. - EXP( - H/hcri_iwm ) ) * hcri_iwm 
    8989      !!              where hcri_iwm is the characteristic length scale of the bottom  
    9090      !!              intensification, ecri_iwm a map of available power, and H the ocean depth. 
    9191      !!                 2. Pycnocline-intensified low-mode dissipation 
    92       !!                     zemx_iwm(z) = ( epyc_iwm / rau0 ) * ( sqrt(rn2(z))^nn_zpyc ) 
     92      !!                     zemx_iwm(z) = ( epyc_iwm / rho0 ) * ( sqrt(rn2(z))^nn_zpyc ) 
    9393      !!                                   / SUM( sqrt(rn2(z))^nn_zpyc * e3w(z) ) 
    9494      !!              where epyc_iwm is a map of available power, and nn_zpyc 
     
    9696      !!              energy dissipation. 
    9797      !!                 3. WKB-height dependent high mode dissipation 
    98       !!                     zemx_iwm(z) = ( ebot_iwm / rau0 ) * rn2(z) * EXP(-z_wkb(z)/hbot_iwm) 
     98      !!                     zemx_iwm(z) = ( ebot_iwm / rho0 ) * rn2(z) * EXP(-z_wkb(z)/hbot_iwm) 
    9999      !!                                   / SUM( rn2(z) * EXP(-z_wkb(z)/hbot_iwm) * e3w(z) ) 
    100100      !!              where hbot_iwm is the characteristic length scale of the WKB bottom  
     
    150150         DO ji = 1, jpi 
    151151            zhdep(ji,jj) = gdepw_0(ji,jj,mbkt(ji,jj)+1)       ! depth of the ocean 
    152             zfact(ji,jj) = rau0 * (  1._wp - EXP( -zhdep(ji,jj) / hcri_iwm(ji,jj) )  ) 
     152            zfact(ji,jj) = rho0 * (  1._wp - EXP( -zhdep(ji,jj) / hcri_iwm(ji,jj) )  ) 
    153153            IF( zfact(ji,jj) /= 0._wp )   zfact(ji,jj) = ecri_iwm(ji,jj) / zfact(ji,jj) 
    154154         END DO 
     
    179179         DO jj = 1, jpj 
    180180            DO ji = 1, jpi 
    181                IF( zfact(ji,jj) /= 0 )   zfact(ji,jj) = epyc_iwm(ji,jj) / ( rau0 * zfact(ji,jj) ) 
     181               IF( zfact(ji,jj) /= 0 )   zfact(ji,jj) = epyc_iwm(ji,jj) / ( rho0 * zfact(ji,jj) ) 
    182182            END DO 
    183183         END DO 
     
    196196         DO jj= 1, jpj 
    197197            DO ji = 1, jpi 
    198                IF( zfact(ji,jj) /= 0 )   zfact(ji,jj) = epyc_iwm(ji,jj) / ( rau0 * zfact(ji,jj) ) 
     198               IF( zfact(ji,jj) /= 0 )   zfact(ji,jj) = epyc_iwm(ji,jj) / ( rho0 * zfact(ji,jj) ) 
    199199            END DO 
    200200         END DO 
     
    246246      DO jj = 1, jpj 
    247247         DO ji = 1, jpi 
    248             IF( zfact(ji,jj) /= 0 )   zfact(ji,jj) = ebot_iwm(ji,jj) / ( rau0 * zfact(ji,jj) ) 
     248            IF( zfact(ji,jj) /= 0 )   zfact(ji,jj) = ebot_iwm(ji,jj) / ( rho0 * zfact(ji,jj) ) 
    249249         END DO 
    250250      END DO 
     
    259259      ! Calculate molecular kinematic viscosity 
    260260      znu_t(:,:,:) = 1.e-4_wp * (  17.91_wp - 0.53810_wp * ts(:,:,:,jp_tem,Kmm) + 0.00694_wp * ts(:,:,:,jp_tem,Kmm) * ts(:,:,:,jp_tem,Kmm)  & 
    261          &                                  + 0.02305_wp * ts(:,:,:,jp_sal,Kmm)  ) * tmask(:,:,:) * r1_rau0 
     261         &                                  + 0.02305_wp * ts(:,:,:,jp_sal,Kmm)  ) * tmask(:,:,:) * r1_rho0 
    262262      DO jk = 2, jpkm1 
    263263         znu_w(:,:,jk) = 0.5_wp * ( znu_t(:,:,jk-1) + znu_t(:,:,jk) ) * wmask(:,:,jk) 
     
    305305         END DO 
    306306         CALL mpp_sum( 'zdfiwm', zztmp ) 
    307          zztmp = rau0 * zztmp ! Global integral of rauo * Kz * N^2 = power contributing to mixing  
     307         zztmp = rho0 * zztmp ! Global integral of rauo * Kz * N^2 = power contributing to mixing  
    308308         ! 
    309309         IF(lwp) THEN 
     
    349349                                    !* output useful diagnostics: Kz*N^2 ,  
    350350!!gm Kz*N2 should take into account the ratio avs/avt if it is used.... (see diaar5) 
    351                                     !  vertical integral of rau0 * Kz * N^2 , energy density (zemx_iwm) 
     351                                    !  vertical integral of rho0 * Kz * N^2 , energy density (zemx_iwm) 
    352352      IF( iom_use("bflx_iwm") .OR. iom_use("pcmap_iwm") ) THEN 
    353353         ALLOCATE( z2d(jpi,jpj) , z3d(jpi,jpj,jpk) ) 
     
    357357            z2d(:,:) = z2d(:,:) + e3w(:,:,jk,Kmm) * z3d(:,:,jk) * wmask(:,:,jk) 
    358358         END DO 
    359          z2d(:,:) = rau0 * z2d(:,:) 
     359         z2d(:,:) = rho0 * z2d(:,:) 
    360360         CALL iom_put( "bflx_iwm", z3d ) 
    361361         CALL iom_put( "pcmap_iwm", z2d ) 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/tests/CANAL/EXPREF/namelist_cfg

    r12406 r12443  
    208208      ln_bt_av      = .true.     ! Time filtering of barotropic variables 
    209209         nn_bt_flt     = 1          ! Time filter choice  = 0 None 
    210          !                          !                     = 1 Boxcar over   nn_baro sub-steps 
    211          !                          !                     = 2 Boxcar over 2*nn_baro  "    " 
     210         !                          !                     = 1 Boxcar over   nn_e sub-steps 
     211         !                          !                     = 2 Boxcar over 2*nn_e  "    " 
    212212      ln_bt_auto    = .false.    ! Number of sub-step defined from: 
    213          nn_baro      =  24         ! =F : the number of sub-step in rn_Dt seconds 
     213         nn_e      =  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

    r12406 r12443  
    169169 
    170170      IF ( iom_use("taubot") ) THEN                ! bottom stress 
    171          zztmp = rau0 * 0.25 
     171         zztmp = rho0 * 0.25 
    172172         z2d(:,:) = 0._wp 
    173173         DO jj = 2, jpjm1 
     
    212212      IF( iom_use('w_masstr') .OR. iom_use('w_masstr2') ) THEN   ! vertical mass transport & its square value 
    213213         ! Caution: in the VVL case, it only correponds to the baroclinic mass transport. 
    214          z2d(:,:) = rau0 * e1e2t(:,:) 
     214         z2d(:,:) = rho0 * e1e2t(:,:) 
    215215         DO jk = 1, jpk 
    216216            z3d(:,:,jk) = ww(:,:,jk) * z2d(:,:) 
     
    272272            END DO 
    273273         END DO 
    274          CALL iom_put( "heatc", rau0_rcp * z2d )   ! vertically integrated heat content (J/m2) 
     274         CALL iom_put( "heatc", rho0_rcp * z2d )   ! vertically integrated heat content (J/m2) 
    275275      ENDIF 
    276276 
     
    284284            END DO 
    285285         END DO 
    286          CALL iom_put( "saltc", rau0 * z2d )          ! vertically integrated salt content (PSU*kg/m2) 
     286         CALL iom_put( "saltc", rho0 * z2d )          ! vertically integrated salt content (PSU*kg/m2) 
    287287      ENDIF 
    288288      ! 
     
    296296            END DO 
    297297         END DO 
    298          CALL iom_put( "salt2c", rau0 * z2d )          ! vertically integrated salt content (PSU*kg/m2) 
     298         CALL iom_put( "salt2c", rho0 * z2d )          ! vertically integrated salt content (PSU*kg/m2) 
    299299      ENDIF 
    300300      ! 
     
    395395         z2d(:,:) = 0.e0 
    396396         DO jk = 1, jpkm1 
    397             z3d(:,:,jk) = rau0 * uu(:,:,jk,Kmm) * e2u(:,:) * e3u(:,:,jk,Kmm) * umask(:,:,jk) 
     397            z3d(:,:,jk) = rho0 * uu(:,:,jk,Kmm) * e2u(:,:) * e3u(:,:,jk,Kmm) * umask(:,:,jk) 
    398398            z2d(:,:) = z2d(:,:) + z3d(:,:,jk) 
    399399         END DO 
     
    432432         z3d(:,:,jpk) = 0.e0 
    433433         DO jk = 1, jpkm1 
    434             z3d(:,:,jk) = rau0 * vv(:,:,jk,Kmm) * e1v(:,:) * e3v(:,:,jk,Kmm) * vmask(:,:,jk) 
     434            z3d(:,:,jk) = rho0 * vv(:,:,jk,Kmm) * e1v(:,:) * e3v(:,:,jk,Kmm) * vmask(:,:,jk) 
    435435         END DO 
    436436         CALL iom_put( "v_masstr", z3d )              ! mass transport in j-direction 
     
    473473         END DO 
    474474         CALL lbc_lnk( 'diawri', z2d, 'T', -1. ) 
    475          CALL iom_put( "tosmint", rau0 * z2d )        ! Vertical integral of temperature 
     475         CALL iom_put( "tosmint", rho0 * z2d )        ! Vertical integral of temperature 
    476476      ENDIF 
    477477      IF( iom_use("somint") ) THEN 
     
    485485         END DO 
    486486         CALL lbc_lnk( 'diawri', z2d, 'T', -1. ) 
    487          CALL iom_put( "somint", rau0 * z2d )         ! Vertical integral of salinity 
     487         CALL iom_put( "somint", rho0 * z2d )         ! Vertical integral of salinity 
    488488      ENDIF 
    489489 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/tests/CANAL/MY_SRC/domvvl.F90

    r12424 r12443  
    609609         ELSE 
    610610            tilde_e3t_b(:,:,:) = tilde_e3t_n(:,:,:) &  
    611             &         + atfp * ( tilde_e3t_b(:,:,:) - 2.0_wp * tilde_e3t_n(:,:,:) + tilde_e3t_a(:,:,:) ) 
     611            &         + rn_atfp * ( tilde_e3t_b(:,:,:) - 2.0_wp * tilde_e3t_n(:,:,:) + tilde_e3t_a(:,:,:) ) 
    612612         ENDIF 
    613613         tilde_e3t_n(:,:,:) = tilde_e3t_a(:,:,:) 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/tests/CANAL/MY_SRC/usrdef_istate.F90

    r10425 r12443  
    218218         ! 
    219219         zr_lambda2 = 1._wp / zlambda**2 
    220          zP0 = rau0 * zf0 * zumax * zlambda * SQRT(EXP(1._wp)/2._wp) 
     220         zP0 = rho0 * zf0 * zumax * zlambda * SQRT(EXP(1._wp)/2._wp) 
    221221         ! 
    222222         DO jj=1, jpj 
     
    225225               zy = gphit(ji,jj) * 1.e3 
    226226               ! Surface pressure: P(x,y,z) = F(z) * Psurf(x,y) 
    227                zpsurf = zP0 * EXP(-(zx**2+zy**2)*zr_lambda2) - rau0 * ff_t(ji,jj) * rn_uzonal * zy 
     227               zpsurf = zP0 * EXP(-(zx**2+zy**2)*zr_lambda2) - rho0 * ff_t(ji,jj) * rn_uzonal * zy 
    228228               ! Sea level: 
    229229               pssh(ji,jj) = 0. 
     
    231231                  zdt = pssh(ji,jj) 
    232232                  zdzF = (1._wp - EXP(zdt-zH)) / (zH - 1._wp + EXP(-zH))   ! F'(z) 
    233                   zrho1 = rau0 * (1._wp + zn2*zdt/grav) - zdzF * zpsurf / grav    ! -1/g Dz(P) = -1/g * F'(z) * Psurf(x,y) 
     233                  zrho1 = rho0 * (1._wp + zn2*zdt/grav) - zdzF * zpsurf / grav    ! -1/g Dz(P) = -1/g * F'(z) * Psurf(x,y) 
    234234                  pssh(ji,jj) = zpsurf / (zrho1*grav) * ptmask(ji,jj,1)   ! ssh = Psurf / (Rho*g) 
    235235               END DO 
     
    237237               DO jk=1,jpk 
    238238                  zdt =  pdept(ji,jj,jk)  
    239                   zrho1 = rau0 * (1._wp + zn2*zdt/grav) 
     239                  zrho1 = rho0 * (1._wp + zn2*zdt/grav) 
    240240                  IF (zdt < zH) THEN 
    241241                     zdzF = (1._wp-EXP(zdt-zH)) / (zH-1._wp + EXP(-zH))   ! F'(z) 
    242242                     zrho1 = zrho1 - zdzF * zpsurf / grav    ! -1/g Dz(P) = -1/g * F'(z) * Psurf(x,y) 
    243243                  ENDIF 
    244                   !               pts(ji,jj,jk,jp_tem) = (20._wp + (rau0-zrho1) / 0.28_wp) * ptmask(ji,jj,jk) 
    245                   pts(ji,jj,jk,jp_tem) = (10._wp + (rau0-zrho1) / 0.28_wp) * ptmask(ji,jj,jk) 
     244                  !               pts(ji,jj,jk,jp_tem) = (20._wp + (rho0-zrho1) / 0.28_wp) * ptmask(ji,jj,jk) 
     245                  pts(ji,jj,jk,jp_tem) = (10._wp + (rho0-zrho1) / 0.28_wp) * ptmask(ji,jj,jk) 
    246246               END DO 
    247247            END DO 
     
    261261                  IF (zdu < zH) THEN 
    262262                     zf = (zH-1._wp-zdu+EXP(zdu-zH)) / (zH-1._wp+EXP(-zH)) 
    263                      zdyPs = - za * zy * EXP(-(zx**2+zy**2)*zr_lambda2) - rau0 * ff_t(ji,jj) * rn_uzonal 
    264                      pu(ji,jj,jk) = - zf / ( rau0 * ff_t(ji,jj) ) * zdyPs * ptmask(ji,jj,jk) * ptmask(ji+1,jj,jk) 
     263                     zdyPs = - za * zy * EXP(-(zx**2+zy**2)*zr_lambda2) - rho0 * ff_t(ji,jj) * rn_uzonal 
     264                     pu(ji,jj,jk) = - zf / ( rho0 * ff_t(ji,jj) ) * zdyPs * ptmask(ji,jj,jk) * ptmask(ji+1,jj,jk) 
    265265                  ELSE 
    266266                     pu(ji,jj,jk) = 0._wp 
     
    279279                     zf = (zH-1._wp-zdv+EXP(zdv-zH)) / (zH-1._wp+EXP(-zH)) 
    280280                     zdxPs = - za * zx * EXP(-(zx**2+zy**2)*zr_lambda2) 
    281                      pv(ji,jj,jk) = zf / ( rau0 * ff_f(ji,jj) ) * zdxPs * ptmask(ji,jj,jk) * ptmask(ji,jj+1,jk) 
     281                     pv(ji,jj,jk) = zf / ( rho0 * ff_f(ji,jj) ) * zdxPs * ptmask(ji,jj,jk) * ptmask(ji,jj+1,jk) 
    282282                  ELSE 
    283283                     pv(ji,jj,jk) = 0._wp 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/tests/ISOMIP+/MY_SRC/eosbn2.F90

    r12353 r12443  
    191191      !!                   ***  ROUTINE eos_insitu  *** 
    192192      !! 
    193       !! ** Purpose :   Compute the in situ density (ratio rho/rau0) from 
     193      !! ** Purpose :   Compute the in situ density (ratio rho/rho0) from 
    194194      !!       potential temperature and salinity using an equation of state 
    195195      !!       selected in the nameos namelist 
    196196      !! 
    197       !! ** Method  :   prd(t,s,z) = ( rho(t,s,z) - rau0 ) / rau0 
     197      !! ** Method  :   prd(t,s,z) = ( rho(t,s,z) - rho0 ) / rho0 
    198198      !!         with   prd    in situ density anomaly      no units 
    199199      !!                t      TEOS10: CT or EOS80: PT      Celsius 
     
    201201      !!                z      depth                        meters 
    202202      !!                rho    in situ density              kg/m^3 
    203       !!                rau0   reference density            kg/m^3 
     203      !!                rho0   reference density            kg/m^3 
    204204      !! 
    205205      !!     ln_teos10 : polynomial TEOS-10 equation of state is used for rho(t,s,z). 
     
    210210      !! 
    211211      !!     ln_seos : simplified equation of state 
    212       !!              prd(t,s,z) = ( -a0*(1+lambda/2*(T-T0)+mu*z+nu*(S-S0))*(T-T0) + b0*(S-S0) ) / rau0 
     212      !!              prd(t,s,z) = ( -a0*(1+lambda/2*(T-T0)+mu*z+nu*(S-S0))*(T-T0) + b0*(S-S0) ) / rho0 
    213213      !!              linear case function of T only: rn_alpha<>0, other coefficients = 0 
    214214      !!              linear eos function of T and S: rn_alpha and rn_beta<>0, other coefficients=0 
     
    216216      !! 
    217217      !!     ln_leos : linear ISOMIP equation of state 
    218       !!              prd(t,s,z) = ( -a0*(T-T0) + b0*(S-S0) ) / rau0 
     218      !!              prd(t,s,z) = ( -a0*(T-T0) + b0*(S-S0) ) / rho0 
    219219      !!              setup for ISOMIP linear eos 
    220220      !! 
     
    273273                  zn  = ( ( zn3 * zh + zn2 ) * zh + zn1 ) * zh + zn0 
    274274                  ! 
    275                   prd(ji,jj,jk) = (  zn * r1_rau0 - 1._wp  ) * ztm  ! density anomaly (masked) 
     275                  prd(ji,jj,jk) = (  zn * r1_rho0 - 1._wp  ) * ztm  ! density anomaly (masked) 
    276276                  ! 
    277277               END DO 
     
    293293                     &  - rn_nu * zt * zs 
    294294                     !                                  
    295                   prd(ji,jj,jk) = zn * r1_rau0 * ztm                ! density anomaly (masked) 
     295                  prd(ji,jj,jk) = zn * r1_rho0 * ztm                ! density anomaly (masked) 
    296296               END DO 
    297297            END DO 
     
    308308                  ztm = tmask(ji,jj,jk) 
    309309                  ! 
    310                   zn =  rau0 * ( - rn_a0 * zt + rn_b0 * zs ) 
     310                  zn =  rho0 * ( - rn_a0 * zt + rn_b0 * zs ) 
    311311                  !                                  
    312                   prd(ji,jj,jk) = zn * r1_rau0 * ztm                ! density anomaly (masked) 
     312                  prd(ji,jj,jk) = zn * r1_rho0 * ztm                ! density anomaly (masked) 
    313313               END DO 
    314314            END DO 
     
    328328      !!                  ***  ROUTINE eos_insitu_pot  *** 
    329329      !! 
    330       !! ** Purpose :   Compute the in situ density (ratio rho/rau0) and the 
     330      !! ** Purpose :   Compute the in situ density (ratio rho/rho0) and the 
    331331      !!      potential volumic mass (Kg/m3) from potential temperature and 
    332332      !!      salinity fields using an equation of state selected in the 
     
    410410                        prhop(ji,jj,jk) = prhop(ji,jj,jk) + zn0_sto(jsmp)                      ! potential density referenced at the surface 
    411411                        ! 
    412                         prd(ji,jj,jk) = prd(ji,jj,jk) + (  zn_sto(jsmp) * r1_rau0 - 1._wp  )   ! density anomaly (masked) 
     412                        prd(ji,jj,jk) = prd(ji,jj,jk) + (  zn_sto(jsmp) * r1_rho0 - 1._wp  )   ! density anomaly (masked) 
    413413                     END DO 
    414414                     prhop(ji,jj,jk) = 0.5_wp * prhop(ji,jj,jk) * ztm / nn_sto_eos 
     
    454454                     prhop(ji,jj,jk) = zn0 * ztm                           ! potential density referenced at the surface 
    455455                     ! 
    456                      prd(ji,jj,jk) = (  zn * r1_rau0 - 1._wp  ) * ztm      ! density anomaly (masked) 
     456                     prd(ji,jj,jk) = (  zn * r1_rho0 - 1._wp  ) * ztm      ! density anomaly (masked) 
    457457                  END DO 
    458458               END DO 
     
    473473                     &  + rn_b0 * ( 1._wp - 0.5_wp*rn_lambda2*zs ) * zs   & 
    474474                     &  - rn_nu * zt * zs 
    475                   prhop(ji,jj,jk) = ( rau0 + zn ) * ztm 
     475                  prhop(ji,jj,jk) = ( rho0 + zn ) * ztm 
    476476                  !                                                     ! density anomaly (masked) 
    477477                  zn = zn - ( rn_a0 * rn_mu1 * zt + rn_b0 * rn_mu2 * zs ) * zh 
    478                   prd(ji,jj,jk) = zn * r1_rau0 * ztm 
     478                  prd(ji,jj,jk) = zn * r1_rho0 * ztm 
    479479                  ! 
    480480               END DO 
     
    492492                  ztm = tmask(ji,jj,jk) 
    493493                  !                                                     ! potential density referenced at the surface 
    494                   zn =  rau0 * ( - rn_a0 * zt + rn_b0 * zs ) 
    495                   prhop(ji,jj,jk) = ( rau0 + zn ) * ztm 
     494                  zn =  rho0 * ( - rn_a0 * zt + rn_b0 * zs ) 
     495                  prhop(ji,jj,jk) = ( rho0 + zn ) * ztm 
    496496                  !                                                     ! density anomaly (masked) 
    497                   prd(ji,jj,jk) = zn * r1_rau0 * ztm 
     497                  prd(ji,jj,jk) = zn * r1_rho0 * ztm 
    498498                  ! 
    499499               END DO 
     
    514514      !!                  ***  ROUTINE eos_insitu_2d  *** 
    515515      !! 
    516       !! ** Purpose :   Compute the in situ density (ratio rho/rau0) from 
     516      !! ** Purpose :   Compute the in situ density (ratio rho/rho0) from 
    517517      !!      potential temperature and salinity using an equation of state 
    518518      !!      selected in the nameos namelist. * 2D field case 
     
    569569               zn  = ( ( zn3 * zh + zn2 ) * zh + zn1 ) * zh + zn0 
    570570               ! 
    571                prd(ji,jj) = zn * r1_rau0 - 1._wp               ! unmasked in situ density anomaly 
     571               prd(ji,jj) = zn * r1_rho0 - 1._wp               ! unmasked in situ density anomaly 
    572572               ! 
    573573            END DO 
     
    589589                  &  - rn_nu * zt * zs 
    590590                  ! 
    591                prd(ji,jj) = zn * r1_rau0               ! unmasked in situ density anomaly 
     591               prd(ji,jj) = zn * r1_rho0               ! unmasked in situ density anomaly 
    592592               ! 
    593593            END DO 
     
    605605               zh    = pdep (ji,jj)                         ! depth at the partial step level 
    606606               ! 
    607                zn =  rau0 * ( - rn_a0 * zt + rn_b0 * zs ) 
    608                   ! 
    609                prd(ji,jj) = zn * r1_rau0               ! unmasked in situ density anomaly 
     607               zn =  rho0 * ( - rn_a0 * zt + rn_b0 * zs ) 
     608                  ! 
     609               prd(ji,jj) = zn * r1_rho0               ! unmasked in situ density anomaly 
    610610               ! 
    611611            END DO 
     
    676676                  zn  = ( ( zn3 * zh + zn2 ) * zh + zn1 ) * zh + zn0 
    677677                  ! 
    678                   pab(ji,jj,jk,jp_tem) = zn * r1_rau0 * ztm 
     678                  pab(ji,jj,jk,jp_tem) = zn * r1_rho0 * ztm 
    679679                  ! 
    680680                  ! beta 
     
    697697                  zn  = ( ( zn3 * zh + zn2 ) * zh + zn1 ) * zh + zn0 
    698698                  ! 
    699                   pab(ji,jj,jk,jp_sal) = zn / zs * r1_rau0 * ztm 
     699                  pab(ji,jj,jk,jp_sal) = zn / zs * r1_rho0 * ztm 
    700700                  ! 
    701701               END DO 
     
    714714                  ! 
    715715                  zn  = rn_a0 * ( 1._wp + rn_lambda1*zt + rn_mu1*zh ) + rn_nu*zs 
    716                   pab(ji,jj,jk,jp_tem) = zn * r1_rau0 * ztm   ! alpha 
     716                  pab(ji,jj,jk,jp_tem) = zn * r1_rho0 * ztm   ! alpha 
    717717                  ! 
    718718                  zn  = rn_b0 * ( 1._wp - rn_lambda2*zs - rn_mu2*zh ) - rn_nu*zt 
    719                   pab(ji,jj,jk,jp_sal) = zn * r1_rau0 * ztm   ! beta 
     719                  pab(ji,jj,jk,jp_sal) = zn * r1_rho0 * ztm   ! beta 
    720720                  ! 
    721721               END DO 
     
    733733                  ztm = tmask(ji,jj,jk)                   ! land/sea bottom mask = surf. mask 
    734734                  ! 
    735                   zn  = rn_a0 * rau0 
    736                   pab(ji,jj,jk,jp_tem) = zn * r1_rau0 * ztm   ! alpha 
    737                   ! 
    738                   zn  = rn_b0 * rau0 
    739                   pab(ji,jj,jk,jp_sal) = zn * r1_rau0 * ztm   ! beta 
     735                  zn  = rn_a0 * rho0 
     736                  pab(ji,jj,jk,jp_tem) = zn * r1_rho0 * ztm   ! alpha 
     737                  ! 
     738                  zn  = rn_b0 * rho0 
     739                  pab(ji,jj,jk,jp_sal) = zn * r1_rho0 * ztm   ! beta 
    740740                  ! 
    741741               END DO 
     
    809809               zn  = ( ( zn3 * zh + zn2 ) * zh + zn1 ) * zh + zn0 
    810810               ! 
    811                pab(ji,jj,jp_tem) = zn * r1_rau0 
     811               pab(ji,jj,jp_tem) = zn * r1_rho0 
    812812               ! 
    813813               ! beta 
     
    830830               zn  = ( ( zn3 * zh + zn2 ) * zh + zn1 ) * zh + zn0 
    831831               ! 
    832                pab(ji,jj,jp_sal) = zn / zs * r1_rau0 
     832               pab(ji,jj,jp_sal) = zn / zs * r1_rho0 
    833833               ! 
    834834               ! 
     
    848848               ! 
    849849               zn  = rn_a0 * ( 1._wp + rn_lambda1*zt + rn_mu1*zh ) + rn_nu*zs 
    850                pab(ji,jj,jp_tem) = zn * r1_rau0   ! alpha 
     850               pab(ji,jj,jp_tem) = zn * r1_rho0   ! alpha 
    851851               ! 
    852852               zn  = rn_b0 * ( 1._wp - rn_lambda2*zs - rn_mu2*zh ) - rn_nu*zt 
    853                pab(ji,jj,jp_sal) = zn * r1_rau0   ! beta 
     853               pab(ji,jj,jp_sal) = zn * r1_rho0   ! beta 
    854854               ! 
    855855            END DO 
     
    867867               zh    = pdep (ji,jj)                   ! depth at the partial step level 
    868868               ! 
    869                zn  = rn_a0 * rau0 
    870                pab(ji,jj,jp_tem) = zn * r1_rau0   ! alpha 
    871                ! 
    872                zn  = rn_b0 * rau0 
    873                pab(ji,jj,jp_sal) = zn * r1_rau0   ! beta 
     869               zn  = rn_a0 * rho0 
     870               pab(ji,jj,jp_tem) = zn * r1_rho0   ! alpha 
     871               ! 
     872               zn  = rn_b0 * rho0 
     873               pab(ji,jj,jp_sal) = zn * r1_rho0   ! beta 
    874874               ! 
    875875            END DO 
     
    941941         zn  = ( ( zn3 * zh + zn2 ) * zh + zn1 ) * zh + zn0 
    942942         ! 
    943          pab(jp_tem) = zn * r1_rau0 
     943         pab(jp_tem) = zn * r1_rho0 
    944944         ! 
    945945         ! beta 
     
    962962         zn  = ( ( zn3 * zh + zn2 ) * zh + zn1 ) * zh + zn0 
    963963         ! 
    964          pab(jp_sal) = zn / zs * r1_rau0 
     964         pab(jp_sal) = zn / zs * r1_rho0 
    965965         ! 
    966966         ! 
     
    973973         ! 
    974974         zn  = rn_a0 * ( 1._wp + rn_lambda1*zt + rn_mu1*zh ) + rn_nu*zs 
    975          pab(jp_tem) = zn * r1_rau0   ! alpha 
     975         pab(jp_tem) = zn * r1_rho0   ! alpha 
    976976         ! 
    977977         zn  = rn_b0 * ( 1._wp - rn_lambda2*zs - rn_mu2*zh ) - rn_nu*zt 
    978          pab(jp_sal) = zn * r1_rau0   ! beta 
     978         pab(jp_sal) = zn * r1_rho0   ! beta 
    979979         ! 
    980980      CASE( np_leos )                  !==  linear ISOMIP EOS  ==! 
     
    984984         zh    = pdep                    ! depth at the partial step level 
    985985         ! 
    986          zn  = rn_a0 * rau0 
    987          pab(jp_tem) = zn * r1_rau0   ! alpha 
    988          ! 
    989          zn  = rn_b0 * rau0 
    990          pab(jp_sal) = zn * r1_rau0   ! beta 
     986         zn  = rn_a0 * rho0 
     987         pab(jp_tem) = zn * r1_rho0   ! alpha 
     988         ! 
     989         zn  = rn_b0 * rho0 
     990         pab(jp_sal) = zn * r1_rho0   ! beta 
    991991         ! 
    992992      CASE DEFAULT 
     
    12141214      !! ** Method  :   PE is defined analytically as the vertical  
    12151215      !!                   primitive of EOS times -g integrated between 0 and z>0. 
    1216       !!                pen is the nonlinear bsq-PE anomaly: pen = ( PE - rau0 gz ) / rau0 gz - rd 
     1216      !!                pen is the nonlinear bsq-PE anomaly: pen = ( PE - rho0 gz ) / rho0 gz - rd 
    12171217      !!                                                      = 1/z * /int_0^z rd dz - rd  
    12181218      !!                                where rd is the density anomaly (see eos_rhd function) 
    12191219      !!                ab_pe are partial derivatives of PE anomaly with respect to T and S: 
    1220       !!                    ab_pe(1) = - 1/(rau0 gz) * dPE/dT + drd/dT = - d(pen)/dT 
    1221       !!                    ab_pe(2) =   1/(rau0 gz) * dPE/dS + drd/dS =   d(pen)/dS 
     1220      !!                    ab_pe(1) = - 1/(rho0 gz) * dPE/dT + drd/dT = - d(pen)/dT 
     1221      !!                    ab_pe(2) =   1/(rho0 gz) * dPE/dS + drd/dS =   d(pen)/dS 
    12221222      !! 
    12231223      !! ** Action  : - pen         : PE anomaly given at T-points 
     
    12671267                  zn  = ( zn2 * zh + zn1 ) * zh + zn0 
    12681268                  ! 
    1269                   ppen(ji,jj,jk)  = zn * zh * r1_rau0 * ztm 
     1269                  ppen(ji,jj,jk)  = zn * zh * r1_rho0 * ztm 
    12701270                  ! 
    12711271                  ! alphaPE non-linear anomaly 
     
    12821282                  zn  = ( zn2 * zh + zn1 ) * zh + zn0 
    12831283                  !                               
    1284                   pab_pe(ji,jj,jk,jp_tem) = zn * zh * r1_rau0 * ztm 
     1284                  pab_pe(ji,jj,jk,jp_tem) = zn * zh * r1_rho0 * ztm 
    12851285                  ! 
    12861286                  ! betaPE non-linear anomaly 
     
    12971297                  zn  = ( zn2 * zh + zn1 ) * zh + zn0 
    12981298                  !                               
    1299                   pab_pe(ji,jj,jk,jp_sal) = zn / zs * zh * r1_rau0 * ztm 
     1299                  pab_pe(ji,jj,jk,jp_sal) = zn / zs * zh * r1_rho0 * ztm 
    13001300                  ! 
    13011301               END DO 
     
    13121312                  zh  = gdept(ji,jj,jk,Kmm)              ! depth in meters  at t-point 
    13131313                  ztm = tmask(ji,jj,jk)                ! tmask 
    1314                   zn  = 0.5_wp * zh * r1_rau0 * ztm 
     1314                  zn  = 0.5_wp * zh * r1_rho0 * ztm 
    13151315                  !                                    ! Potential Energy 
    13161316                  ppen(ji,jj,jk) = ( rn_a0 * rn_mu1 * zt + rn_b0 * rn_mu2 * zs ) * zn 
     
    13321332                  zh  = gdept(ji,jj,jk,Kmm)                ! depth in meters  at t-point 
    13331333                  ztm = tmask(ji,jj,jk)                  ! tmask 
    1334                   zn  = 0.5_wp * zh * r1_rau0 * ztm 
     1334                  zn  = 0.5_wp * zh * r1_rho0 * ztm 
    13351335                  !                                    ! Potential Energy 
    13361336                  ppen(ji,jj,jk) = 0. 
     
    13791379      IF(lwm) WRITE( numond, nameos ) 
    13801380      ! 
    1381       rau0        = 1027.51_wp     !: volumic mass of reference     [kg/m3] 
     1381      rho0        = 1027.51_wp     !: volumic mass of reference     [kg/m3] 
    13821382      rcp         = 3974.00_wp     !: heat capacity     [J/K] 
    13831383      ! 
     
    17931793            WRITE(numout,*) '   ==>>>   use of simplified eos:    ' 
    17941794            WRITE(numout,*) '              rhd(dT=T-10,dS=S-35,Z) = [-a0*(1+lambda1/2*dT+mu1*Z)*dT ' 
    1795             WRITE(numout,*) '                                       + b0*(1+lambda2/2*dT+mu2*Z)*dS - nu*dT*dS] / rau0' 
     1795            WRITE(numout,*) '                                       + b0*(1+lambda2/2*dT+mu2*Z)*dS - nu*dT*dS] / rho0' 
    17961796            WRITE(numout,*) '              with the following coefficients :' 
    17971797            WRITE(numout,*) '                 thermal exp. coef.    rn_a0      = ', rn_a0 
     
    18101810            WRITE(numout,*) 
    18111811            WRITE(numout,*) '          use of linear ISOMIP eos:    rhd(dT=T-(-1),dS=S-(34.2),Z) = ' 
    1812             WRITE(numout,*) '             [ -a0*dT + b0*dS ]/rau0' 
     1812            WRITE(numout,*) '             [ -a0*dT + b0*dS ]/rho0' 
    18131813            WRITE(numout,*) 
    18141814            WRITE(numout,*) '             thermal exp. coef.    rn_a0      = ', rn_a0 
     
    18221822      END SELECT 
    18231823      ! 
    1824       rau0_rcp    = rau0 * rcp  
    1825       r1_rau0     = 1._wp / rau0 
     1824      rho0_rcp    = rho0 * rcp  
     1825      r1_rho0     = 1._wp / rho0 
    18261826      r1_rcp      = 1._wp / rcp 
    1827       r1_rau0_rcp = 1._wp / rau0_rcp  
     1827      r1_rho0_rcp = 1._wp / rho0_rcp  
    18281828      ! 
    18291829      IF(lwp) THEN 
     
    18401840      IF(lwp) WRITE(numout,*) 
    18411841      IF(lwp) WRITE(numout,*) '   Associated physical constant' 
    1842       IF(lwp) WRITE(numout,*) '      volumic mass of reference           rau0  = ', rau0   , ' kg/m^3' 
    1843       IF(lwp) WRITE(numout,*) '      1. / rau0                        r1_rau0  = ', r1_rau0, ' m^3/kg' 
     1842      IF(lwp) WRITE(numout,*) '      volumic mass of reference           rho0  = ', rho0   , ' kg/m^3' 
     1843      IF(lwp) WRITE(numout,*) '      1. / rho0                        r1_rho0  = ', r1_rho0, ' m^3/kg' 
    18441844      IF(lwp) WRITE(numout,*) '      ocean specific heat                 rcp   = ', rcp    , ' J/Kelvin' 
    1845       IF(lwp) WRITE(numout,*) '      rau0 * rcp                       rau0_rcp = ', rau0_rcp 
    1846       IF(lwp) WRITE(numout,*) '      1. / ( rau0 * rcp )           r1_rau0_rcp = ', r1_rau0_rcp 
     1845      IF(lwp) WRITE(numout,*) '      rho0 * rcp                       rho0_rcp = ', rho0_rcp 
     1846      IF(lwp) WRITE(numout,*) '      1. / ( rho0 * rcp )           r1_rho0_rcp = ', r1_rho0_rcp 
    18471847      ! 
    18481848   END SUBROUTINE eos_init 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/tests/ISOMIP+/MY_SRC/sbcfwb.F90

    r12406 r12443  
    122122            ! avoid the model to blow up for large ssh drop (isomip OCEAN3 with melt switch off and uniform T/S) 
    123123            IF (ln_isfcpl .AND. ln_isfcpl_cons) THEN 
    124                z_fwf = z_fwf + glob_sum( 'sbcfwb',  e1e2t(:,:) * risfcpl_cons_ssh(:,:) * rau0 ) 
     124               z_fwf = z_fwf + glob_sum( 'sbcfwb',  e1e2t(:,:) * risfcpl_cons_ssh(:,:) * rho0 ) 
    125125            END IF 
    126126            ! 
     
    155155            a_fwb_b = a_fwb                           ! mean sea level taking into account the ice+snow 
    156156                                                      ! sum over the global domain 
    157             a_fwb   = glob_sum( 'sbcfwb', e1e2t(:,:) * ( ssh(:,:,Kmm) + snwice_mass(:,:) * r1_rau0 ) ) 
     157            a_fwb   = glob_sum( 'sbcfwb', e1e2t(:,:) * ( ssh(:,:,Kmm) + snwice_mass(:,:) * r1_rho0 ) ) 
    158158            a_fwb   = a_fwb * 1.e+3 / ( area * rday * 365. )     ! convert in Kg/m3/s = mm/s 
    159159!!gm        !                                                      !!bug 365d year  
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/tests/LOCK_EXCHANGE/EXPREF/namelist_FCT2_flux_cen2_cfg

    r12406 r12443  
    144144      ln_bt_av      = .true.     ! Time filtering of barotropic variables 
    145145         nn_bt_flt     = 1          ! Time filter choice  = 0 None 
    146          !                          !                     = 1 Boxcar over   nn_baro sub-steps 
    147          !                          !                     = 2 Boxcar over 2*nn_baro  "    " 
     146         !                          !                     = 1 Boxcar over   nn_e sub-steps 
     147         !                          !                     = 2 Boxcar over 2*nn_e  "    " 
    148148      ln_bt_auto    = .true.     ! Number of sub-step defined from: 
    149          nn_baro      =  1         ! =F : the number of sub-step in rn_Dt seconds 
     149         nn_e      =  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

    r12406 r12443  
    216216      ln_bt_av      = .true.     ! Time filtering of barotropic variables 
    217217         nn_bt_flt     = 1          ! Time filter choice  = 0 None 
    218          !                          !                     = 1 Boxcar over   nn_baro sub-steps 
    219          !                          !                     = 2 Boxcar over 2*nn_baro  "    " 
     218         !                          !                     = 1 Boxcar over   nn_e sub-steps 
     219         !                          !                     = 2 Boxcar over 2*nn_e  "    " 
    220220      ln_bt_auto    = .true.     ! Number of sub-step defined from: 
    221          nn_baro      =  1         ! =F : the number of sub-step in rn_Dt seconds 
     221         nn_e      =  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

    r12406 r12443  
    144144      ln_bt_av      = .true.     ! Time filtering of barotropic variables 
    145145         nn_bt_flt     = 1          ! Time filter choice  = 0 None 
    146          !                          !                     = 1 Boxcar over   nn_baro sub-steps 
    147          !                          !                     = 2 Boxcar over 2*nn_baro  "    " 
     146         !                          !                     = 1 Boxcar over   nn_e sub-steps 
     147         !                          !                     = 2 Boxcar over 2*nn_e  "    " 
    148148      ln_bt_auto    = .true.     ! Number of sub-step defined from: 
    149          nn_baro      =  1         ! =F : the number of sub-step in rn_Dt seconds 
     149         nn_e      =  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

    r12406 r12443  
    144144      ln_bt_av      = .true.     ! Time filtering of barotropic variables 
    145145         nn_bt_flt     = 1          ! Time filter choice  = 0 None 
    146          !                          !                     = 1 Boxcar over   nn_baro sub-steps 
    147          !                          !                     = 2 Boxcar over 2*nn_baro  "    " 
     146         !                          !                     = 1 Boxcar over   nn_e sub-steps 
     147         !                          !                     = 2 Boxcar over 2*nn_e  "    " 
    148148      ln_bt_auto    = .true.     ! Number of sub-step defined from: 
    149          nn_baro      =  1         ! =F : the number of sub-step in rn_Dt seconds 
     149         nn_e      =  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

    r12406 r12443  
    144144      ln_bt_av      = .true.     ! Time filtering of barotropic variables 
    145145         nn_bt_flt     = 1          ! Time filter choice  = 0 None 
    146          !                          !                     = 1 Boxcar over   nn_baro sub-steps 
    147          !                          !                     = 2 Boxcar over 2*nn_baro  "    " 
     146         !                          !                     = 1 Boxcar over   nn_e sub-steps 
     147         !                          !                     = 2 Boxcar over 2*nn_e  "    " 
    148148      ln_bt_auto    = .true.     ! Number of sub-step defined from: 
    149          nn_baro      =  1         ! =F : the number of sub-step in rn_Dt seconds 
     149         nn_e      =  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

    r12406 r12443  
    144144      ln_bt_av      = .true.     ! Time filtering of barotropic variables 
    145145         nn_bt_flt     = 1          ! Time filter choice  = 0 None 
    146          !                          !                     = 1 Boxcar over   nn_baro sub-steps 
    147          !                          !                     = 2 Boxcar over 2*nn_baro  "    " 
     146         !                          !                     = 1 Boxcar over   nn_e sub-steps 
     147         !                          !                     = 2 Boxcar over 2*nn_e  "    " 
    148148      ln_bt_auto    = .true.     ! Number of sub-step defined from: 
    149          nn_baro      =  1         ! =F : the number of sub-step in rn_Dt seconds 
     149         nn_e      =  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

    r12406 r12443  
    144144      ln_bt_av      = .true.     ! Time filtering of barotropic variables 
    145145         nn_bt_flt     = 1          ! Time filter choice  = 0 None 
    146          !                          !                     = 1 Boxcar over   nn_baro sub-steps 
    147          !                          !                     = 2 Boxcar over 2*nn_baro  "    " 
     146         !                          !                     = 1 Boxcar over   nn_e sub-steps 
     147         !                          !                     = 2 Boxcar over 2*nn_e  "    " 
    148148      ln_bt_auto    = .true.     ! Number of sub-step defined from: 
    149          nn_baro      =  1         ! =F : the number of sub-step in rn_Dt seconds 
     149         nn_e      =  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

    r12406 r12443  
    144144      ln_bt_av      = .true.     ! Time filtering of barotropic variables 
    145145         nn_bt_flt     = 1          ! Time filter choice  = 0 None 
    146          !                          !                     = 1 Boxcar over   nn_baro sub-steps 
    147          !                          !                     = 2 Boxcar over 2*nn_baro  "    " 
     146         !                          !                     = 1 Boxcar over   nn_e sub-steps 
     147         !                          !                     = 2 Boxcar over 2*nn_e  "    " 
    148148      ln_bt_auto    = .true.     ! Number of sub-step defined from: 
    149          nn_baro      =  1         ! =F : the number of sub-step in rn_Dt seconds 
     149         nn_e      =  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

    r12406 r12443  
    144144      ln_bt_av      = .true.     ! Time filtering of barotropic variables 
    145145         nn_bt_flt     = 1          ! Time filter choice  = 0 None 
    146          !                          !                     = 1 Boxcar over   nn_baro sub-steps 
    147          !                          !                     = 2 Boxcar over 2*nn_baro  "    " 
     146         !                          !                     = 1 Boxcar over   nn_e sub-steps 
     147         !                          !                     = 2 Boxcar over 2*nn_e  "    " 
    148148      ln_bt_auto    = .true.     ! Number of sub-step defined from: 
    149          nn_baro      =  1         ! =F : the number of sub-step in rn_Dt seconds 
     149         nn_e      =  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

    r12406 r12443  
    144144      ln_bt_av      = .true.     ! Time filtering of barotropic variables 
    145145         nn_bt_flt     = 1          ! Time filter choice  = 0 None 
    146          !                          !                     = 1 Boxcar over   nn_baro sub-steps 
    147          !                          !                     = 2 Boxcar over 2*nn_baro  "    " 
     146         !                          !                     = 1 Boxcar over   nn_e sub-steps 
     147         !                          !                     = 2 Boxcar over 2*nn_e  "    " 
    148148      ln_bt_auto    = .true.     ! Number of sub-step defined from: 
    149          nn_baro      =  1         ! =F : the number of sub-step in rn_Dt seconds 
     149         nn_e      =  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

    r12406 r12443  
    144144      ln_bt_av      = .true.     ! Time filtering of barotropic variables 
    145145         nn_bt_flt     = 1          ! Time filter choice  = 0 None 
    146          !                          !                     = 1 Boxcar over   nn_baro sub-steps 
    147          !                          !                     = 2 Boxcar over 2*nn_baro  "    " 
     146         !                          !                     = 1 Boxcar over   nn_e sub-steps 
     147         !                          !                     = 2 Boxcar over 2*nn_e  "    " 
    148148      ln_bt_auto    = .true.     ! Number of sub-step defined from: 
    149          nn_baro      =  1         ! =F : the number of sub-step in rn_Dt seconds 
     149         nn_e      =  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

    r12406 r12443  
    144144      ln_bt_av      = .true.     ! Time filtering of barotropic variables 
    145145         nn_bt_flt     = 1          ! Time filter choice  = 0 None 
    146          !                          !                     = 1 Boxcar over   nn_baro sub-steps 
    147          !                          !                     = 2 Boxcar over 2*nn_baro  "    " 
     146         !                          !                     = 1 Boxcar over   nn_e sub-steps 
     147         !                          !                     = 2 Boxcar over 2*nn_e  "    " 
    148148      ln_bt_auto    = .true.     ! Number of sub-step defined from: 
    149          nn_baro      =  1         ! =F : the number of sub-step in rn_Dt seconds 
     149         nn_e      =  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/MY_SRC/usrdef_istate.F90

    r10074 r12443  
    6060      ! 
    6161      !  rn_a0 =  0.2   !  thermal expension coefficient (nn_eos= 1) 
    62       !  rho = rau0 - rn_a0 * (T-10)  
     62      !  rho = rho0 - rn_a0 * (T-10)  
    6363      !  delta_T = 25 degrees  ==>>  delta_rho = 25 * rn_a0 = 5 kg/m3 
    6464      ! 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/tests/OVERFLOW/EXPREF/namelist_sco_FCT2_flux_cen-ahm1000_cfg

    r12406 r12443  
    155155      ln_bt_av      = .true.     ! Time filtering of barotropic variables 
    156156         nn_bt_flt     = 1          ! Time filter choice  = 0 None 
    157          !                          !                     = 1 Boxcar over   nn_baro sub-steps 
    158          !                          !                     = 2 Boxcar over 2*nn_baro  "    " 
     157         !                          !                     = 1 Boxcar over   nn_e sub-steps 
     158         !                          !                     = 2 Boxcar over 2*nn_e  "    " 
    159159      ln_bt_auto    = .true.     ! Number of sub-step defined from: 
    160          nn_baro      =  1         ! =F : the number of sub-step in rn_Dt seconds 
     160         nn_e      =  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

    r12406 r12443  
    155155      ln_bt_av      = .true.     ! Time filtering of barotropic variables 
    156156         nn_bt_flt     = 1          ! Time filter choice  = 0 None 
    157          !                          !                     = 1 Boxcar over   nn_baro sub-steps 
    158          !                          !                     = 2 Boxcar over 2*nn_baro  "    " 
     157         !                          !                     = 1 Boxcar over   nn_e sub-steps 
     158         !                          !                     = 2 Boxcar over 2*nn_e  "    " 
    159159      ln_bt_auto    = .true.     ! Number of sub-step defined from: 
    160          nn_baro      =  1         ! =F : the number of sub-step in rn_Dt seconds 
     160         nn_e      =  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

    r12406 r12443  
    155155      ln_bt_av      = .true.     ! Time filtering of barotropic variables 
    156156         nn_bt_flt     = 1          ! Time filter choice  = 0 None 
    157          !                          !                     = 1 Boxcar over   nn_baro sub-steps 
    158          !                          !                     = 2 Boxcar over 2*nn_baro  "    " 
     157         !                          !                     = 1 Boxcar over   nn_e sub-steps 
     158         !                          !                     = 2 Boxcar over 2*nn_e  "    " 
    159159      ln_bt_auto    = .true.     ! Number of sub-step defined from: 
    160          nn_baro      =  1         ! =F : the number of sub-step in rn_Dt seconds 
     160         nn_e      =  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

    r12406 r12443  
    155155      ln_bt_av      = .true.     ! Time filtering of barotropic variables 
    156156         nn_bt_flt     = 1          ! Time filter choice  = 0 None 
    157          !                          !                     = 1 Boxcar over   nn_baro sub-steps 
    158          !                          !                     = 2 Boxcar over 2*nn_baro  "    " 
     157         !                          !                     = 1 Boxcar over   nn_e sub-steps 
     158         !                          !                     = 2 Boxcar over 2*nn_e  "    " 
    159159      ln_bt_auto    = .true.     ! Number of sub-step defined from: 
    160          nn_baro      =  1         ! =F : the number of sub-step in rn_Dt seconds 
     160         nn_e      =  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

    r12406 r12443  
    155155      ln_bt_av      = .true.     ! Time filtering of barotropic variables 
    156156         nn_bt_flt     = 1          ! Time filter choice  = 0 None 
    157          !                          !                     = 1 Boxcar over   nn_baro sub-steps 
    158          !                          !                     = 2 Boxcar over 2*nn_baro  "    " 
     157         !                          !                     = 1 Boxcar over   nn_e sub-steps 
     158         !                          !                     = 2 Boxcar over 2*nn_e  "    " 
    159159      ln_bt_auto    = .true.     ! Number of sub-step defined from: 
    160          nn_baro      =  1         ! =F : the number of sub-step in rn_Dt seconds 
     160         nn_e      =  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

    r12406 r12443  
    216216      ln_bt_av      = .true.     ! Time filtering of barotropic variables 
    217217         nn_bt_flt     = 1          ! Time filter choice  = 0 None 
    218          !                          !                     = 1 Boxcar over   nn_baro sub-steps 
    219          !                          !                     = 2 Boxcar over 2*nn_baro  "    " 
     218         !                          !                     = 1 Boxcar over   nn_e sub-steps 
     219         !                          !                     = 2 Boxcar over 2*nn_e  "    " 
    220220      ln_bt_auto    = .true.     ! Number of sub-step defined from: 
    221          nn_baro      =  1         ! =F : the number of sub-step in rn_Dt seconds 
     221         nn_e      =  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

    r12406 r12443  
    155155      ln_bt_av      = .true.     ! Time filtering of barotropic variables 
    156156         nn_bt_flt     = 1          ! Time filter choice  = 0 None 
    157          !                          !                     = 1 Boxcar over   nn_baro sub-steps 
    158          !                          !                     = 2 Boxcar over 2*nn_baro  "    " 
     157         !                          !                     = 1 Boxcar over   nn_e sub-steps 
     158         !                          !                     = 2 Boxcar over 2*nn_e  "    " 
    159159      ln_bt_auto    = .true.     ! Number of sub-step defined from: 
    160          nn_baro      =  1         ! =F : the number of sub-step in rn_Dt seconds 
     160         nn_e      =  1         ! =F : the number of sub-step in rn_Dt seconds 
    161161/ 
    162162!----------------------------------------------------------------------- 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/tests/OVERFLOW/MY_SRC/usrdef_istate.F90

    r10074 r12443  
    6060      ! 
    6161      !  rn_a0 =  0.2   !  thermal expension coefficient (nn_eos= 1) 
    62       !  rho = rau0 - rn_a0 * (T-10)  
     62      !  rho = rho0 - rn_a0 * (T-10)  
    6363      !  delta_T = 10 degrees  ==>>  delta_rho = 10 * rn_a0 = 2 kg/m3 
    6464      ! 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/tests/VORTEX/EXPREF/1_namelist_cfg

    r12406 r12443  
    213213      ln_bt_av      = .true.     ! Time filtering of barotropic variables 
    214214         nn_bt_flt     = 1          ! Time filter choice  = 0 None 
    215          !                          !                     = 1 Boxcar over   nn_baro sub-steps 
    216          !                          !                     = 2 Boxcar over 2*nn_baro  "    " 
     215         !                          !                     = 1 Boxcar over   nn_e sub-steps 
     216         !                          !                     = 2 Boxcar over 2*nn_e  "    " 
    217217      ln_bt_auto    = .false.    ! Number of sub-step defined from: 
    218          nn_baro      =  24         ! =F : the number of sub-step in rn_Dt seconds 
     218         nn_e      =  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

    r12406 r12443  
    204204      ln_bt_av      = .true.     ! Time filtering of barotropic variables 
    205205         nn_bt_flt     = 1          ! Time filter choice  = 0 None 
    206          !                          !                     = 1 Boxcar over   nn_baro sub-steps 
    207          !                          !                     = 2 Boxcar over 2*nn_baro  "    " 
     206         !                          !                     = 1 Boxcar over   nn_e sub-steps 
     207         !                          !                     = 2 Boxcar over 2*nn_e  "    " 
    208208      ln_bt_auto    = .false.    ! Number of sub-step defined from: 
    209          nn_baro      =  24         ! =F : the number of sub-step in rn_Dt seconds 
     209         nn_e      =  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

    r12424 r12443  
    637637         ELSE 
    638638            tilde_e3t_b(:,:,:) = tilde_e3t_n(:,:,:) &  
    639             &         + atfp * ( tilde_e3t_b(:,:,:) - 2.0_wp * tilde_e3t_n(:,:,:) + tilde_e3t_a(:,:,:) ) 
     639            &         + rn_atfp * ( tilde_e3t_b(:,:,:) - 2.0_wp * tilde_e3t_n(:,:,:) + tilde_e3t_a(:,:,:) ) 
    640640         ENDIF 
    641641         tilde_e3t_n(:,:,:) = tilde_e3t_a(:,:,:) 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/tests/VORTEX/MY_SRC/usrdef_istate.F90

    r10425 r12443  
    6969      zH = 0.5_wp * 5000._wp 
    7070      ! 
    71       zP0 = rau0 * zf0 * zumax * zlambda * SQRT(EXP(1._wp)/2._wp) 
     71      zP0 = rho0 * zf0 * zumax * zlambda * SQRT(EXP(1._wp)/2._wp) 
    7272      ! 
    7373      ! Sea level: 
     
    7777            zx = glamt(ji,jj) * 1.e3 
    7878            zy = gphit(ji,jj) * 1.e3 
    79             zrho1 = rau0 + za * EXP(-(zx**2+zy**2)/zlambda**2) 
     79            zrho1 = rho0 + za * EXP(-(zx**2+zy**2)/zlambda**2) 
    8080            pssh(ji,jj) = zP0 * EXP(-(zx**2+zy**2)/zlambda**2)/(zrho1*grav) * ptmask(ji,jj,1) 
    8181         END DO 
     
    8989            DO jk=1,jpk 
    9090               zdt =  pdept(ji,jj,jk)  
    91                zrho1 = rau0 * (1._wp + zn2*zdt/grav) 
     91               zrho1 = rho0 * (1._wp + zn2*zdt/grav) 
    9292               IF (zdt < zH) THEN 
    9393                  zrho1 = zrho1 - zP0 * (1._wp-EXP(zdt-zH)) & 
    9494                          & * EXP(-(zx**2+zy**2)/zlambda**2) / (grav*(zH -1._wp + exp(-zH))); 
    9595               ENDIF 
    96                pts(ji,jj,jk,jp_tem) = (20._wp + (rau0-zrho1) / 0.28_wp) * ptmask(ji,jj,jk) 
     96               pts(ji,jj,jk,jp_tem) = (20._wp + (rho0-zrho1) / 0.28_wp) * ptmask(ji,jj,jk) 
    9797            END DO 
    9898         END DO 
     
    103103      ! 
    104104      ! velocities: 
    105       za = 2._wp * zP0 / (zf0 * rau0 * zlambda**2) 
     105      za = 2._wp * zP0 / (zf0 * rho0 * zlambda**2) 
    106106      DO ji=1, jpim1 
    107107         DO jj=1, jpj 
  • NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/tests/WAD/EXPREF/namelist_cfg

    r12406 r12443  
    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_Dt seconds 
     346   nn_e       =  12     ! =F : the number of sub-step in rn_Dt seconds 
    347347/ 
    348348!----------------------------------------------------------------------- 
Note: See TracChangeset for help on using the changeset viewer.