Changeset 11475 for NEMO/branches/2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps_rewrite_time_filterswap/src/OCE/DYN/dynatf.F90
- Timestamp:
- 2019-08-28T09:22:35+02:00 (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps_rewrite_time_filterswap/src/OCE/DYN/dynatf.F90
r11050 r11475 20 20 !! 3.6 ! 2014-04 (G. Madec) add the diagnostic of the time filter trends 21 21 !! 3.7 ! 2015-11 (J. Chanut) Free surface simplification 22 !! 4.1 ! 2019-0 5 (D. Storkey) Rename dynnxt -> dynatf. Now just does time filtering.22 !! 4.1 ! 2019-08 (A. Coward, D. Storkey) Rename dynnxt.F90 -> dynatf.F90. Now just does time filtering. 23 23 !!------------------------------------------------------------------------- 24 24 … … 69 69 !! *** ROUTINE dyn_atf *** 70 70 !! 71 !! !!!!!!!!!!!!!!!!!! REWRITE HEADER COMMENTS !!!!!!!!!!!!!!!!!72 !!73 71 !! ** Purpose : Finalize after horizontal velocity. Apply the boundary 74 !! condition on the after velocity, achieve the time stepping 75 !! by applying the Asselin filter on now fields and swapping 76 !! the fields. 72 !! condition on the after velocity and apply the Asselin time 73 !! filter to the now fields. 77 74 !! 78 75 !! ** Method : * Ensure after velocities transport matches time splitting … … 84 81 !! at the AGRIF zoom boundaries (lk_agrif=T) 85 82 !! 86 !! * Apply the time filter applied and swap of the dynamics83 !! * Apply the Asselin time filter to the now fields 87 84 !! arrays to start the next time step: 88 !! (puu( :,:,:,Kbb),pvv(:,:,:,Kbb)) = (puu(:,:,:,Kmm),pvv(:,:,:,Kmm)) + atfp [ (puu(:,:,:,Kbb),pvv(:,:,:,Kbb)) + (puu(:,:,:,Kaa),pvv(:,:,:,Kaa)) - 2 (puu(:,:,:,Kmm),pvv(:,:,:,Kmm)) ]89 !! (puu(:,:,:,Kmm),pvv(:,:,:,Kmm)) = (puu(:,:,:,Kaa),pvv(:,:,:,Kaa)).85 !! (puu(Kmm),pvv(Kmm)) = (puu(Kmm),pvv(Kmm)) 86 !! + atfp [ (puu(Kbb),pvv(Kbb)) + (puu(Kaa),pvv(Kaa)) - 2 (puu(Kmm),pvv(Kmm)) ] 90 87 !! Note that with flux form advection and non linear free surface, 91 88 !! the time filter is applied on thickness weighted velocity. 92 !! As a result, dyn_atf MUST be called after tra_ nxt.89 !! As a result, dyn_atf MUST be called after tra_atf. 93 90 !! 94 !! ** Action : uu(Kmm),vv(Kmm) filtered now horizontal velocity91 !! ** Action : puu(Kmm),pvv(Kmm) filtered now horizontal velocity 95 92 !!---------------------------------------------------------------------- 96 93 INTEGER , INTENT(in ) :: kt ! ocean time-step index
Note: See TracChangeset
for help on using the changeset viewer.