Changeset 10877 for NEMO/branches/2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps/src/OCE/DYN/dynadv.F90
- Timestamp:
- 2019-04-16T11:50:59+02:00 (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps/src/OCE/DYN/dynadv.F90
r10874 r10877 53 53 CONTAINS 54 54 55 SUBROUTINE dyn_adv( kt )55 SUBROUTINE dyn_adv( kt, Kbb, Kmm, puu, pvv, Krhs ) 56 56 !!--------------------------------------------------------------------- 57 57 !! *** ROUTINE dyn_adv *** … … 59 59 !! ** Purpose : compute the ocean momentum advection trend. 60 60 !! 61 !! ** Method : - Update ( ua,va) with the advection term following n_dynadv61 !! ** Method : - Update (puu(:,:,:,Krhs),pvv(:,:,:,Krhs)) with the advection term following n_dynadv 62 62 !! 63 63 !! NB: in flux form advection (ln_dynadv_cen2 or ln_dynadv_ubs=T) … … 66 66 !! (see dynvor module). 67 67 !!---------------------------------------------------------------------- 68 INTEGER, INTENT( in ) :: kt ! ocean time-step index 68 INTEGER , INTENT( in ) :: kt ! ocean time-step index 69 INTEGER , INTENT( in ) :: Kbb, Kmm, Krhs ! ocean time level indices 70 REAL(wp), DIMENSION(jpi,jpj,jpk,jpt), INTENT(inout) :: puu, pvv ! ocean velocities and RHS of momentum equation 69 71 !!---------------------------------------------------------------------- 70 72 ! … … 73 75 SELECT CASE( n_dynadv ) !== compute advection trend and add it to general trend ==! 74 76 CASE( np_VEC_c2 ) 75 CALL dyn_keg ( kt, nn_dynkeg ) ! vector form : horizontal gradient of kinetic energy76 CALL dyn_zad ( kt )! vector form : vertical advection77 CALL dyn_keg ( kt, nn_dynkeg, Kmm, puu, pvv, Krhs ) ! vector form : horizontal gradient of kinetic energy 78 CALL dyn_zad ( kt, Kmm, puu, pvv, Krhs ) ! vector form : vertical advection 77 79 CASE( np_FLX_c2 ) 78 CALL dyn_adv_cen2( kt )! 2nd order centered scheme80 CALL dyn_adv_cen2( kt, Kmm, puu, pvv, Krhs ) ! 2nd order centered scheme 79 81 CASE( np_FLX_ubs ) 80 CALL dyn_adv_ubs ( kt )! 3rd order UBS scheme (UP3)82 CALL dyn_adv_ubs ( kt, Kbb, Kmm, puu, pvv, Krhs ) ! 3rd order UBS scheme (UP3) 81 83 END SELECT 82 84 !
Note: See TracChangeset
for help on using the changeset viewer.