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/KERNEL-03_Storkey_Coward_RK3_stage2/tests – NEMO

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/tests
Files:
37 edited

Legend:

Unmodified
Added
Removed

  • 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.