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 3917 for branches/2012/dev_v3_4_STABLE_2012/NEMOGCM/NEMO/NST_SRC/agrif_top_sponge.F90 – NEMO

Ignore:
Timestamp:
2013-06-12T23:04:27+02:00 (11 years ago)
Author:
smasson
Message:

v3_4_STABLE: (kind of) bugfix in agrif, see tricket #1111

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2012/dev_v3_4_STABLE_2012/NEMOGCM/NEMO/NST_SRC/agrif_top_sponge.F90

    r3294 r3917  
    3636      REAL(wp) :: ztra, zabe1, zabe2, zbtr 
    3737      REAL(wp), POINTER, DIMENSION(:,:) :: localviscsponge 
    38       REAL(wp), POINTER, DIMENSION(:,:,:,:) :: trbdiff, ztru, ztrv, ztab 
     38      REAL(wp), POINTER, DIMENSION(:,:,:,:) :: trbdiff, ztru, ztrv, ztab4d 
    3939 
    4040#if defined SPONGE_TOP 
    4141      CALL wrk_alloc( jpi, jpj, localviscsponge ) 
    42       CALL wrk_alloc( jpi, jpj, jpk, jptra, trbdiff, ztru, ztrv, ztab ) 
     42      CALL wrk_alloc( jpi, jpj, jpk, jptra, trbdiff, ztru, ztrv, ztab4d ) 
    4343 
    4444      timecoeff = REAL(Agrif_NbStepint(),wp)/Agrif_rhot() 
     
    4646      Agrif_SpecialValue=0. 
    4747      Agrif_UseSpecialValue = .TRUE. 
    48       ztab = 0.e0 
    49       CALL Agrif_Bc_Variable(ztab, tra_id,calledweight=timecoeff,procname=interptrn) 
     48      ztab4d = 0.e0 
     49      CALL Agrif_Bc_Variable(ztab4d, tra_id,calledweight=timecoeff,procname=interptrn) 
    5050      Agrif_UseSpecialValue = .FALSE. 
    5151 
    52       trbdiff(:,:,:,:) = trb(:,:,:,:) - ztab(:,:,:,:) 
     52      trbdiff(:,:,:,:) = trb(:,:,:,:) - ztab4d(:,:,:,:) 
    5353 
    5454      spongearea = 2 + 2 * Agrif_irhox() 
     
    147147  
    148148      CALL wrk_dealloc( jpi, jpj, localviscsponge ) 
    149       CALL wrk_dealloc( jpi, jpj, jpk, jptra, trbdiff, ztru, ztrv, ztab ) 
     149      CALL wrk_dealloc( jpi, jpj, jpk, jptra, trbdiff, ztru, ztrv, ztab4d ) 
    150150 
    151151#endif 
Note: See TracChangeset for help on using the changeset viewer.