1 | MODULE agrif_oce |
---|
2 | !!====================================================================== |
---|
3 | !! *** MODULE agrif_oce *** |
---|
4 | !! AGRIF : define in memory AGRIF variables |
---|
5 | !!---------------------------------------------------------------------- |
---|
6 | !! History : 2.0 ! 2007-12 (R. Benshila) Original code |
---|
7 | !!---------------------------------------------------------------------- |
---|
8 | #if defined key_agrif |
---|
9 | !!---------------------------------------------------------------------- |
---|
10 | !! 'key_agrif' AGRIF zoom |
---|
11 | !!---------------------------------------------------------------------- |
---|
12 | USE par_oce ! ocean parameters |
---|
13 | USE dom_oce ! domain parameters |
---|
14 | |
---|
15 | IMPLICIT NONE |
---|
16 | PRIVATE |
---|
17 | |
---|
18 | PUBLIC agrif_oce_alloc ! routine called by nemo_init in nemogcm.F90 |
---|
19 | |
---|
20 | ! !!* Namelist namagrif: AGRIF parameters |
---|
21 | LOGICAL , PUBLIC :: ln_spc_dyn = .FALSE. !: |
---|
22 | INTEGER , PUBLIC :: nn_cln_update = 3 !: update frequency |
---|
23 | INTEGER , PUBLIC, PARAMETER :: nn_sponge_len = 2 !: Sponge width (in number of parent grid points) |
---|
24 | REAL(wp), PUBLIC :: rn_sponge_tra = 2800. !: sponge coeff. for tracers |
---|
25 | REAL(wp), PUBLIC :: rn_sponge_dyn = 2800. !: sponge coeff. for dynamics |
---|
26 | LOGICAL , PUBLIC :: ln_chk_bathy = .FALSE. !: check of parent bathymetry |
---|
27 | |
---|
28 | ! !!! OLD namelist names |
---|
29 | INTEGER , PUBLIC :: nbcline = 0 !: update counter |
---|
30 | INTEGER , PUBLIC :: nbclineupdate !: update frequency |
---|
31 | REAL(wp), PUBLIC :: visc_tra !: sponge coeff. for tracers |
---|
32 | REAL(wp), PUBLIC :: visc_dyn !: sponge coeff. for dynamics |
---|
33 | |
---|
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_doupd = .TRUE. !: if true: send update from current grid |
---|
38 | LOGICAL , PUBLIC :: lk_agrif_debug = .FALSE. !: if true: print debugging info |
---|
39 | |
---|
40 | LOGICAL , PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: tabspongedone_tsn |
---|
41 | # if defined key_top |
---|
42 | LOGICAL , PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: tabspongedone_trn |
---|
43 | # endif |
---|
44 | LOGICAL , PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: tabspongedone_u |
---|
45 | LOGICAL , PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: tabspongedone_v |
---|
46 | REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: fsaht_spu, fsaht_spv !: sponge diffusivities |
---|
47 | REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: fsahm_spt, fsahm_spf !: sponge viscosities |
---|
48 | |
---|
49 | ! Barotropic arrays used to store open boundary data during |
---|
50 | ! time-splitting loop: |
---|
51 | REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:) :: ubdy_w, vbdy_w, hbdy_w |
---|
52 | REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:) :: ubdy_e, vbdy_e, hbdy_e |
---|
53 | REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:) :: ubdy_n, vbdy_n, hbdy_n |
---|
54 | REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:) :: ubdy_s, vbdy_s, hbdy_s |
---|
55 | |
---|
56 | INTEGER :: tsn_id ! AGRIF profile for tracers interpolation and update |
---|
57 | INTEGER :: un_interp_id, vn_interp_id ! AGRIF profiles for interpolations |
---|
58 | INTEGER :: un_update_id, vn_update_id ! AGRIF profiles for udpates |
---|
59 | INTEGER :: tsn_sponge_id, un_sponge_id, vn_sponge_id ! AGRIF profiles for sponge layers |
---|
60 | # if defined key_top |
---|
61 | INTEGER :: trn_id, trn_sponge_id |
---|
62 | # endif |
---|
63 | INTEGER :: unb_id, vnb_id, ub2b_interp_id, vb2b_interp_id |
---|
64 | INTEGER :: ub2b_update_id, vb2b_update_id |
---|
65 | INTEGER :: e3t_id, e1u_id, e2v_id, sshn_id |
---|
66 | INTEGER :: scales_t_id |
---|
67 | # if defined key_zdftke |
---|
68 | INTEGER :: avt_id, avm_id, en_id |
---|
69 | # endif |
---|
70 | INTEGER :: umsk_id, vmsk_id |
---|
71 | INTEGER :: kindic_agr |
---|
72 | |
---|
73 | !!---------------------------------------------------------------------- |
---|
74 | !! NEMO/NST 3.3.1 , NEMO Consortium (2011) |
---|
75 | !! $Id$ |
---|
76 | !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) |
---|
77 | !!---------------------------------------------------------------------- |
---|
78 | CONTAINS |
---|
79 | |
---|
80 | INTEGER FUNCTION agrif_oce_alloc() |
---|
81 | !!---------------------------------------------------------------------- |
---|
82 | !! *** FUNCTION agrif_oce_alloc *** |
---|
83 | !!---------------------------------------------------------------------- |
---|
84 | INTEGER, DIMENSION(2) :: ierr |
---|
85 | !!---------------------------------------------------------------------- |
---|
86 | ierr(:) = 0 |
---|
87 | ! |
---|
88 | ALLOCATE( fsaht_spu(jpi,jpj), fsaht_spv(jpi,jpj), & |
---|
89 | & fsahm_spt(jpi,jpj), fsahm_spf(jpi,jpj), & |
---|
90 | & tabspongedone_tsn(jpi,jpj), & |
---|
91 | # if defined key_top |
---|
92 | & tabspongedone_trn(jpi,jpj), & |
---|
93 | # endif |
---|
94 | & tabspongedone_u (jpi,jpj), & |
---|
95 | & tabspongedone_v (jpi,jpj), STAT = ierr(1) ) |
---|
96 | |
---|
97 | ALLOCATE( ubdy_w(jpj), vbdy_w(jpj), hbdy_w(jpj), & |
---|
98 | & ubdy_e(jpj), vbdy_e(jpj), hbdy_e(jpj), & |
---|
99 | & ubdy_n(jpi), vbdy_n(jpi), hbdy_n(jpi), & |
---|
100 | & ubdy_s(jpi), vbdy_s(jpi), hbdy_s(jpi), STAT = ierr(2) ) |
---|
101 | |
---|
102 | agrif_oce_alloc = MAXVAL(ierr) |
---|
103 | ! |
---|
104 | END FUNCTION agrif_oce_alloc |
---|
105 | |
---|
106 | #endif |
---|
107 | !!====================================================================== |
---|
108 | END MODULE agrif_oce |
---|