- Timestamp:
- 2019-08-16T12:32:43+02:00 (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/UKMO/dev_r5518_GO6_package_FOAMv14/NEMOGCM/NEMO/OPA_SRC/ZDF/zdfevd.F90
r7061 r11442 25 25 USE lbclnk ! ocean lateral boundary conditions (or mpp link) 26 26 USE timing ! Timing 27 USE stopack 27 28 28 29 IMPLICIT NONE … … 30 31 31 32 PUBLIC zdf_evd ! called by step.F90 33 REAL(wp), SAVE, ALLOCATABLE, DIMENSION(:,:) :: rn_avevd0 32 34 33 35 !! * Substitutions … … 69 71 IF(lwp) WRITE(numout,*) '~~~~~~~ ' 70 72 IF(lwp) WRITE(numout,*) 73 ALLOCATE ( rn_avevd0(jpi,jpj) ) 74 rn_avevd0(:,:) = rn_avevd 71 75 ENDIF 72 76 73 77 zavt_evd(:,:,:) = avt(:,:,:) ! set avt prior to evd application 78 79 IF( ln_stopack .AND. ( nn_spp_aevd > 0 ) ) THEN 80 rn_avevd0(:,:) = rn_avevd 81 CALL spp_gen(kt, rn_avevd0, nn_spp_aevd, rn_aevd_sd, jk_spp_aevd) 82 ENDIF 74 83 75 84 SELECT CASE ( nn_evdm ) … … 88 97 IF( MIN( rn2(ji,jj,jk), rn2b(ji,jj,jk) ) <= -1.e-12 ) THEN 89 98 #endif 90 avt (ji ,jj ,jk) = rn_avevd * tmask(ji ,jj ,jk)91 avm (ji ,jj ,jk) = rn_avevd * tmask(ji ,jj ,jk)92 avmu(ji ,jj ,jk) = rn_avevd * umask(ji ,jj ,jk)93 avmu(ji-1,jj ,jk) = rn_avevd * umask(ji-1,jj ,jk)94 avmv(ji ,jj ,jk) = rn_avevd * vmask(ji ,jj ,jk)95 avmv(ji ,jj-1,jk) = rn_avevd * vmask(ji ,jj-1,jk)99 avt (ji ,jj ,jk) = rn_avevd0(ji,jj) * tmask(ji ,jj ,jk) 100 avm (ji ,jj ,jk) = rn_avevd0(ji,jj) * tmask(ji ,jj ,jk) 101 avmu(ji ,jj ,jk) = rn_avevd0(ji,jj) * umask(ji ,jj ,jk) 102 avmu(ji-1,jj ,jk) = rn_avevd0(ji,jj) * umask(ji-1,jj ,jk) 103 avmv(ji ,jj ,jk) = rn_avevd0(ji,jj) * vmask(ji ,jj ,jk) 104 avmv(ji ,jj-1,jk) = rn_avevd0(ji,jj) * vmask(ji ,jj-1,jk) 96 105 ENDIF 97 106 END DO … … 115 124 IF( MIN( rn2(ji,jj,jk), rn2b(ji,jj,jk) ) <= -1.e-12 ) & 116 125 #endif 117 avt(ji,jj,jk) = rn_avevd * tmask(ji,jj,jk)126 avt(ji,jj,jk) = rn_avevd0(ji,jj) * tmask(ji,jj,jk) 118 127 END DO 119 128 END DO
Note: See TracChangeset
for help on using the changeset viewer.