Changeset 5989 for branches/2014/dev_r4650_UKMO10_Tidally_Meaned_Diagnostics/NEMOGCM/NEMO/OPA_SRC/ZDF/zdfevd.F90
- Timestamp:
- 2015-12-03T09:10:32+01:00 (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2014/dev_r4650_UKMO10_Tidally_Meaned_Diagnostics/NEMOGCM/NEMO/OPA_SRC/ZDF/zdfevd.F90
r5260 r5989 7 7 !! History : OPA ! 1997-06 (G. Madec, A. Lazar) Original code 8 8 !! NEMO 1.0 ! 2002-06 (G. Madec) F90: Free form and module 9 !! - ! 2005-06 (C. Ethe) KPP parameterization10 9 !! 3.2 ! 2009-03 (M. Leclair, G. Madec, R. Benshila) test on both before & after 11 10 !!---------------------------------------------------------------------- 12 11 13 12 !!---------------------------------------------------------------------- 14 !! zdf_evd : increase the momentum and tracer Kz at the location of15 !! statically unstable portion of the water column (ln_zdfevd=T)13 !! zdf_evd : increase the momentum and tracer Kz at the location of 14 !! statically unstable portion of the water column (ln_zdfevd=T) 16 15 !!---------------------------------------------------------------------- 17 16 USE oce ! ocean dynamics and tracers variables 18 17 USE dom_oce ! ocean space and time domain variables 19 18 USE zdf_oce ! ocean vertical physics variables 20 USE zdfkpp ! KPP vertical mixing19 ! 21 20 USE in_out_manager ! I/O manager 22 21 USE iom ! for iom_put 23 22 USE lbclnk ! ocean lateral boundary conditions (or mpp link) 23 USE wrk_nemo ! work arrays 24 24 USE timing ! Timing 25 25 … … 53 53 !! References : Lazar, A., these de l'universite Paris VI, France, 1997 54 54 !!---------------------------------------------------------------------- 55 USE oce, zavt_evd => ua , zavm_evd => va ! (ua,va) used ua workspace56 !57 55 INTEGER, INTENT( in ) :: kt ! ocean time-step indexocean time step 58 56 ! 59 57 INTEGER :: ji, jj, jk ! dummy loop indices 58 REAL(wp), POINTER, DIMENSION(:,:,:) :: zavt_evd, zavm_evd 60 59 !!---------------------------------------------------------------------- 61 60 ! … … 68 67 IF(lwp) WRITE(numout,*) 69 68 ENDIF 70 69 ! 70 CALL wrk_alloc( jpi,jpj,jpk, zavt_evd, zavm_evd ) 71 ! 71 72 zavt_evd(:,:,:) = avt(:,:,:) ! set avt prior to evd application 72 73 ! 73 74 SELECT CASE ( nn_evdm ) 74 75 ! … … 80 81 DO jj = 2, jpj ! no vector opt. 81 82 DO ji = 2, jpi 82 #if defined key_zdfkpp83 ! no evd mixing in the boundary layer with KPP84 IF( MIN( rn2(ji,jj,jk), rn2b(ji,jj,jk) ) <= -1.e-12 .AND. fsdepw(ji,jj,jk) > hkpp(ji,jj) ) THEN85 #else86 83 IF( MIN( rn2(ji,jj,jk), rn2b(ji,jj,jk) ) <= -1.e-12 ) THEN 87 #endif88 84 avt (ji ,jj ,jk) = rn_avevd * tmask(ji ,jj ,jk) 89 85 avm (ji ,jj ,jk) = rn_avevd * tmask(ji ,jj ,jk) … … 107 103 DO jj = 1, jpj ! loop over the whole domain (no lbc_lnk call) 108 104 DO ji = 1, jpi 109 #if defined key_zdfkpp110 ! no evd mixing in the boundary layer with KPP111 IF( MIN( rn2(ji,jj,jk), rn2b(ji,jj,jk) ) <= -1.e-12 .AND. fsdepw(ji,jj,jk) > hkpp(ji,jj) ) &112 #else113 105 IF( MIN( rn2(ji,jj,jk), rn2b(ji,jj,jk) ) <= -1.e-12 ) & 114 #endif115 106 avt(ji,jj,jk) = rn_avevd * tmask(ji,jj,jk) 116 107 END DO … … 123 114 CALL iom_put( "avt_evd", zavt_evd ) ! output this change 124 115 ! 116 CALL wrk_dealloc( jpi,jpj,jpk, zavt_evd, zavm_evd ) 117 ! 125 118 IF( nn_timing == 1 ) CALL timing_stop('zdf_evd') 126 119 !
Note: See TracChangeset
for help on using the changeset viewer.