Changeset 1681 for trunk/NEMO/OPA_SRC/ZDF/zdfevd.F90
- Timestamp:
- 2009-10-27T18:44:08+01:00 (14 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/OPA_SRC/ZDF/zdfevd.F90
r1546 r1681 53 53 !! References : Lazar, A., these de l'universite Paris VI, France, 1997 54 54 !!---------------------------------------------------------------------- 55 USE oce, zav_evd => ua ! use ua as workspace 55 USE oce, zavt_evd => ua ! use ua as workspace 56 USE oce, zavm_evd => va ! use va as workspace 56 57 !! 57 58 INTEGER, INTENT( in ) :: kt ! ocean time-step indexocean time step … … 67 68 ENDIF 68 69 69 zav _evd(:,:,:) = avt(:,:,:) ! set avt prior to evd application70 zavt_evd(:,:,:) = avt(:,:,:) ! set avt prior to evd application 70 71 71 72 SELECT CASE ( nn_evdm ) 72 73 73 74 CASE ( 1 ) ! enhance vertical eddy viscosity and diffusivity (if rn2<-1.e-12) 75 76 zavm_evd(:,:,:) = avm(:,:,:) ! set avm prior to evd application 77 ! 74 78 DO jk = 1, jpkm1 75 79 #if defined key_vectopt_loop … … 99 103 CALL lbc_lnk( avmu, 'U', 1. ) ; CALL lbc_lnk( avmv, 'V', 1. ) 100 104 ! 105 zavm_evd(:,:,:) = avm(:,:,:) - zavm_evd(:,:,:) ! change in avm due to evd 106 CALL iom_put( "avm_evd", zavm_evd ) ! output this change 107 ! 101 108 CASE DEFAULT ! enhance vertical eddy diffusivity only (if rn2<-1.e-12) 102 109 DO jk = 1, jpkm1 … … 122 129 END SELECT 123 130 124 zav _evd(:,:,:) = avt(:,:,:) - zav_evd(:,:,:)! change in avt due to evd125 CALL iom_put( "avt_evd", zav _evd )! output this change131 zavt_evd(:,:,:) = avt(:,:,:) - zavt_evd(:,:,:) ! change in avt due to evd 132 CALL iom_put( "avt_evd", zavt_evd ) ! output this change 126 133 ! 127 134 END SUBROUTINE zdf_evd
Note: See TracChangeset
for help on using the changeset viewer.