Changeset 2486

Timestamp:

2010-12-19T16:07:25+01:00 (14 years ago)

Author:

rblod

Message:

Correct Agrif inconstency for ssh, nemo_v3_2 version, see ticket 669

Location:

trunk/NEMO

Files:

• NST_SRC/agrif_opa_interp.F90 (modified) (5 diffs)
• NST_SRC/agrif_opa_update.F90 (modified) (1 diff)
• OPA_SRC/DYN/dynnxt.F90 (modified) (2 diffs)
• OPA_SRC/DYN/sshwzv.F90 (modified) (3 diffs)

trunk/NEMO/NST_SRC/agrif_opa_interp.F90

@@ -29 +29 @@
    PRIVATE
     
-   PUBLIC   Agrif_tra, Agrif_dyn, interpu, interpv
+   PUBLIC   Agrif_tra, Agrif_dyn, Agrif_ssh, interpu, interpv
 
 #  include "domzgr_substitute.h90"  
@@ -282 +282 @@
          END DO
 
-         sshn(2,:)=sshn(3,:)
-         sshb(2,:)=sshb(3,:)
-
       ENDIF
 
@@ -363 +360 @@
          END DO
 
-         sshn(nlci-1,:)=sshn(nlci-2,:)
-         sshb(nlci-1,:)=sshb(nlci-2,:)        
       ENDIF
 
@@ -438 +433 @@
          END DO
 
-         sshn(:,2)=sshn(:,3)
-         sshb(:,2)=sshb(:,3)
       ENDIF
 
@@ -514 +507 @@
          END DO
 
-         sshn(:,nlcj-1)=sshn(:,nlcj-2)
-         sshb(:,nlcj-1)=sshb(:,nlcj-2)                
       ENDIF
 
    END SUBROUTINE Agrif_dyn
+
+
+   SUBROUTINE Agrif_ssh( kt )
+      !!----------------------------------------------------------------------
+      !!                  ***  ROUTINE Agrif_ssh  ***
+      !!----------------------------------------------------------------------  
+      INTEGER, INTENT(in) ::   kt
+      !!
+      !!----------------------------------------------------------------------  
+
+      IF( Agrif_Root() )   RETURN
+
+      IF((nbondi == -1).OR.(nbondi == 2)) THEN
+         ssha(2,:)=ssha(3,:)
+         sshn(2,:)=sshn(3,:)
+      ENDIF
+
+      IF((nbondi == 1).OR.(nbondi == 2)) THEN
+         ssha(nlci-1,:)=ssha(nlci-2,:)
+         sshn(nlci-1,:)=sshn(nlci-2,:)        
+      ENDIF
+
+      IF((nbondj == -1).OR.(nbondj == 2)) THEN
+         ssha(:,2)=sshn(:,3)
+         sshn(:,2)=sshb(:,3)
+      ENDIF
+
+      IF((nbondj == 1).OR.(nbondj == 2)) THEN
+         ssha(:,nlcj-1)=ssha(:,nlcj-2)
+         ssha(:,nlcj-1)=sshn(:,nlcj-2)                
+      ENDIF
+
+   END SUBROUTINE Agrif_ssh
 
 

trunk/NEMO/NST_SRC/agrif_opa_update.F90

@@ -93 +93 @@
       !!   *** ROUTINE recompute_diags ***
       !!---------------------------------------------
-      USE divcur
-      USE sshwzv
-      USE cla_div
-
       INTEGER, INTENT(in) :: kt
-
-      ta = hdivb
-      sa = rotb
-      Call div_cur(kt)
-
-      hdivb = ta
-      rotb  = sa
-
-      IF( n_cla == 1 ) CALL div_cla( kt )
-!      CALL wzv( kt )
 
    END SUBROUTINE recompute_diags

trunk/NEMO/OPA_SRC/DYN/dynnxt.F90

@@ -32 +32 @@
    USE bdydta          ! unstructured open boundary conditions
    USE bdydyn          ! unstructured open boundary conditions
-   USE agrif_opa_update
    USE agrif_opa_interp
    USE in_out_manager  ! I/O manager
@@ -259 +258 @@
       ENDIF
 
-#if defined key_agrif
-      ! Update velocity at AGRIF zoom boundaries
-      IF (.NOT.Agrif_Root())    CALL Agrif_Update_Dyn( kt )
-#endif      
-
       IF(ln_ctl)   CALL prt_ctl( tab3d_1=un, clinfo1=' nxt  - Un: ', mask1=umask,   &
          &                       tab3d_2=vn, clinfo2=' Vn: '       , mask2=vmask )

trunk/NEMO/OPA_SRC/DYN/sshwzv.F90

@@ -27 +27 @@
    USE diaar5, ONLY :   lk_diaar5
    USE iom
+   USE agrif_opa_interp
+   USE agrif_opa_update
 
    IMPLICIT NONE
@@ -156 +158 @@
       ssha(:,:) = (  sshb(:,:) - z2dt * ( zraur * emp(:,:) + zhdiv(:,:) )  ) * tmask(:,:,1)
 
+#if defined key_agrif
+      CALL agrif_ssh(kt)
+#endif
 #if defined key_obc
       IF ( Agrif_Root() ) THEN 
@@ -275 +280 @@
       ENDIF
       !
+#if defined key_agrif
+      ! Update velocity at AGRIF zoom boundaries
+      IF (.NOT.Agrif_Root())    CALL Agrif_Update_Dyn( kt )
+#endif
+
       IF(ln_ctl)   CALL prt_ctl(tab2d_1=sshb    , clinfo1=' sshb  - : ', mask1=tmask, ovlap=1 )
       !