| 1 | !!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> |
|---|
| 2 | !! NEMO/TOP1 : 1 - tracer definition (namtrc ) |
|---|
| 3 | !! 2 - tracer data initialisation (namtrc_dta) |
|---|
| 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) |
|---|
| 8 | !! 7 - dynamical tracer trends (namtrc_trd) |
|---|
| 9 | !! 8 - tracer output diagonstics (namtrc_dia) |
|---|
| 10 | !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> |
|---|
| 11 | !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' |
|---|
| 12 | &namtrc ! tracers definition |
|---|
| 13 | !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, |
|---|
| 14 | nn_dttrc = 1 ! time step frequency for passive sn_tracers |
|---|
| 15 | nn_writetrc = 360 ! time step frequency for sn_tracer outputs |
|---|
| 16 | ln_top_euler = .true. ! Euler timestepping |
|---|
| 17 | ln_rsttr = .false. ! start from a restart file (T) or not (F) |
|---|
| 18 | nn_rsttr = 0 ! restart control = 0 initial time step is not compared to the restart file value |
|---|
| 19 | ! = 1 do not use the value in the restart file |
|---|
| 20 | ! = 2 calendar parameters read in the restart file |
|---|
| 21 | cn_trcrst_in = "restart_trc" ! suffix of pass. sn_tracer restart name (input) |
|---|
| 22 | cn_trcrst_out = "restart_trc" ! suffix of pass. sn_tracer restart name (output) |
|---|
| 23 | ln_trcdta = .false. ! Initialisation from data input file (T) or not (F) |
|---|
| 24 | ! |
|---|
| 25 | ! ! name ! title of the field ! initial data ! initial data ! save ! |
|---|
| 26 | ! ! ! ! units ! from file ! or not ! |
|---|
| 27 | ! ! ! ! ! or not ! ! |
|---|
| 28 | sn_tracer(1) = 'DUMMY' , 'Dummy tracer' , 'dummy units', .false. , .true. |
|---|
| 29 | / |
|---|
| 30 | !----------------------------------------------------------------------- |
|---|
| 31 | &namtrc_dta ! Initialisation from data input file |
|---|
| 32 | !----------------------------------------------------------------------- |
|---|
| 33 | ! |
|---|
| 34 | ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! |
|---|
| 35 | ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! |
|---|
| 36 | sn_trcdta(1) = 'data_1m_O2_nomask' , -1 , 'O2' , .false. , .true. , 'yearly' , '' , '' |
|---|
| 37 | sn_trcdta(2) = 'data_1m_PO4_nomask' , -1 , 'PO4' , .false. , .true. , 'yearly' , '' , '' |
|---|
| 38 | sn_trcdta(3) = 'data_1m_NO3_nomask' , -1 , 'NO3' , .false. , .true. , 'yearly' , '' , '' |
|---|
| 39 | sn_trcdta(6) = 'data_1m_Si_nomask' , -1 , 'Si' , .false. , .true. , 'yearly' , '' , '' |
|---|
| 40 | sn_trcdta(43) = 'data_1m_DOC_nomask' , -12 , 'DOC' , .false. , .true. , 'yearly' , '' , '' |
|---|
| 41 | sn_trcdta(49) = 'data_1m_DIC_nomask' , -12 , 'DIC' , .false. , .true. , 'yearly' , '' , '' |
|---|
| 42 | sn_trcdta(50) = 'data_1m_Alkalini_nomask', -12 , 'Alkalini', .false. , .true. , 'yearly' , '' , '' |
|---|
| 43 | ! |
|---|
| 44 | cn_dir = '/users/home/ans017/WORKS/ORCA2_INPUT/' ! root directory for the location of the data files |
|---|
| 45 | rn_trfac(1) = 22.4 ! conversion factor from ml/l to mmol m3 (1 mol O2 = 22.4 l) |
|---|
| 46 | rn_trfac(2) = 1.0 ! - - - - |
|---|
| 47 | rn_trfac(3) = 1.0 ! - - - - |
|---|
| 48 | rn_trfac(6) = 1.0 ! - - - - |
|---|
| 49 | rn_trfac(43) = 12.e6 ! conversion factor from mol/l to mg/m3 |
|---|
| 50 | rn_trfac(49) = 12.0 ! conversion factor from umol/l to mg/m3 |
|---|
| 51 | rn_trfac(50) = 1.0 ! - - - - |
|---|
| 52 | / |
|---|
| 53 | !----------------------------------------------------------------------- |
|---|
| 54 | &namtrc_bc |
|---|
| 55 | !----------------------------------------------------------------------- |
|---|
| 56 | ! ! file name ! frequency (hr) ! variable ! time interp. ! clim !'yearly' or ! weights ! rotation |
|---|
| 57 | ! ! ! (if <0 mon) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing |
|---|
| 58 | sn_trcsbc(3) = 'ndeposition.orca' , -1 , 'ndep' , .false. , .true. , 'yearly' , '' , '' |
|---|
| 59 | |
|---|
| 60 | rn_trsfac(3) = 1.0 |
|---|
| 61 | ! |
|---|
| 62 | cn_dir = '/users/home/ans017/WORKS/ORCA2_INPUT/' ! root directory for the location of the boundary condition files |
|---|
| 63 | / |
|---|
| 64 | !----------------------------------------------------------------------- |
|---|
| 65 | &namtrc_adv ! advection scheme for passive tracer |
|---|
| 66 | !----------------------------------------------------------------------- |
|---|
| 67 | ln_trcadv_cen2 = .false. ! 2nd order centered scheme |
|---|
| 68 | ln_trcadv_tvd = .false. ! TVD scheme |
|---|
| 69 | ln_trcadv_muscl = .true. ! MUSCL scheme |
|---|
| 70 | ln_trcadv_muscl2 = .false. ! MUSCL2 scheme + cen2 at boundaries |
|---|
| 71 | ln_trcadv_ubs = .false. ! UBS scheme |
|---|
| 72 | ln_trcadv_qck = .false. ! QUICKEST scheme |
|---|
| 73 | / |
|---|
| 74 | !----------------------------------------------------------------------- |
|---|
| 75 | &namtrc_ldf ! lateral diffusion scheme for passive tracer |
|---|
| 76 | !----------------------------------------------------------------------- |
|---|
| 77 | ln_trcldf_diff = .true. ! performs lateral diffusion (T) or not (F) |
|---|
| 78 | ! ! Type of the operator : |
|---|
| 79 | ln_trcldf_lap = .true. ! laplacian operator |
|---|
| 80 | ln_trcldf_bilap = .false. ! bilaplacian operator |
|---|
| 81 | ! Direction of action : |
|---|
| 82 | ln_trcldf_level = .false. ! iso-level |
|---|
| 83 | ln_trcldf_hor = .false. ! horizontal (geopotential) (require "key_ldfslp" when ln_sco=T) |
|---|
| 84 | ln_trcldf_iso = .true. ! iso-neutral (require "key_ldfslp") |
|---|
| 85 | ! ! Coefficient |
|---|
| 86 | rn_ahtrc_0 = 2000. ! horizontal eddy diffusivity for tracers [m2/s] |
|---|
| 87 | rn_ahtrb_0 = 0. ! background eddy diffusivity for ldf_iso [m2/s] |
|---|
| 88 | / |
|---|
| 89 | !----------------------------------------------------------------------- |
|---|
| 90 | &namtrc_zdf ! vertical physics |
|---|
| 91 | !----------------------------------------------------------------------- |
|---|
| 92 | ln_trczdf_exp = .false. ! split explicit (T) or implicit (F) time stepping |
|---|
| 93 | nn_trczdf_exp = 3 ! number of sub-timestep for ln_trczdfexp=T |
|---|
| 94 | / |
|---|
| 95 | !----------------------------------------------------------------------- |
|---|
| 96 | &namtrc_rad ! treatment of negative concentrations |
|---|
| 97 | !----------------------------------------------------------------------- |
|---|
| 98 | ln_trcrad = .true. ! artificially correct negative concentrations (T) or not (F) |
|---|
| 99 | / |
|---|
| 100 | !----------------------------------------------------------------------- |
|---|
| 101 | &namtrc_dmp ! passive tracer newtonian damping |
|---|
| 102 | !----------------------------------------------------------------------- |
|---|
| 103 | ln_trcdmp = .false. ! add a damping term (T) or not (F) |
|---|
| 104 | nn_hdmp_tr = -1 ! horizontal shape =-1, damping in Med and Red Seas only |
|---|
| 105 | ! =XX, damping poleward of XX degrees (XX>0) |
|---|
| 106 | ! + F(distance-to-coast) + Red and Med Seas |
|---|
| 107 | nn_zdmp_tr = 1 ! vertical shape =0 damping throughout the water column |
|---|
| 108 | ! =1 no damping in the mixing layer (kz criteria) |
|---|
| 109 | ! =2 no damping in the mixed layer (rho crieria) |
|---|
| 110 | rn_surf_tr = 50. ! surface time scale of damping [days] |
|---|
| 111 | rn_bot_tr = 360. ! bottom time scale of damping [days] |
|---|
| 112 | rn_dep_tr = 800. ! depth of transition between rn_surf and rn_bot [meters] |
|---|
| 113 | nn_file_tr = 0 ! create a damping.coeff NetCDF file (=1) or not (=0) |
|---|
| 114 | / |
|---|
| 115 | !----------------------------------------------------------------------- |
|---|
| 116 | &namtrc_trd ! diagnostics on tracer trends ('key_trdtrc') |
|---|
| 117 | ! or mixed-layer trends ('key_trdmld_trc') |
|---|
| 118 | !---------------------------------------------------------------------- |
|---|
| 119 | nn_trd_trc = 360 ! time step frequency and tracers trends |
|---|
| 120 | nn_ctls_trc = 0 ! control surface type in mixed-layer trends (0,1 or n<jpk) |
|---|
| 121 | rn_ucf_trc = 86400 ! unit conversion factor (=1 -> /seconds ; =86400. -> /day) |
|---|
| 122 | ln_trdmld_trc_restart = .false. ! restart for ML diagnostics |
|---|
| 123 | ln_trdmld_trc_instant = .false. ! flag to diagnose trends of instantantaneous or mean ML T/S |
|---|
| 124 | ln_trdtrc(1) = .true. |
|---|
| 125 | ln_trdtrc(2) = .true. |
|---|
| 126 | ln_trdtrc(23) = .true. |
|---|
| 127 | / |
|---|
| 128 | !----------------------------------------------------------------------- |
|---|
| 129 | &namtrc_dia ! parameters for passive tracer additional diagnostics |
|---|
| 130 | !---------------------------------------------------------------------- |
|---|
| 131 | ln_diatrc = .false. ! save additional diag. (T) or not (F) |
|---|
| 132 | ln_diabio = .false. |
|---|
| 133 | nn_writedia = 5475 ! time step frequency for diagnostics |
|---|
| 134 | nn_writebio = 5475 |
|---|
| 135 | / |
|---|
