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/2014/dev_r4650_UKMO14.12_STAND_ALONE_OBSOPER/NEMOGCM/NEMO/TOP_SRC – NEMO

source: branches/2014/dev_r4650_UKMO14.12_STAND_ALONE_OBSOPER/NEMOGCM/NEMO/TOP_SRC/oce_trc.F90 @ 5600

Last change on this file since 5600 was 5600, checked in by andrewryan, 9 years ago

merged in latest version of trunk alongside changes to SAO_SRC to be compatible with latest OBS

  • Property svn:keywords set to Id
File size: 8.7 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   !!   '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   !: Diurnal Cycle
90   USE sbc_oce , ONLY :   ncpl_qsr_freq   =>   ncpl_qsr_freq   !: qsr coupling frequency per days from atmospher
91   USE sbc_oce , ONLY :   ln_rnf     =>    ln_rnf     !: runoffs / runoff mouths
92   USE sbc_oce , ONLY :   fr_i       =>    fr_i       !: ice fraction (between 0 to 1)
93   USE sbc_oce , ONLY :   nn_ice_embd => nn_ice_embd  !: flag for  levitating/embedding sea-ice in the ocean
94   USE traqsr  , ONLY :   rn_abs     =>    rn_abs     !: fraction absorbed in the very near surface
95   USE traqsr  , ONLY :   rn_si0     =>    rn_si0     !: very near surface depth of extinction
96   USE traqsr  , ONLY :   ln_qsr_bio =>    ln_qsr_bio !: flag to use or not the biological fluxes for light
97   USE sbcrnf  , ONLY :   rnfmsk     =>    rnfmsk     !: mixed adv scheme in runoffs vicinity (hori.)
98   USE sbcrnf  , ONLY :   rnfmsk_z   =>    rnfmsk_z   !: mixed adv scheme in runoffs vicinity (vert.)
99   USE sbcrnf  , ONLY :   h_rnf      =>    h_rnf      !: river runoff   [Kg/m2/s]
100   USE sbcrnf  , ONLY :   nk_rnf     =>    nk_rnf     !: depth of runoff in model level
101
102   USE trc_oce
103
104   !* lateral diffusivity (tracers) *
105   USE ldftra_oce , ONLY :  rldf     =>   rldf        !: multiplicative coef. for lateral diffusivity
106   USE ldftra_oce , ONLY :  rn_aht_0 =>   rn_aht_0    !: horizontal eddy diffusivity for tracers (m2/s)
107   USE ldftra_oce , ONLY :  aht0     =>   aht0        !: horizontal eddy diffusivity for tracers (m2/s)
108   USE ldftra_oce , ONLY :  ahtb0    =>   ahtb0       !: background eddy diffusivity for isopycnal diff. (m2/s)
109   USE ldftra_oce , ONLY :  ahtu     =>   ahtu        !: lateral diffusivity coef. at u-points
110   USE ldftra_oce , ONLY :  ahtv     =>   ahtv        !: lateral diffusivity coef. at v-points
111   USE ldftra_oce , ONLY :  ahtw     =>   ahtw        !: lateral diffusivity coef. at w-points
112   USE ldftra_oce , ONLY :  ahtt     =>   ahtt        !: lateral diffusivity coef. at t-points
113   USE ldftra_oce , ONLY :  aeiv0    =>   aeiv0       !: eddy induced velocity coefficient (m2/s)
114   USE ldftra_oce , ONLY :  aeiu     =>   aeiu        !: eddy induced velocity coef. at u-points (m2/s)   
115   USE ldftra_oce , ONLY :  aeiv     =>   aeiv        !: eddy induced velocity coef. at v-points (m2/s)
116   USE ldftra_oce , ONLY :  aeiw     =>   aeiw        !: eddy induced velocity coef. at w-points (m2/s)
117   USE ldftra_oce , ONLY :  lk_traldf_eiv  =>  lk_traldf_eiv     !: eddy induced velocity flag
118
119   !* vertical diffusion *
120   USE zdf_oce , ONLY :   avt        =>   avt         !: vert. diffusivity coef. at w-point for temp 
121# if defined key_zdfddm
122   USE zdfddm  , ONLY :   avs        =>   avs         !: salinity vertical diffusivity coeff. at w-point
123# endif
124
125   !* mixing & mixed layer depth *
126   USE zdfmxl , ONLY :   nmln        =>   nmln        !: number of level in the mixed layer
127   USE zdfmxl , ONLY :   hmld        =>   hmld        !: mixing layer depth (turbocline)
128   USE zdfmxl , ONLY :   hmlp        =>   hmlp        !: mixed layer depth  (rho=rho0+zdcrit) (m)
129   USE zdfmxl , ONLY :   hmlpt       =>   hmlpt       !: mixed layer depth at t-points (m)
130
131   !* direction of lateral diffusion *
132   USE ldfslp , ONLY :   lk_ldfslp  =>  lk_ldfslp     !: slopes flag
133# if   defined key_ldfslp
134   USE ldfslp , ONLY :   uslp       =>   uslp         !: i-direction slope at u-, w-points
135   USE ldfslp , ONLY :   vslp       =>   vslp         !: j-direction slope at v-, w-points
136   USE ldfslp , ONLY :   wslpi      =>   wslpi        !: i-direction slope at u-, w-points
137   USE ldfslp , ONLY :   wslpj      =>   wslpj        !: j-direction slope at v-, w-points
138# endif
139
140   USE diaar5 , ONLY :   lk_diaar5  =>   lk_diaar5
141#else
142   !!----------------------------------------------------------------------
143   !!  Empty module :                                     No passive tracer
144   !!----------------------------------------------------------------------
145#endif
146
147   !!----------------------------------------------------------------------
148   !! NEMO/TOP 3.3 , NEMO Consortium (2010)
149   !! $Id$
150   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt)
151   !!======================================================================
152END MODULE oce_trc
Note: See TracBrowser for help on using the repository browser.