Ignore:
Timestamp:
2020-12-01T13:44:17+01:00 (5 months ago)
Author:
jchanut
Message:

#2222, added:

  • Divergence conserving interpolation for transport (Balsara's scheme)
  • Possibility to shift barotropic interface inside the domain (currently disabled - needed to match volume within sponge)
  • PPR library as default vertical interpolation scheme
  • Trilinear interpolation over partial bottom cells
File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2020/dev_r13312_AGRIF-03-04_jchanut_vinterp_tstep/src/NST/agrif_oce.F90

    r13565 r13937  
    3131   ! 
    3232   INTEGER , PUBLIC, PARAMETER ::   nn_sponge_len = 2  !: Sponge width (in number of parent grid points) 
     33   INTEGER , PUBLIC, PARAMETER ::   nn_shift_bar = 0   !: nb of coarse grid points by which we shift 2d interface 
    3334 
    3435   LOGICAL , PUBLIC :: spongedoneT = .FALSE.       !: tracer   sponge layer indicator 
     
    3637   LOGICAL , PUBLIC :: lk_agrif_fstep = .TRUE.     !: if true: first step 
    3738   LOGICAL , PUBLIC :: lk_agrif_debug = .FALSE.    !: if true: print debugging info 
    38  
     39   LOGICAL , PUBLIC :: lk_tint2d_notinterp = .FALSE. !: if true, no time interp 
    3940   LOGICAL , PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: tabspongedone_tsn 
    4041# if defined key_top 
     
    4748   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: fspu, fspv !: sponge arrays 
    4849   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: fspt, fspf !:   "      " 
     50   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: fspu_2d,fspv_2d  !: sponge arrays (2d mode) 
     51   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: fspt_2d, fspf_2d !:   "       "     "   " 
    4952 
    5053   ! Barotropic arrays used to store open boundary data during time-splitting loop: 
     
    5255   INTEGER , PUBLIC,              SAVE                 ::  Kbb_a, Kmm_a, Krhs_a   !: AGRIF module-specific copies of time-level indices 
    5356 
    54    REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: ht0_parent, hu0_parent, hv0_parent 
    55    INTEGER,  PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: mbkt_parent, mbku_parent, mbkv_parent 
     57   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:)   :: ht0_parent, hu0_parent, hv0_parent 
     58   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: e3t0_parent, e3u0_parent, e3v0_parent 
     59   INTEGER,  PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:)   :: mbkt_parent, mbku_parent, mbkv_parent 
     60 
    5661 
    5762   INTEGER, PUBLIC :: ts_interp_id, ts_update_id                              ! AGRIF profile for tracers interpolation and update 
    5863   INTEGER, PUBLIC :: un_interp_id, vn_interp_id                              ! AGRIF profiles for interpolations 
    5964   INTEGER, PUBLIC :: un_update_id, vn_update_id                              ! AGRIF profiles for udpates 
    60    INTEGER, PUBLIC :: ts_sponge_id, un_sponge_id, vn_sponge_id                ! AGRIF profiles for sponge layers 
     65   INTEGER, PUBLIC :: ts_sponge_id, un_sponge_id, vn_sponge_id                ! AGRIF profiles for sponge layers (3d) 
     66   INTEGER, PUBLIC :: unb_sponge_id, vnb_sponge_id                            ! AGRIF profiles for sponge layers (2d) 
    6167   INTEGER, PUBLIC :: tsini_id, uini_id, vini_id, sshini_id                   ! AGRIF profile for initialization 
    6268# if defined key_top 
    6369   INTEGER, PUBLIC :: trn_id, trn_sponge_id 
    6470# endif   
    65    INTEGER, PUBLIC :: unb_id, vnb_id, ub2b_interp_id, vb2b_interp_id 
    66    INTEGER, PUBLIC :: ub2b_update_id, vb2b_update_id 
    67    INTEGER, PUBLIC :: e3t_id, e1u_id, e2v_id, sshn_id 
     71   INTEGER, PUBLIC :: unb_interp_id, vnb_interp_id, ub2b_interp_id, vb2b_interp_id 
     72   INTEGER, PUBLIC :: ub2b_update_id, vb2b_update_id, unb_update_id, vnb_update_id 
     73   INTEGER, PUBLIC :: ub2b_cor_id, vb2b_cor_id 
     74   INTEGER, PUBLIC :: e3t_id, sshn_id 
    6875   INTEGER, PUBLIC :: scales_t_id 
    6976   INTEGER, PUBLIC :: avt_id, avm_id, en_id                ! TKE related identificators 
    70    INTEGER, PUBLIC :: mbkt_id, ht0_id 
     77   INTEGER, PUBLIC :: mbkt_id, ht0_id, e3t0_interp_id 
    7178   INTEGER, PUBLIC :: glamt_id, gphit_id 
    7279   INTEGER, PUBLIC :: batupd_id 
     
    98105      ALLOCATE( fspu(jpi,jpj), fspv(jpi,jpj),          & 
    99106         &      fspt(jpi,jpj), fspf(jpi,jpj),               & 
     107         &      fspu_2d(jpi,jpj), fspv_2d(jpi,jpj),         & 
     108         &      fspt_2d(jpi,jpj), fspf_2d(jpi,jpj),         & 
    100109         &      tabspongedone_tsn(jpi,jpj),                 & 
    101110         &      utint_stage(jpi,jpj), vtint_stage(jpi,jpj), & 
Note: See TracChangeset for help on using the changeset viewer.