source: trunk/NEMOGCM/NEMO/TOP_SRC/oce_trc.F90 @ 7646

Last change on this file since 7646 was 7646, checked in by timgraham, 4 years ago

Merge of dev_merge_2016 into trunk. UPDATE TO ARCHFILES NEEDED for XIOS2.
LIM_SRC_s/limrhg.F90 to follow in next commit due to change of kind (I'm unable to do it in this commit).
Merged using the following steps:

1) svn merge —reintegrate svn+ssh://forge.ipsl.jussieu.fr/ipsl/forge/projets/nemo/svn/trunk .
2) Resolve minor conflicts in sette.sh and namelist_cfg for ORCA2LIM3 (due to a change in trunk after branch was created)
3) svn commit
4) svn switch svn+ssh://forge.ipsl.jussieu.fr/ipsl/forge/projets/nemo/svn/trunk
5) svn merge svn+ssh://forge.ipsl.jussieu.fr/ipsl/forge/projets/nemo/svn/branches/2016/dev_merge_2016 .
6) At this stage I checked out a clean copy of the branch to compare against what is about to be committed to the trunk.
6) svn commit #Commit code to the trunk

In this commit I have also reverted a change to Fcheck_archfile.sh which was causing problems on the Paris machine.

  • Property svn:keywords set to Id
