1 | !!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> |
---|
2 | !! NEMO/TOP1 : 1 - tracer definition (namtoptrc) |
---|
3 | !! namelists 2 - dynamical tracer trends (namtoptrd) |
---|
4 | !! 6 - tracer advection (namtopadv) |
---|
5 | !! 7 - tracer bottom boundary (namtopbbl) |
---|
6 | !! 8 - tracer lateral diffusion (namtopldf) |
---|
7 | !! 3 - tracer vertical physics (namtopzdf) |
---|
8 | !! 9 - tracer newtonian damping (namtopdmp) |
---|
9 | !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> |
---|
10 | !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' |
---|
11 | &namtoptrc ! tracers definition |
---|
12 | !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, |
---|
13 | ndttrc = 1 ! time step frequency for passive tracers |
---|
14 | nwritetrc = 5475 ! time step frequency for tracer outputs |
---|
15 | ln_rsttr = .false. ! start from a restart file (T) or not (F) |
---|
16 | nrsttr = 0 ! restart control = 0 initial time step is not compared to the restart file value |
---|
17 | ! = 1 do not use the value in the restart file |
---|
18 | ! = 2 calendar parameters read in the restart file |
---|
19 | cn_trcrst_in = "restart_trc" ! suffix of pass. tracer restart name (input) |
---|
20 | cn_trcrst_out = "restart_trc" ! suffix of pass. tracer restart name (output) |
---|
21 | ! |
---|
22 | ! ! name ! title of the field ! units ! initial data ! save ! |
---|
23 | ! ! ! ! ! from file ! or not ! |
---|
24 | ! ! ! ! ! or not ! ! |
---|
25 | tracer(1) = 'DIC ' , 'Dissolved inorganic Concentration ', 'molC/L' , .true. , .true. |
---|
26 | tracer(2) = 'Alkalini' , 'Total Alkalinity Concentration ', 'eq/L ' , .true. , .true. |
---|
27 | tracer(3) = 'O2 ' , 'Dissolved Oxygen Concentration ', 'molO2/L' , .true. , .true. |
---|
28 | tracer(4) = 'CaCO3 ' , 'Calcite Concentration ', 'molC/L' , .false. , .true. |
---|
29 | tracer(5) = 'PO4 ' , 'Phosphate Concentration ', 'molC/L' , .true. , .true. |
---|
30 | tracer(6) = 'POC ' , 'Small organic carbon Concentration ', 'molC/L' , .false. , .true. |
---|
31 | tracer(7) = 'Si ' , 'Silicate Concentration ', 'molSi/L' , .true. , .true. |
---|
32 | tracer(8) = 'PHY ' , 'Nanophytoplankton Concentration ', 'molC/L' , .false. , .true. |
---|
33 | tracer(9) = 'ZOO ' , 'Microzooplankton Concentration ', 'molC/L' , .false. , .true. |
---|
34 | tracer(10) = 'DOC ' , 'Dissolved organic Concentration ', 'molC/L' , .false. , .true. |
---|
35 | tracer(11) = 'PHY2 ' , 'Diatoms Concentration ', 'molC/L' , .false. , .true. |
---|
36 | tracer(12) = 'ZOO2 ' , 'Mesozooplankton Concentration ', 'molC/L' , .false. , .true. |
---|
37 | tracer(13) = 'BSi ' , 'Diatoms Silicate Concentration ', 'molC/L' , .false. , .true. |
---|
38 | tracer(14) = 'Fer ' , 'Dissolved Iron Concentration ', 'molFe/L' , .true. , .true. |
---|
39 | tracer(15) = 'BFe ' , 'Big iron particles Concentration ', 'molFe/L' , .false. , .true. |
---|
40 | tracer(16) = 'GOC ' , 'Big organic carbon Concentration ', 'molC/L' , .false. , .true. |
---|
41 | tracer(17) = 'SFe ' , 'Small iron particles Concentration ', 'molFe/L' , .false. , .true. |
---|
42 | tracer(18) = 'DFe ' , 'Diatoms iron Concentration ', 'molFe/L' , .false. , .true. |
---|
43 | tracer(19) = 'DSi ' , 'Sinking biogenic Silicate Concentration', 'molC/L' , .false. , .true. |
---|
44 | tracer(20) = 'NFe ' , 'Nano iron Concentration ', 'molFe/L' , .false. , .true. |
---|
45 | tracer(21) = 'NCHL ' , 'Nano chlorophyl Concentration ', 'gChl/L' , .false. , .true. |
---|
46 | tracer(22) = 'DCHL ' , 'Diatoms chlorophyl Concentration ', 'gChl/L' , .false. , .true. |
---|
47 | tracer(23) = 'NO3 ' , 'Nitrates Concentration ', 'molC/L' , .true. , .true. |
---|
48 | tracer(24) = 'NH4 ' , 'Ammonium Concentration ', 'molC/L' , .false. , .true. |
---|
49 | / |
---|
50 | !----------------------------------------------------------------------- |
---|
51 | &namtopadv ! advection scheme for passive tracer |
---|
52 | !----------------------------------------------------------------------- |
---|
53 | ln_trcadv_cen2 = .false. ! 2nd order centered scheme |
---|
54 | ln_trcadv_tvd = .false. ! TVD scheme |
---|
55 | ln_trcadv_muscl = .true. ! MUSCL scheme |
---|
56 | ln_trcadv_muscl2 = .false. ! MUSCL2 scheme + cen2 at boundaries |
---|
57 | ln_trcadv_smolar = .false. ! SMOLAR scheme |
---|
58 | rsc = 1. ! tuning coefficient for Smol-Car. scheme |
---|
59 | ncortrc = 1 ! number of corrective phases for Smol-Car. scheme |
---|
60 | crosster = .false. ! computes Smol-Car crossterms (T) or not (F) |
---|
61 | / |
---|
62 | !----------------------------------------------------------------------- |
---|
63 | &namtopbbl ! bottom boundary layer scheme for passive tracer |
---|
64 | !----------------------------------------------------------------------- |
---|
65 | atrcbbl = 1000. ! lateral mixing coefficient in the bbl [m2/s] |
---|
66 | / |
---|
67 | !----------------------------------------------------------------------- |
---|
68 | &namtopldf ! lateral diffusion scheme for passive tracer |
---|
69 | !----------------------------------------------------------------------- |
---|
70 | ln_trcldf_diff = .true. ! performs lateral diffusion (T) or not (F) |
---|
71 | ! ! Type of the operator : |
---|
72 | ln_trcldf_lap = .true. ! laplacian operator |
---|
73 | ln_trcldf_bilap = .false. ! bilaplacian operator |
---|
74 | ! Direction of action : |
---|
75 | ln_trcldf_level = .false. ! iso-level |
---|
76 | ln_trcldf_hor = .false. ! horizontal (geopotential) (require "key_ldfslp" when ln_sco=T) |
---|
77 | ln_trcldf_iso = .true. ! iso-neutral (require "key_ldfslp") |
---|
78 | ! ! Coefficient |
---|
79 | ahtrc0 = 2000. ! horizontal eddy diffusivity for tracers [m2/s] |
---|
80 | ahtrb0 = 0. ! background eddy diffusivity for ldf_iso [m2/s] |
---|
81 | aeivtr0 = 2000. ! eddy induced velocity coefficient [m2/s] (require "key_trcldf_eiv") |
---|
82 | trcrat = 1. ! ratio betweeen passive and active tracer diffusion coeff |
---|
83 | / |
---|
84 | !----------------------------------------------------------------------- |
---|
85 | &namtopzdf ! vertical physics |
---|
86 | !----------------------------------------------------------------------- |
---|
87 | ln_trczdf_exp = .false. ! split explicit (T) or implicit (F) time stepping |
---|
88 | n_trczdf_exp = 3 ! number of sub-timestep for ln_trczdfexp=T |
---|
89 | / |
---|
90 | !----------------------------------------------------------------------- |
---|
91 | &namtoprad ! treatment of negative concentrations |
---|
92 | !----------------------------------------------------------------------- |
---|
93 | ln_trcrad = .true. ! artificially correct negative concentrations (T) or not (F) |
---|
94 | / |
---|
95 | !----------------------------------------------------------------------- |
---|
96 | &namtopdmp ! passive tracer newtonian damping ('key_trcdmp') |
---|
97 | !----------------------------------------------------------------------- |
---|
98 | ndmptr = 20 ! type of damping in passive tracers |
---|
99 | ! ='latitude', damping poleward of 'ndmp' degrees and function |
---|
100 | ! of the distance-to-coast. Red and Med Seas as ndmptr=-1 |
---|
101 | ! =-1 damping only in Med and Red Seas |
---|
102 | ndmpftr = 0 ! create a damping.coeff NetCDF file (=1) or not (=0) |
---|
103 | nmldmptr = 1 ! type of damping: =0 damping throughout the water column |
---|
104 | ! =1 no damping in the mixed layer defined by avt >5cm2/s ) |
---|
105 | ! =2 no damping in the mixed layer defined rho<rho(surf)+.01 ) |
---|
106 | sdmptr = 50. ! surface time scale for internal damping (days) |
---|
107 | bdmptr = 360. ! bottom time scale for internal damping (days) |
---|
108 | hdmptr = 800. ! depth of transition between sdmptr and bdmptr (meters) |
---|
109 | / |
---|
110 | !----------------------------------------------------------------------- |
---|
111 | &namtoptrd ! diagnostics on tracer trends |
---|
112 | ! or mixed-layer trends ('key_trdmld_trc') |
---|
113 | !---------------------------------------------------------------------- |
---|
114 | ntrd_trc = 5475 ! time step frequency and tracers trends |
---|
115 | nctls_trc = 0 ! control surface type in mixed-layer trends (0,1 or n<jpk) |
---|
116 | ucf_trc = 1 ! unit conversion factor (=1 -> /seconds ; =86400. -> /day) |
---|
117 | ln_trdmld_trc_restart = .false. ! restart for ML diagnostics |
---|
118 | ln_trdmld_trc_instant = .true. ! flag to diagnose trends of instantantaneous or mean ML T/S |
---|
119 | luttrd(1) = .true. |
---|
120 | luttrd(2) = .true. |
---|
121 | luttrd(3) = .false. |
---|
122 | luttrd(4) = .false. |
---|
123 | luttrd(5) = .false. |
---|
124 | luttrd(6) = .false. |
---|
125 | luttrd(7) = .false. |
---|
126 | luttrd(8) = .false. |
---|
127 | luttrd(9) = .false. |
---|
128 | luttrd(10) = .false. |
---|
129 | luttrd(11) = .false. |
---|
130 | luttrd(12) = .false. |
---|
131 | luttrd(13) = .false. |
---|
132 | luttrd(14) = .false. |
---|
133 | luttrd(15) = .false. |
---|
134 | luttrd(16) = .false. |
---|
135 | luttrd(17) = .false. |
---|
136 | luttrd(18) = .false. |
---|
137 | luttrd(19) = .false. |
---|
138 | luttrd(20) = .false. |
---|
139 | luttrd(21) = .false. |
---|
140 | luttrd(22) = .false. |
---|
141 | luttrd(23) = .true. |
---|
142 | luttrd(24) = .false. |
---|
143 | / |
---|