Changeset 1858 for branches/DEV_r1837_mass_heat_salt_fluxes/NEMO/OPA_SRC/DYN/dynzdf_imp.F90

Timestamp:

2010-05-04T10:39:48+02:00 (14 years ago)

Author:

gm

Message:

ticket:#665 : step 1 - heat content of freezing-melting ice

File:

• branches/DEV_r1837_mass_heat_salt_fluxes/NEMO/OPA_SRC/DYN/dynzdf_imp.F90 (modified) (2 diffs)

branches/DEV_r1837_mass_heat_salt_fluxes/NEMO/OPA_SRC/DYN/dynzdf_imp.F90

@@ -59 +59 @@
       !!         ua = ua + dz( avmu dz(u) )
       !!
-      !! ** Action : - Update (ua,va) arrays with the after vertical diffusive
-      !!               mixing trend.
+      !! ** Action : -  (ua,va) updated with the after vertical diffusive mixing trend
       !!---------------------------------------------------------------------
-      !! * Modules used
-      USE oce, ONLY :  zwd   => ta,   &                ! use ta as workspace
-                       zws   => sa                     ! use sa as workspace
-
-      !! * Arguments
-      INTEGER , INTENT( in ) ::   kt                   ! ocean time-step index
-      REAL(wp), INTENT( in ) ::  p2dt                  ! vertical profile of tracer time-step
-
-      !! * Local declarations
-      INTEGER ::   ji, jj, jk                          ! dummy loop indices
-      REAL(wp) ::   zrau0r, z2dtf, zcoef, zzws, zrhs   ! temporary scalars
-      REAL(wp) ::   zzwi                               ! temporary scalars
-      REAL(wp), DIMENSION(jpi,jpj,jpk) ::   zwi        ! temporary workspace arrays
+      USE oce, ONLY :  zwd   => ta   ! use ta as workspace
+      USE oce, ONLY :  zws   => sa   !  -  sa         -
+      !!
+      INTEGER , INTENT( in ) ::   kt    ! ocean time-step index
+      REAL(wp), INTENT( in ) ::  p2dt   ! vertical profile of tracer time-step
+      !!
+      INTEGER ::   ji, jj, jk                     ! dummy loop indices
+      REAL(wp) ::   z2dtf, zcoef, zzws, zrhs      ! temporary scalars
+      REAL(wp) ::   zzwi                          !    -         -
+      REAL(wp), DIMENSION(jpi,jpj,jpk) ::   zwi   ! 3D workspace 
       !!----------------------------------------------------------------------
 
@@ -82 +78 @@
          IF(lwp) WRITE(numout,*) '~~~~~~~~~~~ '
       ENDIF
-
-      ! 0. Local constant initialization
-      ! --------------------------------
-      zrau0r = 1. / rau0      ! inverse of the reference density
 
       ! 1. Vertical diffusion on u