- Timestamp:
- 2016-07-19T10:38:35+02:00 (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/NERC/dev_r5549_BDY_ZEROGRAD/NEMOGCM/NEMO/OPA_SRC/ZDF/zdfevd.F90
r4990 r6808 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 … … 29 29 PUBLIC zdf_evd ! called by step.F90 30 30 31 !! * Substitutions32 # include "domzgr_substitute.h90"33 31 !!---------------------------------------------------------------------- 34 32 !! NEMO/OPA 4.0 , NEMO Consortium (2011) … … 53 51 !! References : Lazar, A., these de l'universite Paris VI, France, 1997 54 52 !!---------------------------------------------------------------------- 55 USE oce, zavt_evd => ua , zavm_evd => va ! (ua,va) used ua workspace56 !57 53 INTEGER, INTENT( in ) :: kt ! ocean time-step indexocean time step 58 54 ! 59 55 INTEGER :: ji, jj, jk ! dummy loop indices 56 REAL(wp), POINTER, DIMENSION(:,:,:) :: zavt_evd, zavm_evd 60 57 !!---------------------------------------------------------------------- 61 58 ! … … 68 65 IF(lwp) WRITE(numout,*) 69 66 ENDIF 70 67 ! 68 CALL wrk_alloc( jpi,jpj,jpk, zavt_evd, zavm_evd ) 69 ! 71 70 zavt_evd(:,:,:) = avt(:,:,:) ! set avt prior to evd application 72 71 ! 73 72 SELECT CASE ( nn_evdm ) 74 73 ! … … 80 79 DO jj = 2, jpj ! no vector opt. 81 80 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 81 IF( MIN( rn2(ji,jj,jk), rn2b(ji,jj,jk) ) <= -1.e-12 ) THEN 87 #endif88 82 avt (ji ,jj ,jk) = rn_avevd * tmask(ji ,jj ,jk) 89 83 avm (ji ,jj ,jk) = rn_avevd * tmask(ji ,jj ,jk) … … 107 101 DO jj = 1, jpj ! loop over the whole domain (no lbc_lnk call) 108 102 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 103 IF( MIN( rn2(ji,jj,jk), rn2b(ji,jj,jk) ) <= -1.e-12 ) & 114 #endif115 104 avt(ji,jj,jk) = rn_avevd * tmask(ji,jj,jk) 116 105 END DO … … 123 112 CALL iom_put( "avt_evd", zavt_evd ) ! output this change 124 113 ! 114 CALL wrk_dealloc( jpi,jpj,jpk, zavt_evd, zavm_evd ) 115 ! 125 116 IF( nn_timing == 1 ) CALL timing_stop('zdf_evd') 126 117 !
Note: See TracChangeset
for help on using the changeset viewer.