File size: 8.5 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 :   jp_tem   =>   jp_tem     !: indice for temperature
23   USE par_oce , ONLY :   jp_sal   =>   jp_sal     !: indice for salinity
24
25   USE in_out_manager                           !* IO manager *
26   USE wrk_nemo                                 !* Memory Allocation *
27   USE timing                                   !* Timing *
28   USE lib_mpp                                  !* MPP library                         
29   USE lib_fortran                              !* Fortran utilities                         
30   USE lbclnk                                   !* Lateral boundary conditions                         
31   USE phycst                                   !* physical constants *
32   USE c1d                                      !* 1D configuration
33
34   USE dom_oce                                  !* model domain *
35
36   USE domvvl, ONLY : un_td, vn_td          !: thickness diffusion transport
37   USE domvvl, ONLY : ln_vvl_ztilde         !: ztilde vertical coordinate
38   USE domvvl, ONLY : ln_vvl_layer          !: level  vertical coordinate
39
40   !* ocean fields: here now and after fields *
41   USE oce , ONLY :   un      =>    un      !: i-horizontal velocity (m s-1)
42   USE oce , ONLY :   vn      =>    vn      !: j-horizontal velocity (m s-1)
43   USE oce , ONLY :   wn      =>    wn      !: vertical velocity (m s-1) 
44   USE oce , ONLY :   tsn     =>    tsn     !: 4D array contaning ( tn, sn )
45   USE oce , ONLY :   tsb     =>    tsb     !: 4D array contaning ( tb, sb )
46   USE oce , ONLY :   tsa     =>    tsa     !: 4D array contaning ( ta, sa )
47   USE oce , ONLY :   rhop    =>    rhop    !: potential volumic mass (kg m-3)
48   USE oce , ONLY :   rhd     =>    rhd     !: in situ density anomalie rhd=(rho-rau0)/rau0 (no units)
49   USE oce , ONLY :   hdivn   =>    hdivn   !: horizontal divergence (1/s)
50   USE oce , ONLY :   sshn    =>    sshn    !: sea surface height at t-point [m]   
51   USE oce , ONLY :   sshb    =>    sshb    !: sea surface height at t-point [m]   
52   USE oce , ONLY :   ssha    =>    ssha    !: sea surface height at t-point [m]   
53   USE oce , ONLY :   rab_n   =>    rab_n   !: local thermal/haline expension ratio at T-points
54
55   !* surface fluxes *
56   USE sbc_oce , ONLY :   utau       =>    utau       !: i-surface stress component
57   USE sbc_oce , ONLY :   vtau       =>    vtau       !: j-surface stress component
58   USE sbc_oce , ONLY :   wndm       =>    wndm       !: 10m wind speed
59   USE sbc_oce , ONLY :   qsr        =>    qsr        !: penetrative solar radiation (w m-2)
60   USE sbc_oce , ONLY :   emp        =>    emp        !: freshwater budget: volume flux               [Kg/m2/s]
61   USE sbc_oce , ONLY :   emp_b      =>    emp_b      !: freshwater budget: volume flux               [Kg/m2/s]
62   USE sbc_oce , ONLY :   fmmflx     =>    fmmflx     !: freshwater budget: volume flux               [Kg/m2/s]
63   USE sbc_oce , ONLY :   rnf        =>    rnf        !: river runoff   [Kg/m2/s]
64   USE sbc_oce , ONLY :   rnf_b      =>    rnf_b      !: river runoff at previus step   [Kg/m2/s]
65   USE sbc_oce , ONLY :   ln_dm2dc   =>    ln_dm2dc   !: Diurnal Cycle
66   USE sbc_oce , ONLY :   ln_cpl     =>    ln_cpl     !: ocean-atmosphere coupled formulation
67   USE sbc_oce , ONLY :   ncpl_qsr_freq   =>   ncpl_qsr_freq   !: qsr coupling frequency per days from atmospher
68   USE sbc_oce , ONLY :   ln_rnf     =>    ln_rnf     !: runoffs / runoff mouths
69   USE sbc_oce , ONLY :   fr_i       =>    fr_i       !: ice fraction (between 0 to 1)
70   USE sbc_oce , ONLY :   nn_ice_embd => nn_ice_embd  !: flag for  levitating/embedding sea-ice in the ocean
71   USE sbc_oce , ONLY :   atm_co2    =>    atm_co2    !  atmospheric pCO2
72   USE traqsr  , ONLY :   rn_abs     =>    rn_abs     !: fraction absorbed in the very near surface
73   USE traqsr  , ONLY :   rn_si0     =>    rn_si0     !: very near surface depth of extinction
74   USE traqsr  , ONLY :   ln_qsr_bio =>    ln_qsr_bio !: flag to use or not the biological fluxes for light
75   USE sbcrnf  , ONLY :   rnfmsk     =>    rnfmsk     !: mixed adv scheme in runoffs vicinity (hori.)
76   USE sbcrnf  , ONLY :   rnfmsk_z   =>    rnfmsk_z   !: mixed adv scheme in runoffs vicinity (vert.)
77   USE sbcrnf  , ONLY :   h_rnf      =>    h_rnf      !: river runoff   [Kg/m2/s]
78   USE sbcrnf  , ONLY :   nk_rnf     =>    nk_rnf     !: depth of runoff in model level
79   USE sbcrnf  , ONLY :   rn_rfact   =>    rn_rfact   !: multiplicative factor for runoff
80
81   USE trc_oce
82
83!!gm : I don't understand this as ldftra (where everything is defined) is used by TRC in all cases (ON/OFF-line)
84!!gm   so the following lines should be removed....   logical should be the one of TRC namelist
85!!gm   In case off coarsening....  the ( ahtu, ahtv, aeiu, aeiv) arrays are needed that's all.
86   !* lateral diffusivity (tracers) *
87   USE ldftra , ONLY :  rn_aht_0     =>   rn_aht_0     !:   laplacian lateral eddy diffusivity [m2/s]
88   USE ldftra , ONLY :  rn_bht_0     =>   rn_bht_0     !: bilaplacian lateral eddy diffusivity [m4/s]
89   USE ldftra , ONLY :  ahtu         =>   ahtu         !: lateral diffusivity coef. at u-points
90   USE ldftra , ONLY :  ahtv         =>   ahtv         !: lateral diffusivity coef. at v-points
91   USE ldftra , ONLY :  rn_aeiv_0    =>   rn_aeiv_0    !: eddy induced velocity coefficient (m2/s)
92   USE ldftra , ONLY :  aeiu         =>   aeiu         !: eddy induced velocity coef. at u-points (m2/s)   
93   USE ldftra , ONLY :  aeiv         =>   aeiv         !: eddy induced velocity coef. at v-points (m2/s)
94   USE ldftra , ONLY :  ln_ldfeiv    =>   ln_ldfeiv    !: eddy induced velocity flag
95     
96!!gm this should be : ln_trcldf_triad (TRC namelist)
97   USE ldfslp , ONLY :  ln_traldf_triad => ln_traldf_triad   !: triad scheme (Griffies et al.)
98
99   !* direction of lateral diffusion *
100   USE ldfslp , ONLY :   l_ldfslp  =>  l_ldfslp       !: slopes flag
101   USE ldfslp , ONLY :   uslp       =>   uslp         !: i-slope at u-point
102   USE ldfslp , ONLY :   vslp       =>   vslp         !: j-slope at v-point
103   USE ldfslp , ONLY :   wslpi      =>   wslpi        !: i-slope at w-point
104   USE ldfslp , ONLY :   wslpj      =>   wslpj        !: j-slope at w-point
105!!gm end
106
107   !* vertical diffusion *
108   USE zdf_oce , ONLY :   avt        =>   avt         !: vert. diffusivity coef. at w-point for temp 
109# if defined key_zdfddm
110   USE zdfddm  , ONLY :   avs        =>   avs         !: salinity vertical diffusivity coeff. at w-point
111# endif
112
113   !* mixing & mixed layer depth *
114   USE zdfmxl , ONLY :   nmln        =>   nmln        !: number of level in the mixed layer
115   USE zdfmxl , ONLY :   hmld        =>   hmld        !: mixing layer depth (turbocline)
116   USE zdfmxl , ONLY :   hmlp        =>   hmlp        !: mixed layer depth  (rho=rho0+zdcrit) (m)
117   USE zdfmxl , ONLY :   hmlpt       =>   hmlpt       !: mixed layer depth at t-points (m)
118
119#else
120   !!----------------------------------------------------------------------
121   !!  Empty module :                                     No passive tracer
122   !!----------------------------------------------------------------------
123#endif
124
125   !!----------------------------------------------------------------------
126   !! NEMO/TOP 3.3 , NEMO Consortium (2010)
127   !! $Id$
128   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt)
129   !!======================================================================
130END MODULE oce_trc
Note: See TracBrowser for help on using the repository browser.