source: NEMO/trunk/src/TOP/oce_trc.F90 @ 10351

Last change on this file since 10351 was 10351, checked in by mathiot, 2 years ago

fix #1685 (replace 5.e-4_wp by avt_c in dmp routine as it is computed in zdfmxl + apply test on avt instead on avs in trcdmp)

  • Property svn:keywords set to Id
File size: 6.6 KB
RevLine 
[186]1MODULE oce_trc
2   !!======================================================================
3   !!                      ***  MODULE  oce_trc  ***
[945]4   !! TOP :   variables shared between ocean and passive tracers
[186]5   !!======================================================================
[945]6   !! History :   1.0  !  2004-03  (C. Ethe)  original code
7   !!             2.0  !  2007-12 (C. Ethe, G. Madec)  rewritting
[186]8   !!----------------------------------------------------------------------
[5836]9   !                                            !* Domain size *
[945]10   USE par_oce , ONLY :   jpi      =>   jpi        !: first  dimension of grid --> i
11   USE par_oce , ONLY :   jpj      =>   jpj        !: second dimension of grid --> j 
12   USE par_oce , ONLY :   jpk      =>   jpk        !: number of levels 
13   USE par_oce , ONLY :   jpim1    =>   jpim1      !: jpi - 1
14   USE par_oce , ONLY :   jpjm1    =>   jpjm1      !: jpj - 1
15   USE par_oce , ONLY :   jpkm1    =>   jpkm1      !: jpk - 1 
16   USE par_oce , ONLY :   jpij     =>   jpij       !: jpi x jpj
[2528]17   USE par_oce , ONLY :   jp_tem   =>   jp_tem     !: indice for temperature
18   USE par_oce , ONLY :   jp_sal   =>   jp_sal     !: indice for salinity
[186]19
[5836]20   USE in_out_manager                           !* IO manager *
21   USE timing                                   !* Timing *
22   USE lib_mpp                                  !* MPP library                         
23   USE lib_fortran                              !* Fortran utilities                         
24   USE lbclnk                                   !* Lateral boundary conditions                         
25   USE phycst                                   !* physical constants *
26   USE c1d                                      !* 1D configuration
[7646]27
[5836]28   USE dom_oce                                  !* model domain *
[2715]29
[4610]30   USE domvvl, ONLY : un_td, vn_td          !: thickness diffusion transport
31   USE domvvl, ONLY : ln_vvl_ztilde         !: ztilde vertical coordinate
32   USE domvvl, ONLY : ln_vvl_layer          !: level  vertical coordinate
[186]33
[945]34   !* ocean fields: here now and after fields *
35   USE oce , ONLY :   un      =>    un      !: i-horizontal velocity (m s-1)
36   USE oce , ONLY :   vn      =>    vn      !: j-horizontal velocity (m s-1)
37   USE oce , ONLY :   wn      =>    wn      !: vertical velocity (m s-1) 
[2528]38   USE oce , ONLY :   tsn     =>    tsn     !: 4D array contaning ( tn, sn )
39   USE oce , ONLY :   tsb     =>    tsb     !: 4D array contaning ( tb, sb )
40   USE oce , ONLY :   tsa     =>    tsa     !: 4D array contaning ( ta, sa )
[945]41   USE oce , ONLY :   rhop    =>    rhop    !: potential volumic mass (kg m-3)
42   USE oce , ONLY :   rhd     =>    rhd     !: in situ density anomalie rhd=(rho-rau0)/rau0 (no units)
43   USE oce , ONLY :   hdivn   =>    hdivn   !: horizontal divergence (1/s)
[3294]44   USE oce , ONLY :   sshn    =>    sshn    !: sea surface height at t-point [m]   
45   USE oce , ONLY :   sshb    =>    sshb    !: sea surface height at t-point [m]   
46   USE oce , ONLY :   ssha    =>    ssha    !: sea surface height at t-point [m]   
[5836]47   USE oce , ONLY :   rab_n   =>    rab_n   !: local thermal/haline expension ratio at T-points
[201]48
[945]49   !* surface fluxes *
50   USE sbc_oce , ONLY :   utau       =>    utau       !: i-surface stress component
51   USE sbc_oce , ONLY :   vtau       =>    vtau       !: j-surface stress component
[1011]52   USE sbc_oce , ONLY :   wndm       =>    wndm       !: 10m wind speed
[4306]53   USE sbc_oce , ONLY :   qsr        =>    qsr        !: penetrative solar radiation (w m-2)
[2528]54   USE sbc_oce , ONLY :   emp        =>    emp        !: freshwater budget: volume flux               [Kg/m2/s]
[3294]55   USE sbc_oce , ONLY :   emp_b      =>    emp_b      !: freshwater budget: volume flux               [Kg/m2/s]
[4148]56   USE sbc_oce , ONLY :   fmmflx     =>    fmmflx     !: freshwater budget: volume flux               [Kg/m2/s]
[2528]57   USE sbc_oce , ONLY :   rnf        =>    rnf        !: river runoff   [Kg/m2/s]
[7646]58   USE sbc_oce , ONLY :   rnf_b      =>    rnf_b      !: river runoff at previus step   [Kg/m2/s]
[5385]59   USE sbc_oce , ONLY :   ln_dm2dc   =>    ln_dm2dc   !: Diurnal Cycle
[7646]60   USE sbc_oce , ONLY :   ln_cpl     =>    ln_cpl     !: ocean-atmosphere coupled formulation
[5385]61   USE sbc_oce , ONLY :   ncpl_qsr_freq   =>   ncpl_qsr_freq   !: qsr coupling frequency per days from atmospher
[2528]62   USE sbc_oce , ONLY :   ln_rnf     =>    ln_rnf     !: runoffs / runoff mouths
[1034]63   USE sbc_oce , ONLY :   fr_i       =>    fr_i       !: ice fraction (between 0 to 1)
[7646]64   USE sbc_oce , ONLY :   atm_co2    =>    atm_co2    !  atmospheric pCO2
[1445]65   USE traqsr  , ONLY :   rn_abs     =>    rn_abs     !: fraction absorbed in the very near surface
66   USE traqsr  , ONLY :   rn_si0     =>    rn_si0     !: very near surface depth of extinction
67   USE traqsr  , ONLY :   ln_qsr_bio =>    ln_qsr_bio !: flag to use or not the biological fluxes for light
[945]68   USE sbcrnf  , ONLY :   rnfmsk     =>    rnfmsk     !: mixed adv scheme in runoffs vicinity (hori.)
69   USE sbcrnf  , ONLY :   rnfmsk_z   =>    rnfmsk_z   !: mixed adv scheme in runoffs vicinity (vert.)
[3294]70   USE sbcrnf  , ONLY :   h_rnf      =>    h_rnf      !: river runoff   [Kg/m2/s]
[5385]71   USE sbcrnf  , ONLY :   nk_rnf     =>    nk_rnf     !: depth of runoff in model level
[7646]72   USE sbcrnf  , ONLY :   rn_rfact   =>    rn_rfact   !: multiplicative factor for runoff
[186]73
[2528]74   USE trc_oce
[5836]75     
[9490]76!!gm  this can be removed if :
77!!gm    in trcadv.F90 and trcsub.F90  we add a USE ldfslp
[5836]78   !* direction of lateral diffusion *
79   USE ldfslp , ONLY :   l_ldfslp  =>  l_ldfslp       !: slopes flag
80   USE ldfslp , ONLY :   uslp       =>   uslp         !: i-slope at u-point
81   USE ldfslp , ONLY :   vslp       =>   vslp         !: j-slope at v-point
82   USE ldfslp , ONLY :   wslpi      =>   wslpi        !: i-slope at w-point
83   USE ldfslp , ONLY :   wslpj      =>   wslpj        !: j-slope at w-point
[10097]84   USE ldfslp , ONLY :   ln_traldf_triad => ln_traldf_triad   !: use of triad scheme
85   USE ldfslp , ONLY :   ln_traldf_iso => ln_traldf_iso   !: use of isopycnal scheme
[9490]86!!gm end
[5836]87
[945]88   !* vertical diffusion *
[10351]89   USE zdf_oce , ONLY :   avs        =>   avs         !: vert. diffusivity coef. for salinity    (w-point)
90   USE zdf_oce , ONLY :   avt        =>   avt         !: vert. diffusivity coef. for temperature (w-point)
[186]91
[945]92   !* mixing & mixed layer depth *
[1177]93   USE zdfmxl , ONLY :   nmln        =>   nmln        !: number of level in the mixed layer
[945]94   USE zdfmxl , ONLY :   hmld        =>   hmld        !: mixing layer depth (turbocline)
95   USE zdfmxl , ONLY :   hmlp        =>   hmlp        !: mixed layer depth  (rho=rho0+zdcrit) (m)
96   USE zdfmxl , ONLY :   hmlpt       =>   hmlpt       !: mixed layer depth at t-points (m)
[10351]97   USE zdfmxl , ONLY :   avt_c       =>   avt_c       !: Kz criterion for the turbocline depth
[281]98
[186]99END MODULE oce_trc
Note: See TracBrowser for help on using the repository browser.