CNRS - NERC - Met OFFICE - MERCATOR-ocean - INGV - CMCC NEMO team Ocean General Circulation Model version 3.4 (2011) AAAAAAAA nemo_ctl: Control prints & Benchmark ~~~~~~~ Namelist namctl run control (for debugging) ln_ctl = F level of print nn_print = 0 Start i indice for SUM control nn_ictls = 0 End i indice for SUM control nn_ictle = 0 Start j indice for SUM control nn_jctls = 0 End j indice for SUM control nn_jctle = 0 number of proc. following i nn_isplt = 1 number of proc. following j nn_jsplt = 1 benchmark parameter (0/1) nn_bench = 0 mpp_init(2) : NO massively parallel processing ~~~~~~~~~~~: nperio = 4 npolj = 4 nimpp = 1 njmpp = 1 file : timing.output open ok unit = 18 status = REPLACE form = FORMATTED access = SEQUENTIAL phy_cst : initialization of ocean parameters and constants ~~~~~~~ Domain info dimension of model Local domain Global domain Data domain jpi : 182 jpiglo : 182 jpidta : 182 jpj : 149 jpjglo : 149 jpjdta : 149 jpk : 31 jpk : 31 jpkdta : 31 jpij : 27118 mpp local domain info (mpp) jpni : 1 jpreci : 1 jpnj : 1 jprecj : 1 jpnij : 1 lateral domain boundary condition type : jperio = 4 Constants mathematical constant rpi = 3.1415926535897931 day rday = 86400.000000000000 s sideral year rsiyea = 31558149.010110732 s sideral day rsiday = 86164.099655911836 s omega omega = 7.29211508304606178E-005 s-1 nb of months per year raamo = 12.000000000000000 months nb of hours per day rjjhh = 24.000000000000000 hours nb of minutes per hour rhhmm = 60.000000000000000 mn nb of seconds per minute rmmss = 60.000000000000000 s earth radius ra = 6371229.0000000000 m gravity grav = 9.8066499999999994 m/s^2 triple point of temperature rtt = 273.16000000000003 K freezing point of water rt0 = 273.14999999999998 K melting point of snow rt0_snow = 273.14999999999998 K melting point of ice rt0_ice = 273.05000000000001 K ocean reference volumic mass rau0 = 1035.0000000000000 kg/m^3 ocean reference specific volume rau0r = 9.66183574879227029E-004 m^3/Kg ocean specific heat rcp = 4000.0000000000000 1. / ( rau0 * rcp ) = ro0cpr = 2.41545893719806766E-007 thermal conductivity of the snow = 0.22000000000000000 J/s/m/K thermal conductivity of the ice = 2.0343960000000001 J/s/m/K density times specific heat for snow = 690690.00000000000 J/m^3/K density times specific heat for ice = 1883700.0000000000 J/m^3/K volumetric latent heat fusion of sea ice = 300330000.00000000 J/m volumetric latent heat fusion of snow = 110121000.00000000 J/m latent heat of sublimation of snow = 2800000.0000000000 J/kg density of sea ice = 900.00000000000000 kg/m^3 density of snow = 330.00000000000000 kg/m^3 emissivity of snow or ice = 0.96999999999999997 salinity of ice = 6.0000000000000000 psu salinity of sea = 34.700000000000003 psu latent heat of evaporation (water) = 2500000.0000000000 J/m^3 correction factor for solar radiation = 0.90000000000000002 von Karman constant = 0.40000000000000002 Stefan-Boltzmann constant = 5.66999999999999982E-008 J/s/m^2/K^4 conversion: degre ==> radian rad = 1.74532925199432955E-002 smallest real computer value rsmall = 1.11022302462515654E-016 eos_init : equation of state ~~~~~~~~ Namelist nameos : set eos parameters flag for eq. of state and N^2 nn_eos = 0 thermal exp. coef. (linear) rn_alpha = 2.00000000000000010E-004 saline exp. coef. (linear) rn_beta = 7.69999999999999960E-004 use of Jackett & McDougall (1994) equation of state and McDougall (1987) Brunt-Vaisala frequency dom_cfg : set the ocean configuration ~~~~~~~ ocean model configuration used : cp_cfg = orca jp_cfg = 2 global domain lateral boundaries jperio= 4, cyclic east-west and north fold with T-point pivot dom_glo : domain: data / local ~~~~~~~ data input domain : jpidta = 182 jpjdta = 149 jpkdta = 31 global or zoom domain: jpiglo = 182 jpjglo = 149 jpk = 31 local domain : jpi = 182 jpj = 149 jpk = 31 south-west indices jpizoom = 1 jpjzoom = 1 conversion local ==> data i-index domain 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 172 173 174 175 176 177 178 179 180 181 182 conversion data ==> local i-index domain starting index 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 172 173 174 175 176 177 178 179 180 181 182 ending index 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 172 173 174 175 176 177 178 179 180 181 182 conversion local ==> data j-index domain 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 conversion data ==> local j-index domain starting index 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 ending index 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 zoom flags : lzoom = F (T = zoom, F = global ) lzoom_e = F (T = forced closed east boundary) lzoom_w = F (T = forced closed west boundary) lzoom_s = F (T = forced closed South boundary) lzoom_n = F (T = forced closed North boundary) ORCA configuration: antarctic/arctic zoom flags : lzoom_arct = F (T= arctic zoom, F=global) lzoom_anta = F (T=antarctic zoom, F=global) dom_init : domain initialization ~~~~~~~~ dom_nam : domain initialization through namelist read ~~~~~~~ Namelist namrun job number nn_no = 0 experiment name for output cn_exp = O2LP_SHORT restart logical ln_rstart = T control of time step nn_rstctl = 2 number of the first time step nn_it000 = 76 number of the last time step nn_itend = 150 initial calendar date aammjj nn_date0 = 10101 leap year calendar (0/1) nn_leapy = 0 initial state output nn_istate = 0 frequency of restart file nn_stock = 75 frequency of output file nn_write = 5475 multi file dimgout ln_dimgnnn = F mask land points ln_mskland = F overwrite an existing file ln_clobber = T NetCDF chunksize (bytes) nn_chunksz = 0 The IOIPSL calendar is "noleap", i.e. no leap year Namelist namdom : space & time domain flag read/compute bathymetry nn_bathy = 1 min depth of the ocean (>0) or rn_hmin = -3.0000000000000000 min number of ocean level (<0) minimum thickness of partial rn_e3zps_min = 20.000000000000000 (m) step level rn_e3zps_rat = 0.10000000000000001 create mesh/mask file(s) nn_msh = 0 = 0 no file created = 1 mesh_mask = 2 mesh and mask = 3 mesh_hgr, msh_zgr and mask ocean time step rn_rdt = 5760.0000000000000 asselin time filter parameter rn_atfp = 0.10000000000000001 time-splitting: nb of sub time-step nn_baro = 64 acceleration of converge nn_acc = 0 nn_acc=1: surface tracer rdt rn_rdtmin = 28800.000000000000 bottom tracer rdt rdtmax = 28800.000000000000 depth of transition rn_rdth = 800.00000000000000 suppression of closed seas (=0) nn_closea = 0 Namelist namcla cross land advection nn_cla = 0 dom_clo : closed seas ~~~~~~~ dom_hgr : define the horizontal mesh from ithe following par_oce parameters ~~~~~~~ type of horizontal mesh jphgr_msh = 0 position of the first row and ppglam0 = 999999.00000000000 column grid-point (degrees) ppgphi0 = 999999.00000000000 zonal grid-spacing (degrees) ppe1_deg = 999999.00000000000 meridional grid-spacing (degrees) ppe2_deg = 999999.00000000000 zonal grid-spacing (meters) ppe1_m = 999999.00000000000 meridional grid-spacing (meters) ppe2_m = 999999.00000000000 curvilinear coordinate on the sphere read in "coordinate" file hgr_read : read the horizontal coordinates ~~~~~~~~ jpiglo = 182 jpjglo = 149 jpk = 31 iom_nf90_open ~~~ open existing file: coordinates.nc in READ mode ---> coordinates.nc OK read glamt (rec: 1) in coordinates.nc ok read glamu (rec: 1) in coordinates.nc ok read glamv (rec: 1) in coordinates.nc ok read glamf (rec: 1) in coordinates.nc ok read gphit (rec: 1) in coordinates.nc ok read gphiu (rec: 1) in coordinates.nc ok read gphiv (rec: 1) in coordinates.nc ok read gphif (rec: 1) in coordinates.nc ok read e1t (rec: 1) in coordinates.nc ok read e1u (rec: 1) in coordinates.nc ok read e1v (rec: 1) in coordinates.nc ok read e1f (rec: 1) in coordinates.nc ok read e2t (rec: 1) in coordinates.nc ok read e2u (rec: 1) in coordinates.nc ok read e2v (rec: 1) in coordinates.nc ok read e2f (rec: 1) in coordinates.nc ok iom_close ~~~ close file: coordinates.nc ok orca_r2: Gibraltar : e2u reduced to 20 km orca_r2: Bab el Mandeb: e2u reduced to 30 km e1v reduced to 18 km orca_r2: Danish Straits : e2u reduced to 10 km dom_zgr : vertical coordinate ~~~~~~~ Namelist namzgr : set vertical coordinate z-coordinate - full steps ln_zco = F z-coordinate - partial steps ln_zps = T s- or hybrid z-s-coordinate ln_sco = F zgr_z : Reference vertical z-coordinates ~~~~~~~ Value of coefficients for vertical mesh: zsur = -4762.9614354630003 za0 = 255.58049070440001 za1 = 245.58132232489999 zkth = 21.433361979379999 zacr = 3.0000000000000000 Reference z-coordinate depth and scale factors: level gdept gdepw e3t e3w 1 5.00 0.00 10.00 10.00 2 15.00 10.00 10.00 10.00 3 25.00 20.00 10.00 10.00 4 35.01 30.00 10.01 10.00 5 45.01 40.01 10.01 10.01 6 55.03 50.02 10.02 10.02 7 65.06 60.04 10.04 10.03 8 75.13 70.09 10.09 10.06 9 85.25 80.18 10.17 10.12 10 95.49 90.35 10.33 10.24 11 105.97 100.69 10.65 10.47 12 116.90 111.36 11.27 10.91 13 128.70 122.65 12.47 11.77 14 142.20 135.16 14.78 13.43 15 158.96 150.03 19.23 16.65 16 181.96 169.42 27.66 22.78 17 216.65 197.37 43.26 34.30 18 272.48 241.13 70.88 55.21 19 364.30 312.74 116.11 90.99 20 511.53 429.72 181.55 146.43 21 732.20 611.89 261.03 220.35 22 1033.22 872.87 339.39 301.42 23 1405.70 1211.59 402.26 373.31 24 1830.89 1612.98 444.87 426.00 25 2289.77 2057.13 470.55 459.47 26 2768.24 2527.22 484.95 478.83 27 3257.48 3011.90 492.70 489.44 28 3752.44 3504.46 496.78 495.07 29 4250.40 4001.16 498.90 498.02 30 4749.91 4500.02 500.00 499.54 31 5250.23 5000.00 500.56 500.33 zgr_bat : defines level and meter bathymetry ~~~~~~~ iom_nf90_open ~~~ open existing file: bathy_meter.nc in READ mode ---> bathy_meter.nc OK read Bathymetry (rec: 1) in bathy_meter.nc ok iom_close ~~~ close file: bathy_meter.nc ok orca_r2: Gibraltar strait open at i= 140 j= 102 orca_r2: Bab el Mandeb strait open at i= 160 j= 88 Minimum ocean depth: 30.003215807563720 minimum number of ocean levels : 3 zgr_zps : z-coordinate with partial steps ~~~~~~~ mbathy is recomputed : bathy_level file is NOT used zgr_bat_ctl : check the bathymetry ~~~~~~~~~~~ suppress isolated ocean grid points ----------------------------------- the number of ocean level at grid-point (i,j) = 61 2 is changed from 21 to 20 the number of ocean level at grid-point (i,j) = 66 4 is changed from 21 to 20 the number of ocean level at grid-point (i,j) = 44 6 is changed from 21 to 20 the number of ocean level at grid-point (i,j) = 129 7 is changed from 20 to 19 the number of ocean level at grid-point (i,j) = 45 8 is changed from 21 to 20 the number of ocean level at grid-point (i,j) = 46 9 is changed from 21 to 20 the number of ocean level at grid-point (i,j) = 49 9 is changed from 20 to 19 the number of ocean level at grid-point (i,j) = 55 19 is changed from 29 to 28 the number of ocean level at grid-point (i,j) = 51 20 is changed from 28 to 27 the number of ocean level at grid-point (i,j) = 53 20 is changed from 28 to 27 the number of ocean level at grid-point (i,j) = 110 26 is changed from 28 to 27 the number of ocean level at grid-point (i,j) = 12 28 is changed from 30 to 29 the number of ocean level at grid-point (i,j) = 120 28 is changed from 28 to 27 the number of ocean level at grid-point (i,j) = 14 29 is changed from 30 to 29 the number of ocean level at grid-point (i,j) = 121 29 is changed from 28 to 27 the number of ocean level at grid-point (i,j) = 114 32 is changed from 29 to 28 the number of ocean level at grid-point (i,j) = 42 33 is changed from 30 to 29 the number of ocean level at grid-point (i,j) = 131 36 is changed from 30 to 29 the number of ocean level at grid-point (i,j) = 130 37 is changed from 30 to 29 the number of ocean level at grid-point (i,j) = 129 38 is changed from 30 to 29 the number of ocean level at grid-point (i,j) = 109 40 is changed from 11 to 10 the number of ocean level at grid-point (i,j) = 144 40 is changed from 30 to 29 the number of ocean level at grid-point (i,j) = 179 40 is changed from 28 to 27 the number of ocean level at grid-point (i,j) = 180 41 is changed from 28 to 27 the number of ocean level at grid-point (i,j) = 110 43 is changed from 8 to 6 the number of ocean level at grid-point (i,j) = 58 44 is changed from 30 to 29 the number of ocean level at grid-point (i,j) = 165 46 is changed from 28 to 27 the number of ocean level at grid-point (i,j) = 11 47 is changed from 30 to 29 the number of ocean level at grid-point (i,j) = 46 47 is changed from 25 to 24 the number of ocean level at grid-point (i,j) = 99 47 is changed from 29 to 28 the number of ocean level at grid-point (i,j) = 45 48 is changed from 26 to 25 the number of ocean level at grid-point (i,j) = 47 49 is changed from 28 to 27 the number of ocean level at grid-point (i,j) = 40 52 is changed from 30 to 29 the number of ocean level at grid-point (i,j) = 17 55 is changed from 30 to 29 the number of ocean level at grid-point (i,j) = 68 55 is changed from 30 to 29 the number of ocean level at grid-point (i,j) = 40 56 is changed from 26 to 25 the number of ocean level at grid-point (i,j) = 46 56 is changed from 28 to 26 the number of ocean level at grid-point (i,j) = 48 56 is changed from 27 to 26 the number of ocean level at grid-point (i,j) = 45 57 is changed from 29 to 28 the number of ocean level at grid-point (i,j) = 43 58 is changed from 30 to 29 the number of ocean level at grid-point (i,j) = 67 58 is changed from 29 to 28 the number of ocean level at grid-point (i,j) = 38 59 is changed from 30 to 29 the number of ocean level at grid-point (i,j) = 55 59 is changed from 29 to 28 the number of ocean level at grid-point (i,j) = 169 60 is changed from 30 to 29 the number of ocean level at grid-point (i,j) = 40 62 is changed from 27 to 26 the number of ocean level at grid-point (i,j) = 25 64 is changed from 29 to 28 the number of ocean level at grid-point (i,j) = 42 64 is changed from 26 to 25 the number of ocean level at grid-point (i,j) = 80 64 is changed from 30 to 29 the number of ocean level at grid-point (i,j) = 172 65 is changed from 29 to 28 the number of ocean level at grid-point (i,j) = 36 69 is changed from 24 to 23 the number of ocean level at grid-point (i,j) = 39 69 is changed from 27 to 26 the number of ocean level at grid-point (i,j) = 16 72 is changed from 5 to 4 the number of ocean level at grid-point (i,j) = 25 72 is changed from 26 to 25 the number of ocean level at grid-point (i,j) = 34 73 is changed from 28 to 27 the number of ocean level at grid-point (i,j) = 24 74 is changed from 26 to 25 the number of ocean level at grid-point (i,j) = 124 74 is changed from 30 to 29 the number of ocean level at grid-point (i,j) = 17 78 is changed from 6 to 5 the number of ocean level at grid-point (i,j) = 18 80 is changed from 10 to 7 the number of ocean level at grid-point (i,j) = 12 82 is changed from 7 to 6 the number of ocean level at grid-point (i,j) = 38 82 is changed from 30 to 29 the number of ocean level at grid-point (i,j) = 102 83 is changed from 27 to 26 the number of ocean level at grid-point (i,j) = 19 84 is changed from 25 to 24 the number of ocean level at grid-point (i,j) = 12 87 is changed from 7 to 6 the number of ocean level at grid-point (i,j) = 19 88 is changed from 29 to 28 the number of ocean level at grid-point (i,j) = 166 88 is changed from 25 to 24 the number of ocean level at grid-point (i,j) = 20 89 is changed from 29 to 28 the number of ocean level at grid-point (i,j) = 104 89 is changed from 29 to 28 the number of ocean level at grid-point (i,j) = 91 90 is changed from 28 to 27 the number of ocean level at grid-point (i,j) = 99 91 is changed from 29 to 28 the number of ocean level at grid-point (i,j) = 173 94 is changed from 27 to 26 the number of ocean level at grid-point (i,j) = 99 95 is changed from 27 to 26 the number of ocean level at grid-point (i,j) = 173 98 is changed from 27 to 26 the number of ocean level at grid-point (i,j) = 24 101 is changed from 10 to 9 the number of ocean level at grid-point (i,j) = 174 101 is changed from 27 to 26 the number of ocean level at grid-point (i,j) = 145 102 is changed from 25 to 24 the number of ocean level at grid-point (i,j) = 30 103 is changed from 25 to 24 the number of ocean level at grid-point (i,j) = 170 105 is changed from 20 to 19 the number of ocean level at grid-point (i,j) = 31 106 is changed from 28 to 27 the number of ocean level at grid-point (i,j) = 108 106 is changed from 16 to 15 the number of ocean level at grid-point (i,j) = 132 108 is changed from 30 to 29 the number of ocean level at grid-point (i,j) = 154 109 is changed from 25 to 24 the number of ocean level at grid-point (i,j) = 33 110 is changed from 21 to 20 the number of ocean level at grid-point (i,j) = 38 110 is changed from 26 to 25 the number of ocean level at grid-point (i,j) = 112 110 is changed from 18 to 16 the number of ocean level at grid-point (i,j) = 159 110 is changed from 5 to 4 the number of ocean level at grid-point (i,j) = 39 111 is changed from 24 to 23 the number of ocean level at grid-point (i,j) = 110 111 is changed from 17 to 15 the number of ocean level at grid-point (i,j) = 139 112 is changed from 13 to 9 the number of ocean level at grid-point (i,j) = 46 113 is changed from 30 to 29 the number of ocean level at grid-point (i,j) = 45 114 is changed from 30 to 29 the number of ocean level at grid-point (i,j) = 64 114 is changed from 30 to 29 the number of ocean level at grid-point (i,j) = 135 114 is changed from 26 to 25 the number of ocean level at grid-point (i,j) = 47 115 is changed from 28 to 27 the number of ocean level at grid-point (i,j) = 141 115 is changed from 12 to 10 the number of ocean level at grid-point (i,j) = 61 116 is changed from 8 to 7 the number of ocean level at grid-point (i,j) = 66 116 is changed from 30 to 29 the number of ocean level at grid-point (i,j) = 64 117 is changed from 14 to 10 the number of ocean level at grid-point (i,j) = 140 117 is changed from 14 to 13 the number of ocean level at grid-point (i,j) = 60 118 is changed from 6 to 5 the number of ocean level at grid-point (i,j) = 68 118 is changed from 29 to 28 the number of ocean level at grid-point (i,j) = 104 118 is changed from 6 to 5 the number of ocean level at grid-point (i,j) = 124 118 is changed from 27 to 26 the number of ocean level at grid-point (i,j) = 141 118 is changed from 14 to 13 the number of ocean level at grid-point (i,j) = 143 118 is changed from 16 to 12 the number of ocean level at grid-point (i,j) = 37 119 is changed from 17 to 16 the number of ocean level at grid-point (i,j) = 110 120 is changed from 15 to 14 the number of ocean level at grid-point (i,j) = 142 120 is changed from 18 to 17 the number of ocean level at grid-point (i,j) = 52 121 is changed from 11 to 10 the number of ocean level at grid-point (i,j) = 111 122 is changed from 20 to 19 the number of ocean level at grid-point (i,j) = 108 125 is changed from 18 to 17 the number of ocean level at grid-point (i,j) = 130 125 is changed from 22 to 21 the number of ocean level at grid-point (i,j) = 53 126 is changed from 8 to 7 the number of ocean level at grid-point (i,j) = 107 126 is changed from 18 to 17 the number of ocean level at grid-point (i,j) = 132 126 is changed from 23 to 22 the number of ocean level at grid-point (i,j) = 149 126 is changed from 14 to 10 the number of ocean level at grid-point (i,j) = 54 127 is changed from 10 to 7 the number of ocean level at grid-point (i,j) = 58 127 is changed from 6 to 5 the number of ocean level at grid-point (i,j) = 47 128 is changed from 7 to 6 the number of ocean level at grid-point (i,j) = 99 128 is changed from 15 to 14 the number of ocean level at grid-point (i,j) = 100 129 is changed from 14 to 12 the number of ocean level at grid-point (i,j) = 109 129 is changed from 8 to 7 the number of ocean level at grid-point (i,j) = 101 130 is changed from 8 to 7 the number of ocean level at grid-point (i,j) = 110 130 is changed from 9 to 7 the number of ocean level at grid-point (i,j) = 106 131 is changed from 13 to 12 the number of ocean level at grid-point (i,j) = 108 133 is changed from 11 to 8 the number of ocean level at grid-point (i,j) = 75 135 is changed from 19 to 18 the number of ocean level at grid-point (i,j) = 155 135 is changed from 9 to 6 the number of ocean level at grid-point (i,j) = 149 137 is changed from 19 to 18 the number of ocean level at grid-point (i,j) = 154 138 is changed from 17 to 16 the number of ocean level at grid-point (i,j) = 68 139 is changed from 20 to 19 the number of ocean level at grid-point (i,j) = 128 141 is changed from 19 to 18 the number of ocean level at grid-point (i,j) = 143 141 is changed from 19 to 18 the number of ocean level at grid-point (i,j) = 152 141 is changed from 18 to 16 the number of ocean level at grid-point (i,j) = 155 141 is changed from 18 to 17 the number of ocean level at grid-point (i,j) = 151 142 is changed from 17 to 16 the number of ocean level at grid-point (i,j) = 150 143 is changed from 18 to 17 the number of ocean level at grid-point (i,j) = 118 145 is changed from 19 to 17 137 ocean grid points suppressed east-west cyclic boundary conditions on mbathy: nperio = 4 zgr_bot_level : ocean bottom k-index of T-, U-, V- and W-levels ~~~~~~~~~~~~~ dommsk : ocean mask ~~~~~~ Namelist namlbc lateral momentum boundary cond. rn_shlat = 2.0000000000000000 consistency with analytical form ln_vorlat = F ocean lateral no-slip dom_stp : time stepping setting ~~~~~~~ synchronous time stepping dynamics and tracer time step = 1.6000000000000001 hours dyn_nept_init : Simplified Neptune module enabled ~~~~~~~~~~~~~ --> Reading namelist namdyn_nept parameters: ln_neptsimp = F ln_smooth_neptvel = F rn_tslse = 12000.000000000000 rn_tslsp = 3000.0000000000000 ln_neptramp = T rn_htrmin = 100.00000000000000 rn_htrmax = 200.00000000000000 --> neptune velocities will not be smoothed --> ln_neptsimp disabled istate_ini : Initialization of the dynamics and tracers ~~~~~~~~~~ dta_tsd_init : Temperature & Salinity data ~~~~~~~~~~~~ Namelist namtsd Initialisation of ocean T & S with T &S input data ln_tsd_init = T damping of ocean T & S toward T &S input data ln_tsd_tradmp = T ===>>> : W A R N I N G =============== dta_tsd_init: ocean restart and T & S data intialisation, we keep the restart T & S values and set ln_tsd_init to FALSE dta_tsd : Temperature & Salinity data ~~~~~~~ namtsd Namelist list of files and frequency (>0: in hours ; <0 in months) root filename: ./data_1m_potential_temperature_nomask variable name: votemper frequency: -1 time interp: T climatology: T weights : pairing : data type: yearly root filename: ./data_1m_salinity_nomask variable name: vosaline frequency: -1 time interp: T climatology: T weights : pairing : data type: yearly rst_read : read oce NetCDF restart file ~~~~~~~~ ===>>> : E R R O R =========== iom_open ~~~ File O2LP_LONG_00000075_restart.nc* not found *** Info used values : date ndastp : 10100 number of elapsed days since the begining of run : 5.0000000000000000 ==============>> 1/2 time step before the start of the run DATE Y/M/D = 0/12/31 nsec_day: 83520 nsec_week:****** ======>> time-step = 76 New day, DATE Y/M/D = 0001/01/01 nday_year = 001 nsec_year = 2880 nsec_month = 2880 nsec_day = 2880 nsec_week = -169920 sbc_init : surface boundary condition setting ~~~~~~~~ Namelist namsbc (partly overwritten with CPP key setting) frequency update of sbc (and ice) nn_fsbc = 5 Type of sbc : analytical formulation ln_ana = F flux formulation ln_flx = F CLIO bulk formulation ln_blk_clio = F CORE bulk formulation ln_blk_core = T MFS bulk formulation ln_blk_mfs = F coupled formulation (T if key_sbc_cpl) ln_cpl = F Misc. options of sbc : Patm gradient added in ocean & ice Eqs. ln_apr_dyn = F ice management in the sbc (=0/1/2/3) nn_ice = 2 daily mean to diurnal cycle qsr ln_dm2dc = F runoff / runoff mouths ln_rnf = T Sea Surface Restoring on SST and/or SSS ln_ssr = T FreshWater Budget control (=0/1/2) nn_fwb = 3 closed sea (=0/1) (set in namdom) nn_closea = 0 CORE bulk formulation zdf_init: vertical physics ~~~~~~~~ Namelist namzdf : set vertical mixing mixing parameters vertical eddy viscosity rn_avm0 = 1.20000000000000003E-004 vertical eddy diffusivity rn_avt0 = 1.20000000000000003E-005 constant background or profile nn_avb = 0 horizontal variation for avtb nn_havtb = 0 time splitting / backward scheme ln_zdfexp = F number of time step nn_zdfexp = 3 enhanced vertical diffusion ln_zdfevd = T applied on momentum (=1/0) nn_evdm = 0 vertical coefficient for evd rn_avevd = 100.00000000000000 non-penetrative convection (npc) ln_zdfnpc = F npc call frequency nn_npc = 1 npc print frequency nn_npcp = 365 vertical mixing option : TKE dependent eddy coefficients convection : use enhanced vertical dif. scheme use the 1.5 turbulent closure zdf_bfr : momentum bottom friction ~~~~~~~ Namelist nam_bfr : set bottom friction parameters linear botton friction friction coef. rn_bfri1 = 4.00000000000000019E-004 implicit bottom friction switch ln_bfrimp = T zdf_tke_init : tke turbulent closure scheme - initialisation ~~~~~~~~~~~~ Namelist namzdf_tke : set tke mixing parameters coef. to compute avt rn_ediff = 0.10000000000000001 Kolmogoroff dissipation coef. rn_ediss = 0.69999999999999996 tke surface input coef. rn_ebb = 67.829999999999998 minimum value of tke rn_emin = 9.99999999999999955E-007 surface minimum value of tke rn_emin0 = 1.00000000000000005E-004 background shear (>0) rn_bshear = 9.99999999999999945E-021 mixing length type nn_mxl = 2 prandl number flag nn_pdl = 1 surface mixing length = F(stress) or not ln_mxl0 = T surface mixing length minimum value rn_mxl0 = 4.00000000000000008E-002 flag to take into acc. Langmuir circ. ln_lc = T coef to compute verticla velocity of LC rn_lc = 0.14999999999999999 test param. to add tke induced by wind nn_etau = 1 flag for computation of exp. tke profile nn_htau = 1 fraction of en which pene. the thermocline rn_efr = 5.00000000000000028E-002 critical Richardson nb with your parameters ri_cri = 0.22222222222222221 use a surface mixing length = F(stress) : set rn_mxl0 = rmxl_min ===>>>> : previous run without tke scheme, en computed by iterative loop zdf_ddm : double diffusive mixing ~~~~~~~ Namelist namzdf_ddm : set dd mixing parameter maximum avs for dd mixing rn_avts = 1.00000000000000005E-004 heat/salt buoyancy flux ratio rn_hsbfr = 1.6000000000000001 ldf_tra_init : lateral tracer physics ~~~~~~~~~~~~ Namelist namtra_ldf : lateral mixing parameters (type, direction, coefficients) laplacian operator ln_traldf_lap = T bilaplacian operator ln_traldf_bilap = F iso-level ln_traldf_level = F horizontal (geopotential) ln_traldf_hor = F iso-neutral ln_traldf_iso = T iso-neutral (Griffies) ln_traldf_grif = F Griffies strmfn diagnostics ln_traldf_gdia = F lateral eddy diffusivity rn_aht_0 = 2000.0000000000000 background hor. diffusivity rn_ahtb_0 = 0.0000000000000000 eddy induced velocity coef. rn_aeiv_0 = 2000.0000000000000 maximum isoppycnal slope rn_slpmax = 1.00000000000000002E-002 pure lateral mixing in ML ln_triad_iso = F lateral mixing on bottom ln_botmix_grif = F tracer mixing coef. = F( latitude, longitude) harmonic tracer diffusion (default) ldf_tra_c2d : 2D eddy diffusivity and eddy ~~~~~~~~~~~ -- induced velocity coefficients ORCA R2 or R4 case Constant values used for eddy diffusivity coefficients Variation lat/lon only for eddy induced velocity coefficients ldf_dyn : lateral momentum physics ~~~~~~~ Namelist nam_dynldf : set lateral mixing parameters laplacian operator ln_dynldf_lap = T bilaplacian operator ln_dynldf_bilap = F iso-level ln_dynldf_level = F horizontal (geopotential) ln_dynldf_hor = T iso-neutral ln_dynldf_iso = F horizontal laplacian eddy viscosity rn_ahm_0_lap = 40000.000000000000 background viscosity rn_ahmb_0 = 0.0000000000000000 horizontal bilaplacian eddy viscosity rn_ahm_0_blp = 0.0000000000000000 momentum mixing coef. = F( latitude, longitude, depth) harmonic momentum diff. (default) ldf_dyn_c3d : 3D lateral eddy viscosity coefficient ~~~~~~~~~~~ laplacian operator: ahm proportional to e1 maximum grid-spacing = 222397.84689709079 maximum value for ahm = 40000.000000000000 ahm profile : jk ahm depth t-level 1 1.0000 4.9999 2 0.9994 15.0003 3 0.9989 25.0018 4 0.9982 35.0054 5 0.9976 45.0133 6 0.9970 55.0295 7 0.9963 65.0618 8 0.9956 75.1255 9 0.9948 85.2504 10 0.9940 95.4943 11 0.9932 105.9699 12 0.9923 116.8962 13 0.9913 128.6979 14 0.9901 142.1952 15 0.9885 158.9606 16 0.9861 181.9628 17 0.9821 216.6479 18 0.9745 272.4767 19 0.9583 364.3030 20 0.9192 511.5348 21 0.8191 732.2009 22 0.6067 1033.2173 23 0.3759 1405.6975 24 0.2766 1830.8850 25 0.2544 2289.7679 26 0.2506 2768.2423 27 0.2501 3257.4789 28 0.2500 3752.4422 29 0.2500 4250.4012 30 0.2500 4749.9133 31 0.2500 5250.2266 ahm profile : jk ahm depth t-level 1 1.0000 4.9999 2 0.9994 15.0003 3 0.9989 25.0018 4 0.9982 35.0054 5 0.9976 45.0133 6 0.9970 55.0295 7 0.9963 65.0618 8 0.9956 75.1255 9 0.9948 85.2504 10 0.9940 95.4943 11 0.9932 105.9699 12 0.9923 116.8962 13 0.9913 128.6979 14 0.9901 142.1952 15 0.9885 158.9606 16 0.9861 181.9628 17 0.9821 216.6479 18 0.9745 272.4767 19 0.9583 364.3030 20 0.9192 511.5348 21 0.8191 732.2009 22 0.6067 1033.2173 23 0.3759 1405.6975 24 0.2766 1830.8850 25 0.2544 2289.7679 26 0.2506 2768.2423 27 0.2501 3257.4789 28 0.2500 3752.4422 29 0.2500 4250.4012 30 0.2500 4749.9133 31 0.2500 5250.2266 ORCA R1, R2 or R4: overwrite the previous definition of ahm ================= ldfdyn_c3d_orca : 3D eddy viscosity coefficient ~~~~~~~~~~~~~~~ orca R1, R2 or R4 configuration: reduced in the surface Eq. strip file : ahmcoef open ok unit = 19 status = OLD form = FORMATTED access = SEQUENTIAL 1D zcoef array ~~~~~~~~~~~~~~ jk zcoef 1 1.000 2 1.000 3 1.000 4 1.000 5 1.000 6 1.000 7 1.000 8 1.000 9 1.000 10 1.000 11 1.000 12 1.000 13 1.000 14 1.000 15 1.000 16 1.000 17 1.000 18 1.000 19 1.000 20 1.000 21 1.000 22 2.000 23 3.000 24 5.000 25 7.000 26 9.000 27 10.000 28 10.000 29 10.000 30 10.000 31 10.000 3D ahm1 array (k=1) 1 21 41 61 81 101 121 141 161 181 149 0.00 40.00 40.00 40.00 40.00 40.00 40.00 40.00 40.00 40.00 129 40.00 40.00 40.00 40.00 40.00 40.00 40.00 40.00 40.00 40.00 109 40.00 40.00 40.00 40.00 40.00 40.00 40.00 40.00 40.00 40.00 89 40.00 40.00 30.21 30.21 30.21 40.00 30.21 40.00 40.00 40.00 69 2.01 40.00 28.22 2.01 2.01 2.01 40.00 2.01 40.00 2.01 49 40.00 40.00 40.00 40.00 40.00 40.00 40.00 40.00 40.00 40.00 29 40.00 40.00 40.00 40.00 40.00 40.00 40.00 40.00 40.00 40.00 9 40.00 40.00 40.00 40.00 40.00 40.00 40.00 40.00 40.00 40.00 3D ahm2 array (k=1) 1 21 41 61 81 101 121 141 161 181 149 40.00 40.00 40.00 40.00 40.00 40.00 40.00 40.00 40.00 40.00 129 40.00 40.00 40.00 40.00 40.00 40.00 40.00 40.00 40.00 40.00 109 40.00 40.00 40.00 40.00 40.00 40.00 40.00 40.00 40.00 40.00 89 40.00 40.00 32.72 32.72 32.72 40.00 32.72 40.00 40.00 40.00 69 2.00 40.00 17.39 2.00 2.00 2.00 40.00 2.00 40.00 2.00 49 40.00 40.00 40.00 40.00 40.00 40.00 40.00 40.00 40.00 40.00 29 40.00 40.00 40.00 40.00 40.00 40.00 40.00 40.00 40.00 40.00 9 40.00 40.00 40.00 40.00 40.00 40.00 40.00 40.00 40.00 40.00 3D ahm2 array (k=jpk) 1 21 41 61 81 101 121 141 161 181 149 40.00 40.00 40.00 40.00 40.00 40.00 40.00 40.00 40.00 40.00 129 40.00 40.00 40.00 40.00 40.00 40.00 40.00 40.00 40.00 40.00 109 40.00 40.00 40.00 40.00 40.00 40.00 40.00 40.00 40.00 40.00 89 40.00 40.00 40.00 40.00 40.00 40.00 40.00 40.00 40.00 40.00 69 20.00 40.00 40.00 20.00 20.00 20.00 40.00 20.00 40.00 20.00 49 40.00 40.00 40.00 40.00 40.00 40.00 40.00 40.00 40.00 40.00 29 40.00 40.00 40.00 40.00 40.00 40.00 40.00 40.00 40.00 40.00 9 40.00 40.00 40.00 40.00 40.00 40.00 40.00 40.00 40.00 40.00 ldf_slp_init : direction of lateral mixing ~~~~~~~~~~~~ Horizontal mixing in s-coordinate: slope = slope of s-surfaces tra_qsr_init : penetration of the surface solar radiation ~~~~~~~~~~~~ Namelist namtra_qsr : set the parameter of penetration Light penetration (T) or not (F) ln_traqsr = T RGB (Red-Green-Blue) light penetration ln_qsr_rgb = T 2 band light penetration ln_qsr_2bd = F bio-model light penetration ln_qsr_bio = F RGB : Chl data (=1) or cst value (=0) nn_chldta = 1 RGB & 2 bands: fraction of light (rn_si1) rn_abs = 0.57999999999999996 RGB & 2 bands: shortess depth of extinction rn_si0 = 0.34999999999999998 2 bands: longest depth of extinction rn_si1 = 23.000000000000000 R-G-B light penetration - Chl data trc_oce_rgb : Initialisation of the optical look-up table ~~~~~~~~~~~ RGB longest depth of extinction r_si2 = 61.804697156983934 jc = 1 Chl = 1.00000000000000002E-002 irgb = 1 jc = 2 Chl = 1.09999999999999994E-002 irgb = 2 jc = 3 Chl = 1.29999999999999994E-002 irgb = 3 jc = 4 Chl = 1.40000000000000003E-002 irgb = 4 jc = 5 Chl = 1.60000000000000003E-002 irgb = 5 jc = 6 Chl = 1.79999999999999986E-002 irgb = 6 jc = 7 Chl = 2.00000000000000004E-002 irgb = 7 jc = 8 Chl = 2.19999999999999987E-002 irgb = 8 jc = 9 Chl = 2.50000000000000014E-002 irgb = 9 jc = 10 Chl = 2.80000000000000006E-002 irgb = 10 jc = 11 Chl = 3.20000000000000007E-002 irgb = 11 jc = 12 Chl = 3.50000000000000033E-002 irgb = 12 jc = 13 Chl = 4.00000000000000008E-002 irgb = 13 jc = 14 Chl = 4.49999999999999983E-002 irgb = 14 jc = 15 Chl = 5.00000000000000028E-002 irgb = 15 jc = 16 Chl = 5.60000000000000012E-002 irgb = 16 jc = 17 Chl = 6.30000000000000004E-002 irgb = 17 jc = 18 Chl = 7.09999999999999937E-002 irgb = 18 jc = 19 Chl = 7.90000000000000008E-002 irgb = 19 jc = 20 Chl = 8.89999999999999958E-002 irgb = 20 jc = 21 Chl = 0.10000000000000001 irgb = 21 jc = 22 Chl = 0.11200000000000000 irgb = 22 jc = 23 Chl = 0.12600000000000000 irgb = 23 jc = 24 Chl = 0.14099999999999999 irgb = 24 jc = 25 Chl = 0.15800000000000000 irgb = 25 jc = 26 Chl = 0.17799999999999999 irgb = 26 jc = 27 Chl = 0.20000000000000001 irgb = 27 jc = 28 Chl = 0.22400000000000000 irgb = 28 jc = 29 Chl = 0.25100000000000000 irgb = 29 jc = 30 Chl = 0.28199999999999997 irgb = 30 jc = 31 Chl = 0.31600000000000000 irgb = 31 jc = 32 Chl = 0.35499999999999998 irgb = 32 jc = 33 Chl = 0.39800000000000002 irgb = 33 jc = 34 Chl = 0.44700000000000001 irgb = 34 jc = 35 Chl = 0.50100000000000000 irgb = 35 jc = 36 Chl = 0.56200000000000006 irgb = 36 jc = 37 Chl = 0.63100000000000001 irgb = 37 jc = 38 Chl = 0.70799999999999996 irgb = 38 jc = 39 Chl = 0.79400000000000004 irgb = 39 jc = 40 Chl = 0.89100000000000001 irgb = 40 jc = 41 Chl = 1.0000000000000000 irgb = 41 jc = 42 Chl = 1.1220000000000001 irgb = 42 jc = 43 Chl = 1.2589999999999999 irgb = 43 jc = 44 Chl = 1.4130000000000000 irgb = 44 jc = 45 Chl = 1.5850000000000000 irgb = 45 jc = 46 Chl = 1.7780000000000000 irgb = 46 jc = 47 Chl = 1.9950000000000001 irgb = 47 jc = 48 Chl = 2.2389999999999999 irgb = 48 jc = 49 Chl = 2.5120000000000000 irgb = 49 jc = 50 Chl = 2.8180000000000001 irgb = 50 jc = 51 Chl = 3.1619999999999999 irgb = 51 jc = 52 Chl = 3.5480000000000000 irgb = 52 jc = 53 Chl = 3.9809999999999999 irgb = 53 jc = 54 Chl = 4.4669999999999996 irgb = 54 jc = 55 Chl = 5.0119999999999996 irgb = 55 jc = 56 Chl = 5.6230000000000002 irgb = 56 jc = 57 Chl = 6.3099999999999996 irgb = 57 jc = 58 Chl = 7.0789999999999997 irgb = 58 jc = 59 Chl = 7.9429999999999996 irgb = 59 jc = 60 Chl = 8.9120000000000008 irgb = 60 jc = 61 Chl = 10.000000000000000 irgb = 61 level of light extinction = 25 ref depth = 2527.2169154512344 m Chlorophyll read in a file tra_qsr_init : Solar penetration function of read chlorophyll ~~~~~~~~~~~~ namtra_qsr Namelist list of files and frequency (>0: in hours ; <0 in months) root filename: ./chlorophyll variable name: CHLA frequency: -1 time interp: T climatology: T weights : pairing : data type: yearly tra_bbc : Bottom Boundary Condition (bbc), apply a Geothermal heating ~~~~~~~ Namelist nambbc : set bbc parameters Apply a geothermal heating at ocean bottom ln_trabbc = T type of geothermal flux nn_geoflx = 2 Constant geothermal flux value rn_geoflx_cst = 8.64000000000000046E-002 *** variable geothermal heat flux iom_nf90_open ~~~ open existing file: geothermal_heating.nc in READ mode ---> geothermal_heating.nc OK read heatflow (rec: 1) in geothermal_heating.nc ok iom_close ~~~ close file: geothermal_heating.nc ok tra_bbl_init : bottom boundary layer initialisation ~~~~~~~~~~~~ Namelist nambbl : set bbl parameters diffusive bbl (=1) or not (=0) nn_bbl_ldf = 1 advective bbl (=1/2) or not (=0) nn_bbl_adv = 0 diffusive bbl coefficient rn_ahtbbl = 1000.0000000000000 m2/s advective bbl coefficient rn_gambbl = 10.000000000000000 s tra_dmp_init : T and S newtonian damping ~~~~~~~ Namelist namtra_dmp : set damping parameter add a damping termn or not ln_tradmp = T T and S damping option nn_hdmp = -1 mixed layer damping option nn_zdmp = 0 (zoom: forced to 0) surface time scale (days) rn_surf = 50.000000000000000 bottom time scale (days) rn_bot = 360.00000000000000 depth of transition (meters) rn_dep = 800.00000000000000 create a damping.coeff file nn_file = 0 tracer damping in the Med & Red seas only tracer damping throughout the water column dtacof : Global domain of ORCA ------------------------------ ORCA configuration: Damping in Med and Red Seas tra_adv_init : choice/control of the tracer advection scheme ~~~~~~~~~~~ Namelist namtra_adv : chose a advection scheme for tracers 2nd order advection scheme ln_traadv_cen2 = F TVD advection scheme ln_traadv_tvd = T MUSCL advection scheme ln_traadv_muscl = F MUSCL2 advection scheme ln_traadv_muscl2 = F UBS advection scheme ln_traadv_ubs = F QUICKEST advection scheme ln_traadv_qck = F TVD scheme is used tra_ldf_init : lateral tracer diffusive operator ~~~~~~~~~~~ Namelist namtra_ldf already read in ldftra module see ldf_tra_init report for lateral mixing parameters Rotated laplacian operator tra:ldf_ano : lateral diffusion acting on the full fields ~~~~~~~~~~~ tra_zdf_init : vertical tracer physics scheme ~~~~~~~~~~~ Implicit (euler backward) scheme dyn_adv_init : choice/control of the momentum advection scheme ~~~~~~~~~~~ Namelist namdyn_adv : chose a advection formulation & scheme for momentum Vector/flux form (T/F) ln_dynadv_vec = T 2nd order centred advection scheme ln_dynadv_cen2 = F 3rd order UBS advection scheme ln_dynadv_ubs = F vector form : keg + zad + vor is used dyn_vor_init : vorticity term : read namelist and control the consistency ~~~~~~~~~~~~ Namelist namdyn_vor : oice of the vorticity term scheme energy conserving scheme ln_dynvor_ene = F enstrophy conserving scheme ln_dynvor_ens = F mixed enstrophy/energy conserving scheme ln_dynvor_mix = F enstrophy and energy conserving scheme ln_dynvor_een = T Vector form advection : vorticity = Coriolis + relative vorticity vorticity scheme : energy and enstrophy conserving scheme dyn_ldf_init : Choice of the lateral diffusive operator on dynamics ~~~~~~~~~~~ Namelist nam_dynldf : set lateral mixing parameters (type, direction, coefficients) laplacian operator ln_dynldf_lap = T bilaplacian operator ln_dynldf_bilap = F iso-level ln_dynldf_level = F horizontal (geopotential) ln_dynldf_hor = T iso-neutral ln_dynldf_iso = F laplacian operator dyn_hpg_init : hydrostatic pressure gradient initialisation ~~~~~~~~~~~~ Namelist namdyn_hpg : choice of hpg scheme z-coord. - full steps ln_hpg_zco = F z-coord. - partial steps (interpolation) ln_hpg_zps = T s-coord. (standard jacobian formulation) ln_hpg_sco = F s-coord. (Density Jacobian: Cubic polynomial) ln_hpg_djc = F s-coord. (Pressure Jacobian: Cubic polynomial) ln_hpg_prj = F time stepping: centered (F) or semi-implicit (T) ln_dynhpg_imp = F dyn_zdf_init : vertical dynamics physics scheme ~~~~~~~~~~~ Implicit (euler backward) scheme dyn_spg_init : choice of the surface pressure gradient scheme ~~~~~~~~~~~ Explicit free surface lk_dynspg_exp = F Free surface with time splitting lk_dynspg_ts = F Filtered free surface cst volume lk_dynspg_flt = T filtered free surface file : solver.stat open ok unit = 20 status = REPLACE form = FORMATTED access = SEQUENTIAL solver_init : solver to compute the surface pressure gradient ~~~~~~~~~~~ Namelist namsol : set solver parameters type of elliptic solver nn_solv = 2 absolute/relative (0/1) precision nn_sol_arp = 0 minimum iterations for solver nn_nmin = 300 maximum iterations for solver nn_nmax = 800 frequency for test nn_nmod = 10 absolute precision of solver rn_eps = 9.99999999999999955E-007 absolute precision for SOR solver rn_resmax = 1.00000000000000004E-010 optimal coefficient of sor rn_sor = 1.9199999999999999 a successive-over-relaxation solver with extra outer halo is used with jpr2di = 0 and jpr2dj = 0 trc_init : initial set up of the passive tracers ~~~~~~~ trc_nam : read the passive tracer namelists ~~~~~~~ Namelist : namtrc time step freq. for passive tracer nn_dttrc = 1 restart for passive tracer ln_rsttr = T control of time step for passive tracer nn_rsttr = 2 first time step for pass. trac. nittrc000 = 76 frequency of outputs for passive tracers nn_writetrc = 5475 Read inputs data from file ln_trcdta = T tracer nb : 1 short name : DIC tracer nb : 2 short name : Alkalini tracer nb : 3 short name : O2 tracer nb : 4 short name : CaCO3 tracer nb : 5 short name : PO4 tracer nb : 6 short name : POC tracer nb : 7 short name : Si tracer nb : 8 short name : PHY tracer nb : 9 short name : ZOO tracer nb : 10 short name : DOC tracer nb : 11 short name : PHY2 tracer nb : 12 short name : ZOO2 tracer nb : 13 short name : DSi tracer nb : 14 short name : Fer tracer nb : 15 short name : BFe tracer nb : 16 short name : GOC tracer nb : 17 short name : SFe tracer nb : 18 short name : DFe tracer nb : 19 short name : GSi tracer nb : 20 short name : NFe tracer nb : 21 short name : NCHL tracer nb : 22 short name : DCHL tracer nb : 23 short name : NO3 tracer nb : 24 short name : NH4 Passive Tracer time step rdttrc = 5760.0000000000000 Namelist : namtrc_dia save additionnal diagnostics arrays ln_diatrc = F save additionnal biology diagnostics arrays ln_diabio = F frequency of outputs for additional arrays nn_writedia = 5475 frequency of outputs for biological trends nn_writebio = 10 trc_nam_trp: read namelist for tracer transport ~~~~~~~~~~~ trc_adv_ctl : choice/control of the tracer advection scheme ~~~~~~~~~~~ Namelist namtrc_adv : chose a advection scheme for tracers 2nd order advection scheme ln_trcadv_cen2 = F TVD advection scheme ln_trcadv_tvd = F MUSCL advection scheme ln_trcadv_muscl = T MUSCL2 advection scheme ln_trcadv_muscl2 = F UBS advection scheme ln_trcadv_ubs = F QUICKEST advection scheme ln_trcadv_qck = F trc:ldf_ctl : lateral tracer diffusive operator ~~~~~~~~~~~ Namelist namtrc_ldf : set lateral mixing parameters (type, direction, coefficients) perform lateral diffusion or not ln_trcldf_diff = T laplacian operator ln_trcldf_lap = T bilaplacian operator ln_trcldf_bilap = F iso-level ln_trcldf_level = F horizontal (geopotential) ln_trcldf_hor = F iso-neutral ln_trcldf_iso = T diffusivity coefficient rn_ahtrc_0 = 2000.0000000000000 background hor. diffusivity rn_ahtrb_0 = 0.0000000000000000 Namelist namtrc_zdf : set vertical diffusion parameters time splitting / backward scheme ln_trczdf_exp = F number of time step nn_trczdf_exp = 3 Namelist namtrc_rad : treatment of negative concentrations correct artificially negative concen. or not ln_trcrad = T ===>>> : W A R N I N G =============== trc_nam: passive tracer restart and data intialisation, we keep the restart values and set ln_trcdta to FALSE read a restart file for passive tracer : O2LP_LONG_00000075_restart_trc LOBSTER not used trc_nam_pisces : read PISCES namelists ~~~~~~~~~~~~~~ Namelist : nampisbio frequence pour la biologie nrdttrc = 1 POC sinking speed wsbio = 2.0000000000000000 half saturation constant for mortality xkmort = 9.99999999999999955E-008 Fe/C in zooplankton ferat3 = 1.00000000000000008E-005 Big particles sinking speed wsbio2 = 30.000000000000000 Namelist : nampisdmp Relaxation of tracer to glodap mean value ln_pisdmp = T Frequency of Relaxation nn_pisdmp = 5475 Restoring of tracer to initial value on closed seas ln_pisclo = F CFC not used C14 not used MY_TRC not used trc_rst_cal : read the TOP restart file for calendar ~~~~~~~~~~~~ ===>>> : E R R O R =========== iom_open ~~~ File O2LP_LONG_00000075_restart_trc.nc* not found *** Info read in restart : previous time-step : 0 *** restart option nn_rsttr = 2 : calendar parameters read in restart ===>>> : E R R O R =========== ===>>>> : problem with nittrc000 for the restart verify the restart file or rerun with nn_rsttr = 0 (namelist) trc_ini_pisces : PISCES biochemical model initialisation ~~~~~~~~~~~~~~ Passive Tracer time step rfact = 5760.0000000000000 rdt = 5760.0000000000000 PISCES Biology time step rfact2 = 5760.0000000000000 trc_oce_rgb : Initialisation of the optical look-up table ~~~~~~~~~~~ RGB longest depth of extinction r_si2 = 61.804697156983934 jc = 1 Chl = 1.00000000000000002E-002 irgb = 1 jc = 2 Chl = 1.09999999999999994E-002 irgb = 2 jc = 3 Chl = 1.29999999999999994E-002 irgb = 3 jc = 4 Chl = 1.40000000000000003E-002 irgb = 4 jc = 5 Chl = 1.60000000000000003E-002 irgb = 5 jc = 6 Chl = 1.79999999999999986E-002 irgb = 6 jc = 7 Chl = 2.00000000000000004E-002 irgb = 7 jc = 8 Chl = 2.19999999999999987E-002 irgb = 8 jc = 9 Chl = 2.50000000000000014E-002 irgb = 9 jc = 10 Chl = 2.80000000000000006E-002 irgb = 10 jc = 11 Chl = 3.20000000000000007E-002 irgb = 11 jc = 12 Chl = 3.50000000000000033E-002 irgb = 12 jc = 13 Chl = 4.00000000000000008E-002 irgb = 13 jc = 14 Chl = 4.49999999999999983E-002 irgb = 14 jc = 15 Chl = 5.00000000000000028E-002 irgb = 15 jc = 16 Chl = 5.60000000000000012E-002 irgb = 16 jc = 17 Chl = 6.30000000000000004E-002 irgb = 17 jc = 18 Chl = 7.09999999999999937E-002 irgb = 18 jc = 19 Chl = 7.90000000000000008E-002 irgb = 19 jc = 20 Chl = 8.89999999999999958E-002 irgb = 20 jc = 21 Chl = 0.10000000000000001 irgb = 21 jc = 22 Chl = 0.11200000000000000 irgb = 22 jc = 23 Chl = 0.12600000000000000 irgb = 23 jc = 24 Chl = 0.14099999999999999 irgb = 24 jc = 25 Chl = 0.15800000000000000 irgb = 25 jc = 26 Chl = 0.17799999999999999 irgb = 26 jc = 27 Chl = 0.20000000000000001 irgb = 27 jc = 28 Chl = 0.22400000000000000 irgb = 28 jc = 29 Chl = 0.25100000000000000 irgb = 29 jc = 30 Chl = 0.28199999999999997 irgb = 30 jc = 31 Chl = 0.31600000000000000 irgb = 31 jc = 32 Chl = 0.35499999999999998 irgb = 32 jc = 33 Chl = 0.39800000000000002 irgb = 33 jc = 34 Chl = 0.44700000000000001 irgb = 34 jc = 35 Chl = 0.50100000000000000 irgb = 35 jc = 36 Chl = 0.56200000000000006 irgb = 36 jc = 37 Chl = 0.63100000000000001 irgb = 37 jc = 38 Chl = 0.70799999999999996 irgb = 38 jc = 39 Chl = 0.79400000000000004 irgb = 39 jc = 40 Chl = 0.89100000000000001 irgb = 40 jc = 41 Chl = 1.0000000000000000 irgb = 41 jc = 42 Chl = 1.1220000000000001 irgb = 42 jc = 43 Chl = 1.2589999999999999 irgb = 43 jc = 44 Chl = 1.4130000000000000 irgb = 44 jc = 45 Chl = 1.5850000000000000 irgb = 45 jc = 46 Chl = 1.7780000000000000 irgb = 46 jc = 47 Chl = 1.9950000000000001 irgb = 47 jc = 48 Chl = 2.2389999999999999 irgb = 48 jc = 49 Chl = 2.5120000000000000 irgb = 49 jc = 50 Chl = 2.8180000000000001 irgb = 50 jc = 51 Chl = 3.1619999999999999 irgb = 51 jc = 52 Chl = 3.5480000000000000 irgb = 52 jc = 53 Chl = 3.9809999999999999 irgb = 53 jc = 54 Chl = 4.4669999999999996 irgb = 54 jc = 55 Chl = 5.0119999999999996 irgb = 55 jc = 56 Chl = 5.6230000000000002 irgb = 56 jc = 57 Chl = 6.3099999999999996 irgb = 57 jc = 58 Chl = 7.0789999999999997 irgb = 58 jc = 59 Chl = 7.9429999999999996 irgb = 59 jc = 60 Chl = 8.9120000000000008 irgb = 60 jc = 61 Chl = 10.000000000000000 irgb = 61 level of light extinction = 25 ref depth = 2527.2169154512344 m Namelist parameters for nutrient limitations, nampislim ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ mean rainratio caco3r = 0.16000000000000000 NO3, PO4 half saturation conc0 = 1.99999999999999991E-006 half saturation constant for Si uptake xksi1 = 1.99999999999999991E-006 half saturation constant for Si/C xksi2 = 3.32999999999999988E-006 2nd half-sat. of DOC remineralization xkdoc = 4.16999999999999999E-004 Phosphate half saturation for diatoms conc1 = 7.99999999999999964E-006 Iron half saturation for phyto conc2 = 2.00000000000000012E-009 Max iron half saturation for phyto conc2m = 4.00000000000000025E-009 Iron half saturation for diatoms conc3 = 2.99999999999999998E-009 Maxi iron half saturation for diatoms conc3m = 8.99999999999999953E-009 Minimum size criteria for diatoms xsizedia = 4.99999999999999977E-007 Minimum size criteria for nanophyto xsizephy = 9.99999999999999955E-007 NH4 half saturation for phyto concnnh4 = 9.99999999999999955E-008 NH4 half saturation for diatoms concdnh4 = 3.99999999999999982E-007 Fe half saturation for bacteria concfebac = 2.99999999999999998E-011 optimal Fe quota for nano. qnfelim = 6.99999999999999989E-006 Optimal Fe quota for diatoms qdfelim = 6.99999999999999989E-006 Namelist parameters for phytoplankton growth, nampisprod ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ Enable new parame. of production (T/F) ln_newprod = F mean Si/C ratio grosip = 0.15100000000000000 P-I slope pislope = 3.0000000000000000 excretion ratio of nanophytoplankton excret = 5.00000000000000028E-002 excretion ratio of diatoms excret2 = 5.00000000000000028E-002 P-I slope for diatoms pislope2 = 2.0000000000000000 Minimum Chl/C in nanophytoplankton chlcnm = 3.30000000000000016E-002 Minimum Chl/C in diatoms chlcdm = 4.00000000000000008E-002 Maximum Fe/C in nanophytoplankton fecnm = 4.00000000000000033E-005 Minimum Fe/C in diatoms fecdm = 4.00000000000000033E-005 Namelist parameters for remineralization, nampisrem ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ remineralisation rate of POC xremip = 2.50000000000000014E-002 remineralization rate of DOC xremik = 0.25000000000000000 remineralization rate of Si xsirem = 3.00000000000000006E-003 fast remineralization rate of Si xsiremlab = 2.50000000000000014E-002 fraction of labile biogenic silica xsilab = 0.31000000000000000 scavenging rate of Iron xlam1 = 5.00000000000000010E-003 NH4 nitrification rate nitrif = 5.00000000000000028E-002 halk saturation constant for anoxia oxymin = 9.99999999999999955E-007 ligand concentration in the ocean ligand = 5.99999999999999996E-010 Namelist parameters for phytoplankton mortality, nampismort ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ quadratic mortality of phytoplankton wchl = 1.00000000000000002E-003 maximum quadratic mortality of diatoms wchld = 2.00000000000000004E-002 phytoplankton mortality rate mprat = 1.00000000000000002E-002 Diatoms mortality rate mprat2 = 1.00000000000000002E-002 Phytoplankton minimum mortality rate mpratm = 1.00000000000000002E-002 Namelist parameters for microzooplankton, nampiszoo ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ part of calcite not dissolved in microzoo guts part = 0.50000000000000000 microzoo preference for POC xpref2c = 0.10000000000000001 microzoo preference for nano xpref2p = 1.0000000000000000 microzoo preference for diatoms xpref2d = 0.59999999999999998 diatoms feeding threshold for microzoo xthreshdia = 1.00000000000000002E-008 nanophyto feeding threshold for microzoo xthreshphy = 1.99999999999999991E-007 poc feeding threshold for microzoo xthreshpoc = 1.00000000000000002E-008 feeding threshold for microzooplankton xthresh = 0.0000000000000000 exsudation rate of microzooplankton resrat = 2.99999999999999989E-002 microzooplankton mortality rate mzrat = 0.0000000000000000 maximal microzoo grazing rate grazrat = 3.0000000000000000 non assimilated fraction of P by microzoo unass = 0.29999999999999999 Efficicency of microzoo growth epsher = 0.33000000000000002 Fraction of microzoo excretion as DOM sigma1 = 0.59999999999999998 half sturation constant for grazing 1 xkgraz = 2.00000000000000016E-005 Namelist parameters for mesozooplankton, nampismes ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ part of calcite not dissolved in mesozoo guts part2 = 0.75000000000000000 mesozoo preference for phyto xprefc = 1.0000000000000000 mesozoo preference for POC xprefp = 0.29999999999999999 mesozoo preference for zoo xprefz = 1.0000000000000000 mesozoo preference for poc xprefpoc = 0.29999999999999999 microzoo feeding threshold for mesozoo xthresh2zoo = 1.00000000000000002E-008 diatoms feeding threshold for mesozoo xthresh2dia = 1.00000000000000002E-008 nanophyto feeding threshold for mesozoo xthresh2phy = 1.99999999999999991E-007 poc feeding threshold for mesozoo xthresh2poc = 1.00000000000000002E-008 feeding threshold for mesozooplankton xthresh2 = 0.0000000000000000 exsudation rate of mesozooplankton resrat2 = 5.00000000000000010E-003 mesozooplankton mortality rate mzrat2 = 2.99999999999999989E-002 maximal mesozoo grazing rate grazrat2 = 0.69999999999999996 mesozoo flux feeding rate grazflux = 3000.0000000000000 non assimilated fraction of P by mesozoo unass2 = 0.29999999999999999 Efficicency of Mesozoo growth epsher2 = 0.33000000000000002 Fraction of mesozoo excretion as DOM sigma2 = 0.59999999999999998 half sturation constant for grazing 2 xkgraz2 = 2.00000000000000016E-005 namelist : nampissed ~~~~~~~~~~~~~~~~~ dust input from the atmosphere ln_dust = F river input of nutrients ln_river = F atmospheric deposition of n ln_ndepo = F fe input from sediments ln_sedinput = F coastal release of iron sedfeinput = 1.00000000000000006E-009 solubility of the dust dustsolub = 2.00000000000000004E-002 sinking speed of the dust wdust = 2.0000000000000000 nitrogen fixation rate nitrfix = 9.99999999999999955E-008 nitrogen fixation sensitivty to light diazolight = 50.000000000000000 fe half-saturation cste for diazotrophs concfediaz = 1.00000000000000004E-010 Total input of elements from river supply ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ N Supply : 0.0000000000000000 TgN/yr Si Supply : 0.0000000000000000 TgSi/yr Alk Supply : 0.0000000000000000 Teq/yr DIC Supply : 0.0000000000000000 TgC/yr Total input of elements from atmospheric supply ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ N Supply : 0.0000000000000000 TgN/yr Namelist parameters for CaCO3 dissolution, nampiscal ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ diss. rate constant calcite (per month) kdca = 6.0000000000000000 order of reaction for calcite dissolution nca = 1.0000000000000000 Namelist parameters for air-sea exchange, nampisext ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ Choice for reading in the atm pCO2 file or constant value, ln_co2int = F Constant Atmospheric pCO2 value atcco2 = 287.00000000000000 Namelist nampisatm : Atmospheric Pressure as external forcing constant atmopsheric pressure (F) or from a file (T) ln_presatm = F Initialization of PISCES tracers done trc_rst_read : read data in the TOP restart file ~~~~~~~~~~~~ trc_rst_read_pisces : Read specific variables from pisces model ~~~~~~~~~~~~~~ *** Total number of passive tracer jptra = 24 *** Total volume of ocean = 1.32653272679519360E+018 *** Total inital content of all tracers tracer nb : 1 name : DIC initial content : 0.0000000000E+00 tracer nb : 2 name : Alkalini initial content : 0.0000000000E+00 tracer nb : 3 name : O2 initial content : 0.0000000000E+00 tracer nb : 4 name : CaCO3 initial content : 0.0000000000E+00 tracer nb : 5 name : PO4 initial content : 0.0000000000E+00 tracer nb : 6 name : POC initial content : 0.0000000000E+00 tracer nb : 7 name : Si initial content : 0.0000000000E+00 tracer nb : 8 name : PHY initial content : 0.0000000000E+00 tracer nb : 9 name : ZOO initial content : 0.0000000000E+00 tracer nb : 10 name : DOC initial content : 0.0000000000E+00 tracer nb : 11 name : PHY2 initial content : 0.0000000000E+00 tracer nb : 12 name : ZOO2 initial content : 0.0000000000E+00 tracer nb : 13 name : DSi initial content : 0.0000000000E+00 tracer nb : 14 name : Fer initial content : 0.0000000000E+00 tracer nb : 15 name : BFe initial content : 0.0000000000E+00 tracer nb : 16 name : GOC initial content : 0.0000000000E+00 tracer nb : 17 name : SFe initial content : 0.0000000000E+00 tracer nb : 18 name : DFe initial content : 0.0000000000E+00 tracer nb : 19 name : GSi initial content : 0.0000000000E+00 tracer nb : 20 name : NFe initial content : 0.0000000000E+00 tracer nb : 21 name : NCHL initial content : 0.0000000000E+00 tracer nb : 22 name : DCHL initial content : 0.0000000000E+00 tracer nb : 23 name : NO3 initial content : 0.0000000000E+00 tracer nb : 24 name : NH4 initial content : 0.0000000000E+00 dia_ptr_init : poleward transport and msf initialization ~~~~~~~~~~~~ Namelist namptr : set ptr parameters Poleward heat & salt transport (T) or not (F) ln_diaptr = F Overturning heat & salt transport ln_ptrcomp = T T & S zonal mean and meridional stream function ln_diaznl = T Global (F) or glo/Atl/Pac/Ind/Indo-Pac basins ln_subbas = T Frequency of computation nn_fptr = 1 Frequency of outputs nn_fwri = 15 dia_hsb_init : check the heat and salt budgets ~~~~~~~~~~~~ Namelist namhsb : set hsb parameters Switch for hsb diagnostic (T) or not (F) ln_diahsb = F Euler time step switch is 1 AAAAAAAA AAAAAAAA ===>>> : E R R O R =========== 3 error have been found