[782] | 1 | MODULE agrif_oce |
---|
[1605] | 2 | !!====================================================================== |
---|
[782] | 3 | !! *** MODULE agrif_oce *** |
---|
[1605] | 4 | !! AGRIF : define in memory AGRIF variables |
---|
[782] | 5 | !!---------------------------------------------------------------------- |
---|
[1605] | 6 | !! History : 2.0 ! 2007-12 (R. Benshila) Original code |
---|
[782] | 7 | !!---------------------------------------------------------------------- |
---|
[1605] | 8 | #if defined key_agrif |
---|
[782] | 9 | !!---------------------------------------------------------------------- |
---|
[1605] | 10 | !! 'key_agrif' AGRIF zoom |
---|
| 11 | !!---------------------------------------------------------------------- |
---|
[782] | 12 | USE par_oce ! ocean parameters |
---|
| 13 | USE dom_oce ! domain parameters |
---|
[5656] | 14 | |
---|
[782] | 15 | IMPLICIT NONE |
---|
[2715] | 16 | PRIVATE |
---|
[782] | 17 | |
---|
[2715] | 18 | PUBLIC agrif_oce_alloc ! routine called by nemo_init in nemogcm.F90 |
---|
[12377] | 19 | |
---|
[1605] | 20 | ! !!* Namelist namagrif: AGRIF parameters |
---|
[13216] | 21 | LOGICAL , PUBLIC :: ln_init_chfrpar = .FALSE. !: set child grids initial state from parent |
---|
[12377] | 22 | LOGICAL , PUBLIC :: ln_agrif_2way = .TRUE. !: activate two way nesting |
---|
| 23 | LOGICAL , PUBLIC :: ln_spc_dyn = .FALSE. !: use zeros (.false.) or not (.true.) in |
---|
| 24 | !: bdys dynamical fields interpolation |
---|
[13351] | 25 | LOGICAL , PUBLIC :: ln_vert_remap = .FALSE. !: use vertical remapping |
---|
[13337] | 26 | REAL(wp), PUBLIC :: rn_sponge_tra = 0.002 !: sponge coeff. for tracers |
---|
| 27 | REAL(wp), PUBLIC :: rn_sponge_dyn = 0.002 !: sponge coeff. for dynamics |
---|
[12377] | 28 | REAL(wp), PUBLIC :: rn_trelax_tra = 0.01 !: time relaxation parameter for tracers |
---|
| 29 | REAL(wp), PUBLIC :: rn_trelax_dyn = 0.01 !: time relaxation parameter for momentum |
---|
[5656] | 30 | LOGICAL , PUBLIC :: ln_chk_bathy = .FALSE. !: check of parent bathymetry |
---|
[12377] | 31 | ! |
---|
| 32 | INTEGER , PUBLIC, PARAMETER :: nn_sponge_len = 2 !: Sponge width (in number of parent grid points) |
---|
[13216] | 33 | |
---|
[5656] | 34 | LOGICAL , PUBLIC :: spongedoneT = .FALSE. !: tracer sponge layer indicator |
---|
| 35 | LOGICAL , PUBLIC :: spongedoneU = .FALSE. !: dynamics sponge layer indicator |
---|
| 36 | LOGICAL , PUBLIC :: lk_agrif_fstep = .TRUE. !: if true: first step |
---|
| 37 | LOGICAL , PUBLIC :: lk_agrif_debug = .FALSE. !: if true: print debugging info |
---|
[1605] | 38 | |
---|
[5656] | 39 | LOGICAL , PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: tabspongedone_tsn |
---|
| 40 | # if defined key_top |
---|
| 41 | LOGICAL , PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: tabspongedone_trn |
---|
| 42 | # endif |
---|
| 43 | LOGICAL , PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: tabspongedone_u |
---|
| 44 | LOGICAL , PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: tabspongedone_v |
---|
[12377] | 45 | INTEGER , PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: utint_stage |
---|
| 46 | INTEGER , PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: vtint_stage |
---|
| 47 | REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: fspu, fspv !: sponge arrays |
---|
| 48 | REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: fspt, fspf !: " " |
---|
[2715] | 49 | |
---|
[9019] | 50 | ! Barotropic arrays used to store open boundary data during time-splitting loop: |
---|
[12377] | 51 | REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: ubdy, vbdy, hbdy |
---|
| 52 | INTEGER , PUBLIC, SAVE :: Kbb_a, Kmm_a, Krhs_a !: AGRIF module-specific copies of time-level indices |
---|
[5656] | 53 | |
---|
[12377] | 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 |
---|
[9019] | 56 | |
---|
[13351] | 57 | INTEGER, PUBLIC :: ts_interp_id, ts_update_id ! AGRIF profile for tracers interpolation and update |
---|
[9019] | 58 | INTEGER, PUBLIC :: un_interp_id, vn_interp_id ! AGRIF profiles for interpolations |
---|
| 59 | INTEGER, PUBLIC :: un_update_id, vn_update_id ! AGRIF profiles for udpates |
---|
[13351] | 60 | INTEGER, PUBLIC :: ts_sponge_id, un_sponge_id, vn_sponge_id ! AGRIF profiles for sponge layers |
---|
[13216] | 61 | INTEGER, PUBLIC :: tsini_id, uini_id, vini_id, sshini_id ! AGRIF profile for initialization |
---|
[5656] | 62 | # if defined key_top |
---|
[9019] | 63 | INTEGER, PUBLIC :: trn_id, trn_sponge_id |
---|
[5656] | 64 | # endif |
---|
[9019] | 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 |
---|
| 68 | INTEGER, PUBLIC :: scales_t_id |
---|
| 69 | INTEGER, PUBLIC :: avt_id, avm_id, en_id ! TKE related identificators |
---|
[12377] | 70 | INTEGER, PUBLIC :: mbkt_id, ht0_id |
---|
[13286] | 71 | INTEGER, PUBLIC :: glamt_id, gphit_id |
---|
[9019] | 72 | INTEGER, PUBLIC :: kindic_agr |
---|
[13216] | 73 | |
---|
| 74 | ! North fold |
---|
| 75 | !$AGRIF_DO_NOT_TREAT |
---|
| 76 | LOGICAL, PUBLIC :: use_sign_north |
---|
[13351] | 77 | REAL, PUBLIC :: sign_north |
---|
[13216] | 78 | LOGICAL, PUBLIC :: l_ini_child = .FALSE. |
---|
| 79 | LOGICAL, PUBLIC :: l_vremap = .FALSE. |
---|
| 80 | !$AGRIF_END_DO_NOT_TREAT |
---|
[9019] | 81 | |
---|
[1605] | 82 | !!---------------------------------------------------------------------- |
---|
[9598] | 83 | !! NEMO/NST 4.0 , NEMO Consortium (2018) |
---|
[1605] | 84 | !! $Id$ |
---|
[10068] | 85 | !! Software governed by the CeCILL license (see ./LICENSE) |
---|
[2715] | 86 | !!---------------------------------------------------------------------- |
---|
| 87 | CONTAINS |
---|
| 88 | |
---|
| 89 | INTEGER FUNCTION agrif_oce_alloc() |
---|
| 90 | !!---------------------------------------------------------------------- |
---|
| 91 | !! *** FUNCTION agrif_oce_alloc *** |
---|
| 92 | !!---------------------------------------------------------------------- |
---|
[5656] | 93 | INTEGER, DIMENSION(2) :: ierr |
---|
| 94 | !!---------------------------------------------------------------------- |
---|
| 95 | ierr(:) = 0 |
---|
| 96 | ! |
---|
[12377] | 97 | ALLOCATE( fspu(jpi,jpj), fspv(jpi,jpj), & |
---|
| 98 | & fspt(jpi,jpj), fspf(jpi,jpj), & |
---|
| 99 | & tabspongedone_tsn(jpi,jpj), & |
---|
| 100 | & utint_stage(jpi,jpj), vtint_stage(jpi,jpj), & |
---|
[5656] | 101 | # if defined key_top |
---|
| 102 | & tabspongedone_trn(jpi,jpj), & |
---|
[12377] | 103 | # endif |
---|
| 104 | & ht0_parent(jpi,jpj), mbkt_parent(jpi,jpj), & |
---|
| 105 | & hu0_parent(jpi,jpj), mbku_parent(jpi,jpj), & |
---|
| 106 | & hv0_parent(jpi,jpj), mbkv_parent(jpi,jpj), & |
---|
[5656] | 107 | & tabspongedone_u (jpi,jpj), & |
---|
| 108 | & tabspongedone_v (jpi,jpj), STAT = ierr(1) ) |
---|
| 109 | |
---|
[12377] | 110 | ALLOCATE( ubdy(jpi,jpj), vbdy(jpi,jpj), hbdy(jpi,jpj), STAT = ierr(2) ) |
---|
[5656] | 111 | |
---|
| 112 | agrif_oce_alloc = MAXVAL(ierr) |
---|
| 113 | ! |
---|
[2715] | 114 | END FUNCTION agrif_oce_alloc |
---|
| 115 | #endif |
---|
[1605] | 116 | !!====================================================================== |
---|
[782] | 117 | END MODULE agrif_oce |
---|