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 15533 – NEMO

Changeset 15533


Ignore:
Timestamp:
2021-11-24T13:07:20+01:00 (6 months ago)
Author:
cdllod
Message:

renamed two constants

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2021/dev_r15388_updated_zdfiwm/src/OCE/ZDF/zdfiwm.F90

    r15526 r15533  
    4242   LOGICAL ::  ln_tsdiff   ! account for differential T/S wave-driven mixing (=T) or not (=F) 
    4343 
    44    REAL(wp)::  z1_6 = 1._wp / 6._wp 
    45    REAL(wp)::  znu  = 1.4e-6_wp   ! molecular kinematic viscosity 
     44   REAL(wp)::  r1_6 = 1._wp / 6._wp 
     45   REAL(wp)::  rnu  = 1.4e-6_wp   ! molecular kinematic viscosity 
    4646 
    4747   REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:) ::   ebot_iwm   ! bottom-intensified dissipation above abyssal hills (W/m2) 
     
    219219      ! Calculate turbulence intensity parameter Reb 
    220220      DO_3D( nn_hls-1, nn_hls-1, nn_hls-1, nn_hls-1, 2, jpkm1 ) 
    221          zReb(ji,jj,jk) = zemx_iwm(ji,jj,jk) / MAX( 1.e-20_wp, znu * rn2(ji,jj,jk) ) 
     221         zReb(ji,jj,jk) = zemx_iwm(ji,jj,jk) / MAX( 1.e-20_wp, rnu * rn2(ji,jj,jk) ) 
    222222      END_3D 
    223223      ! 
    224224      ! Define internal wave-induced diffusivity 
    225225      DO_3D( nn_hls-1, nn_hls-1, nn_hls-1, nn_hls-1, 2, jpkm1 ) 
    226          zav_wave(ji,jj,jk) = zReb(ji,jj,jk) * z1_6 * znu  ! This corresponds to a constant mixing efficiency of 1/6 
     226         zav_wave(ji,jj,jk) = zReb(ji,jj,jk) * r1_6 * rnu  ! This corresponds to a constant mixing efficiency of 1/6 
    227227      END_3D 
    228228      ! 
     
    230230         DO_3D( nn_hls-1, nn_hls-1, nn_hls-1, nn_hls-1, 2, jpkm1 ) ! energetic (Reb > 480) and buoyancy-controlled (Reb <10.224) regimes 
    231231            IF( zReb(ji,jj,jk) > 480.00_wp ) THEN 
    232                zav_wave(ji,jj,jk) = 3.6515_wp * znu * SQRT( zReb(ji,jj,jk) ) 
     232               zav_wave(ji,jj,jk) = 3.6515_wp * rnu * SQRT( zReb(ji,jj,jk) ) 
    233233            ELSEIF( zReb(ji,jj,jk) < 10.224_wp ) THEN 
    234                zav_wave(ji,jj,jk) = 0.052125_wp * znu * zReb(ji,jj,jk) * SQRT( zReb(ji,jj,jk) ) 
     234               zav_wave(ji,jj,jk) = 0.052125_wp * rnu * zReb(ji,jj,jk) * SQRT( zReb(ji,jj,jk) ) 
    235235            ENDIF 
    236236         END_3D 
     
    248248         DO_3D( nn_hls-1, nn_hls-1, nn_hls-1, nn_hls-1, 2, jpkm1 ) ! Calculate S/T diffusivity ratio as a function of Reb (else it is set to 1) 
    249249            zav_ratio(ji,jj,jk) = ( 0.505_wp + & 
    250                &                    0.495_wp * TANH( 0.92_wp * ( LOG10( MAX( 1.e-20, zReb(ji,jj,jk) * 5._wp * z1_6 ) ) - 0.60_wp ) ) & 
     250               &                    0.495_wp * TANH( 0.92_wp * ( LOG10( MAX( 1.e-20, zReb(ji,jj,jk) * 5._wp * r1_6 ) ) - 0.60_wp ) ) & 
    251251               &                  ) * wmask(ji,jj,jk) 
    252252         END_3D 
     
    372372      ! ensures that avt remains larger than its molecular value (=1.4e-7). Therefore, avtb should  
    373373      ! be set here to a very small value, and avmb to its (uniform) molecular value (=1.4e-6). 
    374       avmb(:) = znu              ! molecular value 
     374      avmb(:) = rnu              ! molecular value 
    375375      avtb(:) = 1.e-10_wp        ! very small diffusive minimum (background avt is specified in zdf_iwm)     
    376376      avtb_2d(:,:) = 1._wp       ! uniform  
Note: See TracChangeset for help on using the changeset viewer.