Changeset 2486 for trunk/NEMO/NST_SRC


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

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

Location:
trunk/NEMO/NST_SRC
Files:
2 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 
Note: See TracChangeset for help on using the changeset viewer.