Ignore:
Timestamp:
2019-07-15T15:20:49+02:00 (15 months ago)
Author:
smasson
Message:

dev_r10984_HPC-13 : introduce a logical to force vertical interpolation if the number of vertical levels in bdy dta == jpk, see #2285

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2019/dev_r10984_HPC-13_IRRMANN_BDY_optimization/src/OCE/BDY/bdydta.F90

    r11224 r11268  
    331331      LOGICAL                                ::   ln_full_vel   ! =T => full velocities in 3D boundary data 
    332332      !                                                         ! =F => baroclinic velocities in 3D boundary data 
    333       INTEGER                                ::   ipk,ipl           ! 
     333      LOGICAL                                ::   ln_zinterp    ! =T => requires a vertical interpolation of the bdydta 
     334      INTEGER                                ::   ipk,ipl       ! 
    334335      INTEGER                                ::   idvar         ! variable ID 
    335336      INTEGER                                ::   indims        ! number of dimensions of the variable 
     
    348349      NAMELIST/nambdy_dta/ cn_dir, bn_tem, bn_sal, bn_u3d, bn_v3d, bn_ssh, bn_u2d, bn_v2d  
    349350      NAMELIST/nambdy_dta/ bn_a_i, bn_h_i, bn_h_s 
    350       NAMELIST/nambdy_dta/ ln_full_vel 
     351      NAMELIST/nambdy_dta/ ln_full_vel, ln_zinterp 
    351352      !!--------------------------------------------------------------------------- 
    352353      ! 
     
    361362      ENDIF 
    362363      bf(:,:)%clrootname = 'NOT USED'   ! default definition used as a flag in fld_read to do nothing. 
    363       bf(:,:)%ltotvel = .FALSE.         ! default definition 
     364      bf(:,:)%lzint      = .FALSE.      ! default definition 
     365      bf(:,:)%ltotvel    = .FALSE.      ! default definition 
    364366  
    365367      ! Read namelists 
     
    521523                  CALL fld_fill( bf_alias, bn_alias, cn_dir, 'bdy_dta', cl3//' '//ctmp1, ctmp2 )   ! use namelist info 
    522524                  IF( bf_alias(1)%ln_tint ) ALLOCATE( bf_alias(1)%fdta( iszdim, 1, ipk, 2 ) ) 
    523                   bf_alias(1)%imap => idx_bdy(jbdy)%nbmap(1:iszdim,igrd)      ! associate the mapping used for this bdy 
    524                   bf_alias(1)%igrd = igrd                                     ! used only for vertical integration of 3D arrays 
     525                  bf_alias(1)%imap    => idx_bdy(jbdy)%nbmap(1:iszdim,igrd)   ! associate the mapping used for this bdy 
     526                  bf_alias(1)%igrd    = igrd                                  ! used only for vertical integration of 3D arrays 
    525527                  bf_alias(1)%ltotvel = ln_full_vel                           ! T if u3d is full velocity 
     528                  bf_alias(1)%lzint   = ln_zinterp                            ! T if it requires a vertical interpolation 
    526529               ENDIF 
    527530 
Note: See TracChangeset for help on using the changeset viewer.