New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
oce_trc.F90 in branches/2015/dev_r5003_MERCATOR6_CRS/NEMOGCM/NEMO/TOP_SRC – NEMO

source: branches/2015/dev_r5003_MERCATOR6_CRS/NEMOGCM/NEMO/TOP_SRC/oce_trc.F90 @ 5601

Last change on this file since 5601 was 5601, checked in by cbricaud, 9 years ago

commit changes/bugfix/... for crs ; ok with time-splitting/fixed volume

  • Property svn:keywords set to Id
File size: 22.6 KB
Line 
1MODULE 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#if defined key_crs
12
13  !* Domain size *
14   USE par_oce , ONLY :   jpi      =>   jpi        !: first  dimension of grid --> i
15   USE par_oce , ONLY :   jpj      =>   jpj        !: second dimension of grid --> j 
16   USE par_oce , ONLY :   jpk      =>   jpk        !: number of levels 
17   USE par_oce , ONLY :   jpim1    =>   jpim1      !: jpi - 1
18   USE par_oce , ONLY :   jpjm1    =>   jpjm1      !: jpj - 1
19   USE par_oce , ONLY :   jpkm1    =>   jpkm1      !: jpk - 1 
20   USE par_oce , ONLY :   jpij     =>   jpij       !: jpi x jpj
21   USE par_oce , ONLY :   lk_esopa =>   lk_esopa   !: flag to activate the all option
22   USE par_oce , ONLY :   jp_tem   =>   jp_tem     !: indice for temperature
23   USE par_oce , ONLY :   jp_sal   =>   jp_sal     !: indice for salinity
24
25   !* IO manager *
26   USE in_out_manager
27
28   !* Memory Allocation *
29   USE wrk_nemo
30
31   !* Timing *
32   USE timing, ONLY : timing_start , timing_stop
33
34   !* MPP library                         
35   USE lib_mpp
36
37   !* Fortran utilities                         
38   USE lib_fortran
39
40   !* Lateral boundary conditions                         
41   USE lbclnk
42
43   !* physical constants *
44   USE phycst
45
46   !* 1D configuration
47   USE c1d
48
49   !* model domain *
50   USE dom_oce , ONLY : narea => narea
51   USE dom_oce , ONLY : nproc => nproc
52   USE dom_oce , ONLY : nimpp => nimpp
53   USE dom_oce , ONLY : njmpp => njmpp
54   USE dom_oce , ONLY : nreci => nreci
55   USE dom_oce , ONLY : nrecj => nrecj
56   USE dom_oce , ONLY : nlci  => nlci
57   USE dom_oce , ONLY : nldi  => nldi
58   USE dom_oce , ONLY : nlei  => nlei
59   USE dom_oce , ONLY : nlcj  => nlcj
60   USE dom_oce , ONLY : nldj  => nldj
61   USE dom_oce , ONLY : nlej  => nlej
62   USE dom_oce , ONLY : nlcit  => nlcit
63   USE dom_oce , ONLY : nldit  => nldit
64   USE dom_oce , ONLY : nleit  => nleit
65   USE dom_oce , ONLY : nlcjt  => nlcjt
66   USE dom_oce , ONLY : nldjt  => nldjt
67   USE dom_oce , ONLY : nlejt  => nlejt
68   USE dom_oce , ONLY : nimppt => nimppt
69   USE dom_oce , ONLY : njmppt => njmppt
70   USE dom_oce , ONLY : ibonit => ibonit
71   USE dom_oce , ONLY : ibonjt => ibonjt
72   USE dom_oce , ONLY : lk_vvl => lk_vvl
73   USE dom_oce , ONLY : rdt => rdt
74   USE dom_oce , ONLY : ln_zco => ln_zco
75   USE dom_oce , ONLY : ln_zps => ln_zps
76   USE dom_oce , ONLY : ln_sco => ln_sco
77   USE dom_oce , ONLY : neuler => neuler
78
79   USE crs,  ONLY : mi0 => mi0 
80   USE crs,  ONLY : mi1 => mi1 
81   USE crs,  ONLY : mj0 => mj0 
82   USE crs,  ONLY : mj1 => mj1 
83
84   USE dom_oce , ONLY :  lzoom => lzoom 
85   !USE dom_oce , ONLY :  =>
86
87   !* horizontal mesh *
88   USE crs , ONLY :   glamt      =>   glamt_crs      !: longitude of t-point (degre) 
89   USE crs , ONLY :   glamu      =>   glamu_crs      !: longitude of t-point (degre) 
90   USE crs , ONLY :   glamv      =>   glamv_crs      !: longitude of t-point (degre) 
91   USE crs , ONLY :   glamf      =>   glamf_crs      !: longitude of t-point (degre) 
92   USE crs , ONLY :   gphit      =>   gphit_crs      !: latitude  of t-point (degre)   
93   USE crs , ONLY :   gphiu      =>   gphiu_crs      !: latitude  of t-point (degre)   
94   USE crs , ONLY :   gphiv      =>   gphiv_crs      !: latitude  of t-point (degre)   
95   USE crs , ONLY :   gphif      =>   gphif_crs      !: latitude  of t-point (degre)   
96   USE crs , ONLY :   e1t        =>   e1t_crs        !: horizontal scale factors at t-point (m) 
97   USE crs , ONLY :   e2t        =>   e2t_crs        !: horizontal scale factors at t-point (m)   
98   USE crs , ONLY :   e1e2t      =>   e1e2t_crs      !: cell surface at t-point (m2)
99   USE crs , ONLY :   e1u        =>   e1u_crs        !: horizontal scale factors at u-point (m)
100   USE crs , ONLY :   e2u        =>   e2u_crs        !: horizontal scale factors at u-point (m)
101   USE crs , ONLY :   e1v        =>   e1v_crs        !: horizontal scale factors at v-point (m)
102   USE crs , ONLY :   e2v        =>   e2v_crs        !: horizontal scale factors at v-point (m) 
103   USE crs , ONLY :   e3t        =>  e3t_crs         !: vertical scale factors at t-
104   USE crs , ONLY :   e3t_0      =>  e3t_crs         !: vertical scale factors at t-
105   USE crs , ONLY :   fse3t      =>  e3t_crs
106   USE crs , ONLY :   fse3t_b      =>  e3t_crs
107   USE crs , ONLY :   fse3t_a      =>  e3t_crs
108   USE crs , ONLY :   fse3w      =>  e3w_crs
109   USE crs , ONLY :   e3u        =>  e3u_crs         !: vertical scale factors at u-
110   USE crs , ONLY :   e3u_0      =>  e3u_crs         !: vertical scale factors at u-
111   USE crs , ONLY :   e3v        =>  e3v_crs         !: vertical scale factors v-
112   USE crs , ONLY :   e3v_0      =>  e3v_crs         !: vertical scale factors v-
113   USE crs , ONLY :   e3w        =>  e3w_crs         !: w-points (m)
114   USE crs , ONLY :   e3w_0      =>  e3w_crs         !: w-points (m)
115   USE crs , ONLY :   e3f        =>  e3f_crs         !: f-points (m)
116   USE crs , ONLY :   ff         =>  ff_crs         !: f-points (m)
117
118   USE crs , ONLY :   gdept_0    =>  gdept_crs       !: depth of t-points (m)
119   USE dom_oce , ONLY :   gdept_1d   =>  gdept_1d      !: depth of t-points (m)
120#if defined key_zco
121   USE crs , ONLY :   gdept      =>  gdept_crs       !: depth of t-points (m)
122   USE crs , ONLY :   gdepw      =>  gdepw_crs       !: depth of t-points (m)
123#endif
124  !* masks, bathymetry *
125   USE crs , ONLY :   mbkt       =>   mbkt_crs       !: vertical index of the bottom last T- ocean level
126   USE crs , ONLY :   mbku       =>   mbku_crs       !: vertical index of the bottom last U- ocean level
127   USE crs , ONLY :   mbkv       =>   mbkv_crs       !: vertical index of the bottom last V- ocean level
128   USE crs , ONLY :   tmask_i    =>   tmask_i_crs    !: Interior mask at t-points
129   USE crs , ONLY :   tmask      =>   tmask_crs      !: land/ocean mask at t-points
130   USE crs , ONLY :   umask      =>   umask_crs      !: land/ocean mask at u-points   
131   USE crs , ONLY :   vmask      =>   vmask_crs      !: land/ocean mask at v-points
132   USE crs , ONLY :   fmask      =>   fmask_crs      !: land/ocean mask at f-points
133
134 !* ocean fields: here now and after fields *
135!cbr?   USE crs , ONLY :   ua      =>    ua_crs      !: i-horizontal velocity (m s-1)
136!cbr?   USE crs , ONLY :   va      =>    va_crs      !: j-horizontal velocity (m s-1)
137   USE crs , ONLY :   un      =>    un_crs      !: i-horizontal velocity (m s-1)
138   USE crs , ONLY :   vn      =>    vn_crs      !: j-horizontal velocity (m s-1)
139   USE crs , ONLY :   wn      =>    wn_crs      !: vertical velocity (m s-1) 
140   USE crs , ONLY :   tsn     =>    tsn_crs     !: 4D array contaning ( tn, sn )
141   USE oce , ONLY :   tsb     =>    tsb     !: 4D array contaning ( tb, sb )
142   USE oce , ONLY :   tsa     =>    tsa     !: 4D array contaning ( ta, sa )
143   USE oce , ONLY :   rhop    =>    rhop    !: potential volumic mass (kg m-3)
144   USE crs , ONLY :   rhd     =>    rhd_crs    !: in situ density anomalie rhd=(rho-rau0)/rau0 (no units)
145   USE crs , ONLY :   rn2b    =>    rb2_crs     !: in situ density anomalie rhd=(rho-rau0)/rau0 (no units)
146   USE crs , ONLY :   rab_n   =>    rab_crs_n     !: in situ density anomalie rhd=(rho-rau0)/rau0 (no units)
147   USE crs , ONLY :   hdivn   =>    hdivn_crs   !: horizontal divergence (1/s)
148   USE crs , ONLY :   hdivb   =>    hdivb_crs   !: horizontal divergence (1/s)
149   USE crs , ONLY :   sshb    =>    sshb_crs    !: sea surface height at t-point [m]   
150   USE crs , ONLY :   sshn    =>    sshn_crs    !: sea surface height at t-point [m]   
151   USE crs , ONLY :   ssha    =>    ssha_crs    !: sea surface height at t-point [m]   
152
153   !* surface fluxes *
154   USE crs , ONLY :   utau       =>    utau_crs       !: i-surface stress component
155   USE crs , ONLY :   vtau       =>    vtau_crs       !: j-surface stress component
156   USE crs , ONLY :   wndm       =>    wndm_crs       !: 10m wind speed
157   USE crs , ONLY :   qsr        =>    qsr_crs        !: penetrative solar radiation (w m-2) 
158   USE crs , ONLY :   emp        =>    emp_crs        !: freshwater budget: volume flux               [Kg/m2/s]
159   USE crs , ONLY :   emp_b      =>    emp_b_crs      !: freshwater budget: volume flux               [Kg/m2/s]
160   USE crs , ONLY :   sfx        =>    sfx_crs        !: freshwater budget: concentration/dillution   [Kg/m2/s]
161   USE crs , ONLY :   fmmflx     =>    fmmflx_crs     !: freshwater budget: volume flux               [Kg/m2/s]
162   USE crs , ONLY :   rnf        =>    rnf_crs        !: river runoff   [Kg/m2/s]
163   USE crs , ONLY :   fr_i       =>    fr_i_crs       !: ice fraction (between 0 to 1)
164   USE trcnam_trp , ONLY :  aht0     =>   rn_ahtrc_0        !: horizontal eddy diffusivity for tracers (m2/s)
165   USE crs , ONLY :  ahtu     =>   ahtu_crs        !: lateral diffusivity coef. at u-points
166   USE crs , ONLY :  ahtv     =>   ahtv_crs        !: lateral diffusivity coef. at v-points
167   USE crs , ONLY :  ahtw     =>   ahtw_crs        !: lateral diffusivity coef. at w-points
168   USE crs , ONLY :  ahtt     =>   ahtt_crs        !: lateral diffusivity coef. at t-points
169   USE ldftra_oce , ONLY :  rldf     =>   rldf
170
171   USE crs , ONLY :   avt        =>   avt_crs         !: vert. diffusivity coef. at w-point for temp 
172#if defined key_zdfddm
173   USE crs , ONLY :   avs        =>   avs_crs         !: salinity vertical diffusivity coeff. at w-point
174#endif
175
176!cbr   USE trc_oce
177   USE trc_oce, ONLY : lk_offline
178   USE trc_oce, ONLY : nn_dttrc
179
180   USE crs , ONLY :   nmln        =>   nmln_crs        !: number of level in the mixed layer
181   USE crs , ONLY :   hmld        =>   hmld_crs        !: mixing layer depth (turbocline)
182   USE crs , ONLY :   hmlp        =>   hmlp_crs        !: mixed layer depth  (rho=rho0+zdcrit) (m)
183   USE crs , ONLY :   hmlpt       =>   hmlpt_crs       !: mixed layer depth at t-points (m)
184
185  !* direction of lateral diffusion *
186#if   defined key_ldfslp
187   USE crs , ONLY :   uslp       =>   uslp_crs         !: i-direction slope at u-, w-points
188   USE crs , ONLY :   vslp       =>   vslp_crs         !: j-direction slope at v-, w-points
189   USE crs , ONLY :   wslpi      =>   wslpi_crs        !: i-direction slope at u-, w-points
190   USE crs , ONLY :   wslpj      =>   wslpj_crs        !: j-direction slope at v-, w-points
191#endif
192
193#else
194
195   !!----------------------------------------------------------------------
196   !!   'key_top'                                                TOP models
197   !!----------------------------------------------------------------------
198
199   !* Domain size *
200   USE par_oce , ONLY :   jpi      =>   jpi        !: first  dimension of grid --> i
201   USE par_oce , ONLY :   jpj      =>   jpj        !: second dimension of grid --> j 
202   USE par_oce , ONLY :   jpk      =>   jpk        !: number of levels 
203   USE par_oce , ONLY :   jpim1    =>   jpim1      !: jpi - 1
204   USE par_oce , ONLY :   jpjm1    =>   jpjm1      !: jpj - 1
205   USE par_oce , ONLY :   jpkm1    =>   jpkm1      !: jpk - 1 
206   USE par_oce , ONLY :   jpij     =>   jpij       !: jpi x jpj
207   USE par_oce , ONLY :   lk_esopa =>   lk_esopa   !: flag to activate the all option
208   USE par_oce , ONLY :   jp_tem   =>   jp_tem     !: indice for temperature
209   USE par_oce , ONLY :   jp_sal   =>   jp_sal     !: indice for salinity
210
211  !* model domain *
212   USE dom_oce , ONLY : narea => narea
213   USE dom_oce , ONLY : nproc => nproc
214   USE dom_oce , ONLY : nimpp => nimpp
215   USE dom_oce , ONLY : njmpp => njmpp
216   USE dom_oce , ONLY : nreci => nreci
217   USE dom_oce , ONLY : nrecj => nrecj
218   USE dom_oce , ONLY : nlci  => nlci
219   USE dom_oce , ONLY : nldi  => nldi
220   USE dom_oce , ONLY : nlei  => nlei
221   USE dom_oce , ONLY : nlcj  => nlcj
222   USE dom_oce , ONLY : nldj  => nldj
223   USE dom_oce , ONLY : nlej  => nlej
224   USE dom_oce , ONLY : nlcit  => nlcit
225   USE dom_oce , ONLY : nldit  => nldit
226   USE dom_oce , ONLY : nleit  => nleit
227   USE dom_oce , ONLY : nlcjt  => nlcjt
228   USE dom_oce , ONLY : nldjt  => nldjt
229   USE dom_oce , ONLY : nlejt  => nlejt
230   USE dom_oce , ONLY : nimppt => nimppt
231   USE dom_oce , ONLY : njmppt => njmppt
232   USE dom_oce , ONLY : ibonit => ibonit
233   USE dom_oce , ONLY : ibonjt => ibonjt
234   USE dom_oce , ONLY : lk_vvl => lk_vvl
235   USE dom_oce , ONLY : rdt => rdt
236   USE dom_oce , ONLY : ln_zco => ln_zco
237   USE dom_oce , ONLY : ln_zps => ln_zps
238   USE dom_oce , ONLY : ln_sco => ln_sco
239   USE dom_oce , ONLY : neuler => neuler
240
241   USE dom_oce,  ONLY : mi0 => mi0
242   USE dom_oce,  ONLY : mi1 => mi1
243   USE dom_oce,  ONLY : mj0 => mj0
244   USE dom_oce,  ONLY : mj1 => mj1
245
246   USE dom_oce , ONLY :   glamt      =>   glamt      !: longitude of t-point (degre) 
247   USE dom_oce , ONLY :   glamu      =>   glamu      !: longitude of t-point (degre) 
248   USE dom_oce , ONLY :   glamv      =>   glamv      !: longitude of t-point (degre) 
249   USE dom_oce , ONLY :   glamf      =>   glamf      !: longitude of t-point (degre) 
250   USE dom_oce , ONLY :   gphit      =>   gphit      !: latitude  of t-point (degre)   
251   USE dom_oce , ONLY :   gphiu      =>   gphiu      !: latitude  of t-point (degre)   
252   USE dom_oce , ONLY :   gphiv      =>   gphiv      !: latitude  of t-point (degre)   
253   USE dom_oce , ONLY :   gphif      =>   gphif     !: latitude  of t-point (degre)   
254   USE dom_oce , ONLY :   e1t        =>   e1t        !: horizontal scale factors at t-point (m) 
255   USE dom_oce , ONLY :   e2t        =>   e2t        !: horizontal scale factors at t-point (m)   
256   USE dom_oce , ONLY :   e1e2t      =>   e1e2t      !: cell surface at t-point (m2)
257   USE dom_oce , ONLY :   e1u        =>   e1u        !: horizontal scale factors at u-point (m)
258   USE dom_oce , ONLY :   e2u        =>   e2u        !: horizontal scale factors at u-point (m)
259   USE dom_oce , ONLY :   e1v        =>   e1v        !: horizontal scale factors at v-point (m)
260   USE dom_oce , ONLY :   e2v        =>   e2v        !: horizontal scale factors at v-point (m) 
261   USE dom_oce , ONLY :   e3t        =>  e3t_0         !: vertical scale factors at t-
262   USE dom_oce , ONLY :   e3t_0      =>  e3t_0         !: vertical scale factors at t-
263   USE dom_oce , ONLY :   fse3t      =>  e3t_0
264   USE dom_oce , ONLY :   fse3t_b      =>  e3t_0
265   USE dom_oce , ONLY :   fse3t_a      =>  e3t_0
266   USE dom_oce , ONLY :   fse3w      =>  e3w_0
267   USE dom_oce , ONLY :   e3u        =>  e3u_0         !: vertical scale factors at u-
268   USE dom_oce , ONLY :   e3u_0      =>  e3u_0         !: vertical scale factors at u-
269   USE dom_oce , ONLY :   e3v        =>  e3v_0         !: vertical scale factors v-
270   USE dom_oce , ONLY :   e3v_0      =>  e3v_0         !: vertical scale factors v-
271   USE dom_oce , ONLY :   e3w        =>  e3w_0         !: w-points (m)
272   USE dom_oce , ONLY :   e3w_0      =>  e3w_0         !: w-points (m)
273   USE dom_oce , ONLY :   e3f        =>  e3f_0         !: f-points (m)
274   USE dom_oce , ONLY :   ff         =>  ff         !: f-points (m)
275   USE dom_oce , ONLY :   gdept_0    =>  gdept_0         !: f-points (m)
276   USE dom_oce , ONLY :   gdept_1d   => gdept_1d          !: f-points (m)
277   USE dom_oce , ONLY :   tmask      => tmask          !: f-points (m)
278   USE dom_oce , ONLY :   umask      => umask          !: f-points (m)
279   USE dom_oce , ONLY :   vmask      => vmask          !: f-points (m)
280   USE dom_oce , ONLY :   tmask_i      => tmask_i          !: f-points (m)
281   USE dom_oce , ONLY :   mbkt      => mbkt          !: f-points (m)
282   USE dom_oce , ONLY :   mbku      => mbku          !: f-points (m)
283   USE dom_oce , ONLY :   mbkv      => mbkv          !: f-points (m)
284
285   !* IO manager *
286   USE in_out_manager   
287 
288   !* Memory Allocation *
289   USE wrk_nemo     
290 
291   !* Timing *
292   USE timing   
293 
294   !* MPP library                         
295   USE lib_mpp 
296
297   !* Fortran utilities                         
298   USE lib_fortran
299
300   !* Lateral boundary conditions                         
301   USE lbclnk
302
303   !* physical constants *
304   USE phycst           
305
306   !* 1D configuration
307   USE c1d                                         
308
309   !* model domain *
310   !cbr USE dom_oce , ONLY : e3w_0
311   USE dom_oce , ONLY :  lzoom => lzoom 
312
313   USE domvvl, ONLY : un_td, vn_td          !: thickness diffusion transport
314   USE domvvl, ONLY : ln_vvl_ztilde         !: ztilde vertical coordinate
315   USE domvvl, ONLY : ln_vvl_layer          !: level  vertical coordinate
316
317   !* ocean fields: here now and after fields *
318   USE oce , ONLY :   ua      =>    ua      !: i-horizontal velocity (m s-1)
319   USE oce , ONLY :   va      =>    va      !: j-horizontal velocity (m s-1)
320   USE oce , ONLY :   un      =>    un      !: i-horizontal velocity (m s-1)
321   USE oce , ONLY :   vn      =>    vn      !: j-horizontal velocity (m s-1)
322   USE oce , ONLY :   wn      =>    wn      !: vertical velocity (m s-1) 
323   USE oce , ONLY :   tsn     =>    tsn     !: 4D array contaning ( tn, sn )
324   USE oce , ONLY :   tsb     =>    tsb     !: 4D array contaning ( tb, sb )
325   USE oce , ONLY :   tsa     =>    tsa     !: 4D array contaning ( ta, sa )
326   USE oce , ONLY :   rhop    =>    rhop    !: potential volumic mass (kg m-3)
327   USE oce , ONLY :   rhd     =>    rhd     !: in situ density anomalie rhd=(rho-rau0)/rau0 (no units)
328   USE oce , ONLY :   rab_n     =>    rab_n   !: in situ density anomalie rhd=(rho-rau0)/rau0 (no units)
329   USE oce , ONLY :   rn2b    =>    rn2b    !: in situ density anomalie rhd=(rho-rau0)/rau0 (no units)
330   USE oce , ONLY :   hdivn   =>    hdivn   !: horizontal divergence (1/s)
331   USE oce , ONLY :   rotn    =>    rotn    !: relative vorticity    [s-1]
332   USE oce , ONLY :   hdivb   =>    hdivb   !: horizontal divergence (1/s)
333   USE oce , ONLY :   rotb    =>    rotb    !: relative vorticity    [s-1]
334   USE oce , ONLY :   sshn    =>    sshn    !: sea surface height at t-point [m]   
335   USE oce , ONLY :   sshb    =>    sshb    !: sea surface height at t-point [m]   
336   USE oce , ONLY :   ssha    =>    ssha    !: sea surface height at t-point [m]   
337   USE oce , ONLY :   l_traldf_rot => l_traldf_rot  !: rotated laplacian operator for lateral diffusion
338
339   !* surface fluxes *
340   USE sbc_oce , ONLY :   utau       =>    utau       !: i-surface stress component
341   USE sbc_oce , ONLY :   vtau       =>    vtau       !: j-surface stress component
342   USE sbc_oce , ONLY :   wndm       =>    wndm       !: 10m wind speed
343   USE sbc_oce , ONLY :   qsr        =>    qsr        !: penetrative solar radiation (w m-2)
344   USE sbc_oce , ONLY :   emp        =>    emp        !: freshwater budget: volume flux               [Kg/m2/s]
345   USE sbc_oce , ONLY :   emp_b      =>    emp_b      !: freshwater budget: volume flux               [Kg/m2/s]
346   USE sbc_oce , ONLY :   fmmflx     =>    fmmflx     !: freshwater budget: volume flux               [Kg/m2/s]
347   USE sbc_oce , ONLY :   rnf        =>    rnf        !: river runoff   [Kg/m2/s]
348   USE sbc_oce , ONLY :   ln_dm2dc   =>    ln_dm2dc   !: Daily mean to Diurnal Cycle short wave (qsr)
349   USE sbc_oce , ONLY :   ln_rnf     =>    ln_rnf     !: runoffs / runoff mouths
350   USE sbc_oce , ONLY :   fr_i       =>    fr_i       !: ice fraction (between 0 to 1)
351   USE traqsr  , ONLY :   rn_abs     =>    rn_abs     !: fraction absorbed in the very near surface
352   USE traqsr  , ONLY :   rn_si0     =>    rn_si0     !: very near surface depth of extinction
353   USE traqsr  , ONLY :   ln_qsr_bio =>    ln_qsr_bio !: flag to use or not the biological fluxes for light
354   USE sbcrnf  , ONLY :   rnfmsk     =>    rnfmsk     !: mixed adv scheme in runoffs vicinity (hori.)
355   USE sbcrnf  , ONLY :   rnfmsk_z   =>    rnfmsk_z   !: mixed adv scheme in runoffs vicinity (vert.)
356   USE sbcrnf  , ONLY :   h_rnf      =>    h_rnf      !: river runoff   [Kg/m2/s]
357
358   USE trc_oce
359
360   !* lateral diffusivity (tracers) *
361   USE ldftra_oce , ONLY :  rldf     =>   rldf        !: multiplicative coef. for lateral diffusivity
362   USE ldftra_oce , ONLY :  rn_aht_0 =>   rn_aht_0    !: horizontal eddy diffusivity for tracers (m2/s)
363   USE ldftra_oce , ONLY :  aht0     =>   aht0        !: horizontal eddy diffusivity for tracers (m2/s)
364   USE ldftra_oce , ONLY :  ahtb0    =>   ahtb0       !: background eddy diffusivity for isopycnal diff. (m2/s)
365   USE ldftra_oce , ONLY :  ahtu     =>   ahtu        !: lateral diffusivity coef. at u-points
366   USE ldftra_oce , ONLY :  ahtv     =>   ahtv        !: lateral diffusivity coef. at v-points
367   USE ldftra_oce , ONLY :  ahtw     =>   ahtw        !: lateral diffusivity coef. at w-points
368   USE ldftra_oce , ONLY :  ahtt     =>   ahtt        !: lateral diffusivity coef. at t-points
369   USE ldftra_oce , ONLY :  aeiv0    =>   aeiv0       !: eddy induced velocity coefficient (m2/s)
370   USE ldftra_oce , ONLY :  aeiu     =>   aeiu        !: eddy induced velocity coef. at u-points (m2/s)   
371   USE ldftra_oce , ONLY :  aeiv     =>   aeiv        !: eddy induced velocity coef. at v-points (m2/s)
372   USE ldftra_oce , ONLY :  aeiw     =>   aeiw        !: eddy induced velocity coef. at w-points (m2/s)
373   USE ldftra_oce , ONLY :  lk_traldf_eiv  =>  lk_traldf_eiv     !: eddy induced velocity flag
374
375   !* vertical diffusion *
376   USE zdf_oce , ONLY :   avt        =>   avt         !: vert. diffusivity coef. at w-point for temp 
377# if defined key_zdfddm
378   USE zdfddm  , ONLY :   avs        =>   avs         !: salinity vertical diffusivity coeff. at w-point
379# endif
380
381   !* mixing & mixed layer depth *
382   USE zdfmxl , ONLY :   nmln        =>   nmln        !: number of level in the mixed layer
383   USE zdfmxl , ONLY :   hmld        =>   hmld        !: mixing layer depth (turbocline)
384   USE zdfmxl , ONLY :   hmlp        =>   hmlp        !: mixed layer depth  (rho=rho0+zdcrit) (m)
385   USE zdfmxl , ONLY :   hmlpt       =>   hmlpt       !: mixed layer depth at t-points (m)
386
387   !* direction of lateral diffusion *
388   USE ldfslp , ONLY :   lk_ldfslp  =>  lk_ldfslp     !: slopes flag
389# if   defined key_ldfslp
390   USE ldfslp , ONLY :   uslp       =>   uslp         !: i-direction slope at u-, w-points
391   USE ldfslp , ONLY :   vslp       =>   vslp         !: j-direction slope at v-, w-points
392   USE ldfslp , ONLY :   wslpi      =>   wslpi        !: i-direction slope at u-, w-points
393   USE ldfslp , ONLY :   wslpj      =>   wslpj        !: j-direction slope at v-, w-points
394# endif
395
396#endif
397#else
398   !!----------------------------------------------------------------------
399   !!  Empty module :                                     No passive tracer
400   !!----------------------------------------------------------------------
401#endif
402
403   !!----------------------------------------------------------------------
404   !! NEMO/TOP 3.3 , NEMO Consortium (2010)
405   !! $Id$
406   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt)
407   !!======================================================================
408END MODULE oce_trc
Note: See TracBrowser for help on using the repository browser.