Changeset 14053 for NEMO/trunk/src/NST/agrif_oce_update.F90

Timestamp:

2020-12-03T14:48:38+01:00 (3 years ago)

Author:

techene

Message:

#2385 added to the trunk

File:

• NEMO/trunk/src/NST/agrif_oce_update.F90 (modified) (9 diffs)

NEMO/trunk/src/NST/agrif_oce_update.F90

@@ -27 +27 @@
    USE vremap         ! Vertical remapping
    USE lbclnk 
-
+#if defined key_qco
+   USE domqco
+#endif
    IMPLICIT NONE
    PRIVATE
@@ -34 +36 @@
    PUBLIC   Update_Scales
 
+   !! * Substitutions
+#  include "domzgr_substitute.h90"
    !!----------------------------------------------------------------------
    !! NEMO/NST 4.0 , NEMO Consortium (2018)
@@ -191 +195 @@
    END SUBROUTINE Agrif_Update_Tke
 
-
    SUBROUTINE Agrif_Update_vvl( )
       !!---------------------------------------------
@@ -201 +204 @@
       IF (lwp.AND.lk_agrif_debug) Write(*,*) 'Update e3 from grid Number',Agrif_Fixed(), 'Step', Agrif_Nb_Step()
       !
+#if ! defined key_qco
       Agrif_UseSpecialValueInUpdate = .TRUE.
       Agrif_SpecialValueFineGrid = 0.
@@ -213 +217 @@
       CALL dom_vvl_update_UVF
       CALL Agrif_ParentGrid_To_ChildGrid()
+#else
+      CALL Agrif_ChildGrid_To_ParentGrid()
+      CALL Agrif_Update_qco
+      CALL Agrif_ParentGrid_To_ChildGrid()
+#endif
       !
    END SUBROUTINE Agrif_Update_vvl
 
+
+#if defined key_qco
+   SUBROUTINE Agrif_Update_qco
+      !!---------------------------------------------
+      !!       *** ROUTINE dom_Update_qco ***
+      !!---------------------------------------------
+      !
+      ! Save arrays prior update (needed for asselin correction)
+      r3t(:,:,Krhs_a) = r3t(:,:,Kmm_a)
+      r3u(:,:,Krhs_a) = r3u(:,:,Kmm_a)
+      r3v(:,:,Krhs_a) = r3v(:,:,Kmm_a)
+
+      ! Update r3x arrays from updated ssh
+      CALL dom_qco_zgr( Kbb_a, Kmm_a )
+      !
+   END SUBROUTINE Agrif_Update_qco
+#endif
+
+
+#if ! defined key_qco
    SUBROUTINE dom_vvl_update_UVF
       !!---------------------------------------------
@@ -224 +253 @@
       REAL(wp):: zcoef
       !!---------------------------------------------
-
       IF (lwp.AND.lk_agrif_debug) Write(*,*) 'Finalize e3 on grid Number', &
                   & Agrif_Fixed(), 'Step', Agrif_Nb_Step()
@@ -290 +318 @@
       !
    END SUBROUTINE dom_vvl_update_UVF
+#endif
 
 #if defined key_vertical
@@ -1332 +1361 @@
    END SUBROUTINE updateAVM
 
+#if ! defined key_qco
    SUBROUTINE updatee3t(ptab_dum, i1, i2, j1, j2, k1, k2, before )
       !!---------------------------------------------
@@ -1443 +1473 @@
       !
    END SUBROUTINE updatee3t
+#endif
 
 #else