New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 2486 for trunk – NEMO

Changeset 2486 for trunk


Ignore:
Timestamp:
2010-12-19T16:07:25+01:00 (13 years ago)
Author:
rblod
Message:

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

Location:
trunk/NEMO
Files:
4 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMO/NST_SRC/agrif_opa_interp.F90

    r1605 r2486  
    2929   PRIVATE 
    3030     
    31    PUBLIC   Agrif_tra, Agrif_dyn, interpu, interpv 
     31   PUBLIC   Agrif_tra, Agrif_dyn, Agrif_ssh, interpu, interpv 
    3232 
    3333#  include "domzgr_substitute.h90"   
     
    282282         END DO 
    283283 
    284          sshn(2,:)=sshn(3,:) 
    285          sshb(2,:)=sshb(3,:) 
    286  
    287284      ENDIF 
    288285 
     
    363360         END DO 
    364361 
    365          sshn(nlci-1,:)=sshn(nlci-2,:) 
    366          sshb(nlci-1,:)=sshb(nlci-2,:)         
    367362      ENDIF 
    368363 
     
    438433         END DO 
    439434 
    440          sshn(:,2)=sshn(:,3) 
    441          sshb(:,2)=sshb(:,3) 
    442435      ENDIF 
    443436 
     
    514507         END DO 
    515508 
    516          sshn(:,nlcj-1)=sshn(:,nlcj-2) 
    517          sshb(:,nlcj-1)=sshb(:,nlcj-2)                 
    518509      ENDIF 
    519510 
    520511   END SUBROUTINE Agrif_dyn 
     512 
     513 
     514   SUBROUTINE Agrif_ssh( kt ) 
     515      !!---------------------------------------------------------------------- 
     516      !!                  ***  ROUTINE Agrif_ssh  *** 
     517      !!----------------------------------------------------------------------   
     518      INTEGER, INTENT(in) ::   kt 
     519      !! 
     520      !!----------------------------------------------------------------------   
     521 
     522      IF( Agrif_Root() )   RETURN 
     523 
     524      IF((nbondi == -1).OR.(nbondi == 2)) THEN 
     525         ssha(2,:)=ssha(3,:) 
     526         sshn(2,:)=sshn(3,:) 
     527      ENDIF 
     528 
     529      IF((nbondi == 1).OR.(nbondi == 2)) THEN 
     530         ssha(nlci-1,:)=ssha(nlci-2,:) 
     531         sshn(nlci-1,:)=sshn(nlci-2,:)         
     532      ENDIF 
     533 
     534      IF((nbondj == -1).OR.(nbondj == 2)) THEN 
     535         ssha(:,2)=sshn(:,3) 
     536         sshn(:,2)=sshb(:,3) 
     537      ENDIF 
     538 
     539      IF((nbondj == 1).OR.(nbondj == 2)) THEN 
     540         ssha(:,nlcj-1)=ssha(:,nlcj-2) 
     541         ssha(:,nlcj-1)=sshn(:,nlcj-2)                 
     542      ENDIF 
     543 
     544   END SUBROUTINE Agrif_ssh 
    521545 
    522546 
  • trunk/NEMO/NST_SRC/agrif_opa_update.F90

    r1587 r2486  
    9393      !!   *** ROUTINE recompute_diags *** 
    9494      !!--------------------------------------------- 
    95       USE divcur 
    96       USE sshwzv 
    97       USE cla_div 
    98  
    9995      INTEGER, INTENT(in) :: kt 
    100  
    101       ta = hdivb 
    102       sa = rotb 
    103       Call div_cur(kt) 
    104  
    105       hdivb = ta 
    106       rotb  = sa 
    107  
    108       IF( n_cla == 1 ) CALL div_cla( kt ) 
    109 !      CALL wzv( kt ) 
    11096 
    11197   END SUBROUTINE recompute_diags 
  • trunk/NEMO/OPA_SRC/DYN/dynnxt.F90

    r1876 r2486  
    3232   USE bdydta          ! unstructured open boundary conditions 
    3333   USE bdydyn          ! unstructured open boundary conditions 
    34    USE agrif_opa_update 
    3534   USE agrif_opa_interp 
    3635   USE in_out_manager  ! I/O manager 
     
    259258      ENDIF 
    260259 
    261 #if defined key_agrif 
    262       ! Update velocity at AGRIF zoom boundaries 
    263       IF (.NOT.Agrif_Root())    CALL Agrif_Update_Dyn( kt ) 
    264 #endif       
    265  
    266260      IF(ln_ctl)   CALL prt_ctl( tab3d_1=un, clinfo1=' nxt  - Un: ', mask1=umask,   & 
    267261         &                       tab3d_2=vn, clinfo2=' Vn: '       , mask2=vmask ) 
  • trunk/NEMO/OPA_SRC/DYN/sshwzv.F90

    r1792 r2486  
    2727   USE diaar5, ONLY :   lk_diaar5 
    2828   USE iom 
     29   USE agrif_opa_interp 
     30   USE agrif_opa_update 
    2931 
    3032   IMPLICIT NONE 
     
    156158      ssha(:,:) = (  sshb(:,:) - z2dt * ( zraur * emp(:,:) + zhdiv(:,:) )  ) * tmask(:,:,1) 
    157159 
     160#if defined key_agrif 
     161      CALL agrif_ssh(kt) 
     162#endif 
    158163#if defined key_obc 
    159164      IF ( Agrif_Root() ) THEN  
     
    275280      ENDIF 
    276281      ! 
     282#if defined key_agrif 
     283      ! Update velocity at AGRIF zoom boundaries 
     284      IF (.NOT.Agrif_Root())    CALL Agrif_Update_Dyn( kt ) 
     285#endif 
     286 
    277287      IF(ln_ctl)   CALL prt_ctl(tab2d_1=sshb    , clinfo1=' sshb  - : ', mask1=tmask, ovlap=1 ) 
    278288      ! 
Note: See TracChangeset for help on using the changeset viewer.