Changeset 7123 for branches/UKMO/r5518_bdy_sponge_temp/NEMOGCM/NEMO/OPA_SRC/TRA/traldf_lap.F90

Timestamp:

2016-10-26T19:32:29+02:00 (8 years ago)

Author:

jcastill

Message:

Changes as in UKMO/dev_r5518_bdy_sponge_temp@5878

File:

• branches/UKMO/r5518_bdy_sponge_temp/NEMOGCM/NEMO/OPA_SRC/TRA/traldf_lap.F90 (modified) (4 diffs)

branches/UKMO/r5518_bdy_sponge_temp/NEMOGCM/NEMO/OPA_SRC/TRA/traldf_lap.F90

@@ -30 +30 @@
    USE lib_mpp         ! distribued memory computing library
    USE timing          ! Timing
+#if defined key_bdy  
+   USE bdy_oce  
+#endif
    USE wrk_nemo        ! Memory allocation
 
@@ -81 +84 @@
       REAL(wp) ::   zsign            ! local scalars
       REAL(wp), POINTER, DIMENSION(:,:,:) ::  ztu, ztv, zaheeu, zaheev
+      REAL(wp), DIMENSION(jpi,jpj)        ::  zfactor  ! multiplier for diffusion
       !!----------------------------------------------------------------------
       !
@@ -90 +94 @@
          WRITE(numout,*) '~~~~~~~~~~~ '
       ENDIF
+#if defined key_bdy 
+     zfactor(:,:) = sponge_factor(:,:) 
+#else 
+     zfactor(:,:) = 1.0 
+#endif 
       !
       CALL wrk_alloc( jpi,jpj,jpk,   ztu, ztv, zaheeu, zaheev ) 
@@ -100 +109 @@
          DO jj = 1, jpjm1
             DO ji = 1, fs_jpim1   ! vector opt.
-               zaheeu(ji,jj,jk) = zsign * pahu(ji,jj,jk) * e2_e1u(ji,jj) * fse3u_n(ji,jj,jk)   !!gm   * umask(ji,jj,jk) pah masked!
-               zaheev(ji,jj,jk) = zsign * pahv(ji,jj,jk) * e1_e2v(ji,jj) * fse3v_n(ji,jj,jk)   !!gm   * vmask(ji,jj,jk)
+               zaheeu(ji,jj,jk) = zfactor(ji,jj) * zsign * pahu(ji,jj,jk) * e2_e1u(ji,jj) * fse3u_n(ji,jj,jk)   !!gm   * umask(ji,jj,jk) pah masked!
+               zaheev(ji,jj,jk) = zfactor(ji,jj) * zsign * pahv(ji,jj,jk) * e1_e2v(ji,jj) * fse3v_n(ji,jj,jk)   !!gm   * vmask(ji,jj,jk)
             END DO
          END DO