Changeset 5836 for trunk/NEMOGCM/NEMO/TOP_SRC/oce_trc.F90
- Timestamp:
- 2015-10-26T15:49:40+01:00 (9 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMOGCM/NEMO/TOP_SRC/oce_trc.F90
r5385 r5836 11 11 !! 'key_top' TOP models 12 12 !!---------------------------------------------------------------------- 13 14 ! * Domain size *13 ! 14 ! !* Domain size * 15 15 USE par_oce , ONLY : jpi => jpi !: first dimension of grid --> i 16 16 USE par_oce , ONLY : jpj => jpj !: second dimension of grid --> j … … 20 20 USE par_oce , ONLY : jpkm1 => jpkm1 !: jpk - 1 21 21 USE par_oce , ONLY : jpij => jpij !: jpi x jpj 22 USE par_oce , ONLY : lk_esopa => lk_esopa !: flag to activate the all option23 22 USE par_oce , ONLY : jp_tem => jp_tem !: indice for temperature 24 23 USE par_oce , ONLY : jp_sal => jp_sal !: indice for salinity 25 24 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 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 USE dom_oce !* model domain * 52 34 53 35 USE domvvl, ONLY : un_td, vn_td !: thickness diffusion transport … … 56 38 57 39 !* 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 40 USE oce , ONLY : un => un !: i-horizontal velocity (m s-1) 61 41 USE oce , ONLY : vn => vn !: j-horizontal velocity (m s-1) … … 66 46 USE oce , ONLY : rhop => rhop !: potential volumic mass (kg m-3) 67 47 USE oce , ONLY : rhd => rhd !: in situ density anomalie rhd=(rho-rau0)/rau0 (no units) 48 USE oce , ONLY : hdivn => hdivn !: horizontal divergence (1/s) 49 USE oce , ONLY : sshn => sshn !: sea surface height at t-point [m] 50 USE oce , ONLY : sshb => sshb !: sea surface height at t-point [m] 51 USE oce , ONLY : ssha => ssha !: sea surface height at t-point [m] 68 52 #if defined key_offline 69 53 USE oce , ONLY : rab_n => rab_n !: local thermal/haline expension ratio at T-points 70 54 #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 diffusion79 55 80 56 !* surface fluxes * … … 102 78 USE trc_oce 103 79 80 !!gm : I don't understand this as ldftra (where everything is defined) is used by TRC in all cases (ON/OFF-line) 81 !!gm so the following lines should be removed.... logical should be the one of TRC namelist 82 !!gm In case off coarsening.... the ( ahtu, ahtv, aeiu, aeiv) arrays are needed that's all. 104 83 !* 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 84 USE ldftra , ONLY : rn_aht_0 => rn_aht_0 !: laplacian lateral eddy diffusivity [m2/s] 85 USE ldftra , ONLY : rn_bht_0 => rn_bht_0 !: bilaplacian lateral eddy diffusivity [m4/s] 86 USE ldftra , ONLY : ahtu => ahtu !: lateral diffusivity coef. at u-points 87 USE ldftra , ONLY : ahtv => ahtv !: lateral diffusivity coef. at v-points 88 USE ldftra , ONLY : rn_aeiv_0 => rn_aeiv_0 !: eddy induced velocity coefficient (m2/s) 89 USE ldftra , ONLY : aeiu => aeiu !: eddy induced velocity coef. at u-points (m2/s) 90 USE ldftra , ONLY : aeiv => aeiv !: eddy induced velocity coef. at v-points (m2/s) 91 USE ldftra , ONLY : ln_ldfeiv => ln_ldfeiv !: eddy induced velocity flag 92 93 !!gm this should be : ln_trcldf_triad (TRC namelist) 94 USE ldfslp , ONLY : ln_traldf_triad => ln_traldf_triad !: triad scheme (Griffies et al.) 95 96 !* direction of lateral diffusion * 97 USE ldfslp , ONLY : l_ldfslp => l_ldfslp !: slopes flag 98 USE ldfslp , ONLY : uslp => uslp !: i-slope at u-point 99 USE ldfslp , ONLY : vslp => vslp !: j-slope at v-point 100 USE ldfslp , ONLY : wslpi => wslpi !: i-slope at w-point 101 USE ldfslp , ONLY : wslpj => wslpj !: j-slope at w-point 102 !!gm end 118 103 119 104 !* vertical diffusion * … … 128 113 USE zdfmxl , ONLY : hmlp => hmlp !: mixed layer depth (rho=rho0+zdcrit) (m) 129 114 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 flag133 # if defined key_ldfslp134 USE ldfslp , ONLY : uslp => uslp !: i-direction slope at u-, w-points135 USE ldfslp , ONLY : vslp => vslp !: j-direction slope at v-, w-points136 USE ldfslp , ONLY : wslpi => wslpi !: i-direction slope at u-, w-points137 USE ldfslp , ONLY : wslpj => wslpj !: j-direction slope at v-, w-points138 # endif139 115 140 116 USE diaar5 , ONLY : lk_diaar5 => lk_diaar5
Note: See TracChangeset
for help on using the changeset viewer.