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)

Author:

davestorkey

Message:

2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps_rewrite_time_filterswap : Update header comments etc
and rename dom_vvl_sf_swp -> dom_vvl_sf_update.

File:

• NEMO/branches/2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps_rewrite_time_filterswap/src/OCE/DYN/dynatf.F90 (modified) (3 diffs)

NEMO/branches/2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps_rewrite_time_filterswap/src/OCE/DYN/dynatf.F90

@@ -20 +20 @@
    !!            3.6  !  2014-04  (G. Madec) add the diagnostic of the time filter trends
    !!            3.7  !  2015-11  (J. Chanut) Free surface simplification
-   !!            4.1  !  2019-05  (D. Storkey) Rename dynnxt -> dynatf. Now just does time filtering.
+   !!            4.1  !  2019-08  (A. Coward, D. Storkey) Rename dynnxt.F90 -> dynatf.F90. Now just does time filtering.
    !!-------------------------------------------------------------------------
   
@@ -69 +69 @@
       !!                  ***  ROUTINE dyn_atf  ***
       !!                   
-      !!  !!!!!!!!!!!!!!!!!!  REWRITE HEADER COMMENTS !!!!!!!!!!!!!!!!!
-      !!
       !! ** Purpose :   Finalize after horizontal velocity. Apply the boundary 
-      !!             condition on the after velocity, achieve the time stepping 
-      !!             by applying the Asselin filter on now fields and swapping 
-      !!             the fields.
+      !!             condition on the after velocity and apply the Asselin time
+      !!             filter to the now fields.
       !!
       !! ** Method  : * Ensure after velocities transport matches time splitting
@@ -84 +81 @@
       !!             at the AGRIF zoom   boundaries (lk_agrif=T)
       !!
-      !!              * Apply the time filter applied and swap of the dynamics
+      !!              * Apply the Asselin time filter to the now fields
       !!             arrays to start the next time step:
-      !!                (puu(:,:,:,Kbb),pvv(:,:,:,Kbb)) = (puu(:,:,:,Kmm),pvv(:,:,:,Kmm)) + atfp [ (puu(:,:,:,Kbb),pvv(:,:,:,Kbb)) + (puu(:,:,:,Kaa),pvv(:,:,:,Kaa)) - 2 (puu(:,:,:,Kmm),pvv(:,:,:,Kmm)) ]
-      !!                (puu(:,:,:,Kmm),pvv(:,:,:,Kmm)) = (puu(:,:,:,Kaa),pvv(:,:,:,Kaa)).
+      !!                (puu(Kmm),pvv(Kmm)) = (puu(Kmm),pvv(Kmm)) 
+      !!                                    + atfp [ (puu(Kbb),pvv(Kbb)) + (puu(Kaa),pvv(Kaa)) - 2 (puu(Kmm),pvv(Kmm)) ]
       !!             Note that with flux form advection and non linear free surface,
       !!             the time filter is applied on thickness weighted velocity.
-      !!             As a result, dyn_atf MUST be called after tra_nxt.
+      !!             As a result, dyn_atf MUST be called after tra_atf.
       !!
-      !! ** Action :   uu(Kmm),vv(Kmm)   filtered now horizontal velocity 
+      !! ** Action :   puu(Kmm),pvv(Kmm)   filtered now horizontal velocity 
       !!----------------------------------------------------------------------
       INTEGER                             , INTENT(in   ) :: kt               ! ocean time-step index