1 | MODULE oce_trc |
---|
2 | !!====================================================================== |
---|
3 | !! *** MODULE oce_trc *** |
---|
4 | !! TOP : variables shared between ocean and passive tracers |
---|
5 | !!====================================================================== |
---|
6 | !! History : 1.0 ! 2004-03 (C. Ethe) original code |
---|
7 | !! 2.0 ! 2007-12 (C. Ethe, G. Madec) rewritting |
---|
8 | !!---------------------------------------------------------------------- |
---|
9 | #if defined key_top |
---|
10 | !!---------------------------------------------------------------------- |
---|
11 | !! 'key_top' TOP models |
---|
12 | !!---------------------------------------------------------------------- |
---|
13 | |
---|
14 | !* Domain size * |
---|
15 | USE par_oce , ONLY : jpi => jpi !: first dimension of grid --> i |
---|
16 | USE par_oce , ONLY : jpj => jpj !: second dimension of grid --> j |
---|
17 | USE par_oce , ONLY : jpk => jpk !: number of levels |
---|
18 | USE par_oce , ONLY : jpim1 => jpim1 !: jpi - 1 |
---|
19 | USE par_oce , ONLY : jpjm1 => jpjm1 !: jpj - 1 |
---|
20 | USE par_oce , ONLY : jpkm1 => jpkm1 !: jpk - 1 |
---|
21 | USE par_oce , ONLY : jpij => jpij !: jpi x jpj |
---|
22 | USE par_oce , ONLY : lk_esopa => lk_esopa !: flag to activate the all option |
---|
23 | USE par_oce , ONLY : jp_tem => jp_tem !: indice for temperature |
---|
24 | USE par_oce , ONLY : jp_sal => jp_sal !: indice for salinity |
---|
25 | |
---|
26 | !* IO manager * |
---|
27 | USE in_out_manager |
---|
28 | |
---|
29 | !* Memory Allocation * |
---|
30 | USE wrk_nemo |
---|
31 | |
---|
32 | !* Timing * |
---|
33 | USE timing |
---|
34 | |
---|
35 | !* MPP library |
---|
36 | USE lib_mpp |
---|
37 | |
---|
38 | !* Fortran utilities |
---|
39 | USE lib_fortran |
---|
40 | |
---|
41 | !* Lateral boundary conditions |
---|
42 | USE lbclnk |
---|
43 | |
---|
44 | !* physical constants * |
---|
45 | USE phycst |
---|
46 | |
---|
47 | !* 1D configuration |
---|
48 | USE c1d |
---|
49 | |
---|
50 | !* model domain * |
---|
51 | USE dom_oce |
---|
52 | |
---|
53 | USE domvvl, ONLY : un_td, vn_td !: thickness diffusion transport |
---|
54 | USE domvvl, ONLY : ln_vvl_ztilde !: ztilde vertical coordinate |
---|
55 | USE domvvl, ONLY : ln_vvl_layer !: level vertical coordinate |
---|
56 | |
---|
57 | !* ocean fields: here now and after fields * |
---|
58 | USE oce , ONLY : ua => ua !: i-horizontal velocity (m s-1) |
---|
59 | USE oce , ONLY : va => va !: j-horizontal velocity (m s-1) |
---|
60 | USE oce , ONLY : un => un !: i-horizontal velocity (m s-1) |
---|
61 | USE oce , ONLY : vn => vn !: j-horizontal velocity (m s-1) |
---|
62 | USE oce , ONLY : wn => wn !: vertical velocity (m s-1) |
---|
63 | USE oce , ONLY : tsn => tsn !: 4D array contaning ( tn, sn ) |
---|
64 | USE oce , ONLY : tsb => tsb !: 4D array contaning ( tb, sb ) |
---|
65 | USE oce , ONLY : tsa => tsa !: 4D array contaning ( ta, sa ) |
---|
66 | USE oce , ONLY : rhop => rhop !: potential volumic mass (kg m-3) |
---|
67 | USE oce , ONLY : rhd => rhd !: in situ density anomalie rhd=(rho-rau0)/rau0 (no units) |
---|
68 | #if defined key_offline |
---|
69 | USE oce , ONLY : rab_n => rab_n !: local thermal/haline expension ratio at T-points |
---|
70 | #endif |
---|
71 | USE oce , ONLY : hdivn => hdivn !: horizontal divergence (1/s) |
---|
72 | USE oce , ONLY : rotn => rotn !: relative vorticity [s-1] |
---|
73 | USE oce , ONLY : hdivb => hdivb !: horizontal divergence (1/s) |
---|
74 | USE oce , ONLY : rotb => rotb !: relative vorticity [s-1] |
---|
75 | USE oce , ONLY : sshn => sshn !: sea surface height at t-point [m] |
---|
76 | USE oce , ONLY : sshb => sshb !: sea surface height at t-point [m] |
---|
77 | USE oce , ONLY : ssha => ssha !: sea surface height at t-point [m] |
---|
78 | USE oce , ONLY : l_traldf_rot => l_traldf_rot !: rotated laplacian operator for lateral diffusion |
---|
79 | |
---|
80 | !* surface fluxes * |
---|
81 | USE sbc_oce , ONLY : utau => utau !: i-surface stress component |
---|
82 | USE sbc_oce , ONLY : vtau => vtau !: j-surface stress component |
---|
83 | USE sbc_oce , ONLY : wndm => wndm !: 10m wind speed |
---|
84 | USE sbc_oce , ONLY : qsr => qsr !: penetrative solar radiation (w m-2) |
---|
85 | USE sbc_oce , ONLY : emp => emp !: freshwater budget: volume flux [Kg/m2/s] |
---|
86 | USE sbc_oce , ONLY : emp_b => emp_b !: freshwater budget: volume flux [Kg/m2/s] |
---|
87 | USE sbc_oce , ONLY : fmmflx => fmmflx !: freshwater budget: volume flux [Kg/m2/s] |
---|
88 | USE sbc_oce , ONLY : rnf => rnf !: river runoff [Kg/m2/s] |
---|
89 | USE sbc_oce , ONLY : ln_dm2dc => ln_dm2dc !: Daily mean to Diurnal Cycle short wave (qsr) |
---|
90 | USE sbc_oce , ONLY : ln_rnf => ln_rnf !: runoffs / runoff mouths |
---|
91 | USE sbc_oce , ONLY : fr_i => fr_i !: ice fraction (between 0 to 1) |
---|
92 | USE traqsr , ONLY : rn_abs => rn_abs !: fraction absorbed in the very near surface |
---|
93 | USE traqsr , ONLY : rn_si0 => rn_si0 !: very near surface depth of extinction |
---|
94 | USE traqsr , ONLY : ln_qsr_bio => ln_qsr_bio !: flag to use or not the biological fluxes for light |
---|
95 | USE sbcrnf , ONLY : rnfmsk => rnfmsk !: mixed adv scheme in runoffs vicinity (hori.) |
---|
96 | USE sbcrnf , ONLY : rnfmsk_z => rnfmsk_z !: mixed adv scheme in runoffs vicinity (vert.) |
---|
97 | USE sbcrnf , ONLY : h_rnf => h_rnf !: river runoff [Kg/m2/s] |
---|
98 | |
---|
99 | USE trc_oce |
---|
100 | |
---|
101 | !* lateral diffusivity (tracers) * |
---|
102 | USE ldftra_oce , ONLY : rldf => rldf !: multiplicative coef. for lateral diffusivity |
---|
103 | USE ldftra_oce , ONLY : rn_aht_0 => rn_aht_0 !: horizontal eddy diffusivity for tracers (m2/s) |
---|
104 | USE ldftra_oce , ONLY : aht0 => aht0 !: horizontal eddy diffusivity for tracers (m2/s) |
---|
105 | USE ldftra_oce , ONLY : ahtb0 => ahtb0 !: background eddy diffusivity for isopycnal diff. (m2/s) |
---|
106 | USE ldftra_oce , ONLY : ahtu => ahtu !: lateral diffusivity coef. at u-points |
---|
107 | USE ldftra_oce , ONLY : ahtv => ahtv !: lateral diffusivity coef. at v-points |
---|
108 | USE ldftra_oce , ONLY : ahtw => ahtw !: lateral diffusivity coef. at w-points |
---|
109 | USE ldftra_oce , ONLY : ahtt => ahtt !: lateral diffusivity coef. at t-points |
---|
110 | USE ldftra_oce , ONLY : aeiv0 => aeiv0 !: eddy induced velocity coefficient (m2/s) |
---|
111 | USE ldftra_oce , ONLY : aeiu => aeiu !: eddy induced velocity coef. at u-points (m2/s) |
---|
112 | USE ldftra_oce , ONLY : aeiv => aeiv !: eddy induced velocity coef. at v-points (m2/s) |
---|
113 | USE ldftra_oce , ONLY : aeiw => aeiw !: eddy induced velocity coef. at w-points (m2/s) |
---|
114 | USE ldftra_oce , ONLY : lk_traldf_eiv => lk_traldf_eiv !: eddy induced velocity flag |
---|
115 | |
---|
116 | !* vertical diffusion * |
---|
117 | USE zdf_oce , ONLY : avt => avt !: vert. diffusivity coef. at w-point for temp |
---|
118 | # if defined key_zdfddm |
---|
119 | USE zdfddm , ONLY : avs => avs !: salinity vertical diffusivity coeff. at w-point |
---|
120 | # endif |
---|
121 | |
---|
122 | !* mixing & mixed layer depth * |
---|
123 | USE zdfmxl , ONLY : nmln => nmln !: number of level in the mixed layer |
---|
124 | USE zdfmxl , ONLY : hmld => hmld !: mixing layer depth (turbocline) |
---|
125 | USE zdfmxl , ONLY : hmlp => hmlp !: mixed layer depth (rho=rho0+zdcrit) (m) |
---|
126 | USE zdfmxl , ONLY : hmlpt => hmlpt !: mixed layer depth at t-points (m) |
---|
127 | |
---|
128 | !* direction of lateral diffusion * |
---|
129 | USE ldfslp , ONLY : lk_ldfslp => lk_ldfslp !: slopes flag |
---|
130 | # if defined key_ldfslp |
---|
131 | USE ldfslp , ONLY : uslp => uslp !: i-direction slope at u-, w-points |
---|
132 | USE ldfslp , ONLY : vslp => vslp !: j-direction slope at v-, w-points |
---|
133 | USE ldfslp , ONLY : wslpi => wslpi !: i-direction slope at u-, w-points |
---|
134 | USE ldfslp , ONLY : wslpj => wslpj !: j-direction slope at v-, w-points |
---|
135 | # endif |
---|
136 | |
---|
137 | #else |
---|
138 | !!---------------------------------------------------------------------- |
---|
139 | !! Empty module : No passive tracer |
---|
140 | !!---------------------------------------------------------------------- |
---|
141 | #endif |
---|
142 | |
---|
143 | !!---------------------------------------------------------------------- |
---|
144 | !! NEMO/TOP 3.3 , NEMO Consortium (2010) |
---|
145 | !! $Id$ |
---|
146 | !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) |
---|
147 | !!====================================================================== |
---|
148 | END MODULE oce_trc |
---|