Changeset 9250
- Timestamp:
- 2018-01-17T08:25:53+01:00 (7 years ago)
- Location:
- branches/2017/dev_merge_2017/NEMOGCM/NEMO/OPA_SRC
- Files:
-
- 5 edited
-
DYN/dynbfr.F90 (modified) (1 diff)
-
DYN/dynzdf.F90 (modified) (1 diff)
-
ZDF/zdfdrg.F90 (modified) (4 diffs)
-
step.F90 (modified) (2 diffs)
-
step_oce.F90 (modified) (1 diff)
Legend:
- Unmodified
- Added
- Removed
-
branches/2017/dev_merge_2017/NEMOGCM/NEMO/OPA_SRC/DYN/dynbfr.F90
r9124 r9250 62 62 ALLOCATE( ztrdu(jpi,jpj,jpk) , ztrdv(jpi,jpj,jpk) ) 63 63 ztrdu(:,:,:) = ua(:,:,:) 64 64 ztrdv(:,:,:) = va(:,:,:) 65 65 ENDIF 66 66 -
branches/2017/dev_merge_2017/NEMOGCM/NEMO/OPA_SRC/DYN/dynzdf.F90
r9019 r9250 91 91 ELSEIF( kt <= nit000 + 1 ) THEN ; r2dt = 2. * rdt ! = 2 rdt (leapfrog) 92 92 ENDIF 93 93 ! 94 ! !* explicit top/bottom drag case 95 IF( .NOT.ln_drgimp ) CALL zdf_drg_exp( kt, ub, vb, ua, va ) ! add top/bottom friction trend to (ua,va) 96 ! 97 ! 94 98 IF( l_trddyn ) THEN !* temporary save of ta and sa trends 95 99 ALLOCATE( ztrdu(jpi,jpj,jpk), ztrdv(jpi,jpj,jpk) ) -
branches/2017/dev_merge_2017/NEMOGCM/NEMO/OPA_SRC/ZDF/zdfdrg.F90
r9209 r9250 16 16 !!---------------------------------------------------------------------- 17 17 !! zdf_drg : update bottom friction coefficient (non-linear bottom friction only) 18 !! zdf_drg_exp : compute the top & bottom friction in explicit case 18 19 !! zdf_drg_init : read in namdrg namelist and control the bottom friction parameters. 19 20 !! drg_init : … … 23 24 USE dom_oce ! ocean space and time domain variables 24 25 USE zdf_oce ! ocean vertical physics variables 26 USE trd_oce ! trends: ocean variables 27 USE trddyn ! trend manager: dynamics 25 28 ! 26 29 USE in_out_manager ! I/O manager … … 34 37 35 38 PUBLIC zdf_drg ! called by zdf_phy 39 PUBLIC zdf_drg_exp ! called by dyn_zdf 36 40 PUBLIC zdf_drg_init ! called by zdf_phy_init 37 41 … … 140 144 141 145 146 SUBROUTINE zdf_drg_exp( kt, pub, pvb, pua, pva ) 147 !!---------------------------------------------------------------------- 148 !! *** ROUTINE zdf_drg_exp *** 149 !! 150 !! ** Purpose : compute and add the explicit top and bottom frictions. 151 !! 152 !! ** Method : in explicit case, 153 !! 154 !! NB: in implicit case the calculation is performed in dynzdf.F90 155 !! 156 !! ** Action : (pua,pva) momentum trend increased by top & bottom friction trend 157 !!--------------------------------------------------------------------- 158 INTEGER , INTENT(in ) :: kt ! ocean time-step index 159 REAL(wp), DIMENSION(jpi,jpj,jpk), INTENT(inout) :: pub, pvb ! the two components of the before velocity 160 REAL(wp), DIMENSION(jpi,jpj,jpk), INTENT(inout) :: pua, pva ! the two components of the velocity tendency 161 !! 162 INTEGER :: ji, jj ! dummy loop indexes 163 INTEGER :: ikbu, ikbv ! local integers 164 REAL(wp) :: zm1_2dt ! local scalar 165 REAL(wp) :: zCdu, zCdv ! - - 166 REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: ztrdu, ztrdv 167 !!--------------------------------------------------------------------- 168 ! 169 !!gm bug : time step is only rdt (not 2 rdt if euler start !) 170 zm1_2dt = - 1._wp / ( 2._wp * rdt ) 171 172 IF( l_trddyn ) THEN ! trends: store the input trends 173 ALLOCATE( ztrdu(jpi,jpj,jpk) , ztrdv(jpi,jpj,jpk) ) 174 ztrdu(:,:,:) = pua(:,:,:) 175 ztrdv(:,:,:) = pva(:,:,:) 176 ENDIF 177 178 DO jj = 2, jpjm1 179 DO ji = 2, jpim1 180 ikbu = mbku(ji,jj) ! deepest wet ocean u- & v-levels 181 ikbv = mbkv(ji,jj) 182 ! 183 ! Apply stability criteria on absolute value : abs(bfr/e3) < 1/(2dt) => bfr/e3 > -1/(2dt) 184 zCdu = 0.5*( rCdU_bot(ji+1,jj)+rCdU_bot(ji,jj) ) / e3u_n(ji,jj,ikbu) 185 zCdv = 0.5*( rCdU_bot(ji,jj+1)+rCdU_bot(ji,jj) ) / e3v_n(ji,jj,ikbv) 186 ! 187 pua(ji,jj,ikbu) = pua(ji,jj,ikbu) + MAX( zCdu , zm1_2dt ) * pub(ji,jj,ikbu) 188 pva(ji,jj,ikbv) = pva(ji,jj,ikbv) + MAX( zCdv , zm1_2dt ) * pvb(ji,jj,ikbv) 189 END DO 190 END DO 191 ! 192 IF( ln_isfcav ) THEN ! ocean cavities 193 DO jj = 2, jpjm1 194 DO ji = 2, jpim1 195 ikbu = miku(ji,jj) ! first wet ocean u- & v-levels 196 ikbv = mikv(ji,jj) 197 ! 198 ! Apply stability criteria on absolute value : abs(bfr/e3) < 1/(2dt) => bfr/e3 > -1/(2dt) 199 zCdu = 0.5*( rCdU_top(ji+1,jj)+rCdU_top(ji,jj) ) / e3u_n(ji,jj,ikbu) ! NB: Cdtop masked 200 zCdv = 0.5*( rCdU_top(ji,jj+1)+rCdU_top(ji,jj) ) / e3v_n(ji,jj,ikbv) 201 ! 202 pua(ji,jj,ikbu) = pua(ji,jj,ikbu) + MAX( zCdu , zm1_2dt ) * pub(ji,jj,ikbu) 203 pva(ji,jj,ikbv) = pva(ji,jj,ikbv) + MAX( zCdv , zm1_2dt ) * pvb(ji,jj,ikbv) 204 END DO 205 END DO 206 ENDIF 207 ! 208 IF( l_trddyn ) THEN ! trends: send trends to trddyn for further diagnostics 209 ztrdu(:,:,:) = pua(:,:,:) - ztrdu(:,:,:) 210 ztrdv(:,:,:) = pva(:,:,:) - ztrdv(:,:,:) 211 CALL trd_dyn( ztrdu(:,:,:), ztrdv(:,:,:), jpdyn_bfr, kt ) 212 DEALLOCATE( ztrdu, ztrdv ) 213 ENDIF 214 ! ! print mean trends (used for debugging) 215 IF(ln_ctl) CALL prt_ctl( tab3d_1=pua, clinfo1=' bfr - Ua: ', mask1=umask, & 216 & tab3d_2=pva, clinfo2= ' Va: ', mask2=vmask, clinfo3='dyn' ) 217 ! 218 END SUBROUTINE zdf_drg_exp 219 220 142 221 SUBROUTINE zdf_drg_init 143 222 !!---------------------------------------------------------------------- -
branches/2017/dev_merge_2017/NEMOGCM/NEMO/OPA_SRC/step.F90
r9210 r9250 199 199 CALL wzv ( kstp ) ! now cross-level velocity 200 200 ENDIF 201 202 IF( .NOT.ln_drgimp) CALL dyn_bfr ( kstp ) ! bottom friction203 201 204 202 CALL dyn_zdf ( kstp ) ! vertical diffusion … … 276 274 !!jc2: dynnxt must be the latest call. e3t_b are indeed updated in that routine 277 275 CALL tra_nxt ( kstp ) ! finalize (bcs) tracer fields at next time step and swap 278 CALL dyn_nxt ( kstp ) ! finalize (bcs) velocities at next time step and swap 276 CALL dyn_nxt ( kstp ) ! finalize (bcs) velocities at next time step and swap (always called after tra_nxt) 279 277 CALL ssh_swp ( kstp ) ! swap of sea surface height 280 278 IF(.NOT.ln_linssh) CALL dom_vvl_sf_swp( kstp ) ! swap of vertical scale factors -
branches/2017/dev_merge_2017/NEMOGCM/NEMO/OPA_SRC/step_oce.F90
r9213 r9250 37 37 USE divhor ! horizontal divergence (div_hor routine) 38 38 USE dynadv ! advection (dyn_adv routine) 39 USE dynbfr ! Bottom friction terms (dyn_bfr routine)40 39 USE dynvor ! vorticity term (dyn_vor routine) 41 40 USE dynhpg ! hydrostatic pressure grad. (dyn_hpg routine)
Note: See TracChangeset
for help on using the changeset viewer.
