Changeset 9939 for NEMO/branches/2018/dev_r9838_ENHANCE04_RK3/src/OCE/TRD/trdken.F90

Timestamp:

2018-07-13T09:28:50+02:00 (6 years ago)

Author:

gm

Message:

#1911 (ENHANCE-04): RK3 branche phased with MLF@9937 branche

Location:

NEMO/branches/2018/dev_r9838_ENHANCE04_RK3

Files:

• . (copied) (copied from NEMO/trunk)
• src/OCE/TRD/trdken.F90 (modified) (5 diffs)

NEMO/branches/2018/dev_r9838_ENHANCE04_RK3/src/OCE/TRD/trdken.F90

@@ -103 +103 @@
          DO jj = 2, jpj
             DO ji = 2, jpi
-               zke(ji,jj,jk) = 0.5_wp * rau0 *( un(ji  ,jj,jk) * putrd(ji  ,jj,jk) * bu(ji  ,jj,jk)  &
+               zke(ji,jj,jk) = 0.5_wp * rho0 *( un(ji  ,jj,jk) * putrd(ji  ,jj,jk) * bu(ji  ,jj,jk)  &
                   &                           + un(ji-1,jj,jk) * putrd(ji-1,jj,jk) * bu(ji-1,jj,jk)  &
                   &                           + vn(ji,jj  ,jk) * pvtrd(ji,jj  ,jk) * bv(ji,jj  ,jk)  &
@@ -127 +127 @@
                            DO jj = 2, jpj
                               DO ji = 2, jpi
-                                 zke2d(ji,jj) = r1_rau0 * 0.5_wp * (   z2dx(ji,jj) + z2dx(ji-1,jj)   &
+                                 zke2d(ji,jj) = r1_rho0 * 0.5_wp * (   z2dx(ji,jj) + z2dx(ji-1,jj)   &
                                  &                                   + z2dy(ji,jj) + z2dy(ji,jj-1)   ) * r1_bt(ji,jj,1)
                               END DO
@@ -184 +184 @@
                     !
                     CALL ken_p2k( kt , zke )
-                      CALL iom_put( "ketrd_convP2K", zke )     ! conversion -rau*g*w
+                      CALL iom_put( "ketrd_convP2K", zke )     ! conversion -rho*g*w
          !
       END SELECT
@@ -197 +197 @@
       !! ** Purpose :   compute rate of conversion from potential to kinetic energy
       !!
-      !! ** Method  : - compute conv defined as -rau*g*w on T-grid points
+      !! ** Method  : - compute conv defined as -rho*g*w on T-grid points
       !! 
       !! ** Work only for full steps and partial steps (ln_hpg_zco or ln_hpg_zps)
@@ -211 +211 @@
       !
       ! Local constant initialization 
-      zcoef = - rau0 * grav * 0.5_wp      
+      zcoef = - rho0 * grav * 0.5_wp      
       
       !  Surface value (also valid in partial step case)