Changeset 113 for trunk/NEMO/OPA_SRC/LDF

Timestamp:

2004-06-28T16:43:28+02:00 (20 years ago)

Author:

opalod

Message:

CT : UPDATE072 : ahtt is proportional to the cube of the maximum of the gridspacing in the to horizontal direction as in ldfdyn_c2d.h90

File:

• trunk/NEMO/OPA_SRC/LDF/ldftra_c2d.h90 (modified) (3 diffs)

trunk/NEMO/OPA_SRC/LDF/ldftra_c2d.h90

@@ -35 +35 @@
       INTEGER :: i1, i2, j1, j2
 # endif
+      REAL(wp) ::   za00, zdx_max
+      
       !!----------------------------------------------------------------------
 
@@ -49 +51 @@
       ENDIF
 
-
+      zdx_max = MAXVAL( e1t(:,:) )
+      IF( lk_mpp )   CALL mpp_max( zdx_max )   ! max over the global domain
+      
       ! biharmonic operator : (T-point)
       ! ====================
+      ! Here: ahtt is proportional to the cube of the maximum of the gridspacing
+      !       in the to horizontal direction
+      
+      za00 = aht0 / ( zdx_max * zdx_max * zdx_max )
 
-      ahtt(:,:) = aht0                ! set ahtt at T-point (here no space variation)
+      ahtt(:,:) = za00 * e1t(:,:) * e1t(:,:) *e1t(:,:)      ! set ahtt at T-point 
 
       CALL lbc_lnk( ahtt, 'T', 1. )   ! Lateral boundary conditions on ( ahtt )
@@ -69 +77 @@
       ! ==================
 
-      ahtu(:,:) = aht0                ! set ahtu = ahtv at u- and v-points,
-      ahtv(:,:) = aht0                ! and ahtw at w-point (idem T-point)
-      ahtw(:,:) = aht0                ! (here no space variation)
+      za00 = aht0 / zdx_max
+
+      ahtu(:,:) = za00 * e1u(:,:)               ! set ahtu = ahtv at u- and v-points,
+      ahtv(:,:) = za00 * e1f(:,:)               ! and ahtw at w-point (idem T-point)
+      ahtw(:,:) = za00 * e1t(:,:)               ! 
 
       CALL lbc_lnk( ahtu, 'U', 1. )   ! Lateral boundary conditions