[1121] | 1 | !!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> |
---|
[2528] | 2 | !! NEMO/TOP1 : 1 - tracer definition (namtrc ) |
---|
[3294] | 3 | !! 2 - tracer data initialisation (namtrc_dta) |
---|
[2528] | 4 | !! 3 - tracer advection (namtrc_adv) |
---|
| 5 | !! 4 - tracer lateral diffusion (namtrc_ldf) |
---|
| 6 | !! 5 - tracer vertical physics (namtrc_zdf) |
---|
| 7 | !! 6 - tracer newtonian damping (namtrc_dmp) |
---|
[3294] | 8 | !! 7 - dynamical tracer trends (namtrc_trd) |
---|
| 9 | !! 8 - tracer output diagonstics (namtrc_dia) |
---|
[1086] | 10 | !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> |
---|
[1121] | 11 | !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' |
---|
[2528] | 12 | &namtrc ! tracers definition |
---|
[1121] | 13 | !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, |
---|
[2528] | 14 | nn_dttrc = 1 ! time step frequency for passive sn_tracers |
---|
| 15 | nn_writetrc = 1460 ! time step frequency for sn_tracer outputs |
---|
| 16 | ln_rsttr = .false. ! start from a restart file (T) or not (F) |
---|
| 17 | nn_rsttr = 0 ! restart control = 0 initial time step is not compared to the restart file value |
---|
[3319] | 18 | ! = 1 do not use the value in the restart file |
---|
| 19 | ! = 2 calendar parameters read in the restart file |
---|
[2528] | 20 | cn_trcrst_in = "restart_trc" ! suffix of pass. sn_tracer restart name (input) |
---|
| 21 | cn_trcrst_out = "restart_trc" ! suffix of pass. sn_tracer restart name (output) |
---|
[3319] | 22 | ln_trcdta = .true. ! Initialisation from data input file (T) or not (F) |
---|
| 23 | ln_trcdmp = .false. ! add a damping termn (T) or not (F) |
---|
[1086] | 24 | ! |
---|
[3295] | 25 | ! ! name ! title of the field ! initial data ! initial data ! save ! |
---|
| 26 | ! ! ! ! units ! from file ! or not ! |
---|
| 27 | ! ! ! ! ! or not ! ! |
---|
[2528] | 28 | sn_tracer(1) = 'DIC ' , 'Dissolved inorganic Concentration ', 'mol-C/L' , .true. , .true. |
---|
| 29 | sn_tracer(2) = 'Alkalini' , 'Total Alkalinity Concentration ', 'eq/L ' , .true. , .true. |
---|
| 30 | sn_tracer(3) = 'O2 ' , 'Dissolved Oxygen Concentration ', 'mol-C/L' , .true. , .true. |
---|
| 31 | sn_tracer(4) = 'CaCO3 ' , 'Calcite Concentration ', 'mol-C/L' , .false. , .true. |
---|
| 32 | sn_tracer(5) = 'PO4 ' , 'Phosphate Concentration ', 'mol-C/L' , .true. , .true. |
---|
| 33 | sn_tracer(6) = 'POC ' , 'Small organic carbon Concentration ', 'mol-C/L' , .false. , .true. |
---|
| 34 | sn_tracer(7) = 'Si ' , 'Silicate Concentration ', 'mol-C/L' , .true. , .true. |
---|
| 35 | sn_tracer(8) = 'PHY ' , 'Nanophytoplankton Concentration ', 'mol-C/L' , .false. , .true. |
---|
| 36 | sn_tracer(9) = 'ZOO ' , 'Microzooplankton Concentration ', 'mol-C/L' , .false. , .true. |
---|
[3319] | 37 | sn_tracer(10) = 'DOC ' , 'Dissolved organic Concentration ', 'mol-C/L' , .true. , .true. |
---|
[2528] | 38 | sn_tracer(11) = 'PHY2 ' , 'Diatoms Concentration ', 'mol-C/L' , .false. , .true. |
---|
| 39 | sn_tracer(12) = 'ZOO2 ' , 'Mesozooplankton Concentration ', 'mol-C/L' , .false. , .true. |
---|
[3295] | 40 | sn_tracer(13) = 'DSi ' , 'Diatoms Silicate Concentration ', 'mol-C/L' , .false. , .true. |
---|
[2528] | 41 | sn_tracer(14) = 'Fer ' , 'Dissolved Iron Concentration ', 'mol-C/L' , .true. , .true. |
---|
| 42 | sn_tracer(15) = 'BFe ' , 'Big iron particles Concentration ', 'mol-C/L' , .false. , .true. |
---|
| 43 | sn_tracer(16) = 'GOC ' , 'Big organic carbon Concentration ', 'mol-C/L' , .false. , .true. |
---|
| 44 | sn_tracer(17) = 'SFe ' , 'Small iron particles Concentration ', 'mol-C/L' , .false. , .true. |
---|
| 45 | sn_tracer(18) = 'DFe ' , 'Diatoms iron Concentration ', 'mol-C/L' , .false. , .true. |
---|
[3295] | 46 | sn_tracer(19) = 'GSi ' , 'Sinking biogenic Silicate Concentration', 'mol-C/L' , .false. , .true. |
---|
[2528] | 47 | sn_tracer(20) = 'NFe ' , 'Nano iron Concentration ', 'mol-C/L' , .false. , .true. |
---|
| 48 | sn_tracer(21) = 'NCHL ' , 'Nano chlorophyl Concentration ', 'mol-C/L' , .false. , .true. |
---|
| 49 | sn_tracer(22) = 'DCHL ' , 'Diatoms chlorophyl Concentration ', 'mol-C/L' , .false. , .true. |
---|
| 50 | sn_tracer(23) = 'NO3 ' , 'Nitrates Concentration ', 'mol-C/L' , .true. , .true. |
---|
| 51 | sn_tracer(24) = 'NH4 ' , 'Ammonium Concentration ', 'mol-C/L' , .false. , .true. |
---|
[1086] | 52 | / |
---|
| 53 | !----------------------------------------------------------------------- |
---|
[3294] | 54 | &namtrc_dta ! Initialisation from data input file |
---|
| 55 | !----------------------------------------------------------------------- |
---|
| 56 | ! |
---|
| 57 | ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! |
---|
| 58 | ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! |
---|
| 59 | sn_trcdta(1) = 'data_DIC_nomask' , -12 , 'DIC' , .false. , .true. , 'yearly' , '' , '' |
---|
| 60 | sn_trcdta(2) = 'data_Alkalini_nomask' , -12 , 'Alkalini', .false. , .true. , 'yearly' , '' , '' |
---|
| 61 | sn_trcdta(3) = 'data_O2_nomask' , -1 , 'O2' , .true. , .true. , 'yearly' , '' , '' |
---|
| 62 | sn_trcdta(5) = 'data_PO4_nomask' , -1 , 'PO4' , .true. , .true. , 'yearly' , '' , '' |
---|
| 63 | sn_trcdta(7) = 'data_Si_nomask' , -1 , 'Si' , .true. , .true. , 'yearly' , '' , '' |
---|
| 64 | sn_trcdta(10) = 'data_DOC_nomask' , -12 , 'DOC' , .false. , .true. , 'yearly' , '' , '' |
---|
| 65 | sn_trcdta(14) = 'data_Fer_nomask' , -12 , 'Fer' , .false. , .true. , 'yearly' , '' , '' |
---|
| 66 | sn_trcdta(23) = 'data_NO3_nomask' , -1 , 'NO3' , .true. , .true. , 'yearly' , '' , '' |
---|
| 67 | ! |
---|
| 68 | cn_dir = './' ! root directory for the location of the data files |
---|
| 69 | rn_trfac(1) = 1.0e-06 ! multiplicative factor |
---|
| 70 | rn_trfac(2) = 1.0e-06 ! - - - - |
---|
| 71 | rn_trfac(3) = 44.6e-06 ! - - - - |
---|
| 72 | rn_trfac(5) = 122.0e-06 ! - - - - |
---|
| 73 | rn_trfac(7) = 1.0e-06 ! - - - - |
---|
| 74 | rn_trfac(10) = 1.0 ! - - - - |
---|
| 75 | rn_trfac(14) = 1.0 ! - - - - |
---|
| 76 | rn_trfac(23) = 7.6e-06 ! - - - - |
---|
| 77 | / |
---|
| 78 | !----------------------------------------------------------------------- |
---|
[2528] | 79 | &namtrc_adv ! advection scheme for passive tracer |
---|
[1086] | 80 | !----------------------------------------------------------------------- |
---|
[1121] | 81 | ln_trcadv_cen2 = .false. ! 2nd order centered scheme |
---|
[2528] | 82 | ln_trcadv_tvd = .false. ! TVD scheme |
---|
[1121] | 83 | ln_trcadv_muscl = .true. ! MUSCL scheme |
---|
| 84 | ln_trcadv_muscl2 = .false. ! MUSCL2 scheme + cen2 at boundaries |
---|
[2528] | 85 | ln_trcadv_ubs = .false. ! UBS scheme |
---|
| 86 | ln_trcadv_qck = .false. ! QUICKEST scheme |
---|
[1086] | 87 | / |
---|
| 88 | !----------------------------------------------------------------------- |
---|
[2528] | 89 | &namtrc_ldf ! lateral diffusion scheme for passive tracer |
---|
[1086] | 90 | !----------------------------------------------------------------------- |
---|
[1121] | 91 | ln_trcldf_diff = .true. ! performs lateral diffusion (T) or not (F) |
---|
| 92 | ! ! Type of the operator : |
---|
| 93 | ln_trcldf_lap = .true. ! laplacian operator |
---|
| 94 | ln_trcldf_bilap = .false. ! bilaplacian operator |
---|
| 95 | ! Direction of action : |
---|
| 96 | ln_trcldf_level = .false. ! iso-level |
---|
| 97 | ln_trcldf_hor = .false. ! horizontal (geopotential) (require "key_ldfslp" when ln_sco=T) |
---|
| 98 | ln_trcldf_iso = .true. ! iso-neutral (require "key_ldfslp") |
---|
| 99 | ! ! Coefficient |
---|
[3294] | 100 | rn_ahtrc_0 = 2000. ! horizontal eddy diffusivity for tracers [m2/s] |
---|
[2528] | 101 | rn_ahtrb_0 = 0. ! background eddy diffusivity for ldf_iso [m2/s] |
---|
[1086] | 102 | / |
---|
| 103 | !----------------------------------------------------------------------- |
---|
[2528] | 104 | &namtrc_zdf ! vertical physics |
---|
[1086] | 105 | !----------------------------------------------------------------------- |
---|
[1121] | 106 | ln_trczdf_exp = .false. ! split explicit (T) or implicit (F) time stepping |
---|
[2528] | 107 | nn_trczdf_exp = 3 ! number of sub-timestep for ln_trczdfexp=T |
---|
[1086] | 108 | / |
---|
| 109 | !----------------------------------------------------------------------- |
---|
[2528] | 110 | &namtrc_rad ! treatment of negative concentrations |
---|
[1259] | 111 | !----------------------------------------------------------------------- |
---|
| 112 | ln_trcrad = .true. ! artificially correct negative concentrations (T) or not (F) |
---|
| 113 | / |
---|
| 114 | !----------------------------------------------------------------------- |
---|
[3294] | 115 | &namtrc_dmp ! passive tracer newtonian damping |
---|
[1086] | 116 | !----------------------------------------------------------------------- |
---|
[2528] | 117 | nn_hdmp_tr = -1 ! horizontal shape =-1, damping in Med and Red Seas only |
---|
| 118 | ! =XX, damping poleward of XX degrees (XX>0) |
---|
| 119 | ! + F(distance-to-coast) + Red and Med Seas |
---|
| 120 | nn_zdmp_tr = 1 ! vertical shape =0 damping throughout the water column |
---|
| 121 | ! =1 no damping in the mixing layer (kz criteria) |
---|
| 122 | ! =2 no damping in the mixed layer (rho crieria) |
---|
| 123 | rn_surf_tr = 50. ! surface time scale of damping [days] |
---|
| 124 | rn_bot_tr = 360. ! bottom time scale of damping [days] |
---|
| 125 | rn_dep_tr = 800. ! depth of transition between rn_surf and rn_bot [meters] |
---|
| 126 | nn_file_tr = 0 ! create a damping.coeff NetCDF file (=1) or not (=0) |
---|
[1086] | 127 | / |
---|
[1183] | 128 | !----------------------------------------------------------------------- |
---|
[2528] | 129 | &namtrc_trd ! diagnostics on tracer trends ('key_trdtrc') |
---|
[1183] | 130 | ! or mixed-layer trends ('key_trdmld_trc') |
---|
| 131 | !---------------------------------------------------------------------- |
---|
[3294] | 132 | nn_trd_trc = 1460 ! time step frequency and tracers trends |
---|
[2528] | 133 | nn_ctls_trc = 0 ! control surface type in mixed-layer trends (0,1 or n<jpk) |
---|
| 134 | rn_ucf_trc = 1 ! unit conversion factor (=1 -> /seconds ; =86400. -> /day) |
---|
[1183] | 135 | ln_trdmld_trc_restart = .false. ! restart for ML diagnostics |
---|
| 136 | ln_trdmld_trc_instant = .true. ! flag to diagnose trends of instantantaneous or mean ML T/S |
---|
[2528] | 137 | ln_trdtrc(1) = .true. |
---|
| 138 | ln_trdtrc(2) = .true. |
---|
| 139 | ln_trdtrc(23) = .true. |
---|
[1183] | 140 | / |
---|
[3294] | 141 | !----------------------------------------------------------------------- |
---|
| 142 | &namtrc_dia ! parameters for passive tracer additional diagnostics |
---|
| 143 | !---------------------------------------------------------------------- |
---|
| 144 | ln_diatrc = .true. ! save additional diag. (T) or not (F) |
---|
| 145 | nn_writedia = 1460 ! time step frequency for diagnostics |
---|
| 146 | / |
---|