About to run experiment CM2.1p1 with model CM2M at Mon Jul 18 16:59:57 JST 2022. The command is: mpiexec.hydra -n 45 /home/t23598/dogarm/MOM5_Copy/MOM5/exec/hu/CM2M/fms_CM2M.x >fms.out NOTE from PE 0: MPP_DOMAINS_SET_STACK_SIZE: stack size set to 32768. &MPP_IO_NML HEADER_BUFFER_VAL = 16384, GLOBAL_FIELD_ON_ROOT_PE = T, IO_CLOCKS_ON = F, SHUFFLE = 0, DEFLATE_LEVEL = -1, CF_COMPLIANCE = F / NOTE from PE 0: MPP_IO_SET_STACK_SIZE: stack size set to 131072. NOTE from PE 0: MPP_DOMAINS_SET_STACK_SIZE: stack size set to 130296. 20220718 170004.037: Memuse(MB) at fms_init= 8.714E+01 9.712E+01 2.980E+00 9.114E+01 Entering coupler_init at 20220718 170004.038 Starting initializing ensemble_manager at 20220718 170004.042 Finished initializing ensemble_manager at 20220718 170004.043 NOTE from PE 0: coupler_init: Atmos PE range: 0 14 ens_01 NOTE from PE 0: coupler_init: Ocean PE range: 15 44 ens_01 NOTE from PE 0: coupler_init: Land PE range: 0 14 ens_01 NOTE from PE 0: coupler_init: Ice PE range: 0 14 ens_01 NOTE from PE 0: coupler_init: Running with CONCURRENT coupling. NOTE from PE 0: coupler_init: Sending LAG fluxes to ocean. Starting to initialize diag_manager at 20220718 170004.054 Starting to initialize diag_manager at 20220718 170004.054 NOTE from PE 0: diag_manager_mod::diag_manager_init: prepend_date only supported when diag_manager_init is called with time_init present. NOTE from PE 15: diag_manager_mod::diag_manager_init: prepend_date only supported when diag_manager_init is called with time_init present. 20220718 170004.074: Memuse(MB) at diag_manager_init= 9.209E+01 1.035E+02 3.434E+00 9.723E+01 Finished initializing diag_manager at 20220718 170004.074 20220718 170004.081: Memuse(MB) at diag_manager_init= 9.299E+01 1.002E+02 2.322E+00 9.577E+01 Finished initializing diag_manager at 20220718 170004.081 WARNING from PE 0: MPP_OPEN: File time_stamp.out opened WRONLY already exists! WARNING from PE 0: MPP_OPEN: File time_stamp.out opened WRONLY already exists! Starting to initialize tracer_manager at 20220718 170004.084 Starting to initialize ocean model at 20220718 170004.091 MOM_COMMIT_HASH= &OCEAN_PARAMETERS_NML CP_OCEAN = 3992.10322329649 , CP_LIQUID_RUNOFF = 4218.00000000000 , CP_SOLID_RUNOFF = 2106.00000000000 , RHO0 = 1035.00000000000 , TFREEZE = 273.150000000000 , OMEGA_EARTH = 7.292100000000000E-005, GRAV = 9.80000000000000 / Finished initializing tracer_manager at 20220718 170004.093 Starting to initialize coupler_types at 20220718 170004.093 ======== STARTING MOM INITIALIZATION ======== ==>NOTE: Using dynamically allocated array option in MOM Dumping coupler_mod/types tree types/ air_sea_gas_flux_generic/ implementation/ ocmip2/ num_parameters = 2 num_flags = 0 use_atm_pressure = T use_10m_wind_speed = T pass_through_ice = F atm/ name[1] = 'pcair' name[2] = 'u10' name[3] = 'psurf' long_name[1] = 'Atmospheric concentration' long_name[2] = 'Wind speed at 10 m' long_name[3] = 'Surface atmospheric pressure' units[1] = 'mol/mol' units[2] = 'm/s' units[3] = 'Pa' ice/ name[1] = 'alpha' name[2] = 'csurf' name[3] = 'sc_no' long_name[1] = 'Solubility w.r.t. atmosphere' long_name[2] = 'Ocean concentration' long_name[3] = 'Schmidt number' units[1] = 'mol/m^3/atm' units[2] = 'mol/m^3' units[3] = 'dimensionless' flux/ name[1] = 'flux' name[2] = 'deltap' name[3] = 'kw' long_name[1] = 'Surface flux' long_name[2] = 'Ocean-air delta pressure' long_name[3] = 'Piston velocity' units[1] = 'mol/m^2/s' units[2] = 'uatm' units[3] = 'm/s' air_sea_gas_flux/ implementation/ ocmip2/ num_parameters = 2 ocmip2_data/ num_parameters = 2 linear/ num_parameters = 3 num_flags = 0 use_atm_pressure = T use_10m_wind_speed = T pass_through_ice = F atm/ name[1] = 'pcair' name[2] = 'u10' name[3] = 'psurf' long_name[1] = 'Atmospheric concentration' long_name[2] = 'Wind speed at 10 m' long_name[3] = 'Surface atmospheric pressure' units[1] = 'mol/mol' units[2] = 'm/s' units[3] = 'Pa' ice/ name[1] = 'alpha' name[2] = 'csurf' long_name[1] = 'Solubility from atmosphere times Schmidt number term' long_name[2] = 'Ocean concentration times Schmidt number term' units[1] = 'mol/m^3/atm' units[2] = 'mol/m^3' flux/ name = 'flux' long_name = 'Surface flux' units = 'mol/m^2/s' air_sea_deposition/ implementation/ dry/ num_parameters = 1 wet/ num_parameters = 1 num_flags = 0 use_atm_pressure = F use_10m_wind_speed = F pass_through_ice = T atm/ name = 'deposition' long_name = 'Atmospheric deposition' units = 'kg/m^2/s' ice/ name = NULL long_name = NULL units = NULL flux/ name = 'flux' long_name = 'Surface deposition' units = 'mol/m^2/s' land_sea_runoff/ implementation/ river/ num_parameters = 1 num_flags = 0 use_atm_pressure = F use_10m_wind_speed = F pass_through_ice = T atm/ name = 'runoff' long_name = 'Concentration in land runoff' units = 'mol/m^3' ice/ name = NULL long_name = NULL units = NULL flux/ name = 'flux' long_name = 'Concentration in land runoff' units = 'mol/m^3' Finished initializing coupler_types at 20220718 170004.094 Beginning to initialize component models at 20220718 170004.094 Starting to initialize atmospheric model at 20220718 170004.094 ncnst= 4 num_prog= 4 pnats= 0 num_family= 0 &OCEAN_MODEL_NML TIME_TENDENCY = twolevel , IMPOSE_INIT_FROM_RESTART = F, REINITIALIZE_THICKNESS = F, BAROCLINIC_SPLIT = 1, BAROTROPIC_SPLIT = 80, SURFACE_HEIGHT_SPLIT = 1, LAYOUT = 6, 5, IO_LAYOUT = 2*0, DEBUG = F, VERTICAL_COORDINATE = geopotential , DT_OCEAN = 7200, CMIP_UNITS = F, CMIP_VERSION = 5, HORIZONTAL_GRID = bgrid , USE_BLOBS = F, USE_VELOCITY_OVERRIDE = F, MASK_TABLE = INPUT/ocean_mask_table , INTRODUCE_BLOBS = F, BETA_TXTY = 0.000000000000000E+000, BETA_TF = 0.000000000000000E+000, BETA_QF = 0.000000000000000E+000, BETA_LWSW = 0.000000000000000E+000, DO_WAVE = F / ==>Note: Running MOM using 30 computer processors. ==>Note: Running MOM with staggered twotime level scheme to compute time tenden cies. This is the default. Mass/volume and tracer are conserved with this sc heme. ==>Note: Time%Time_init = time stamp at very start of the MOM experiment is given by yyyy/mm/dd hh:mm:ss = 1/ 1/ 1 0: 0: 0 ==>Note: Time%model_time = time stamp at start of this leg of the MOM experiment is yyyy/mm/dd hh:mm:ss = 1/ 1/ 1 0: 0: 0 ==>Note: Time%init=.true. =>MOM will start from user specified initial conditions. ==> Note: time steps (seconds) used for MOM dtts (tracer) = 7200.00 dtuv (baroclinic) = 7200.00 dteta (surface height or bottom pressure) = 7200.00 dtbt (barotropic) = 90.00 ==> Note: The velocity equations will be split into baroclinic and barotropic pieces. ==> Note: Using MOM with B-grid horizontal layout of variables on a grid. ==> Note: Using MOM with geopotential vertical coordinate. Beware of vanishing top model grid cells. The equations are Boussinesq, and so conserve volume rather than mass. Use one of the pressure-like coordinates to get non-Boussinesq effects. FV Core Arrays. domain decomposition whalo = 0, ehalo = 0, shalo = 3, nhalo = 3 X-AXIS = 144 Y-AXIS = 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 MPP_PSET_CREATE creating PSETs... npset= 1 MPP_PSET_CREATE: setting stacksize= 1000000 &OCEAN_GRIDS_NML DEBUG_THIS_MODULE = T, VERBOSE_INIT = T, READ_RHO0_PROFILE = F, DO_BITWISE_EXACT_SUM = F, WRITE_GRID = F / MPI-1 is used within mod_comm ==>Note from ocean_grids_mod(set_ocean_grid_size): read grid from mosaic versio n grid Two-level Parallelism: Domains x Threads= 15 1 SMP pseudo domain decomposition in E-W = 1 Cold_start= F Using n_split from the namelist: 5 Total number of sponge layers= 1 Mean cell width (km)= 221.574005894595 Initializing polar filter ......... Starting latitude for algebraic pft= 36.0027517282874 Tracer sphum initialized with surface value of 0.300000E-05 and vertical multiplier of 1.000000 NOTE from PE 15: ==>Note from ocean_grids_mod(set_ocean_grid_size): x_boundary_type is cyclic NOTE from PE 15: ==>Note from ocean_grids_mod(set_ocean_grid_size): y_boundary_type is fold_north_edge ==> Note: Energy conversion errors are nontrivial when using tripolar=.true. The cause is related to the need to update redundantly computed information across the Arctic bipolar fold in a bit-wise exact manner for terms contributing to the energy conversion analysis. The extra code and mpp calls have not been implemented. Restart days= 36500 sec= 0 ---------------------------------- Model Top Pressure (pa)= 100.000000000000 ---------------------------------- Checking Initial condition ... &OCEAN_DOMAINS_NML HALO = 1, MAX_TRACERS = 10, X_CYCLIC_OFFSET = 0, Y_CYCLIC_OFFSET = 0 / mom_domain domain decomposition whalo = 1, ehalo = 1, shalo = 1, nhalo = 1 X-AXIS = 60 60 60 60 60 60 Y-AXIS = 40 40 40 40 40 NOTE from PE 15: mpp_domains_define.inc(mpp_define_io_domain): io domain will not be defined for mom_domain when one or both entry of io_layout is not positive NOTE from PE 15: MPP_DOMAINS_SET_STACK_SIZE: stack size set to 2604000. ocean_super_grid domain decomposition whalo = 2, ehalo = 2, shalo = 2, nhalo = 2 X-AXIS = 120 120 120 120 120 120 Y-AXIS = 80 80 80 80 80 ZS max = 5359.05515463059 min = 0.000000000000000E+000 U max = 102.548035456680 min = -43.4099281670030 V max = 62.7091559302077 min = -55.9461084689003 T max = 311.248683399278 min = 186.069798388561 Q max = 2.314804662529603E-002 min = 9.080161437986342E-008 Q max = 1.067320255585952E-003 min = -4.340937542639236E-021 Q max = 2.072445831417558E-003 min = -3.096888368834090E-021 Q max = 1.00070428169025 min = -3.952117007028043E-014 Total surface pressure (mb) = 984.988152396148 mean dry surface pressure = 982.879999999965 TPW-vapor (kg/m**2) = 21.5117591447274 TPW-total (kg/m**2) = 21.6002167279660 u_srf max = 19.6880672679790 min = -16.1918077199760 v_srf max = 18.8751037121480 min = -22.4035141396967 Check sums for FV input data: chksum (t) = 230A73054496730E chksum (u) = 2BA728B05B3995E2 chksum (v) = 3EF1219D3B3CFD76 chksum (r) = 1E389B5DD9017895 NOTE from PE 0: read_fv_rst: FV netcdf restart file read. Check_eta: allowable min_ps (mb)= 385.413553452405 NOTE from PE 0: strat_cloud_mod: Reading netCDF formatted restart file: INPUT/strat_cloud.res.nc NOTE from PE 15: ==>Note from ocean_grids_mod (set_ocean_hgrid_arrays): altering U-grid arrays at j=0 [chksum] xt 745864707818127360 [chksum] xu 745820727353016320 [chksum] yt -8722290178805985338 [chksum] yu -8723404134886878374 [chksum] dxt -7205131407146550824 [chksum] dxu -7203630101784009518 [chksum] dyt -7133734828288020206 [chksum] dyu -7133846536504224530 [chksum] dat -4971524625182553254 [chksum] dau -4969384579357392896 [chksum] dxtn -7203630101784009518 [chksum] dytn -7133846536504224530 [chksum] dxte -7205131407146550824 [chksum] dyte -7133734828288020206 [chksum] dxun -7202138806629264056 [chksum] dyun -7133972903803456088 [chksum] dxue -7203630101784009518 [chksum] dyue -7133846536504224530 [chksum] dtn+dts -7133734828288020206 [chksum] dun+dus -7133846536504224530 [chksum] dte+dtw -7205131407146550824 [chksum] due+duw -7203630101784009518 [chksum] dte -7340239395967665704 [chksum] dtw -7340239395967665704 [chksum] due -7338738090605124398 [chksum] duw -7338738090605124398 [chksum] sin_rot 0 [chksum] cos_rot 9088264048033660928 &OCEAN_TOPOG_NML FLAT_BOTTOM = F, FLAT_BOTTOM_KMT = 50, FLAT_BOTTOM_HT = 5500.00000000000 , WRITE_TOPOG = F, MIN_THICKNESS = 5.00000000000000 , KMT_RECOMPUTE = F, KMT_RECOMPUTE_OFFSET = 0, DEBUG_THIS_MODULE = T / The shallowest wet ocean model grid cell has depth (meters) 40.00000 and this occurs at (i,j,k) = ( 59, 16, 4) which has (long,lat,depth) = ( -221.5000, -66.5000, 40.00000) Beware that shallow regions (e.g., those shallower than 50m) may be subject to numerical problems if strong surface forcing is not mixed vertically. Such problems may occur especially in shallow regions with kmt==2. Current speeds and/or tracer deviations may become large due to the deposition of wind and/or buoyancy over just a small upper ocean region. Such problems can be resolved by adding sufficient vertical mixing in these regions. Such happens in Nature due to tides and breaking surface waves. Topography checksums [chksum] ht 792435832758655498 [chksum] hu 788729280302577333 [chksum] htr -945180253257382874 [chksum] kmu 750 [chksum] kmt 756 &OCEAN_OBC_NML NOBC = 0, DIRECTION = , NAME = test_obc none none none , IS = 4*-999, IE = 4*-999, JS = 4*-999, JE = 4*-999, IERS = 4*-999, IERE = 4*-999, JERS = 4*-999, JERE = 4*-999, ITRS = 4*-999, ITRE = 4*-999, JTRS = 4*-999, JTRE = 4*-999, OBC_NOR = NOGRAD NOGRAD NOGRAD NOGRAD , OBC_TAN = NOGRAD NOGRAD NOGRAD NOGRAD , OBC_ETA = NOTHIN NOTHIN NOTHIN NOTHIN , OBC_UD = NOGRAD NOGRAD NOGRAD NOGRAD , OBC_TRA = NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NO GRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD NOGRAD , OBC_MIX = NOGRAD NOGRAD NOGRAD NOGRAD , REL_COEF_ETA_IN = 4*0.000000000000000E+000 , REL_COEF_ETA_OUT = 4*0.000000000000000E+000 , REL_ETA_PNTS = 4*1, REL_CLIN_PNTS = 160*1, CTROP_MAX = 4*1.50000000000000 , CTROP_MIN = 4*0.100000000000000 , CTROP_INC = 4*0.000000000000000E+000 , CTROP_SMOOTH = 4*0.700000000000000 , FILENAME_ETA = obc_eta_t.nc none none none , FIELDNAME_ETA = eta_t none none none , FILENAME_UD = obc_ud.nc none none none , FIELDNAME_UD = ud none none none , OBC_CONSIDER_CONVU = 4*F, OBC_ADJUST_FORCING_BT = 4*F, OBC_VERT_ADVEL_T = 4*F, OBC_VERT_ADVEL_U = 4*F, OBC_ENHANCE_VISC_BACK = NONE NONE NONE NONE , OBC_ENHANCE_DIFF_BACK = NONE NONE NONE NONE , ENH_PNTS = 4*1, ENH_FAC_V = 4*0.900000000000000 , ENH_FAC_D = 4*1.00000000000000 , OBC_RELAX_TRACER = 160*F, OBC_FLOW_RELAX = 160*1, OBC_CONSIDER_SOURCES = 160*F, OBC_TRACER_NO_INFLOW = 160*F, REL_COEF_TRACER_IN = 160*0.000000000000000E+000 , REL_COEF_TRACER_OUT = 160*0.000000000000000E+000 , FILENAME_TRACER = INPUT/obc_tr.nc INPUT/obc_tr.nc INPUT/obc_tr.nc INPUT/obc_tr.nc INPUT/obc_tr.nc INPUT/obc_tr.nc INPUT/obc_tr.nc INPUT/obc_tr.nc none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none no ne none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none , FIELDNAME_TRACER = temp_obc temp_obc temp_obc temp_obc salt_obc salt_obc salt_obc salt_obc none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none none no ne none none none none none none none none none none none none none none none none none none none , DEBUG_PHASE_SPEED = F, DEBUG_THIS_MODULE = F, OBC_DAMP_NEWTON = 4*F, DAMP_FACTOR = 4*1.00000000000000 / Number of wet ocean tracer points = 2050088 Number of wet ocean velocity points = 1956338 Number of computed ocean tracer points = 3600000 Wet ocean volume with eta_t=0.0 (T-cells) = 1.32536324265417318E+18 m^3 (not bit reproducible) Ocean surface area (T-cells) = 3.59027326205200875E+14 m^2 (not bit reproducible) Wet ocean volume with eta_u=0.0 (U-cells) = 1.23840706957569946E+18 m^3 (not bit reproducible) Wet ocean surface area (U-cells) = 3.45161542346050062E+14 m^2 (not bit reproducible) Wet ocean + masked-out (land) surface area (T-cells) = 5.07591895230173312E+14 m^2 (not bit reproducible) Wet ocean + masked-out (land) surface area (U-cells) = 5.08013684969460000E+14 m^2 (not bit reproducible) &OCEAN_TRACER_UTIL_NML REBIN_ONTO_RHO_ALL_VALUES = T, DEBUG_DIAGNOSE_MASS_OF_LAYER = F, EPSLN_DIAGNOSE_MASS_OF_LAYER = 1.000000000000000E-005 / Note: ocean_tracer_util: rebin_onto_rho will include density values outside bounds range. &OCEAN_CORIOLIS_NML DEBUG_THIS_MODULE = F, USE_THIS_MODULE = T, ACOR = 0.500000000000000 / ==>Note: Coriolis on Bgrid computed semi-implicitly to remove inertial time step constraint. NOTE from PE 0: topography_mod: Reading NetCDF formatted input data file: INPUT/navy_topography.data.nc NOTE from PE 0: MPP_IO_SET_STACK_SIZE: stack size set to 2332800. ==> Note: 2*pi timesteps/(min inertial period) implies a maximum dtuv for time-explicit Coriolis = 6856. sec. &OCEAN_VELOCITY_NML DEBUG_THIS_MODULE = F, WRITE_A_RESTART = T, MAX_CGINT = 2.00000000000000 , ZERO_TENDENCY = F, ZERO_TENDENCY_EXPLICIT_A = F, ZERO_TENDENCY_EXPLICIT_B = F, ZERO_TENDENCY_IMPLICIT = F, TRUNCATE_VELOCITY = F, TRUNCATE_VERBOSE = T, TRUNCATE_VELOCITY_LAT = 0.000000000000000E+000, TRUNCATE_VELOCITY_VALUE = 2.00000000000000 , ADAMS_BASHFORTH_THIRD = T, ADAMS_BASHFORTH_EPSILON = 0.600000000000000 , USE_CONSTANT_VELOCITY = F, CONSTANT_U = 0.000000000000000E+000, CONSTANT_V = 0.000000000000000E+000, UPDATE_VELOCITY_VIA_UPRIME = T / ==>Note from ocean_velocity_mod: use of twolevel time_tendency necessitates an Adams-Bashforth treatment of velocity advection. Using 3rd order Adams-Bashforth for velocity advection. This is the MOM default. NOTE from PE 15: ==>update_velocity_via_uprime=.true., so keep udrho from external mode solver. Baroclinic time step stability most nearly violated at U-cell (i,j) = ( 90, 200), (lon,lat) = ( -100.00, 90.00). The number of kmu-levels at this point is 45 The dxu grid distance (m) at this point is 0.246515E+05 The dyu grid distance (m) at this point is 0.111202E+06 Due to a specified maximum baroclinic gravity wave speed of 2.00 m/s. "dtuv" must be less than 12034. sec. "dtuv" = 7200. sec. NOTE from PE 0: aerosol_mod: inconsistent nml settings -- not using aerosol timeseries but requesting interannual variation of aerosol amount for so4_anthro -- this aerosol will NOT exhibit interannual variation NOTE from PE 0: aerosol_mod: inconsistent nml settings -- not using aerosol timeseries but requesting interannual variation of aerosol amount for so4_natural -- this aerosol will NOT exhibit interannual variation NOTE from PE 0: aerosol_mod: inconsistent nml settings -- not using aerosol timeseries but requesting interannual variation of aerosol amount for organic_carbon -- this aerosol will NOT exhibit interannual variation NOTE from PE 0: aerosol_mod: inconsistent nml settings -- not using aerosol timeseries but requesting interannual variation of aerosol amount for black_carbon -- this aerosol will NOT exhibit interannual variation NOTE from PE 0: aerosol_mod: inconsistent nml settings -- not using aerosol timeseries but requesting interannual variation of aerosol amount for sea_salt -- this aerosol will NOT exhibit interannual variation NOTE from PE 0: aerosol_mod: inconsistent nml settings -- not using aerosol timeseries but requesting interannual variation of aerosol amount for anthro_dust_0.1 -- this aerosol will NOT exhibit interannual variation NOTE from PE 0: aerosol_mod: inconsistent nml settings -- not using aerosol timeseries but requesting interannual variation of aerosol amount for anthro_dust_0.2 -- this aerosol will NOT exhibit interannual variation NOTE from PE 0: aerosol_mod: inconsistent nml settings -- not using aerosol timeseries but requesting interannual variation of aerosol amount for anthro_dust_0.4 -- this aerosol will NOT exhibit interannual variation NOTE from PE 0: aerosol_mod: inconsistent nml settings -- not using aerosol timeseries but requesting interannual variation of aerosol amount for anthro_dust_0.8 -- this aerosol will NOT exhibit interannual variation NOTE from PE 0: aerosol_mod: inconsistent nml settings -- not using aerosol timeseries but requesting interannual variation of aerosol amount for anthro_dust_1.0 -- this aerosol will NOT exhibit interannual variation NOTE from PE 0: aerosol_mod: inconsistent nml settings -- not using aerosol timeseries but requesting interannual variation of aerosol amount for anthro_dust_2.0 -- this aerosol will NOT exhibit interannual variation NOTE from PE 0: aerosol_mod: inconsistent nml settings -- not using aerosol timeseries but requesting interannual variation of aerosol amount for anthro_dust_4.0 -- this aerosol will NOT exhibit interannual variation NOTE from PE 0: aerosol_mod: inconsistent nml settings -- not using aerosol timeseries but requesting interannual variation of aerosol amount for anthro_dust_8.0 -- this aerosol will NOT exhibit interannual variation NOTE from PE 0: aerosol_mod: inconsistent nml settings -- not using aerosol timeseries but requesting interannual variation of aerosol amount for natural_dust_0.1 -- this aerosol will NOT exhibit interannual variation NOTE from PE 0: aerosol_mod: inconsistent nml settings -- not using aerosol timeseries but requesting interannual variation of aerosol amount for natural_dust_0.2 -- this aerosol will NOT exhibit interannual variation NOTE from PE 0: aerosol_mod: inconsistent nml settings -- not using aerosol timeseries but requesting interannual variation of aerosol amount for natural_dust_0.4 -- this aerosol will NOT exhibit interannual variation NOTE from PE 0: aerosol_mod: inconsistent nml settings -- not using aerosol timeseries but requesting interannual variation of aerosol amount for natural_dust_0.8 -- this aerosol will NOT exhibit interannual variation NOTE from PE 0: aerosol_mod: inconsistent nml settings -- not using aerosol timeseries but requesting interannual variation of aerosol amount for natural_dust_1.0 -- this aerosol will NOT exhibit interannual variation NOTE from PE 0: aerosol_mod: inconsistent nml settings -- not using aerosol timeseries but requesting interannual variation of aerosol amount for natural_dust_2.0 -- this aerosol will NOT exhibit interannual variation NOTE from PE 0: aerosol_mod: inconsistent nml settings -- not using aerosol timeseries but requesting interannual variation of aerosol amount for natural_dust_4.0 -- this aerosol will NOT exhibit interannual variation NOTE from PE 0: aerosol_mod: inconsistent nml settings -- not using aerosol timeseries but requesting interannual variation of aerosol amount for natural_dust_8.0 -- this aerosol will NOT exhibit interannual variation NOTE from PE 0: aerosol_mod: Aerosol data is defined from a single annual cycle for so4_anthro - no interannual variation Aerosol data for so4_anthro obtained from aerosol timeseries for year: 1990 NOTE from PE 0: aerosol_mod: Aerosol data is defined from a single annual cycle for so4_natural - no interannual variation Aerosol data for so4_natural obtained from aerosol timeseries for year: 1990 NOTE from PE 0: aerosol_mod: Aerosol data is defined from a single annual cycle for organic_carbon - no interannual variation Aerosol data for organic_carbon obtained from aerosol timeseries for year: 1990 NOTE from PE 0: aerosol_mod: Aerosol data is defined from a single annual cycle for black_carbon - no interannual variation Aerosol data for black_carbon obtained from aerosol timeseries for year: 1990 NOTE from PE 0: aerosol_mod: Aerosol data is defined from a single annual cycle for sea_salt - no interannual variation Aerosol data for sea_salt obtained from aerosol timeseries for year: 1990 NOTE from PE 0: aerosol_mod: Aerosol data is defined from a single annual cycle for anthro_dust_0.1 - no interannual variation Aerosol data for anthro_dust_0.1 obtained from aerosol timeseries for year: 1990 NOTE from PE 0: aerosol_mod: Aerosol data is defined from a single annual cycle for anthro_dust_0.2 - no interannual variation Aerosol data for anthro_dust_0.2 obtained from aerosol timeseries for year: 1990 NOTE from PE 0: aerosol_mod: Aerosol data is defined from a single annual cycle for anthro_dust_0.4 - no interannual variation Aerosol data for anthro_dust_0.4 obtained from aerosol timeseries for year: 1990 NOTE from PE 0: aerosol_mod: Aerosol data is defined from a single annual cycle for anthro_dust_0.8 - no interannual variation Aerosol data for anthro_dust_0.8 obtained from aerosol timeseries for year: 1990 NOTE from PE 0: aerosol_mod: Aerosol data is defined from a single annual cycle for anthro_dust_1.0 - no interannual variation Aerosol data for anthro_dust_1.0 obtained from aerosol timeseries for year: 1990 NOTE from PE 0: aerosol_mod: Aerosol data is defined from a single annual cycle for anthro_dust_2.0 - no interannual variation Aerosol data for anthro_dust_2.0 obtained from aerosol timeseries for year: 1990 NOTE from PE 0: aerosol_mod: Aerosol data is defined from a single annual cycle for anthro_dust_4.0 - no interannual variation Aerosol data for anthro_dust_4.0 obtained from aerosol timeseries for year: 1990 NOTE from PE 0: aerosol_mod: Aerosol data is defined from a single annual cycle for anthro_dust_8.0 - no interannual variation Aerosol data for anthro_dust_8.0 obtained from aerosol timeseries for year: 1990 NOTE from PE 0: aerosol_mod: Aerosol data is defined from a single annual cycle for natural_dust_0.1 - no interannual variation Aerosol data for natural_dust_0.1 obtained from aerosol timeseries for year: 1990 NOTE from PE 0: aerosol_mod: Aerosol data is defined from a single annual cycle for natural_dust_0.2 - no interannual variation Aerosol data for natural_dust_0.2 obtained from aerosol timeseries for year: 1990 NOTE from PE 0: aerosol_mod: Aerosol data is defined from a single annual cycle for natural_dust_0.4 - no interannual variation Aerosol data for natural_dust_0.4 obtained from aerosol timeseries for year: 1990 NOTE from PE 0: aerosol_mod: Aerosol data is defined from a single annual cycle for natural_dust_0.8 - no interannual variation Aerosol data for natural_dust_0.8 obtained from aerosol timeseries for year: 1990 NOTE from PE 0: aerosol_mod: Aerosol data is defined from a single annual cycle for natural_dust_1.0 - no interannual variation Aerosol data for natural_dust_1.0 obtained from aerosol timeseries for year: 1990 NOTE from PE 0: aerosol_mod: Aerosol data is defined from a single annual cycle for natural_dust_2.0 - no interannual variation Aerosol data for natural_dust_2.0 obtained from aerosol timeseries for year: 1990 NOTE from PE 0: aerosol_mod: Aerosol data is defined from a single annual cycle for natural_dust_4.0 - no interannual variation Aerosol data for natural_dust_4.0 obtained from aerosol timeseries for year: 1990 NOTE from PE 0: aerosol_mod: Aerosol data is defined from a single annual cycle for natural_dust_8.0 - no interannual variation Aerosol data for natural_dust_8.0 obtained from aerosol timeseries for year: 1990 NOTE from PE 0: GLOBAL ATT too long - not reading this metadata Reading TWO_LEVEL restart for velocity from INPUT/ocean_velocity.res.nc Expecting only one time record for each restart field. Finished reading restart for velocity field ===Initial velocity checksums ==> yyyy/mm/dd hh:mm:ss = 1/ 1/ 1 0: 0: 0 From ocean_velocity_mod: initial velocity chksum (taup1) [chksum] Zonal velocity 7256256170631703796 [chksum] Meridional velocity -340983864059611130 [chksum] Advection of u 0 [chksum] Advection of v 0 &OCEAN_BAROTROPIC_NML WRITE_A_RESTART = T, ZERO_TENDENCY = F, ZERO_ETA_IC = F, ZERO_ETA_T = F, ZERO_ETA_U = F, ZERO_CORIOLIS_BT = F, ZERO_NONLINEAR_FORCING_BT = F, ZERO_FORCING_BT = F, ZERO_ETA_TENDENCY = F, BAROTROPIC_TIME_STEPPING_A = T, BAROTROPIC_TIME_STEPPING_B = F, TIDAL_FORCING_M2 = F, TIDAL_FORCING_8 = F, TIDAL_FORCING_IDEAL = F, GEOID_FORCING = F, ALPHAT = 0.948000000000000 , PRED_CORR_GAMMA = 0.200000000000000 , SMOOTH_ETA_T_BT_LAPLACIAN = F, SMOOTH_ETA_T_BT_BIHARMONIC = F, SMOOTH_ETA_T_LAPLACIAN = F, SMOOTH_ETA_T_BIHARMONIC = T, SMOOTH_ANOMPB_BT_LAPLACIAN = F, SMOOTH_ANOMPB_BT_BIHARMONIC = F, SMOOTH_PBOT_T_LAPLACIAN = F, SMOOTH_PBOT_T_BIHARMONIC = T, SMOOTH_PBOT_T_BIHARMONIC_LEGACY = F, SMOOTH_ETA_DIAG_LAPLACIAN = T, SMOOTH_ETA_DIAG_BIHARMONIC = F, VEL_MICOM_LAP = 5.000000000000000E-002, VEL_MICOM_LAP_DIAG = 0.200000000000000 , VEL_MICOM_BIH = 1.000000000000000E-002, VEL_MICOM_BIH_DIAG = 0.100000000000000 , TRUNCATE_ETA = F, VERBOSE_TRUNCATE = T, ETA_MAX = 8.00000000000000 , FRAC_CRIT_CELL_HEIGHT = 0.200000000000000 , VERBOSE_INIT = T, DEBUG_THIS_MODULE = F, DIAG_STEP = 12, ETA_OFFSET = 1.000000000000000E-012, PBOT_OFFSET = 1.000000000000000E-012, INITSUM_WITH_BAR_MOM4P0 = F, INITSUM_WITH_BAR_MOM4P1 = T, IDEAL_INITIAL_ETA = F, IDEAL_INITIAL_ETA_AMPLITUDE = 5.00000000000000 , IDEAL_INITIAL_ETA_XWIDTH = 100000.000000000 , IDEAL_INITIAL_ETA_YWIDTH = 100000.000000000 , UDRHO_BT_LAP = F, UDRHO_BT_BIH = F, UDRHO_LAP = F, UDRHO_BIH = F, UDRHO_LAP_VEL_MICOM = 5.000000000000000E-002, UDRHO_BIH_VEL_MICOM = 1.000000000000000E-002, BAROTROPIC_HALO = 1, DO_BITWISE_EXACT_SUM = F, USE_LEGACY_BAROTROPIC_HALOS = T / horz diff flux domain decomposition whalo = 1, ehalo = 1, shalo = 1, nhalo = 1 X-AXIS = 60 60 60 60 60 60 Y-AXIS = 40 40 40 40 40 NOTE from PE 15: mpp_domains_define.inc(mpp_define_io_domain): io domain will not be defined for horz diff flux when one or both entry of io_layout is not positive Reading TWO_LEVEL restart from INPUT/ocean_barotropic.res.nc Expecting only one time record for each restart field. NOTE from PE 15: Using barotropic predictor-corrector for integrating barotropic dynamics. Predictor-Corrector time filter on barotropic dynamics has value= 0.20 Updating eta_t or pbot_t using a big time step as in MOM4.0. Not recommended for OBC applications. Initialise sum of barotropic sea level with eta_t or pbot_t. This is the default. ==> Note: The barotropic dynamics integrate 160 timesteps for every one baroclinic timestep. ==>Using smooth_eta_t_biharmonic to smooth eta_t. Can produce extrema. Not recommended. NOTE from PE 15: ==>Note from tidal_forcing_init: tidal_forcing=false, so not adding tidal forcing to external mode. Barotropic stability most nearly violated at T-cell (i,j) = ( 90, 200), (lon,lat) = ( -202.50, 89.49). The number of kmt-levels at this point is 47 The dxt grid spacing (m) at this point is 0.246606E+05 The dyt grid spacing (m) at this point is 0.111200E+06 where the barotropic gravity wave speed is ~232.0 m/s. "dtbt" must be less than 115.000 sec. dtbt = 90.000 sec. &OCEAN_THICKNESS_NML DEBUG_THIS_MODULE = F, DEBUG_THIS_MODULE_DETAIL = F, WRITE_A_RESTART = T, FULL_STEP_TOPOGRAPHY = F, INITIALIZE_ZERO_ETA = F, ENFORCE_POSITIVE_DZT = F, DEPTH_MIN_FOR_SIGMA = 1.000000000000000E-002, THICKNESS_METHOD = energetic , READ_RHO0_PROFILE = F, THICKNESS_DZT_MIN = 2.00000000000000 , THICKNESS_DZT_MIN_INIT = 10.0000000000000 , RESCALE_MASS_TO_GET_HT_MOD = F, PBOT0_SIMPLE = F, EPSILON_INIT_THICKNESS = 1.000000000000000E-005, READ_RESCALE_RHO0_MASK = F, RESCALE_RHO0_MASK_GFDL = F, RESCALE_RHO0_BASIN_LABEL = -1.00000000000000 , RESCALE_RHO0_VALUE = 1.00000000000000 , LINEAR_FREE_SURFACE = F, MAX_NUM_BAD_PRINT = 25, UPDATE_DZWU_K0 = F / ==>Note: running ocean_thickness with thickness_method=energetic. Reading TWO_LEVEL restart for rho_dzt from INPUT/ocean_thickness.res.nc Expect only one time record for rho_dzt. &OCEAN_OPERATORS_NML USE_LEGACY_DIV_UD = F / horz diff flux domain decomposition whalo = 1, ehalo = 1, shalo = 1, nhalo = 1 X-AXIS = 60 60 60 60 60 60 Y-AXIS = 40 40 40 40 40 NOTE from PE 15: mpp_domains_define.inc(mpp_define_io_domain): io domain will not be defined for horz diff flux when one or both entry of io_layout is not positive &OCEAN_TRACER_NML DEBUG_THIS_MODULE = F, ZERO_TENDENCY = F, ZERO_TRACER_SOURCE = F, WRITE_A_RESTART = T, OCEAN_TPM_DEBUG = F, TMASK_LIMIT_TS_SAME = T, REMAP_DEPTH_TO_S_INIT = F, INFLOW_NBOUNDARY = F, INTERPOLATE_TPROG_TO_PBOTT = T, INTERPOLATE_TDIAG_TO_PBOTT = F, LIMIT_AGE_TRACER = T, AGE_TRACER_MAX_INIT = 0.000000000000000E+000, FRAZIL_HEATING_BEFORE_VPHYSICS = T, FRAZIL_HEATING_AFTER_VPHYSICS = F, COMPUTE_TMASK_LIMIT_ON = T, USE_TEMPSALT_CHECK_RANGE = F / ==>Note: frazil heating called before vertical physics and before boundary fluxes. This method is retained for legacy purposes: it is NOT recommended for new runs. NOTE from PE 15: ==>Note from ocean_tracer_mod(ocean_prog_tracer_init): Synchronous timesteps have been specified (dtts = dtuv). ==>Note from ocean_tpm_util_mod(otpm_set_tracer_package)[ocean_tracer_mod(ocean _prog_tracer_init)]: Processing tracer package required &OCEAN_TEMPSALT_NML DEBUG_THIS_MODULE = F, TEMPERATURE_VARIABLE = potential_temp , POTTEMP_2ND_ITERATION = T, POTTEMP_EQUAL_CONTEMP = F, T_MIN = -5.00000000000000 , T_MAX = 55.0000000000000 , S_MIN = -1.00000000000000 , S_MAX = 55.0000000000000 , T_MIN_LIMIT = -2.00000000000000 , T_MAX_LIMIT = 32.0000000000000 , S_MIN_LIMIT = 0.000000000000000E+000, S_MAX_LIMIT = 42.0000000000000 , REINIT_TS_WITH_IDEAL = F, REINIT_TS_WITH_IDEAL_EFOLD = 1000.00000000000 , REINIT_TS_WITH_IDEAL_TVALUE = 10.0000000000000 , REINIT_TS_WITH_IDEAL_SVALUE = 30.0000000000000 , TEOS10 = F, DO_FAFMIP_HEAT = F / ==>Note from ocean_tempsalt_mod: MOM prognostic temp = potential temperature. MOM diagnostic temp = conservative temperature. ==>Note from ocean_tpm_util_mod(otpm_set_prog_tracer)[ocean_tempsalt_mod/ocean_ tempsalt_init]: Processing prog tracer temp ==>Note from ocean_tpm_util_mod(otpm_set_prog_tracer)[ocean_tempsalt_mod/ocean_ tempsalt_init]: Processing prog tracer salt ==>Note from ocean_tpm_util_mod(otpm_set_diag_tracer)[ocean_tempsalt_mod/ocean_ tempsalt_init]: Processing diag tracer con_temp ==>Note from ocean_tempsalt_mod: NOT Initializing FAFMIP temperature fields. &OCEAN_FRAZIL_NML USE_THIS_MODULE = T, DEBUG_THIS_MODULE = F, FREEZING_TEMP_SIMPLE = T, FREEZING_TEMP_PRETEOS10 = F, FREEZING_TEMP_TEOS10 = F, FRAZIL_FACTOR = 1.00000000000000 , AIR_SATURATED_WATER = T, FRAZIL_ONLY_IN_SURFACE = T / NOTE from PE 15: ==>Note from ocean_frazil_mod: USING frazil heating. ==>Note: Using simple equation for seawater freezing temperature. Assuming that frazil forms only in the surface(k=1) ocean grid cell. Setting gauge pressure to zero when computing seawater freezing temperature. ==>Note from ocean_tpm_util_mod(otpm_set_diag_tracer)[ocean_frazil_mod/ocean_fr azil_init]: Processing diag tracer frazil NOT Initializing frazil_redist heat diagnostic tracer for FAFMIP purposes. ==>Note from ocean_frazil_mod: using frazil_factor= 1.00 ==>Note from ocean_tpm_util_mod(otpm_set_tracer_package)[ocean_passive_mod(ocea n_passive_init)]: Processing tracer package ocean_passive ==>Note from ocean_passive_mod(ocean_passive_init): No instances of passive tracers in field_table. ==>Note: NOT running with idealized passive tracers. NOTE from PE 15: ==>Note: ocean_passive_mod: NOT using idealized passive tracer module. ==>Note from ocean_tpm_util_mod(otpm_set_tracer_package)[ocean_residency_mod(oc ean_residency_init)]: Processing tracer package ocean_residency ==>Note from ocean_residency_mod(ocean_residency_init): No instances ==>Note from ocean_tpm_util_mod(otpm_set_tracer_package)[ocean_age_tracer_mod(o cean_age_tracer_init)]: Processing tracer package ocean_age_tracer ==>Note from ocean_age_tracer_mod(ocean_age_tracer_init): 1 instances ==>Note from ocean_tpm_util_mod(otpm_set_prog_tracer)[ocean_age_tracer_mod(ocea n_age_tracer_init)]: Processing prog tracer age_global ==>Note from fm_util_mod(fm_util_start_namelist)[ocean_age_tracer_mod(ocean_age _tracer_init)]: Processing namelist ocean_age_tracer/global ==>Note from ocean_tpm_util_mod(otpm_set_tracer_package)[ocmip2_cfc_mod(ocmip2_ cfc_init)]: Processing tracer package ocmip2_cfc ==>Note from ocmip2_cfc_mod(ocmip2_cfc_init): No instances ==>Note from ocean_tpm_util_mod(otpm_set_tracer_package)[ocmip2_he_mod(ocmip2_h e_init)]: Processing tracer package ocmip2_he ==>Note from ocmip2_he_mod(ocmip2_he_init): No instances ==>Note from ocean_tpm_util_mod(otpm_set_tracer_package)[ocean_pert_co2_mod(oce an_pert_co2_init)]: Processing tracer package ocean_pert_co2 ==>Note from ocean_pert_co2_mod(ocean_pert_co2_init): No instances ==>Note from ocean_tpm_util_mod(otpm_set_tracer_package)[ocmip2_abiotic_mod(ocm ip2_abiotic_init)]: Processing tracer package ocmip2_abiotic ==>Note from ocmip2_abiotic_mod(ocmip2_abiotic_init): No instances ==>Note from ocean_tpm_util_mod(otpm_set_tracer_package)[ocmip2_biotic_mod(ocmi p2_biotic_init)]: Processing tracer package ocmip2_biotic ==>Note from ocmip2_biotic_mod(ocmip2_biotic_init): No instances ==>Note from ocean_tpm_util_mod(otpm_set_tracer_package)[ocean_bgc_restore_mod( ocean_bgc_restore_init)]: Processing tracer package ocean_bgc_restore ==>Note from ocean_bgc_restore_mod(ocean_bgc_restore_init): No instances ==>Note from ocean_tpm_util_mod(otpm_set_tracer_package)[ocean_po4_pre_mod(ocea n_po4_pre_init)]: Processing tracer package ocean_po4_pre ==>Note from ocean_po4_pre_mod(ocean_po4_pre_init): No instances ==>Note from ocean_tpm_util_mod(otpm_set_tracer_package)[ocean_ibgc_mod(ocean_i bgc_init)]: Processing tracer package ocean_ibgc ==>Note from ocean_ibgc_mod(ocean_ibgc_init): No instances &GENERIC_TRACER_NML DO_GENERIC_TRACER = T, DO_GENERIC_CFC = T, DO_GENERIC_TOPAZ = F, DO_GENERIC_ERGOM = F, DO_GENERIC_BLING = F, DO_GENERIC_MINIBLING = F, DO_GENERIC_COBALT = F / ==>Note from ocean_tpm_util_mod(otpm_set_tracer_package): Processing tracer package generic_cfc ==>Note from ocean_tpm_util_mod(otpm_set_prog_tracer): Processing prog tracer cfc_11 ==>Note from ocean_tpm_util_mod(otpm_set_prog_tracer): Processing prog tracer cfc_12 ==>Note from ocean_tpm_util_mod(otpm_set_tracer_package)[transport_matrix_mod(t ransport_matrix_init)]: Processing tracer package transport_matrix ==>Note from transport_matrix_mod(transport_matrix_init): No instances Dumping tracer_packages tracer tree tracer_packages/ ocean_age_tracer/ names = 'global' horizontal-advection-scheme = 'mdfl_sweby' vertical-advection-scheme = 'mdfl_sweby' restart_file = 'ocean_age.res.nc' min_tracer_limit = 0.000000000000000E+000 units = 'yr' min_range = 0.000000000000000E+000 max_range = 1.000000000000000E+020 flux_units = 'm' max_tracer_limit = 1.000000000000000E+020 generic_cfc/ names = '_' horizontal-advection-scheme = 'mdfl_sweby' vertical-advection-scheme = 'mdfl_sweby' min_tracer_limit = 0.000000000000000E+000 ppm_hlimiter = 2 ppm_vlimiter = 2 const_init_tracer = T restart_file = 'ocmip2_cfc.res.nc' required/ names = NULL ocean_passive/ names = NULL units = 'dimensionless' conversion = 1.00000000000000 offset = 0.000000000000000E+000 min_tracer = -1000000.00000000 max_tracer = 1000000.00000000 min_range = -10.0000000000000 max_range = 100.000000000000 restart_file = 'ocean_passive.res.nc' psom_limit = T flux_units = 'dimensionless' min_flux_range = -1.000000000000000E+016 max_flux_range = 1.000000000000000E+016 min_tracer_limit = -0.100000000000000 max_tracer_limit = 1.00000000000000 vertical-advection-scheme = 'mdppm' horizontal-advection-scheme = 'mdppm' ocean_residency/ names = NULL units = 'yr' restart_file = 'ocean_residency.res.nc' flux_units = 'm' min_tracer_limit = 0.000000000000000E+000 max_tracer_limit = 1.000000000000000E+020 ocmip2_cfc/ names = NULL restart_file = 'ocmip2_cfc.res.nc' ocmip2_he/ names = NULL restart_file = 'ocmip2_he.res.nc' ocean_pert_co2/ names = NULL restart_file = 'ocean_pert_co2.res.nc' ocmip2_abiotic/ names = NULL restart_file = 'ocmip2_abiotic.res.nc' ocmip2_biotic/ names = NULL restart_file = 'ocmip2_biotic.res.nc' ocean_bgc_restore/ names = NULL restart_file = 'ocean_bgc_restore.res.nc' ocean_po4_pre/ names = NULL restart_file = 'ocean_po4_pre.res.nc' ocean_ibgc/ names = NULL restart_file = 'ocean_ibgc.res.nc' transport_matrix/ names = NULL restart_file = 'transport_matrix.res.nc' Dumping prog_tracers tracer tree prog_tracers/ temp/ horizontal-advection-scheme = 'mdfl_sweby' vertical-advection-scheme = 'mdfl_sweby' restart_file = 'ocean_temp_salt.res.nc' longname = 'Potential temperature' units = 'deg_C' type = '' conversion = 3992.10322329649 offset = 273.150000000000 min_tracer = -5.00000000000000 max_tracer = 55.0000000000000 min_range = -10.0000000000000 max_range = 500.000000000000 use_only_advection = F const_init_tracer = F const_init_value = 0.000000000000000E+000 psom_limit = F ppm_hlimiter = 2 ppm_vlimiter = 2 mdt_scheme = 1 flux_units = 'Watts/m^2' min_flux_range = -1.000000000000000E+016 max_flux_range = 1.000000000000000E+016 min_tracer_limit = -2.00000000000000 max_tracer_limit = 32.0000000000000 salt/ horizontal-advection-scheme = 'mdfl_sweby' vertical-advection-scheme = 'mdfl_sweby' restart_file = 'ocean_temp_salt.res.nc' longname = 'Practical Salinity' units = 'psu' type = '' conversion = 1.000000000000000E-003 offset = 0.000000000000000E+000 min_tracer = -1.00000000000000 max_tracer = 55.0000000000000 min_range = -10.0000000000000 max_range = 100.000000000000 use_only_advection = F const_init_tracer = F const_init_value = 0.000000000000000E+000 psom_limit = F ppm_hlimiter = 2 ppm_vlimiter = 2 mdt_scheme = 1 flux_units = 'kg/(sec*m^2)' min_flux_range = -100000.000000000 max_flux_range = 100000.000000000 min_tracer_limit = 0.000000000000000E+000 max_tracer_limit = 42.0000000000000 age_global/ longname = 'Age (global)' units = 'yr' type = '' conversion = 1.00000000000000 offset = 0.000000000000000E+000 min_tracer = -1.000000000000000E+020 max_tracer = 1.000000000000000E+020 min_range = 0.000000000000000E+000 max_range = 1.000000000000000E+020 use_only_advection = F restart_file = 'ocean_age.res.nc' const_init_tracer = F const_init_value = 0.000000000000000E+000 psom_limit = F ppm_hlimiter = 2 ppm_vlimiter = 2 mdt_scheme = 1 flux_units = 'm' min_flux_range = 1.00000000000000 max_flux_range = 0.000000000000000E+000 min_tracer_limit = 0.000000000000000E+000 max_tracer_limit = 1.000000000000000E+020 vertical-advection-scheme = 'mdfl_sweby' horizontal-advection-scheme = 'mdfl_sweby' cfc_11/ longname = 'cfc_11 Concentration' units = 'mol/kg' type = 'generic' conversion = 1.00000000000000 offset = 0.000000000000000E+000 min_tracer = -1.000000000000000E+020 max_tracer = 1.000000000000000E+020 min_range = 1.00000000000000 max_range = 0.000000000000000E+000 use_only_advection = F restart_file = 'ocmip2_cfc.res.nc' const_init_tracer = T const_init_value = 0.000000000000000E+000 psom_limit = F ppm_hlimiter = 2 ppm_vlimiter = 2 mdt_scheme = 1 flux_units = '' min_flux_range = 1.00000000000000 max_flux_range = 0.000000000000000E+000 min_tracer_limit = 0.000000000000000E+000 max_tracer_limit = 1.000000000000000E+020 vertical-advection-scheme = 'mdfl_sweby' horizontal-advection-scheme = 'mdfl_sweby' cfc_12/ longname = 'cfc_12 Concentration' units = 'mol/kg' type = 'generic' conversion = 1.00000000000000 offset = 0.000000000000000E+000 min_tracer = -1.000000000000000E+020 max_tracer = 1.000000000000000E+020 min_range = 1.00000000000000 max_range = 0.000000000000000E+000 use_only_advection = F restart_file = 'ocmip2_cfc.res.nc' const_init_tracer = T const_init_value = 0.000000000000000E+000 psom_limit = F ppm_hlimiter = 2 ppm_vlimiter = 2 mdt_scheme = 1 flux_units = '' min_flux_range = 1.00000000000000 max_flux_range = 0.000000000000000E+000 min_tracer_limit = 0.000000000000000E+000 max_tracer_limit = 1.000000000000000E+020 vertical-advection-scheme = 'mdfl_sweby' horizontal-advection-scheme = 'mdfl_sweby' Dumping namelists tracer tree namelists/ ocean_age_tracer/ global/ slat = -90.0000000000000 nlat = 90.0000000000000 wlon = 0.000000000000000E+000 elon = 360.000000000000 coastal_only = F t_mask[1] = T t_mask[2] = T t_mask[3] = T t_mask[4] = T t_mask[5] = T t_mask[6] = T t_mask[7] = T t_mask[8] = T t_mask[9] = T t_mask[10] = T t_mask[11] = T t_mask[12] = T age_tracer_type = 'not used' generic_cfc/ ==> Note from ocean_tracer_mod: prognostic temperature = potential temperature. diagnostic temperature = conservative temperature. NOTE from PE 15: g_tracer_set_pointer_4D: Deallocating generic tracer cfc_11 % field NOTE from PE 15: g_tracer_set_pointer_3D: Deallocating generic tracer cfc_11 % tendency NOTE from PE 15: g_tracer_set_pointer_4D: Deallocating generic tracer cfc_12 % field NOTE from PE 15: g_tracer_set_pointer_3D: Deallocating generic tracer cfc_12 % tendency ==>Note from ocean_tracer_mod(ocean_prog_tracer_init): Reading prognostic tracer initial conditions or restarts Initializing tracer number 1 at time level tau. This tracer is called temp Reading restart for prog tracer temp from file ocean_temp_salt.res.nc NOTE from PE 0: interpolator_mod: Using noleap model calendar and julian file calendar for file aerosol.climatology.nc; calendar conversion needed NOTE from PE 0: interpolator_mod :aerosol.climatology.nc is a timeseries file for file aerosol.climatology.nc, the first time slice is mapped to : 1859 Jan 16 00:00:00 for file aerosol.climatology.nc, the last time slice is mapped to: 2010 Dec 16 00:00:00 Initializing src field : so4_anthro Initializing src field : so4_natural Initializing src field : organic_carbon Initializing src field : black_carbon Initializing src field : sea_salt Initializing src field : anthro_dust_0.1 Initializing src field : anthro_dust_0.2 Initializing src field : anthro_dust_0.4 Initializing src field : anthro_dust_0.8 Initializing src field : anthro_dust_1.0 Initializing src field : anthro_dust_2.0 Initializing src field : anthro_dust_4.0 Initializing src field : anthro_dust_8.0 Initializing src field : natural_dust_0.1 Initializing src field : natural_dust_0.2 Initializing src field : natural_dust_0.4 Initializing src field : natural_dust_0.8 Initializing src field : natural_dust_1.0 Initializing src field : natural_dust_2.0 Initializing src field : natural_dust_4.0 Initializing src field : natural_dust_8.0 NOTE from PE 0: radiative_gases_mod: Reading NetCDF formatted restart file: INPUT/radiative_gases.res.nc NOTE from PE 0: radiative_gases_mod: PROCESSING TIMESERIES FOR ch4 Gas value is taken from timeseries at time: 1990 Jan 01 00:00:00 ch4 value is 1.688625000000000E-006 NOTE from PE 0: radiative_gases_mod: PROCESSING TIMESERIES FOR n2o Gas value is taken from timeseries at time: 1990 Jan 01 00:00:00 n2o value is 3.084500000000000E-007 NOTE from PE 0: radiative_gases_mod: PROCESSING TIMESERIES FOR f11 Gas value is taken from timeseries at time: 1990 Jan 01 00:00:00 f11 value is 2.590000000000000E-010 NOTE from PE 0: radiative_gases_mod: PROCESSING TIMESERIES FOR f12 Gas value is taken from timeseries at time: 1990 Jan 01 00:00:00 f12 value is 4.663750000000001E-010 NOTE from PE 0: radiative_gases_mod: PROCESSING TIMESERIES FOR f113 Gas value is taken from timeseries at time: 1990 Jan 01 00:00:00 f113 value is 7.137500000000000E-011 NOTE from PE 0: radiative_gases_mod: PROCESSING TIMESERIES FOR f22 Gas value is taken from timeseries at time: 1990 Jan 01 00:00:00 f22 value is 8.925000000000000E-011 NOTE from PE 0: radiative_gases_mod: PROCESSING TIMESERIES FOR co2 Gas value is taken from timeseries at time: 1990 Jan 01 00:00:00 co2 value is 3.527219000000000E-004 NOTE from PE 0: ozone_mod: Ozone data is defined from a single annual cycle - no interannual variation Ozone data obtained from ozone timeseries for year: 1990 NOTE from PE 0: GLOBAL ATT too long - not reading this metadata After reading ic, linearly interpolate temp to partial cell bottom. Completed initialization of tracer temp at time level tau Initializing tracer number 2 at time level tau. This tracer is called salt Reading restart for prog tracer salt from file ocean_temp_salt.res.nc After reading ic, linearly interpolate salt to partial cell bottom. Completed initialization of tracer salt at time level tau Initializing tracer number 3 at time level tau. This tracer is called age_global Reading restart for prog tracer age_global from file ocean_age.res.nc After reading ic, linearly interpolate age_global to partial cell bottom. Completed initialization of tracer age_global at time level tau Initializing tracer number 4 at time level tau. This tracer is called cfc_11 Initializing the tracer cfc_11 to the constant 0.000000000000000E+000 Completed initialization of tracer cfc_11 at time level tau Initializing tracer number 5 at time level tau. This tracer is called cfc_12 Initializing the tracer cfc_12 to the constant 0.000000000000000E+000 Completed initialization of tracer cfc_12 at time level tau ==>Note from ocean_tracer_mod(ocean_prog_tracer_init): finished reading prognostic tracer restarts. ==>Note from ocean_tpm_util_mod(otpm_set_diag_tracer)[ocean_shortwave_mod(ocean _irradiance_init)]: Processing diag tracer irr ==>Note from ocean_tracer_mod(ocean_diag_tracer_init): 3 diagnostic tracers requested. Dumping ocean diag field tree after reading diag tracer tree diag_tracers/ frazil/ restart_file = 'ocean_frazil.res.nc' longname = 'frazil heating' units = 'J/m^2' type = '' conversion = 1.00000000000000 offset = 0.000000000000000E+000 min_tracer = 0.000000000000000E+000 max_tracer = 1.000000000000000E+020 min_range = -10.0000000000000 max_range = 100.000000000000 const_init_tracer = T const_init_value = 0.000000000000000E+000 con_temp/ longname = 'Conservative temperature' units = 'deg_C' type = '' conversion = 3992.10322329649 offset = 273.150000000000 min_tracer = -5.00000000000000 max_tracer = 55.0000000000000 min_range = -10.0000000000000 max_range = 500.000000000000 restart_file = 'ocean_con_temp.res.nc' const_init_tracer = T const_init_value = 0.000000000000000E+000 irr/ longname = 'Irradiance' units = 'Watts/m^2' type = '' conversion = 1.00000000000000 offset = 0.000000000000000E+000 min_tracer = -1.000000000000000E+020 max_tracer = 1.000000000000000E+020 min_range = 1.00000000000000 max_range = 0.000000000000000E+000 restart_file = '' const_init_tracer = T const_init_value = 0.000000000000000E+000 ==>Note from ocean_tracer_mod(ocean_diag_tracer_init): Reading diagnostic tracer initial conditions and/or restarts Initializing tracer number 1 at time level tau. This tracer is called frazil Initializing diagnostic tracer frazil to constant 0.000000000000000E+000 Initializing tracer number 2 at time level tau. This tracer is called con_temp Initializing diagnostic tracer con_temp to constant 0.000000000000000E+000 Initializing tracer number 3 at time level tau. This tracer is called irr Skipping tracer irr ==>Note from ocean_tracer_mod(ocean_diag_tracer_init): Finished reading diagnostic tracer restarts. &OCEAN_ADVECTION_VELOCITY_NML DEBUG_THIS_MODULE = F, MAX_ADVECTION_VELOCITY = 0.500000000000000 , INFLOW_NBOUNDARY = F, READ_ADVECTION_VELOCITY = F, READ_ADVECTION_TRANSPORT = F, CONSTANT_ADVECTION_VELOCITY = F / Note: Advection stability most nearly violated at T-cell (i,j) = ( 84, 5), (lon,lat) = (-196.500, -77.500). Note: Advection stability most nearly violated at T-cell (i,j) = ( 121, 5), (lon,lat) = (-159.500, -77.500). Assuming a maximum advection velocity of 0.50 m/s, Assuming a maximum advection velocity of 0.50 m/s, Linear stability requires max(dtuv,dtts) be less than 12033.00 sec. Model is now using (dtuv,dtts) = ( 7200.00, 7200.00) sec. Linear stability requires max(dtuv,dtts) be less than 12033.00 sec. Model is now using (dtuv,dtts) = ( 7200.00, 7200.00) sec. &OCEAN_DENSITY_NML S_TEST = 20.0000000000000 , T_TEST = 20.0000000000000 , P_TEST = 1000.00000000000 , PRESS_STANDARD = 0.000000000000000E+000, SN_TEST = 35.0000000000000 , TN_TEST = 20.0000000000000 , EOS_LINEAR = F, ALPHA_LINEAR_EOS = 0.255000000000000 , BETA_LINEAR_EOS = 0.000000000000000E+000, EOS_PRETEOS10 = T, EOS_TEOS10 = F, POTRHO_PRESS = 2000.00000000000 , POTRHO_MIN = 1028.00000000000 , POTRHO_MAX = 1038.00000000000 , NEUTRALRHO_MIN = 1020.00000000000 , NEUTRALRHO_MAX = 1030.00000000000 , LAYER_NK = 80, THETA_MIN = -2.00000000000000 , THETA_MAX = 30.0000000000000 , DEBUG_THIS_MODULE = F, WRITE_A_RESTART = T, RHO0_DENSITY = F, DENSITY_EQUAL_POTRHO = F, BUOYFREQ_SMOOTH_VERT = T, NUM_121_PASSES = 1, EPSLN_DRHODZ = 1.000000000000000E-010, MASK_DOMAIN_RESTART = F, DO_BITWISE_EXACT_SUM = F, DRHODZ_DIAG_STABLE = T, EPSLN_DRHODZ_DIAG = 1.000000000000000E-010, GRAD_NRHO_LRPOTRHO_COMPUTE = F, GRAD_NRHO_LRPOTRHO_MAX = 10.0000000000000 , GRAD_NRHO_LRPOTRHO_MIN = 1.00000000000000 , NEUTRAL_DENSITY_OMEGA = F, NEUTRAL_DENSITY_POTRHO = T, NEUTRAL_DENSITY_THETA = F, SMOOTH_STRATIFICATION_FACTOR = F, UPDATE_DIAGNOSTIC_FACTORS = F, SMAX_DIAG = -1.00000000000000 , SMAX_MIN_IN_COLUMN = F / ==> Note: USING preTEOS10 EOS, as relevant for realistic ocean climate simulations. Subtracting standard atmosphere of 0.000000 dbar for EOS calculation. ==> Note: Computing EOS assuming prognostic temp = potential temperature. ==> Note: Computing EOS assuming prognostic salinity = practical salinity. ==>Note: The Boussinesq rho0 density has a value of (kg/m^3) 1035.00000000000 ==> Note: Computing diagnostic neutral_rho as potential density referenced to pressure potrho_press. ==> Note: Enforcing drhodz_diag < 0, so to use stable stratification for certain diagnostic purposes. ==> Note: Diagnostic factors are NOT computed. So if enable watermass diagnostics, they will be corrupted. preTEOS10 EQUATION OF STATE TEST VALUES s_test(psu) = 20.00, t_test(C) = 20.00, p_test(dbar) = 1000.00 rho ( 20.00, 20.00, 1000.00) = 0.1017728868019642E+04 kg/m^3 diff from JMFWG = 0.0000000000000000E+00 kg/m^3 alpha( 20.00, 20.00, 1000.00) = 0.2525481286927142E-03 1/C diff from JMFWG = 0.8673617379884035E-18 1/C beta ( 20.00, 20.00, 1000.00) = 0.7379638527217574E-03 1/psu squared sound speed ( 20.00, 20.00, 1000.00) = 0.2316107367172913E+03 (m/s)^2 diff from JMFWG = -.1084202172485504E-18 1/psu NEUTRAL DENSITY EQUATION TEST VALUES sn_test(psu) = 35.00, tn_test(C) = 20.00 rho ( 35.00, 20.00) = 0.1024594167511967E+04 kg/m^3 diff from Klocker and McDougall test = -.2728484105318785E-11 kg/m^3 Reading restart for density from INPUT/ocean_density.res.nc NOTE from PE 0: interpolator_mod: Using noleap model calendar and julian file calendar for file o3.climatology.nc; calendar conversion needed NOTE from PE 0: interpolator_mod :o3.climatology.nc is a timeseries file for file o3.climatology.nc, the first time slice is mapped to : 1859 Jan 16 00:00:00 for file o3.climatology.nc, the last time slice is mapped to: 2004 Dec 16 00:00:00 Initializing src field : ozone ==>Note: ocean_density_mod: did not read density derivatives from restart. Initialising salinity for use in density calculation NOTE from PE 0: radiation_driver_mod: Solar data is fixed in time Data used in this experiment is from solar timeseries at time: 1990 Jan 01 00:00:00 NOTE from PE 0: radiation_driver_mod: radiation to be calculated on first step: user asserts that this is a scheduled radiation step; if it is not, restart seamlessness will be lost ==>Note: ocean_density_mod: did not read drhodz_zt from restart. ==>Note: From ocean_density_mod: Boussinesq reference density rho0(kg/m3) = 0.103500000000E+04 Initial rho_average(kg/m3) = 0.103500000275E+04 Since rho0 .ne. rho_average, consider changing rho0 in ocean_parameters.F90 to be equal to rho_average for better accuracy. From ocean_density_mod: density chksums from ocean_density_init yyyy/mm/dd hh:mm:ss = 1/ 1/ 1 0: 0: 0 [chksum] rho(taup1) 9096130062845791024 [chksum] pressure_at_depth -6334309842238097397 [chksum] denominator_r -1533836596824791085 [chksum] drhodT 0 [chksum] drhodS 0 [chksum] drhodz_zt 2778077585614187437 &OCEAN_BLOB_NML BLOB_SMALL_MASS = 1000.00000000000 , DEBUG_THIS_MODULE = F, REALLY_DEBUG = F, DO_BITWISE_EXACT_SUM = F, BITWISE_REPRODUCTION = F, MAX_PROP_THICKNESS = 0.700000000000000 / ==>Note: NOT using the Lagrangian buoyancy blobs scheme. &OCEAN_PRESSURE_NML DEBUG_THIS_MODULE = F, ZERO_PRESSURE_FORCE = F, ZERO_CORRECTION_TERM_GRAD = F, ZERO_DIAGONAL_PRESS_GRAD = F, ZERO_ETA_OVER_H_ZSTAR_PRESSURE = F / ==>NOTE: Running MOM with finite difference formulation of pressure force. &OCEAN_VERT_MIX_NML DEBUG_THIS_MODULE = F, VERT_MIX_SCHEME = kpp_mom4p0, VERBOSE_INIT = T, AIDIF = 1.00000000000000 , VERT_DIFF_BACK_VIA_MAX = T, USE_EXPLICIT_VERT_DIFFUSE = T, USE_DIFF_CBT_TABLE = F, LINEAR_TAPER_DIFF_CBT_TABLE = T, BRYAN_LEWIS_DIFFUSIVITY = T, BRYAN_LEWIS_LAT_DEPEND = T, BRYAN_LEWIS_LAT_TRANSITION = 35.0000000000000 , AFKPH_90 = 0.750000000000000 , DFKPH_90 = 0.950000000000000 , SFKPH_90 = 4.500000000000000E-005, ZFKPH_90 = 250000.000000000 , AFKPH_00 = 0.650000000000000 , DFKPH_00 = 1.15000000000000 , SFKPH_00 = 4.500000000000000E-005, ZFKPH_00 = 250000.000000000 , VERT_VISC_BACK = F, VISC_CBU_BACK_MAX = 1.000000000000000E-002, VISC_CBU_BACK_MIN = 1.000000000000000E-003, VISC_CBU_BACK_ZMID = 50.0000000000000 , VISC_CBU_BACK_ZWID = 30.0000000000000 , HWF_DIFFUSIVITY = F, HWF_DEPTH_TRANSITION = 25000000.0000000 , HWF_MIN_DIFFUSIVITY = 2.000000000000000E-006, HWF_30_DIFFUSIVITY = 2.000000000000000E-005, HWF_N0_2OMEGA = 20.0000000000000 , HWF_DIFFUSIVITY_3D = F, DIFF_CBT_TANH = F, DIFF_CBT_TANH_MAX = 1.000000000000000E-003, DIFF_CBT_TANH_MIN = 2.000000000000000E-005, DIFF_CBT_TANH_ZMID = 150.000000000000 , DIFF_CBT_TANH_ZWID = 30.0000000000000 , J09_DIFFUSIVITY = F, J09_BGMIN = 1.000000000000000E-006, J09_BGMAX = 1.000000000000000E-005, J09_LAT = 20.0000000000000 , QUEBEC_2009_10_BUG = F, VMIX_RESCALE_NONBOUSS = F, VMIX_SET_MIN_DISSIPATION = F, VMIX_MIN_DISS_CONST = 1.000000000000000E-007, VMIX_MIN_DISS_BVFREQ_SCALE = 5.999999999999999E-004, VMIX_MIN_DISS_FLUX_RI_MAX = 0.200000000000000 , SMOOTH_RHO_N2 = T, NUM_121_PASSES = 1, READ_DIFF_CBT_FILE = F / &OCEAN_VERT_UTIL_NML DEBUG_THIS_MODULE = F, SMOOTH_N2 = T, SMOOTH_RI_NUMBER = T, NUM_N2_SMOOTH = 1, NUM_RI_SMOOTH = 1 / NOTE from PE 15: ==>Warning: ocean_vert_mix_init found diff_cbt_table, yet use_diff_cbt_table=.false. =>Note: USING Bryan Lewis vertical background tracer diffusivity. If using diffusivity from vert_tidal_mod, then should turn off Bryan Le wis. =>Note: USING bryan_lewis_lat_depend. This will modify the background latitudi nally, with a transition latitude at bryan_lewis_lat_transition = 35.0000000000000 Bryan-Lewis diffusivities diff_bryan_lewis_00(m^2/sec)( 1) = 0.1075827E-04 diff_bryan_lewis_00(m^2/sec)( 2) = 0.1077134E-04 diff_bryan_lewis_00(m^2/sec)( 3) = 0.1078451E-04 diff_bryan_lewis_00(m^2/sec)( 4) = 0.1079779E-04 diff_bryan_lewis_00(m^2/sec)( 5) = 0.1081118E-04 diff_bryan_lewis_00(m^2/sec)( 6) = 0.1082467E-04 diff_bryan_lewis_00(m^2/sec)( 7) = 0.1083828E-04 diff_bryan_lewis_00(m^2/sec)( 8) = 0.1085200E-04 diff_bryan_lewis_00(m^2/sec)( 9) = 0.1086583E-04 diff_bryan_lewis_00(m^2/sec)( 10) = 0.1087977E-04 diff_bryan_lewis_00(m^2/sec)( 11) = 0.1089383E-04 diff_bryan_lewis_00(m^2/sec)( 12) = 0.1090801E-04 diff_bryan_lewis_00(m^2/sec)( 13) = 0.1092231E-04 diff_bryan_lewis_00(m^2/sec)( 14) = 0.1093673E-04 diff_bryan_lewis_00(m^2/sec)( 15) = 0.1095126E-04 diff_bryan_lewis_00(m^2/sec)( 16) = 0.1096592E-04 diff_bryan_lewis_00(m^2/sec)( 17) = 0.1098071E-04 diff_bryan_lewis_00(m^2/sec)( 18) = 0.1099562E-04 diff_bryan_lewis_00(m^2/sec)( 19) = 0.1101066E-04 diff_bryan_lewis_00(m^2/sec)( 20) = 0.1102583E-04 diff_bryan_lewis_00(m^2/sec)( 21) = 0.1104113E-04 diff_bryan_lewis_00(m^2/sec)( 22) = 0.1105657E-04 diff_bryan_lewis_00(m^2/sec)( 23) = 0.1107301E-04 diff_bryan_lewis_00(m^2/sec)( 24) = 0.1109315E-04 diff_bryan_lewis_00(m^2/sec)( 25) = 0.1112065E-04 diff_bryan_lewis_00(m^2/sec)( 26) = 0.1115937E-04 diff_bryan_lewis_00(m^2/sec)( 27) = 0.1121353E-04 diff_bryan_lewis_00(m^2/sec)( 28) = 0.1128791E-04 diff_bryan_lewis_00(m^2/sec)( 29) = 0.1138826E-04 diff_bryan_lewis_00(m^2/sec)( 30) = 0.1152186E-04 diff_bryan_lewis_00(m^2/sec)( 31) = 0.1169839E-04 diff_bryan_lewis_00(m^2/sec)( 32) = 0.1193134E-04 diff_bryan_lewis_00(m^2/sec)( 33) = 0.1224049E-04 diff_bryan_lewis_00(m^2/sec)( 34) = 0.1265628E-04 diff_bryan_lewis_00(m^2/sec)( 35) = 0.1322832E-04 diff_bryan_lewis_00(m^2/sec)( 36) = 0.1404353E-04 diff_bryan_lewis_00(m^2/sec)( 37) = 0.1526883E-04 diff_bryan_lewis_00(m^2/sec)( 38) = 0.1726850E-04 diff_bryan_lewis_00(m^2/sec)( 39) = 0.2100099E-04 diff_bryan_lewis_00(m^2/sec)( 40) = 0.2981891E-04 diff_bryan_lewis_00(m^2/sec)( 41) = 0.5795053E-04 diff_bryan_lewis_00(m^2/sec)( 42) = 0.9578484E-04 diff_bryan_lewis_00(m^2/sec)( 43) = 0.1085765E-03 diff_bryan_lewis_00(m^2/sec)( 44) = 0.1134187E-03 diff_bryan_lewis_00(m^2/sec)( 45) = 0.1158581E-03 diff_bryan_lewis_00(m^2/sec)( 46) = 0.1173069E-03 diff_bryan_lewis_00(m^2/sec)( 47) = 0.1182583E-03 diff_bryan_lewis_00(m^2/sec)( 48) = 0.1189264E-03 diff_bryan_lewis_00(m^2/sec)( 49) = 0.1194183E-03 diff_bryan_lewis_00(m^2/sec)( 50) = 0.1197934E-03 Bryan-Lewis diffusivities at pole diff_bryan_lewis_90(m^2/sec)( 1) = 0.3019161E-04 diff_bryan_lewis_90(m^2/sec)( 2) = 0.3020241E-04 diff_bryan_lewis_90(m^2/sec)( 3) = 0.3021329E-04 diff_bryan_lewis_90(m^2/sec)( 4) = 0.3022426E-04 diff_bryan_lewis_90(m^2/sec)( 5) = 0.3023532E-04 diff_bryan_lewis_90(m^2/sec)( 6) = 0.3024647E-04 diff_bryan_lewis_90(m^2/sec)( 7) = 0.3025771E-04 diff_bryan_lewis_90(m^2/sec)( 8) = 0.3026904E-04 diff_bryan_lewis_90(m^2/sec)( 9) = 0.3028047E-04 diff_bryan_lewis_90(m^2/sec)( 10) = 0.3029199E-04 diff_bryan_lewis_90(m^2/sec)( 11) = 0.3030360E-04 diff_bryan_lewis_90(m^2/sec)( 12) = 0.3031531E-04 diff_bryan_lewis_90(m^2/sec)( 13) = 0.3032712E-04 diff_bryan_lewis_90(m^2/sec)( 14) = 0.3033903E-04 diff_bryan_lewis_90(m^2/sec)( 15) = 0.3035104E-04 diff_bryan_lewis_90(m^2/sec)( 16) = 0.3036315E-04 diff_bryan_lewis_90(m^2/sec)( 17) = 0.3037537E-04 diff_bryan_lewis_90(m^2/sec)( 18) = 0.3038769E-04 diff_bryan_lewis_90(m^2/sec)( 19) = 0.3040011E-04 diff_bryan_lewis_90(m^2/sec)( 20) = 0.3041264E-04 diff_bryan_lewis_90(m^2/sec)( 21) = 0.3042528E-04 diff_bryan_lewis_90(m^2/sec)( 22) = 0.3043803E-04 diff_bryan_lewis_90(m^2/sec)( 23) = 0.3045162E-04 diff_bryan_lewis_90(m^2/sec)( 24) = 0.3046826E-04 diff_bryan_lewis_90(m^2/sec)( 25) = 0.3049097E-04 diff_bryan_lewis_90(m^2/sec)( 26) = 0.3052296E-04 diff_bryan_lewis_90(m^2/sec)( 27) = 0.3056769E-04 diff_bryan_lewis_90(m^2/sec)( 28) = 0.3062914E-04 diff_bryan_lewis_90(m^2/sec)( 29) = 0.3071204E-04 diff_bryan_lewis_90(m^2/sec)( 30) = 0.3082241E-04 diff_bryan_lewis_90(m^2/sec)( 31) = 0.3096823E-04 diff_bryan_lewis_90(m^2/sec)( 32) = 0.3116067E-04 diff_bryan_lewis_90(m^2/sec)( 33) = 0.3141606E-04 diff_bryan_lewis_90(m^2/sec)( 34) = 0.3175953E-04 diff_bryan_lewis_90(m^2/sec)( 35) = 0.3223209E-04 diff_bryan_lewis_90(m^2/sec)( 36) = 0.3290553E-04 diff_bryan_lewis_90(m^2/sec)( 37) = 0.3391773E-04 diff_bryan_lewis_90(m^2/sec)( 38) = 0.3556963E-04 diff_bryan_lewis_90(m^2/sec)( 39) = 0.3865299E-04 diff_bryan_lewis_90(m^2/sec)( 40) = 0.4593736E-04 diff_bryan_lewis_90(m^2/sec)( 41) = 0.6917653E-04 diff_bryan_lewis_90(m^2/sec)( 42) = 0.1004310E-03 diff_bryan_lewis_90(m^2/sec)( 43) = 0.1109980E-03 diff_bryan_lewis_90(m^2/sec)( 44) = 0.1149980E-03 diff_bryan_lewis_90(m^2/sec)( 45) = 0.1170132E-03 diff_bryan_lewis_90(m^2/sec)( 46) = 0.1182100E-03 diff_bryan_lewis_90(m^2/sec)( 47) = 0.1189960E-03 diff_bryan_lewis_90(m^2/sec)( 48) = 0.1195479E-03 diff_bryan_lewis_90(m^2/sec)( 49) = 0.1199543E-03 diff_bryan_lewis_90(m^2/sec)( 50) = 0.1202641E-03 &OCEAN_VERT_TIDAL_NML USE_THIS_MODULE = F, USE_LEGACY_METHODS = F, DEBUG_THIS_MODULE = F, USE_WAVE_DISSIPATION = F, USE_DRAG_DISSIPATION = F, READ_ROUGHNESS = T, READ_TIDE_SPEED = T, DEFAULT_ROUGHNESS_LENGTH = 25.0000000000000 , DEFAULT_TIDE_SPEED = 1.000000000000000E-002, SHELF_DEPTH_CUTOFF = -1000.00000000000 , DECAY_SCALE = 500.000000000000 , ROUGHNESS_SCALE = 85000.0000000000 , TIDAL_DISS_EFFICIENCY = 0.333330000000000 , MIXING_EFFICIENCY = 0.200000000000000 , MIXING_EFFICIENCY_N2DEPEND = F, MUNK_ANDERSON_P = 0.250000000000000 , MUNK_ANDERSON_SIGMA = 3.00000000000000 , DRAG_DISSIPATION_EFOLD = T, DRAG_DISSIPATION_TIDE_PERIOD = 43200.0000000000 , DRAG_MASK_DEEP = T, DRAG_MASK_DEEP_RATIO = 0.100000000000000 , BOTTOM_DRAG_CD = 2.400000000000000E-003, DRHODZ_MIN = 1.000000000000000E-010, SPEED_MIN = 5.000000000000000E-003, BACKGROUND_DIFFUSIVITY = 1.000000000000000E-005, BACKGROUND_VISCOSITY = 1.000000000000000E-005, MAX_WAVE_DIFFUSIVITY = 5.000000000000000E-003, MAX_DRAG_DIFFUSIVITY = 5.000000000000000E-003, SMOOTH_BVFREQ_BOTTOM = T, VEL_MICOM_SMOOTH = 0.200000000000000 , SMOOTH_RHO_N2 = T, NUM_121_PASSES = 1, WAVE_DIFFUSIVITY_MONOTONIC = T, TIDE_SPEED_DATA_ON_T_GRID = T, READING_ROUGHNESS_AMP = F, READING_ROUGHNESS_LENGTH = F, READ_WAVE_DISSIPATION = F, FIXED_WAVE_DISSIPATION = F, WAVE_ENERGY_FLUX_MAX = 0.100000000000000 , USE_LEEWAVE_DISSIPATION = F, READ_LEEWAVE_DISSIPATION = F, DRAG_DISSIPATION_USE_CDBOT = F / NOTE from PE 15: ==>Note: NOT using ocean_vert_tidal_mod ==>Note from ocean_vert_mix: KPP_mom4p0 for vert diffusivity, viscosity, nonlocal, and barotropic tide drag. &OCEAN_VERT_KPP_MOM4P0_NML USE_THIS_MODULE = T, SHEAR_INSTABILITY = T, DOUBLE_DIFFUSION = T, DIFF_CBT_IW = 1.340000000000000E-007, VISC_CBU_IW = 1.000000000000000E-004, VISC_CBU_LIMIT = 5.000000000000000E-003, DIFF_CBT_LIMIT = 5.000000000000000E-003, VISC_CON_LIMIT = 0.100000000000000 , DIFF_CON_LIMIT = 0.100000000000000 , CONCV = 1.80000000000000 , RICR = 0.300000000000000 , NON_LOCAL_KPP = T, SMOOTH_BLMC = T, LGAM = 1.04000000000000 , CW_0 = 0.150000000000000 , L_SMYTH = 2.00000000000000 , LTMAX = 5.00000000000000 , WSTFAC = 0.600000000000000 , COASTAL_TIDAL_MIX = T, P_TIDE = -0.250000000000000 , SIGMA_TIDE = 3.00000000000000 , INT_TIDAL_MIX = F, INT_TIDE_ZETA1 = 300.000000000000 , INT_TIDE_ZETA2 = 1800.00000000000 , INT_TIDE_MIN_DEPTH = 100.000000000000 , INT_TIDE_Q = 0.333330000000000 , INT_TIDE_GAMMA = 0.200000000000000 , WSFC_COMBINE_RUNOFF_CALVE = T, DO_LANGMUIR = F / ==> NOTE: USING KPP_mom4p0 vertical mixing scheme. This scheme is hard-wired for GEOPOTENTIAL coordinates. It is not generally recommended for use, other than for legacy. ==> NOTE: KPP is typically run with penetrative shortwave heating. ==> NOTE: KPP is typically run with a seasonal and/or diurnal cycle. ==>From ocean_vert_kpp_mom4p0_mod: time step for vert-frict of (secs) 7200.00 ==>From ocean_vert_kpp_mom4p0_mod: time step for vert-diff of (secs) 7200.00 Computing vertical mixing from shear instability in KPP module. Computing vertical mixing from barotropic tides within KPP_mom4p0 module. Completed read of tidal velocity amplitude &OCEAN_BIH_TRACER_NML USE_THIS_MODULE = F, ABIH = 0.000000000000000E+000, TRACER_MIX_MICOM = F, VEL_MICOM = 0.000000000000000E+000, READ_DIFFUSIVITY_MASK = F, HORZ_Z_DIFFUSE = F, HORZ_S_DIFFUSE = T / ==>Note from ocean_bih_tracer_mod: NOT using this module. &OCEAN_LAP_TRACER_NML USE_THIS_MODULE = F, ALAP = 0.000000000000000E+000, TRACER_MIX_MICOM = F, VEL_MICOM = 0.000000000000000E+000, VERBOSE_INIT = T, READ_DIFFUSIVITY_MASK = F, HORZ_Z_DIFFUSE = F, HORZ_S_DIFFUSE = T / ==>Note from ocean_lap_tracer_mod: NOT using this module. &OCEAN_SIGMA_TRANSPORT_NML USE_THIS_MODULE = T, DEBUG_THIS_MODULE = F, TMASK_SIGMA_ON = F, SIGMA_DIFFUSION_ON = T, SIGMA_ADVECTION_ON = F, SIGMA_ADVECTION_SGS_ONLY = F, SIGMA_ADVECTION_CHECK = T, THICKNESS_SIGMA_LAYER = 100.000000000000 , THICKNESS_SIGMA_MAX = 100.000000000000 , THICKNESS_SIGMA_MIN = 100.000000000000 , SIGMA_DIFFUSIVITY = 1000.00000000000 , SIGMA_DIFFUSIVITY_RATIO = 1.000000000000000E-006, TRACER_MIX_MICOM = T, VEL_MICOM = 0.500000000000000 , VERBOSE_INIT = T, SIGMA_JUST_IN_BOTTOM_CELL = T, CAMPINGOOSE_MU = 1.000000000000000E-004, CAMPINGOOSE_DELTA = 0.333300000000000 , SIGMA_UMAX = 1.000000000000000E-002, WRITE_A_RESTART = T, SMOOTH_SIGMA_THICKNESS = T, SMOOTH_SIGMA_VELOCITY = T, SMOOTH_VELMICOM = 0.200000000000000 / NOTE from PE 15: ==>Note from ocean_sigma_transport_mod: USING ocean_sigma_transport_mod. ==>Note: ocean_sigma_transport_mod: using forward time step of (secs) 7200.00 ==>Note: ocean_sigma_transport_mod: sigma_diffusion_on=.true. ==>Note: ocean_sigma_transport_mod: sigma_advection_on=.false. ==>Note from ocean_sigma_transport_mod: using sigma_just_in_bottom=.true., as in mom4p0. flux dom sigma domain decomposition whalo = 1, ehalo = 1, shalo = 1, nhalo = 1 X-AXIS = 60 60 60 60 60 60 Y-AXIS = 40 40 40 40 40 NOTE from PE 15: mpp_domains_define.inc(mpp_define_io_domain): io domain will not be defined for flux dom sigma when one or both entry of io_layout is not positive ==>Note: ocean_sigma_transport_mod: initial thickness of sigma layer (m) = 100.0000 ==>Note: ocean_sigma_transport_mod: smooth_sigma_thickness=.true. => diffuse sigma_thickness. Laplacian diffusivity in sigma layer at (isc, 1) = 0.0000000E+00 m^2/s Laplacian diffusivity in sigma layer at (isc, 2) = 0.0000000E+00 m^2/s Laplacian diffusivity in sigma layer at (isc, 3) = 0.0000000E+00 m^2/s Laplacian diffusivity in sigma layer at (isc, 4) = 0.0000000E+00 m^2/s Laplacian diffusivity in sigma layer at (isc, 5) = 0.0000000E+00 m^2/s Laplacian diffusivity in sigma layer at (isc, 6) = 0.0000000E+00 m^2/s Laplacian diffusivity in sigma layer at (isc, 7) = 0.0000000E+00 m^2/s Laplacian diffusivity in sigma layer at (isc, 8) = 0.0000000E+00 m^2/s Laplacian diffusivity in sigma layer at (isc, 9) = 0.0000000E+00 m^2/s Laplacian diffusivity in sigma layer at (isc, 10) = 0.0000000E+00 m^2/s Laplacian diffusivity in sigma layer at (isc, 11) = 0.0000000E+00 m^2/s Laplacian diffusivity in sigma layer at (isc, 12) = 0.0000000E+00 m^2/s Laplacian diffusivity in sigma layer at (isc, 13) = 0.0000000E+00 m^2/s Laplacian diffusivity in sigma layer at (isc, 14) = 0.0000000E+00 m^2/s Laplacian diffusivity in sigma layer at (isc, 15) = 0.3077527E+05 m^2/s Laplacian diffusivity in sigma layer at (isc, 16) = 0.3169895E+05 m^2/s Laplacian diffusivity in sigma layer at (isc, 17) = 0.3259490E+05 m^2/s Laplacian diffusivity in sigma layer at (isc, 18) = 0.3346402E+05 m^2/s Laplacian diffusivity in sigma layer at (isc, 19) = 0.3430716E+05 m^2/s Laplacian diffusivity in sigma layer at (isc, 20) = 0.3512513E+05 m^2/s Laplacian diffusivity in sigma layer at (isc, 21) = 0.3591871E+05 m^2/s Laplacian diffusivity in sigma layer at (isc, 22) = 0.3668865E+05 m^2/s Laplacian diffusivity in sigma layer at (isc, 23) = 0.3743566E+05 m^2/s Laplacian diffusivity in sigma layer at (isc, 24) = 0.3816042E+05 m^2/s Laplacian diffusivity in sigma layer at (isc, 25) = 0.3886360E+05 m^2/s Laplacian diffusivity in sigma layer at (isc, 26) = 0.3954580E+05 m^2/s Laplacian diffusivity in sigma layer at (isc, 27) = 0.4020764E+05 m^2/s Laplacian diffusivity in sigma layer at (isc, 28) = 0.4084969E+05 m^2/s Laplacian diffusivity in sigma layer at (isc, 29) = 0.4147249E+05 m^2/s Laplacian diffusivity in sigma layer at (isc, 30) = 0.4207658E+05 m^2/s Laplacian diffusivity in sigma layer at (isc, 31) = 0.4266246E+05 m^2/s Laplacian diffusivity in sigma layer at (isc, 32) = 0.4323062E+05 m^2/s Laplacian diffusivity in sigma layer at (isc, 33) = 0.4378151E+05 m^2/s Laplacian diffusivity in sigma layer at (isc, 34) = 0.4431559E+05 m^2/s Laplacian diffusivity in sigma layer at (isc, 35) = 0.4483328E+05 m^2/s Laplacian diffusivity in sigma layer at (isc, 36) = 0.4533499E+05 m^2/s Laplacian diffusivity in sigma layer at (isc, 37) = 0.4582111E+05 m^2/s Laplacian diffusivity in sigma layer at (isc, 38) = 0.4629202E+05 m^2/s Laplacian diffusivity in sigma layer at (isc, 39) = 0.4674809E+05 m^2/s Laplacian diffusivity in sigma layer at (isc, 40) = 0.4718964E+05 m^2/s From ocean_sigma_transport_mod: initial thickness_sigma chksum yyyy/mm/dd hh:mm:ss = 1/ 1/ 1 0: 0: 0 [chksum] thickness_sigma 1591740993298759680 &OCEAN_NPHYSICS_NML USE_THIS_MODULE = T, DEBUG_THIS_MODULE = F, WRITE_A_RESTART = T, USE_NPHYSICSA = T, USE_NPHYSICSB = F, USE_NPHYSICSC = F / [Note] ocean_nphysics.F90: USING ocean_nphysics_mod. &OCEAN_NPHYSICS_UTIL_NML DEBUG_THIS_MODULE = F, NPHYSICS_UTIL_ZERO_INIT = F, SMAX = 2.000000000000000E-003, SWIDTH = 2.000000000000000E-003, EPSLN_DRHODZ = 1.000000000000000E-030, DRHODZ_MOM4P1 = F, DRHODZ_SMOOTH_HORZ = F, DRHODZ_SMOOTH_VERT = F, NUM_121_PASSES = 1, AREDI = 600.000000000000 , AGM = 100.000000000000 , AREDI_EQUAL_AGM = F, TRACER_MIX_MICOM = F, VEL_MICOM = 0.000000000000000E+000, BRYAN_LEWIS_AREDI = F, AHS = 0.000000000000000E+000, AHB = 0.000000000000000E+000, NEUTRAL_HORZ_MIX_BDY = F, VEL_MICOM_BDY = 0.000000000000000E+000, AH_BDY = 0.000000000000000E+000, AGM_LAT_BANDS = F, AGM_LAT_BANDS_BOUNDARY = -999.000000000000 , AGM_LAT_BANDS_RATIO = 1.00000000000000 , ROSSBY_RADIUS_MAX = 100000.000000000 , ROSSBY_RADIUS_MIN = 15000.0000000000 , AGM_READ_RESTART = F, AGM_CLOSURE = T, AGM_CLOSURE_SCALING = 7.000000000000001E-002, AGM_CLOSURE_MAX = 600.000000000000 , AGM_CLOSURE_MIN = 100.000000000000 , AGM_CLOSURE_GROWTH_SCALE = 0.500000000000000 , AGM_CLOSURE_LENGTH_FIXED = F, AGM_CLOSURE_LENGTH = 50000.0000000000 , AGM_CLOSURE_LENGTH_ROSSBY = F, AGM_CLOSURE_LENGTH_BCZONE = F, BCZONE_MAX_PTS = 10, AGM_CLOSURE_BCZONE_CRIT_RATE = 1.400000000000000E-006, AGM_CLOSURE_EDEN_GREATBATCH = F, AGM_CLOSURE_EDEN_GAMMA = 200.000000000000 , AGM_CLOSURE_EDEN_LENGTH_CONST = F, AGM_CLOSURE_EDEN_LENGTH = 10000.0000000000 , AGM_CLOSURE_EADY_SMOOTH_VERT = F, AGM_CLOSURE_EADY_SMOOTH_HORZ = F, AGM_CLOSURE_EADY_AVE_MIXED = F, AGM_CLOSURE_EADY_CAP = F, AGM_CLOSURE_BAROCLINIC = T, AGM_CLOSURE_BUOY_FREQ = 4.000000000000000E-003, AGM_CLOSURE_UPPER_DEPTH = 100.000000000000 , AGM_CLOSURE_LOWER_DEPTH = 2000.00000000000 , AGM_CLOSURE_LENGTH_CAP = F, AGM_CLOSURE_LENGTH_MAX = 50000.0000000000 , AGM_SMOOTH_SPACE = F, VEL_MICOM_SMOOTH = 0.200000000000000 , AGM_SMOOTH_TIME = F, AGM_DAMPING_TIME = 10.0000000000000 , AGM_CLOSURE_GRID_SCALING = F, AGM_CLOSURE_GRID_SCALING_POWER = 2.00000000000000 , AREDI_DIFFUSIVITY_GRID_SCALING = F, AGM_CLOSURE_N2_SCALE = F, AGM_CLOSURE_N2_SCALE_COEFF = 1000.00000000000 , AGM_CLOSURE_N2_SCALE_NREF_CST = F, SMAX_GRAD_GAMMA_SCALAR = 1.000000000000000E-002, EPSLN_DRHODZ_DIAGNOSTICS = 1.000000000000000E-007, WDIANEUTRAL_SMOOTH = T, SMOOTH_ETA_TEND_GM90 = F / bczone domain decomposition whalo = 10, ehalo = 10, shalo = 10, nhalo = 10 X-AXIS = 60 60 60 60 60 60 Y-AXIS = 40 40 40 40 40 NOTE from PE 15: mpp_domains_define.inc(mpp_define_io_domain): io domain will not be defined for bczone when one or both entry of io_layout is not positive NOTE from PE 0: atmos_tracer_utilities_init: Dry dep units set to none kg/(m2 s) in atmos_tracer_utilities for cld_amt [Note] ocean_nphysics_util.F90: Computing 2d flow-dependent tracer diffusivity with agm_closure_baroclinic. The maximum allowable diffusivity (m^2/s) is given by .60000E+03 The minimum allowable diffusivity (m^2/s) is given by .10000E+03 Depths (m) between which compute eady growth and baroclinicity = .10000E+03 .20000E+04 [Note] ocean_nphysics_util.F90: Length and time scales set by vertically averaged baroclinicity |grad(rho)|, as well as the constant buoyancy freq(sec^-1) = .40000E-02 and the constant length scale (m) = .50000E+05 [Note] ocean_nphysics.F90: USING ocean_nphysicsA. &OCEAN_NPHYSICSA_NML USE_THIS_MODULE = T, DEBUG_THIS_MODULE = F, NEUTRAL_PHYSICS_LIMIT = T, USE_GM_SKEW = T, DM_TAPER = T, GKW_TAPER = F, TMASK_NEUTRAL_ON = T, NEUTRAL_PHYSICS_SIMPLE = F, NEUTRAL_TAPER_DIAGONAL = F, NEUTRAL_LINEAR_GM_TAPER = T, NEUTRAL_SINE_TAPER = T, DIFFUSION_ALL_EXPLICIT = F, NEUTRAL_BLAYER_DIAGNOSE = F, TURB_BLAYER_MIN = 0.000000000000000E+000 / [Note] ocean_nphysicsA.F90: USING ocean_nphysicsA_mod. ==> Note from ocean_nphysicsA_mod: using forward time step of (secs) 7200.00 [Note] ocean_nphysicsA.F90: use_gm_skew=.true. so will use GM-skewsion. [Note] ocean_nphysicsA.F90: neutral_physics_limit=.true. Will revert to horizontal diffusion for points where tracer is outside specified range. [Note] ocean_nphysicsA.F90: dm_taper=.true. Will use the tanh scheme of Danabasoglu and McWilliams to taper neutral physics in steep sloped regions [Note] ocean_nphysicsA.F90: neutral_linear_gm_taper=.true., so will linearly taper GM towards the surface when reaching steep neutral slopes in surface bdy. [Note] ocean_nphysicsA.F90: Running with neutral_sine_taper=.true., and so will use sine-taper on fluxes where eddy penetration depth and/or KPP hblt exceeds grid depth. [Note] ocean_nphysicsA.F90: Running with a nontrivial GM transport in steep neutral slope regions. [Note] ocean_nphysicsA.F90: Running w/ neutral_taper_diagonal=.false. so taper_diagonal = 0.0 Diagonal pieces of horizontal flux components are untapered regardless the neutral slope. flux dom neutral domain decomposition whalo = 1, ehalo = 1, shalo = 1, nhalo = 1 X-AXIS = 60 60 60 60 60 60 Y-AXIS = 40 40 40 40 40 NOTE from PE 15: mpp_domains_define.inc(mpp_define_io_domain): io domain will not be defined for flux dom neutral when one or both entry of io_layout is not positive BEGIN CHECKSUM(physics_driver_init):: CHECKSUM:: diff_cu_mo = 4CF3D79D55852ED0 CHECKSUM:: pbltop = A365705D75BEC0F8 CHECKSUM:: cush = D600000000000000 CHECKSUM:: cbmf = 0 CHECKSUM:: diff_t = 31C70AA6CB58456 CHECKSUM:: diff_m = 9744FE919D5136E6 CHECKSUM:: r_convect = 2D0000000000000 CHECKSUM:: radturbten = 0 CHECKSUM:: lw_tendency = DACF5C19915F1AC2 CHECKSUM:: lsc_cloud_area = 3740000000000000 CHECKSUM:: lsc_liquid = 3740000000000000 CHECKSUM:: lsc_ice = 3740000000000000 CHECKSUM:: lsc_droplet_number = 3740000000000000 CHECKSUM:: lsc_ice_number = 3740000000000000 CHECKSUM:: lsc_snow = 0 CHECKSUM:: lsc_rain = 0 CHECKSUM:: lsc_snow_size = 0 CHECKSUM:: lsc_rain_size = 0 Starting Threads : 1 NOTE from PE 0: fv_physics_init: num_phys_windows, numthreads 1 1 NOTE from PE 0: atmos_model_mod: Reading netCDF formatted restart file: INPUT/atmos_coupled.res.nc Finished initializing atmospheric model at 20220718 170004.979 20220718 170004.979: Memuse(MB) at atmos_model_init= 2.203E+02 2.305E+02 2.831E+00 2.239E+02 Starting to initialize land model at 20220718 170005.000 agm_array being read from restart. [chksum] checksum start agm_array 6517150202458889048 aredi_array being read from restart. [chksum] checksum start aredi_array 1904459692424298496 rossby_radius being read from restart. [chksum] checksum start rossby_radius 3342765363438361939 rossby_radius_raw being read from restart. [chksum] checksum start rossby_radius_ra 0 bczone_radius being read from restart. [chksum] checksum start bczone_radius 0 aredi_equal_agm=.false. and agm_read_restart=.false. => aredi_array set to static profiles. [Note] ocean_nphysics_util.F90: aredi_equal_agm=.false. allows aredi_array to differ from agm_array LAND MODEL domain decomposition whalo = 0, ehalo = 0, shalo = 0, nhalo = 0 X-AXIS = 29 29 28 29 29 Y-AXIS = 30 30 30 ---Neutral direction slope check I for linear stability of neutral diffusion--- With a neutral physics time step (secs) of 0.7200000E+04 the most stringent linear stability constraint was found at the following ocean cell: long( 232, 5) = -0.4850000E+02 lat ( 232, 5) = -0.7750000E+02 thick( 232, 5, 1) = 0.7353387E+01 aredi( 232, 5, 1) = 0.6000000E+03 delta_iso = 0.2048309E-01 is the maximum neutral direction slope available for linear stability of the neutral diffusion scheme. The namelist parameter smax should conservatively be <= delta_iso. ---Neutral direction slope check II for linear stability of neutral diffusion--- Assuming maximum Redi neutral diffusion slope of 0.2000000E-02 and neutral physics time step (secs) of 0.7200000E+04 the most stringent linear stability constraint was found at the following ocean cell: long( 1, 1) = -0.2795000E+03 lat ( 1, 1) = -0.8150000E+02 thick( 1, 1, 1)= 0.1000000E+02 A_max = 0.5706825E+04 (m^2/sec) is the maximum neutral diffusivity available for linear stability of the neutral diffusion scheme. Conservatively, neutral diffusivities used in the model should be less than A_max. -------------------------------------------------------------------------------- NOTE from PE 0: soil_mod: Reading NetCDF formatted restart file: INPUT/soil.res.nc &OCEAN_NPHYSICS_NEW_NML USE_THIS_MODULE = F, DRHODZ_SMOOTH_VERT = F, DRHODZ_SMOOTH_HORZ = F, SMAX = 1.000000000000000E-002, VEL_MICOM_SMOOTH = 0.200000000000000 / NOTE from PE 0: land_properties_init: Using static ASCII cover type dataset &OCEAN_NPHYSICS_UTIL_NEW_NML NUM_121_PASSES = 1 / [Note] ocean_nphysics_new.F90: NOT using ocean_nphysics_new. &OCEAN_SUBMESOSCALE_NML USE_THIS_MODULE = F, DEBUG_THIS_MODULE = F, DIAG_STEP = 1200, USE_HBLT_CONSTANT = F, USE_HBLT_EQUAL_MLD = T, SMOOTH_HBLT = T, SMOOTH_HBLT_NUM = 2, CONSTANT_HBLT = 100.000000000000 , COEFFICIENT_CE = 7.000000000000001E-002, TIME_CONSTANT = 86400.0000000000 , FRONT_LENGTH_CONST = 10000.0000000000 , MIN_KBLT = 4, MINIMUM_HBLT = 0.000000000000000E+000, SMOOTH_PSI = T, SMOOTH_PSI_NUM = 2, FRONT_LENGTH_DEFORM_RADIUS = T, LIMIT_PSI = T, USE_PSI_LEGACY = F, LIMIT_PSI_VELOCITY_SCALE = 0.500000000000000 , SUBMESO_LIMIT_FLUX = T, SMOOTH_ADVECT_TRANSPORT = T, SMOOTH_ADVECT_TRANSPORT_NUM = 2, SUBMESO_SKEW_FLUX = T, SUBMESO_ADVECT_FLUX = F, SUBMESO_ADVECT_UPWIND = T, SUBMESO_ADVECT_SWEBY = F, SUBMESO_ADVECT_LIMIT = F, SUBMESO_ADVECT_ZERO_BDY = F, SUBMESO_DIFFUSION = F, SUBMESO_DIFFUSION_SCALE = 1.00000000000000 , SUBMESO_DIFFUSION_BIHARMONIC = F / NOTE from PE 15: ==>Note: NOT using ocean_submesoscale_mod &OCEAN_LAP_FRICTION_NML LAP_FRICTION_SCHEME = general , DEBUG_THIS_MODULE = F, WRITE_A_RESTART = T / ==>Note from ocean_lap_friction_init: general laplacian friction scheme for B-grid is used. &OCEAN_LAPGEN_FRICTION_NML USE_THIS_MODULE = T, DEBUG_THIS_MODULE = F, BOTTOM_5POINT = T, K_SMAG_ISO = 2.00000000000000 , K_SMAG_ANISO = 0.000000000000000E+000, VEL_MICOM_ISO = 0.100000000000000 , VEL_MICOM_ANISO = 0.000000000000000E+000, EQ_VEL_MICOM_ISO = 0.000000000000000E+000, EQ_VEL_MICOM_ANISO = 0.000000000000000E+000, EQ_LAT_MICOM = 0.000000000000000E+000, EQUATORIAL_ZONAL = T, EQUATORIAL_ZONAL_LAT = 20.0000000000000 , EQUATORIAL_NO_SMAG = T, VISCOSITY_NCAR = T, VISCOSITY_NCAR_2000 = T, VISCOSITY_NCAR_2007 = F, NCAR_ISOTROPIC_OFF_EQUATOR = F, NCAR_ONLY_EQUATORIAL = T, VCONST_1 = 8000000.00000000 , VCONST_2 = 0.000000000000000E+000, VCONST_3 = 0.800000000000000 , VCONST_4 = 5.000000000000000E-009, VCONST_5 = 3, VCONST_6 = 300000000.000000 , VCONST_7 = 100.000000000000 , VCONST_8 = 45.0000000000000 , DEBUG_NCAR_A = F, DEBUG_NCAR_B = F, VISC_VEL_SCALE_LENGTH = 150000.000000000 , NEPTUNE = F, NEPTUNE_LENGTH_EQ = 1200.00000000000 , NEPTUNE_LENGTH_POLE = 3000.00000000000 , NEPTUNE_DEPTH_MIN = 100.000000000000 , NEPTUNE_SMOOTH = T, NEPTUNE_SMOOTH_NUM = 1, RESTRICT_POLAR_VISC = T, RESTRICT_POLAR_VISC_LAT = 60.0000000000000 , RESTRICT_POLAR_VISC_RATIO = 0.350000000000000 , NCAR_ISOTROPIC_AT_DEPTH = F, NCAR_ISOTROPIC_DEPTH = 4000.00000000000 , NCAR_ISOTROPIC_AT_DEPTH_VISC = 10000.0000000000 , DIVERGENCE_DAMP = F, DIVERGENCE_DAMP_VEL_MICOM = 0.000000000000000E+000, VISCOSITY_SCALE_BY_ROSSBY = F, VISCOSITY_SCALE_BY_ROSSBY_POWER = 2.00000000000000 , ASYNC_DOMAIN_UPDATE = F, BLOCKSIZE = 10, USE_SIDE_DRAG_FRICTION = F, SIDE_DRAG_FRICTION_SCALING = 1.00000000000000 , SIDE_DRAG_FRICTION_UVMAG_MAX = 10.0000000000000 , SIDE_DRAG_FRICTION_MAX = 1.00000000000000 / NOTE from PE 15: ==> NOTE: USING ocean_lapgen_friction_mod. ==>Note: not using asynchronous domain update in the vertical loop. This may be slow. ==> Note from ocean_lapgen_friction_mod: using forward time step of (secs) 7200.00 ==> Note: Will reduce horizontal friction to a 5point Laplacian on the bottom This helps to alleviate numerical problems with thin bottom partial cells. ==> NOTE: USING background horz viscosities according to NCAR CCSM2.0 algorithm. ==> NOTE: USING NCAR viscosity as formulated in 2000. NCAR vconst_1 (cm^2/sec) = 0.8000E+07 NCAR vconst_2 = 0.000000000000000E+000 NCAR vconst_3 = 0.800000000000000 NCAR vconst_4 (1/cm) = 5.000000000000000E-009 NCAR vconst_5 = 3 NCAR vconst_6 (cm^2/sec) = 0.3000E+09 NCAR vconst_7 (cm/sec) = 100.000000000000 NCAR vconst_8 (degrees) = 45.0000000000000 ==>ncar_only_equatorial=.true. =>ncar scheme only in band +/- lat 20.0000 ==> NOTE: USING horz isotropic viscosity via Smagorinsky. ==> NOTE: Setting horz anisotropic Smagorinsky viscosity to zero. If using anisotropic friction, zonally orient the friction within a latitudinal band of width 20.00000 degrees. Using restrict_polar_visc to lower visc_crit poleward of (deg) 60.00000 by an amount given by the fraction 0.35000 This approach is useful when coupling to ice, where effective (ocn+ice) visc > ocn visc. NOTE from PE 0: topography_mod: Reading NetCDF formatted input data file: INPUT/navy_topography.data.nc lap_viscosity being read from restart. [chksum] start lap_viscosity -4185879036935019400 &OCEAN_BIH_FRICTION_NML BIH_FRICTION_SCHEME = general , DEBUG_THIS_MODULE = F, WRITE_A_RESTART = T / ==>Note from ocean_bih_friction_init: general biharmonic friction scheme for B-grid is used. &OCEAN_BIHGEN_FRICTION_NML USE_THIS_MODULE = F, DEBUG_THIS_MODULE = F, K_SMAG_ISO = 2.00000000000000 , K_SMAG_ANISO = 0.000000000000000E+000, VEL_MICOM_ISO = 8.000000000000000E-002, VEL_MICOM_ANISO = 0.000000000000000E+000, EQ_VEL_MICOM_ISO = 0.000000000000000E+000, EQ_VEL_MICOM_ANISO = 0.000000000000000E+000, EQ_LAT_MICOM = 0.000000000000000E+000, VEL_MICOM_BOTTOM = 1.000000000000000E-002, BOTTOM_5POINT = T, EQUATORIAL_ZONAL = F, EQUATORIAL_ZONAL_LAT = 20.0000000000000 , VISC_CRIT_SCALE = 0.500000000000000 , READ_AISO_BIH_BACK = F, NCAR_BOUNDARY_SCALING = F, NCAR_RESCALE_POWER = 1, NCAR_VCONST_4 = 2.000000000000000E-008, NCAR_VCONST_5 = 3, NCAR_BOUNDARY_SCALING_READ = F, NEPTUNE = F, NEPTUNE_LENGTH_EQ = 4200.00000000000 , NEPTUNE_LENGTH_POLE = 17000.0000000000 , NEPTUNE_DEPTH_MIN = 100.000000000000 , NEPTUNE_SCALING = 1.00000000000000 , NEPTUNE_SMOOTH = T, NEPTUNE_SMOOTH_NUM = 1, VISC_DIVERGE_SCALING = 0.000000000000000E+000, USE_SIDE_DRAG_FRICTION = F, SIDE_DRAG_FRICTION_SCALING = 1.00000000000000 , SIDE_DRAG_FRICTION_UVMAG_MAX = 10.0000000000000 , SIDE_DRAG_FRICTION_MAX = 1.00000000000000 / NOTE from PE 15: ==> NOTE: NOT using ocean_bihgen_friction_mod. NOTE from PE 0: climap_albedo_mod: Reading NetCDF formatted input data file: INPUT/albedo.data.nc &OCEAN_MOMENTUM_SOURCE_NML VERBOSE_INIT = T, USE_THIS_MODULE = F, DEBUG_THIS_MODULE = F, USE_RAYLEIGH_DAMP_TABLE = F, RAYLEIGH_DAMP_EXP_FROM_BOTTOM = F, RAYLEIGH_DAMP_EXP_SCALE = 100.000000000000 , RAYLEIGH_DAMP_EXP_TIME = 864000.000000000 / NOTE from PE 15: ==>Note from ocean_momentum_source_mod: NOT USING this module &OCEAN_FORM_DRAG_NML VERBOSE_INIT = T, USE_THIS_MODULE = F, DEBUG_THIS_MODULE = F, USE_FORM_DRAG_AIKI = F, CPRIME_AIKI = 0.300000000000000 , FORM_DRAG_AIKI_BOTTOM_LAYER = F, FORM_DRAG_AIKI_BOTTOM_KLEVELS = 3, FORM_DRAG_AIKI_SCALE_BY_GRADH = F, FORM_DRAG_AIKI_GRADH_MAX = 5.000000000000000E-002, FORM_DRAG_AIKI_GRADH_POWER = 1.00000000000000 , FORM_DRAG_AIKI_SCALE_BY_GM = F, USE_FORM_DRAG_GBATCH = F, VISC_CBU_FORM_DRAG_MAX = 1.00000000000000 , VEL_FORM_DRAG_MAX = 1.00000000000000 , N_SQUARED_MIN = 1.000000000000000E-010, AGM_FORM_DRAG = 600.000000000000 , FORM_DRAG_GBATCH_SURF_LAYER = F, KSURF_BLAYER_MIN = 3, FORM_DRAG_GBATCH_ALPHA_F2 = F, FORM_DRAG_GBATCH_ALPHA = 300000000.000000 , FORM_DRAG_GBATCH_F2OVERN2 = F, FORM_DRAG_GBATCH_F2OVERNB2 = F, FORM_DRAG_GBATCH_SMOOTH_N2 = F, NUM_121_PASSES = 1, FORM_DRAG_GBATCH_F2OVERNO2 = F, FORM_DRAG_GBATCH_NO = 5.000000000000000E-003 / NOTE from PE 15: ==>Note from ocean_form_drag_mod: NOT USING this module &OCEAN_TRACER_ADVECT_NML DEBUG_THIS_MODULE = F, LIMIT_WITH_UPWIND = F, ADVECT_SWEBY_ALL = T, ZERO_TRACER_ADVECT_HORZ = F, ZERO_TRACER_ADVECT_VERT = F, WRITE_A_RESTART = T, PSOM_LIMIT_PRATHER = F, READ_BASIN_MASK = F, ASYNC_DOMAIN_UPDATE = F / ==>ocean_tracer_advect_mod: advect_sweby_all=.true. so all tracers advected with mdfl_sweby, regardless the settings in field_table. This method exploits mpp_update_domain capabilities and is faster in some cases. flux domain decomposition whalo = 1, ehalo = 1, shalo = 1, nhalo = 1 X-AXIS = 60 60 60 60 60 60 Y-AXIS = 40 40 40 40 40 fourth domain decomposition whalo = 2, ehalo = 2, shalo = 2, nhalo = 2 X-AXIS = 60 60 60 60 60 60 Y-AXIS = 40 40 40 40 40 NOTE from PE 15: mpp_domains_define.inc(mpp_define_io_domain): io domain will not be defined for fourth when one or both entry of io_layout is not positive sixth domain decomposition whalo = 3, ehalo = 3, shalo = 3, nhalo = 3 X-AXIS = 60 60 60 60 60 60 Y-AXIS = 40 40 40 40 40 NOTE from PE 15: mpp_domains_define.inc(mpp_define_io_domain): io domain will not be defined for sixth when one or both entry of io_layout is not positive quicker domain decomposition whalo = 2, ehalo = 2, shalo = 2, nhalo = 2 X-AXIS = 60 60 60 60 60 60 Y-AXIS = 40 40 40 40 40 NOTE from PE 15: mpp_domains_define.inc(mpp_define_io_domain): io domain will not be defined for quicker when one or both entry of io_layout is not positive mdfl domain decomposition whalo = 2, ehalo = 2, shalo = 2, nhalo = 2 X-AXIS = 60 60 60 60 60 60 Y-AXIS = 40 40 40 40 40 NOTE from PE 15: mpp_domains_define.inc(mpp_define_io_domain): io domain will not be defined for mdfl when one or both entry of io_layout is not positive mdppm domain decomposition whalo = 4, ehalo = 4, shalo = 4, nhalo = 4 X-AXIS = 60 60 60 60 60 60 Y-AXIS = 40 40 40 40 40 NOTE from PE 15: mpp_domains_define.inc(mpp_define_io_domain): io domain will not be defined for mdppm when one or both entry of io_layout is not positive mdmdt domain decomposition whalo = 4, ehalo = 4, shalo = 4, nhalo = 4 X-AXIS = 60 60 60 60 60 60 Y-AXIS = 40 40 40 40 40 NOTE from PE 15: mpp_domains_define.inc(mpp_define_io_domain): io domain will not be defined for mdmdt when one or both entry of io_layout is not positive &OCEAN_VELOCITY_ADVECT_NML DEBUG_THIS_MODULE = F, ZERO_VELOCITY_ADVECT_HORZ = F, ZERO_VELOCITY_ADVECT_VERT = F, VELOCITY_ADVECT_CENTERED = T, VELOCITY_ADVECT_UPWIND = F / ==>Note: running MOM with traditional second order centred advection of linear momentum. &OCEAN_CONVECT_NML USE_THIS_MODULE = T, CONVECT_NCON = F, NCON = 7, CONVECT_FULL_SCALAR = T, CONVECT_FULL_VECTOR = F / ==>Note: Using convective adjustment to stabilize gravitationally unstable water columns. ==>Note: Using the Rahmstorf full convection scheme as implemented on scalar machines at GFDL. Using preTEOS-10 implementation of convective adjustment. &OCEAN_SBC_NML TEMP_RESTORE_TSCALE = -10.0000000000000 , SALT_RESTORE_TSCALE = -60.0000000000000 , SALT_RESTORE_UNDER_ICE = F, SALT_RESTORE_AS_SALT_FLUX = F, ETA_RESTORE_TSCALE = -30.0000000000000 , ZERO_NET_PME_ETA_RESTORE = F, ROTATE_WINDS = F, TAUX_SINX = F, TAUY_SINY = F, USE_WATERFLUX = T, WATERFLUX_TAVG = F, MAX_ICE_THICKNESS = 4.00000000000000 , RUNOFFSPREAD = F, CALVINGSPREAD = F, USE_WATERFLUX_OVERRIDE_CALVING = F, USE_WATERFLUX_OVERRIDE_EVAP = F, USE_WATERFLUX_OVERRIDE_FPREC = F, SALINITY_REF = 35.0000000000000 , ZERO_NET_SALT_RESTORE = F, ZERO_NET_WATER_RESTORE = F, ZERO_NET_WATER_COUPLER = F, ZERO_NET_WATER_COUPLE_RESTORE = F, ZERO_NET_SALT_CORRECTION = F, ZERO_NET_WATER_CORRECTION = F, DEBUG_WATER_FLUXES = F, ZERO_WATER_FLUXES = F, ZERO_CALVING_FLUXES = F, ZERO_PME_FLUXES = F, ZERO_RUNOFF_FLUXES = F, ZERO_RIVER_FLUXES = F, CONVERT_RIVER_TO_PME = F, ZERO_HEAT_FLUXES = F, ZERO_SURFACE_STRESS = F, AVG_SFC_VELOCITY = T, AVG_SFC_TEMP_SALT_ETA = T, ICE_SALT_CONCENTRATION = 5.000000000000000E-003, OCEAN_ICE_SALT_LIMIT = 0.000000000000000E+000, RUNOFF_SALINITY = 0.000000000000000E+000, RUNOFF_TEMP_MIN = 0.000000000000000E+000, READ_RESTORE_MASK = F, RESTORE_MASK_GFDL = F, LAND_MODEL_HEAT_FLUXES = F, USE_FULL_PATM_FOR_SEA_LEVEL = F, MAX_DELTA_SALINITY_RESTORE = -0.500000000000000 , DO_FLUX_CORRECTION = F, SALINITY_RESTORE_LIMIT_LOWER = 0.000000000000000E+000, SALINITY_RESTORE_LIMIT_UPPER = 100.000000000000 , TEMP_CORRECTION_SCALE = 0.000000000000000E+000, SALT_CORRECTION_SCALE = 0.000000000000000E+000, TAU_X_CORRECTION_SCALE = 0.000000000000000E+000, TAU_Y_CORRECTION_SCALE = 0.000000000000000E+000, DO_BITWISE_EXACT_SUM = T, SBC_HEAT_FLUXES_CONST = F, SBC_HEAT_FLUXES_CONST_VALUE = 0.000000000000000E+000, SBC_HEAT_FLUXES_CONST_SEASONAL = F, USE_CONSTANT_SSS_FOR_RESTORE = F, CONSTANT_SSS_FOR_RESTORE = 35.0000000000000 , USE_CONSTANT_SST_FOR_RESTORE = F, CONSTANT_SST_FOR_RESTORE = 12.0000000000000 , USE_IDEAL_CALVING = F, USE_IDEAL_RUNOFF = F, CONSTANT_HLF = T, CONSTANT_HLV = T, READ_STOKES_DRIFT = F, DO_LANGMUIR = F, DO_USTAR_CORRECTION = T, DO_FRAZIL_REDIST = T / &OCEAN_SBC_OFAM_NML RESTORE_MASK_OFAM = F, RIVER_TEMP_OFAM = F / sbc domain decomposition whalo = 0, ehalo = 0, shalo = 0, nhalo = 0 X-AXIS = 60 60 60 60 60 60 Y-AXIS = 40 40 40 40 40 NOTE from PE 15: mpp_domains_define.inc(mpp_define_io_domain): io domain will not be defined for sbc when one or both entry of io_layout is not positive ==>Note from ocean_sbc_mod: if inputting river water, enable rivermix_mod to get river tracers into ocean. ==>Note from ocean_sbc_mod: temp_restore_tscale < 0. no surface restoring for t emp ==>Note from ocean_sbc_mod: salt_restore_tscale < 0. no surface restoring for s alt or pme. ==>Note from ocean_sbc_mod: eta_restore_tscale < 0. no surface restoring for et a ==>If coupling, then avg_sfc_velocity=.true. means will pass averaged ocean vel ocity to ice model. ==>If coupling, then avg_sfc_temp_salt_eta=.true. means will pass averaged sst, sss, eta to ice model. ==>Note from ocean_sbc_mod: Using constant latent heat of fusion at ocean surfa ce. ==>Note from ocean_sbc_mod: Using constant latent heat of evaporation at ocean surface. ==>Note that diagnose_sea_level_forcing==.false., so no "eta_tend_" fields will be diagnosed. &OCEAN_BBC_NML BMF_IMPLICIT = F, CDBOT = 1.000000000000000E-003, URESIDUAL = 5.000000000000000E-002, CDBOT_LAW_OF_WALL = F, LAW_OF_WALL_ROUGH_LENGTH = 1.000000000000000E-002, CDBOT_ROUGHNESS_LENGTH = F, USE_GEOTHERMAL_HEATING = F, CONVERT_GEOTHERMAL = 1.000000000000000E-003, CDBOT_HI = 3.000000000000000E-003, CDBOT_LO = 1.000000000000000E-003, CDBOT_GAMMA = 40.0000000000000 , UVMAG_MAX = 10.0000000000000 , BMF_MAX = 1.00000000000000 , DEBUG_THIS_MODULE = F, CDBOT_ROUGHNESS_UAMP = F, CDBOT_HH = 1100.00000000000 , CDBOT_UU = 1.00000000000000 , CDBOT_WAVE = F / &OCEAN_BBC_OFAM_NML READ_TIDE_SPEED = F, URESIDUAL2_MAX = 5.000000000000000E-002 / ==>Note: NOT reading tide_speed for ocean_vert_tidal_mod. NOTE from PE 15: ==>ocean_vert_tidal_mod: Setting tide_speed to default value. &OCEAN_SHORTWAVE_NML USE_THIS_MODULE = T, USE_SHORTWAVE_GFDL = T, USE_SHORTWAVE_CSIRO = F, USE_SHORTWAVE_JERLOV = F, USE_SHORTWAVE_EXT = F / NOTE from PE 15: ==>Note: USING shortwave_mod. &OCEAN_SHORTWAVE_GFDL_NML USE_THIS_MODULE = T, READ_CHL = T, CHL_DEFAULT = 8.000000000000000E-002, ZMAX_PEN = 100.000000000000 , SW_FRAC_TOP = 0.000000000000000E+000, DEBUG_THIS_MODULE = F, ENFORCE_SW_FRAC = T, OVERRIDE_F_VIS = F, SW_MOREL_FIXED_DEPTHS = T, OPTICS_FOR_UNIFORM_CHL = F, OPTICS_MOREL_ANTOINE = T, OPTICS_MANIZZA = F / NOTE from PE 15: ==>Note: USING shortwave_gfdl_mod. =>Note: Using shortwave penetration with GFDL formulaton & Morel-Antoine optics. ==>Note: enforce_sw_frac=.true. enforcing monotonic decrease of sw_frac with depth. ==>Warning: sw_morel_fixed_depths=.true. is unsuitable for time varying thicknesses. Time varying thicknesses are the norm in MOM, so recommend setting sw_morel_fixed_depths=.false. However, to reproduce MOM4.0 algorithm, then set sw_morel_fixed_depths=.true. NOTE from PE 15: ==>Note: Reading in chlorophyll-a from data file for shortwave penetration. =>Note: computing solar shortwave penetration. Assume stf has sw-radiation field included. Hence, solar shortwave penetration effects placed in sw_source will subtract out the effects of shortwave at k=1 to avoid double-counting. ==>Note: Setting optical model coefficients assuming nonuniform chl distribution. &OCEAN_SPONGES_TRACER_NML USE_THIS_MODULE = F, DAMP_COEFF_3D = F / &OCEAN_SPONGES_TRACER_OFAM_NML USE_ADAPTIVE_RESTORE = F, USE_SPONGE_AFTER_INIT = F, USE_NORMALISING = F, USE_HARD_THUMP = F, ATHRESH = 0.500000000000000 , TAUMIN = 720.000000000000 , LAMBDA = 8.300000000000000E-003, NPOWER = 1.00000000000000 , DAYS_TO_RESTORE = 1, SECS_TO_RESTORE = 0, DEFLATE = F, DEFLATE_FRACTION = 0.600000000000000 , LIMIT_TEMP = F, LIMIT_TEMP_MIN = -1.80000000000000 , LIMIT_TEMP_RESTORE = 10800.0000000000 , LIMIT_SALT = F, LIMIT_SALT_MIN = 1.000000000000000E-002, LIMIT_SALT_RESTORE = 3600.00000000000 / ==>Note from ocean_sponges_tracer_mod: NOT using ocean tracer sponges. &OCEAN_SPONGES_VELOCITY_NML USE_THIS_MODULE = F, DAMP_COEFF_3D = F / &OCEAN_SPONGES_VELOCITY_OFAM_NML USE_ADAPTIVE_RESTORE = F, USE_SPONGE_AFTER_INIT = F, USE_NORMALISING = F, USE_HARD_THUMP = F, ATHRESH = 0.500000000000000 , TAUMIN = 720.000000000000 , LAMBDA = 8.300000000000000E-003, NPOWER = 1.00000000000000 , DAYS_TO_RESTORE = 1, SECS_TO_RESTORE = 0 / ==>Note from ocean_sponges_velocity_mod: NOT using this module: no velocity spo nges. &OCEAN_SPONGES_ETA_NML USE_THIS_MODULE = F / &OCEAN_SPONGES_ETA_OFAM_NML USE_ADAPTIVE_RESTORE = F, USE_SPONGE_AFTER_INIT = F, USE_NORMALISING = F, USE_HARD_THUMP = F, ATHRESH = 0.500000000000000 , TAUMIN = 720.000000000000 , LAMBDA = 8.300000000000000E-003, NPOWER = 1.00000000000000 , DAYS_TO_RESTORE = 1, SECS_TO_RESTORE = 0 / ==>Note from ocean_sponges_eta_mod: NOT using this module. &OCEAN_XLANDMIX_NML VERBOSE_INIT = T, USE_THIS_MODULE = T, XLANDMIX_KMT = T / NOTE from PE 15: ==>Note from ocean_xlandmix_mod (ocean_xlandmix_init): USING this module ==>Note: Using xlandmix to connect tracers and mass between non-local ocean cel ls. ==> Note from ocean_xlandmix_mod: using forward time step of (secs) 7200.00 ==> Note from ocean_xlandmix_mod: allowing xlandmix to occur at k=kmt. ==>Warning: xlandmix has not been implemented to connect points across the tripolar fold. Defining extra tracer arrays for xland_domain over k-levels 1 through 28 The model local computational domain has a x,y halo = 1 1 This is smaller than the halo required for xlandmix. For xlandmix, will define a new domain type with halo = 4 xlandmix domain decomposition whalo = 4, ehalo = 4, shalo = 4, nhalo = 4 X-AXIS = 60 60 60 60 60 60 Y-AXIS = 40 40 40 40 40 NOTE from PE 15: mpp_domains_define.inc(mpp_define_io_domain): io domain will not be defined for xlandmix when one or both entry of io_layout is not positive ===== from ocean_xland_init ===== for crossland sea connection pair number 1 mix (i,j) gridpoint ( 274, 146) [long= -6.500 lat= 35.500] with (i,j) gridpoint ( 276, 146) [long= -4.500 lat= 35.500] from level 1 to 28 [depths of 0.000 to 360.155m] ===== from xlandvchk ===== for crossland sea communication pair number 1 mix I,J gridpoints ( 274, 146) and ( 276, 146) from level 1 to 28 (a depth range of 0.360155E+03 m) column volumes =0.362531E+13 and 0.362531E+13 m^3 simulated flow in = flow out = 0.550000E+06 m^3/sec, so mix 0.151711E-06 fraction of 1st column with 0.151711E-06 of 2nd column per sec. ===== from ocean_xland_init ===== for crossland sea connection pair number 2 mix (i,j) gridpoint ( 274, 147) [long= -6.500 lat= 36.500] with (i,j) gridpoint ( 276, 147) [long= -4.500 lat= 36.500] from level 1 to 28 [depths of 0.000 to 360.155m] ===== from xlandvchk ===== for crossland sea communication pair number 2 mix I,J gridpoints ( 274, 147) and ( 276, 147) from level 1 to 28 (a depth range of 0.360155E+03 m) column volumes =0.357963E+13 and 0.357963E+13 m^3 simulated flow in = flow out = 0.550000E+06 m^3/sec, ===== from ocean_xland_init ===== for crossland sea connection pair number 7 mix (i,j) gridpoint ( 290, 168) [long= 9.500 lat= 57.500] with (i,j) gridpoint ( 291, 166) [long= 10.500 lat= 55.500] from level 1 to 2 [depths of 0.000 to 20.000m] ===== from xlandvchk ===== for crossland sea communication pair number 7 mix I,J gridpoints ( 290, 168) and ( 291, 166) from level 1 to 2 (a depth range of 0.200000E+02 m) column volumes =0.132867E+12 and 0.140064E+12 m^3 simulated flow in = flow out = 0.150000E+05 m^3/sec, so mix 0.112895E-06 fraction of 1st column with 0.107094E-06 of 2nd column per sec. ===== from ocean_xland_init ===== for crossland sea connection pair number 8 mix (i,j) gridpoint ( 290, 167) [long= 9.500 lat= 56.500] with (i,j) gridpoint ( 291, 165) [long= 10.500 lat= 54.500] from level 1 to 2 [depths of 0.000 to 20.000m] ===== from xlandvchk ===== for crossland sea communication pair number 8 mix I,J gridpoints ( 290, 167) and ( 291, 165) from level 1 to 2 (a depth range of 0.200000E+02 m) column volumes =0.136486E+12 and 0.143600E+12 m^3 simulated flow in = flow out = 0.150000E+05 m^3/sec, so mix 0.109901E-06 fraction of 1st column with 0.104457E-06 of 2nd column per sec. ===== from ocean_xland_init ===== for crossland sea connection pair number 9 mix (i,j) gridpoint ( 300, 174) [long= 19.500 lat= 63.500] with (i,j) gridpoint ( 302, 174) [long= 21.500 lat= 63.500] from level 1 to 6 [depths of 0.000 to 60.000m] ===== from xlandvchk ===== for crossland sea communication pair number 9 mix I,J gridpoints ( 300, 174) and ( 302, 174) from level 1 to 6 (a depth range of 0.600000E+02 m) column volumes =0.331016E+12 and 0.331016E+12 m^3 simulated flow in = flow out = 0.300000E+04 m^3/sec, so mix 0.906301E-08 fraction of 1st column with 0.906301E-08 of 2nd column per sec. ===== from ocean_xland_init ===== for crossland sea connection pair number 9 mix (i,j) gridpoint ( 300, 174) [long= 19.500 lat= 63.500] with (i,j) gridpoint ( 302, 174) [long= 21.500 lat= 63.500] from level 1 to 6 [depths of 0.000 to 60.000m] ===== from xlandvchk ===== for crossland sea communication pair number 9 mix I,J gridpoints ( 300, 174) and ( 302, 174) from level 1 to 6 (a depth range of 0.600000E+02 m) column volumes =0.331016E+12 and 0.331016E+12 m^3 simulated flow in = flow out = 0.300000E+04 m^3/sec, so mix 0.906301E-08 fraction of 1st column with 0.906301E-08 of 2nd column per sec. ===== from ocean_xland_init ===== for crossland sea connection pair number 10 mix (i,j) gridpoint ( 301, 174) [long= 20.500 lat= 63.500] with (i,j) gridpoint ( 303, 174) [long= 22.500 lat= 63.500] from level 1 to 6 [depths of 0.000 to 60.000m] ===== from xlandvchk ===== for crossland sea communication pair number 10 mix I,J gridpoints ( 301, 174) and ( 303, 174) from level 1 to 6 (a depth range of 0.600000E+02 m) column volumes =0.331016E+12 and 0.331016E+12 m^3 simulated flow in = flow out = 0.300000E+04 m^3/sec, so mix 0.906301E-08 fraction of 1st column with 0.906301E-08 of 2nd column per sec. ===== from ocean_xland_init ===== for crossland sea connection pair number 5 mix (i,j) gridpoint ( 322, 123) [long= 41.500 lat= 14.459] with (i,j) gridpoint ( 324, 120) [long= 43.500 lat= 12.353] from level 1 to 20 [depths of 0.000 to 200.000m] ===== from xlandvchk ===== for crossland sea communication pair number 5 mix I,J gridpoints ( 322, 123) and ( 324, 120) so mix 0.153647E-06 fraction of 1st column with 0.153647E-06 of 2nd column per sec. ===== from ocean_xland_init ===== for crossland sea connection pair number 3 mix (i,j) gridpoint ( 305, 151) [long= 24.500 lat= 40.500] with (i,j) gridpoint ( 309, 152) [long= 28.500 lat= 41.500] from level 1 to 6 [depths of 0.000 to 60.000m] ===== from xlandvchk ===== for crossland sea communication pair number 3 mix I,J gridpoints ( 305, 151) and ( 309, 152) from level 1 to 6 (a depth range of 0.600000E+02 m) column volumes =0.564114E+12 and 0.555619E+12 m^3 simulated flow in = flow out = 0.100000E+05 m^3/sec, so mix 0.177269E-07 fraction of 1st column with 0.179979E-07 of 2nd column per sec. ===== from ocean_xland_init ===== for crossland sea connection pair number 4 mix (i,j) gridpoint ( 306, 151) [long= 25.500 lat= 40.500] with (i,j) gridpoint ( 309, 153) [long= 28.500 lat= 42.500] from level 1 to 6 [depths of 0.000 to 60.000m] ===== from xlandvchk ===== for crossland sea communication pair number 4 mix I,J gridpoints ( 306, 151) and ( 309, 153) from level 1 to 6 (a depth range of 0.600000E+02 m) column volumes =0.564114E+12 and 0.546956E+12 m^3 simulated flow in = flow out = 0.100000E+05 m^3/sec, so mix 0.177269E-07 fraction of 1st column with 0.182830E-07 of 2nd column per sec. ===== from ocean_xland_init ===== for crossland sea connection pair number 5 mix (i,j) gridpoint ( 322, 123) [long= 41.500 lat= 14.459] with (i,j) gridpoint ( 324, 120) [long= 43.500 lat= 12.353] from level 1 to 20 [depths of 0.000 to 200.000m] ===== from xlandvchk ===== for crossland sea communication pair number 5 mix I,J gridpoints ( 322, 123) and ( 324, 120) from level 1 to 20 (a depth range of 0.200000E+03 m) column volumes =0.173288E+13 and 0.165270E+13 m^3 simulated flow in = flow out = 0.180000E+06 m^3/sec, so mix 0.103874E-06 fraction of 1st column with 0.108913E-06 of 2nd column per sec. ===== from ocean_xland_init ===== for crossland sea connection pair number 6 mix (i,j) gridpoint ( 323, 123) [long= 42.500 lat= 14.459] with (i,j) gridpoint ( 324, 119) [long= 43.500 lat= 11.674] from level 1 to 20 [depths of 0.000 to 200.000m] ===== from xlandvchk ===== for crossland sea communication pair number 6 mix I,J gridpoints ( 323, 123) and ( 324, 119) from level 1 to 20 (a depth range of 0.200000E+03 m) column volumes =0.173288E+13 and 0.163671E+13 m^3 simulated flow in = flow out = 0.180000E+06 m^3/sec, so mix 0.103874E-06 fraction of 1st column with 0.109976E-06 of 2nd column per sec. ===== from ocean_xland_init ===== for crossland sea connection pair number 15 mix (i,j) gridpoint ( 336, 137) [long= 55.500 lat= 26.520] with (i,j) gridpoint ( 338, 136) [long= 57.500 lat= 25.542] from level 1 to 4 [depths of 0.000 to 40.000m] ===== from xlandvchk ===== for crossland sea communication pair number 15 mix I,J gridpoints ( 336, 137) and ( 338, 136) from level 1 to 4 (a depth range of 0.400000E+02 m) column volumes =0.434786E+12 and 0.433571E+12 m^3 simulated flow in = flow out = 0.500000E+05 m^3/sec, so mix 0.114999E-06 fraction of 1st column with 0.115321E-06 of 2nd column per sec. ===== from ocean_xland_init ===== for crossland sea connection pair number 16 mix (i,j) gridpoint ( 336, 136) [long= 55.500 lat= 25.542] with (i,j) gridpoint ( 338, 135) [long= 57.500 lat= 24.577] from level 1 to 4 [depths of 0.000 to 40.000m] ===== from xlandvchk ===== for crossland sea communication pair number 16 mix I,J gridpoints ( 336, 136) and ( 338, 135) from level 1 to 4 (a depth range of 0.400000E+02 m) column volumes =0.433571E+12 and 0.431044E+12 m^3 simulated flow in = flow out = 0.500000E+05 m^3/sec, so mix 0.115321E-06 fraction of 1st column with 0.115997E-06 of 2nd column per sec. ===== from ocean_xland_init ===== for crossland sea connection pair number 11 mix (i,j) gridpoint ( 211, 173) [long= -69.500 lat= 62.500] with (i,j) gridpoint ( 212, 171) [long= -68.500 lat= 60.500] from level 1 to 12 [depths of 0.000 to 120.000m] ===== from xlandvchk ===== for crossland sea communication pair number 11 mix I,J gridpoints ( 211, 173) and ( 212, 171) from level 1 to 12 (a depth range of 0.120000E+03 m) column volumes =0.685104E+12 and 0.730617E+12 m^3 simulated flow in = flow out = 0.450000E+06 m^3/sec, so mix 0.656834E-06 fraction of 1st column with 0.615917E-06 of 2nd column per sec. ===== from ocean_xland_init ===== for crossland sea connection pair number 12 mix (i,j) gridpoint ( 211, 172) [long= -69.500 lat= 61.500] with (i,j) gridpoint ( 212, 170) [long= -68.500 lat= 59.500] from level 1 to 12 [depths of 0.000 to 120.000m] ===== from xlandvchk ===== for crossland sea communication pair number 12 mix I,J gridpoints ( 211, 172) and ( 212, 170) from level 1 to 12 (a depth range of 0.120000E+03 m) column volumes =0.707969E+12 and 0.753043E+12 m^3 simulated flow in = flow out = 0.450000E+06 m^3/sec, so mix 0.635621E-06 fraction of 1st column with 0.597575E-06 of 2nd column per sec. ===== from ocean_xland_init ===== for crossland sea connection pair number 13 mix (i,j) gridpoint ( 199, 188) [long= -85.119 lat= 69.329] with (i,j) gridpoint ( 203, 187) [long= -81.008 lat= 69.749] from level 1 to 4 [depths of 0.000 to 40.000m] ===== from xlandvchk ===== for crossland sea communication pair number 13 mix I,J gridpoints ( 199, 188) and ( 203, 187) from level 1 to 4 (a depth range of 0.400000E+02 m) column volumes =0.743800E+11 and 0.843998E+11 m^3 simulated flow in = flow out = 0.350000E+05 m^3/sec, so mix 0.470556E-06 fraction of 1st column with 0.414693E-06 of 2nd column per sec. ===== from ocean_xland_init ===== for crossland sea connection pair number 14 mix (i,j) gridpoint ( 198, 188) [long= -85.942 lat= 69.117] with (i,j) gridpoint ( 202, 187) [long= -81.876 lat= 69.562] from level 1 to 4 [depths of 0.000 to 40.000m] ===== from xlandvchk ===== for crossland sea communication pair number 14 mix I,J gridpoints ( 198, 188) and ( 202, 187) from level 1 to 4 (a depth range of 0.400000E+02 m) column volumes =0.715971E+11 and 0.819460E+11 m^3 simulated flow in = flow out = 0.350000E+05 m^3/sec, so mix 0.488847E-06 fraction of 1st column with 0.427110E-06 of 2nd column per sec. from level 1 to 20 (a depth range of 0.200000E+03 m) column volumes =0.173288E+13 and 0.165270E+13 m^3 simulated flow in = flow out = 0.180000E+06 m^3/sec, so mix 0.103874E-06 fraction of 1st column with 0.108913E-06 of 2nd column per sec. ===== from ocean_xland_init ===== for crossland sea connection pair number 6 mix (i,j) gridpoint ( 323, 123) [long= 42.500 lat= 14.459] with (i,j) gridpoint ( 324, 119) [long= 43.500 lat= 11.674] from level 1 to 20 [depths of 0.000 to 200.000m] ===== from xlandvchk ===== for crossland sea communication pair number 6 mix I,J gridpoints ( 323, 123) and ( 324, 119) from level 1 to 20 (a depth range of 0.200000E+03 m) column volumes =0.173288E+13 and 0.163671E+13 m^3 simulated flow in = flow out = 0.180000E+06 m^3/sec, so mix 0.103874E-06 fraction of 1st column with 0.109976E-06 of 2nd column per sec. &OCEAN_XLANDINSERT_NML VERBOSE_INIT = T, USE_THIS_MODULE = T, DEBUG_THIS_MODULE = F / NOTE from PE 15: ==>Note from ocean_xlandinsert_mod (ocean_xlandinsert_init): USING this module ==>NOte: Using xlandinsert to connect tracers and mass between non-local ocean cells. ==> Note from ocean_xlandinsert_mod: using time step of (secs) 7200.00 ==>Warning: xlandinsert has not been implemented to connect points across the tripolar fold. Using xlandinsert to connect tracers and surface height between non-local ocean cells. Defining extra tracer arrays for xland_domain over k-levels 1 through 20 The model local computational domain has a x,y halo = 1 1 This is smaller than the halo required for xlandinsert. For xlandinsert, define a new domain type with halo = 4 xlandinsert domain decomposition whalo = 4, ehalo = 4, shalo = 4, nhalo = 4 X-AXIS = 60 60 60 60 60 60 Y-AXIS = 40 40 40 40 40 NOTE from PE 15: mpp_domains_define.inc(mpp_define_io_domain): io domain will not be defined for xlandinsert when one or both entry of io_layout is not positive ===== from ocean_xlandinsert_init ===== for crossland sea connection pair number 11 mix (i,j) gridpoint ( 211, 173) [long= -69.500 lat= 62.500] with (i,j) gridpoint ( 212, 171) [long= -68.500 lat= 60.500] from level 1 to 12 [depths of 0.000 to 120.000m] ===== from ocean_xlandinsert_init ===== for crossland sea connection pair number 7 mix (i,j) gridpoint ( 290, 168) [long= 9.500 lat= 57.500] with (i,j) gridpoint ( 291, 166) [long= 10.500 lat= 55.500] from level 1 to 2 [depths of 0.000 to 20.000m] Time scale for insertion via upwind advection is (sec) 86400.000000 ===== from ocean_xlandinsert_init ===== for crossland sea connection pair number 8 mix (i,j) gridpoint ( 290, 167) [long= 9.500 lat= 56.500] with (i,j) gridpoint ( 291, 165) [long= 10.500 lat= 54.500] from level 1 to 2 [depths of 0.000 to 20.000m] Time scale for insertion via upwind advection is (sec) 86400.000000 ===== from ocean_xlandinsert_init ===== for crossland sea connection pair number 9 mix (i,j) gridpoint ( 300, 174) [long= 19.500 lat= 63.500] with (i,j) gridpoint ( 302, 174) [long= 21.500 lat= 63.500] from level 1 to 6 [depths of 0.000 to 60.000m] Time scale for insertion via upwind advection is (sec) 86400.000000 ===== from ocean_xlandinsert_init ===== for crossland sea connection pair number 10 mix (i,j) gridpoint ( 301, 174) [long= 20.500 lat= 63.500] with (i,j) gridpoint ( 303, 174) [long= 22.500 lat= 63.500] from level 1 to 6 [depths of 0.000 to 60.000m] Time scale for insertion via upwind advection is (sec) 86400.000000 ===== from ocean_xlandinsert_init ===== for crossland sea connection pair number 5 mix (i,j) gridpoint ( 322, 123) [long= 41.500 lat= 14.459] with (i,j) gridpoint ( 324, 120) [long= 43.500 lat= 12.353] from level 1 to 20 [depths of 0.000 to 200.000m] Time scale for insertion via upwind advection is (sec) 86400.000000 ===== from ocean_xlandinsert_init ===== for crossland sea connection pair number 6 mix (i,j) gridpoint ( 323, 123) [long= 42.500 lat= 14.459] with (i,j) gridpoint ( 324, 119) [long= 43.500 lat= 11.674] from level 1 to 20 [depths of 0.000 to 200.000m] Time scale for insertion via upwind advection is (sec) 86400.000000 ===== from ocean_xlandinsert_init ===== for crossland sea connection pair number 1 mix (i,j) gridpoint ( 274, 146) [long= -6.500 lat= 35.500] with (i,j) gridpoint ( 276, 146) [long= -4.500 lat= 35.500] from level 1 to 18 [depths of 0.000 to 180.000m] Time scale for insertion via upwind advection is (sec) 86400.000000 ===== from ocean_xlandinsert_init ===== for crossland sea connection pair number 2 mix (i,j) gridpoint ( 274, 147) [long= -6.500 lat= 36.500] with (i,j) gridpoint ( 276, 147) [long= -4.500 lat= 36.500] from level 1 to 18 [depths of 0.000 to 180.000m] Time scale for insertion via upwind advection is (sec) 86400.000000 ===== from ocean_xlandinsert_init ===== for crossland sea connection pair number 3 mix (i,j) gridpoint ( 305, 151) [long= 24.500 lat= 40.500] with (i,j) gridpoint ( 309, 152) [long= 28.500 lat= 41.500] from level 1 to 6 [depths of 0.000 to 60.000m] Time scale for insertion via upwind advection is (sec) 86400.000000 ===== from ocean_xlandinsert_init ===== for crossland sea connection pair number 4 mix (i,j) gridpoint ( 306, 151) [long= 25.500 lat= 40.500] with (i,j) gridpoint ( 309, 153) [long= 28.500 lat= 42.500] from level 1 to 6 [depths of 0.000 to 60.000m] Time scale for insertion via upwind advection is (sec) 86400.000000 ===== from ocean_xlandinsert_init ===== for crossland sea connection pair number 5 mix (i,j) gridpoint ( 322, 123) [long= 41.500 lat= 14.459] with (i,j) gridpoint ( 324, 120) [long= 43.500 lat= 12.353] from level 1 to 20 [depths of 0.000 to 200.000m] Time scale for insertion via upwind advection is (sec) 86400.000000 ===== from ocean_xlandinsert_init ===== for crossland sea connection pair number 6 mix (i,j) gridpoint ( 323, 123) [long= 42.500 lat= 14.459] with (i,j) gridpoint ( 324, 119) [long= 43.500 lat= 11.674] from level 1 to 20 [depths of 0.000 to 200.000m] Time scale for insertion via upwind advection is (sec) 86400.000000 ===== from ocean_xlandinsert_init ===== for crossland sea connection pair number 15 mix (i,j) gridpoint ( 336, 137) [long= 55.500 lat= 26.520] with (i,j) gridpoint ( 338, 136) [long= 57.500 lat= 25.542] from level 1 to 4 [depths of 0.000 to 40.000m] Time scale for insertion via upwind advection is (sec) 86400.000000 ===== from ocean_xlandinsert_init ===== for crossland sea connection pair number 16 mix (i,j) gridpoint ( 336, 136) [long= 55.500 lat= 25.542] with (i,j) gridpoint ( 338, 135) [long= 57.500 lat= 24.577] from level 1 to 4 [depths of 0.000 to 40.000m] Time scale for insertion via upwind advection is (sec) 86400.000000 Time scale for insertion via upwind advection is (sec) 86400.000000 ===== from ocean_xlandinsert_init ===== for crossland sea connection pair number 12 mix (i,j) gridpoint ( 211, 172) [long= -69.500 lat= 61.500] with (i,j) gridpoint ( 212, 170) [long= -68.500 lat= 59.500] from level 1 to 12 [depths of 0.000 to 120.000m] Time scale for insertion via upwind advection is (sec) 86400.000000 ===== from ocean_xlandinsert_init ===== for crossland sea connection pair number 13 mix (i,j) gridpoint ( 199, 188) [long= -85.119 lat= 69.329] with (i,j) gridpoint ( 203, 187) [long= -81.008 lat= 69.749] from level 1 to 4 [depths of 0.000 to 40.000m] Time scale for insertion via upwind advection is (sec) 86400.000000 ===== from ocean_xlandinsert_init ===== for crossland sea connection pair number 14 mix (i,j) gridpoint ( 198, 188) [long= -85.942 lat= 69.117] with (i,j) gridpoint ( 202, 187) [long= -81.876 lat= 69.562] from level 1 to 4 [depths of 0.000 to 40.000m] Time scale for insertion via upwind advection is (sec) 86400.000000 ===== from ocean_xlandinsert_init ===== for crossland sea connection pair number 9 mix (i,j) gridpoint ( 300, 174) [long= 19.500 lat= 63.500] with (i,j) gridpoint ( 302, 174) [long= 21.500 lat= 63.500] from level 1 to 6 [depths of 0.000 to 60.000m] Time scale for insertion via upwind advection is (sec) 86400.000000 ==>Warning: ocean_riverspread_init: n<1 for riverspread table. Will NOT use tab le for ocean_riverspread. &OCEAN_RIVERSPREAD_NML USE_THIS_MODULE = F, DEBUG_THIS_MODULE = F, RIVERSPREAD_DIFFUSION = F, RIVERSPREAD_DIFFUSION_PASSES = 0, VEL_MICOM_SMOOTH = 0.200000000000000 / NOTE from PE 15: ==>From ocean_riverspread_mod: NOT using riverspread module, yet there are n > 0 riverspread points. &OCEAN_RIVERMIX_NML USE_THIS_MODULE = T, DEBUG_THIS_MODULE = F, DEBUG_ALL_IN_TOP_CELL = F, DEBUG_THIS_MODULE_HEAT = F, RIVER_DIFFUSE_TEMP = F, RIVER_DIFFUSE_SALT = F, RIVER_DIFFUSION_THICKNESS = 0.000000000000000E+000, RIVER_DIFFUSIVITY = 0.000000000000000E+000, DISCHARGE_COMBINE_RUNOFF_CALVE = T, RIVER_INSERTION_THICKNESS = 40.0000000000000 , RUNOFF_INSERTION_THICKNESS = 0.000000000000000E+000, CALVING_INSERTION_THICKNESS = 0.000000000000000E+000, DO_BITWISE_EXACT_SUM = F / NOTE from PE 15: ==>From ocean_rivermix_mod: Using rivermix module to mix liquid and/or solid runoff into the ocean. ==>Note: discharging calving+runoff together. The alternative is to separately discharge. ==>Note: if using waterflux and rivers, then will discharge river tracer over 4 grid points in vertical &OCEAN_OVEREXCHANGE_NML USE_THIS_MODULE = F, DEBUG_THIS_MODULE = F, OVEREXCH_NPTS = 4, OVEREXCH_WIDTH = 1, OVEREXCH_WEIGHT_FAR = F, OVERFLOW_MU = 1.000000000000000E-004, OVERFLOW_DELTA = 0.333300000000000 , OVERFLOW_UMAX = 5.00000000000000 , DO_BITWISE_EXACT_SUM = F, OVEREXCH_STABILITY = 0.250000000000000 , OVEREXCH_MIN_THICKNESS = 4.00000000000000 , OVEREXCH_CHECK_EXTREMA = F / NOTE from PE 15: ==>From ocean_overexchange_mod: NOT using overflow exchange scheme. &OCEAN_MIXDOWNSLOPE_NML USE_THIS_MODULE = F, DEBUG_THIS_MODULE = F, MIXDOWNSLOPE_NPTS = 4, MIXDOWNSLOPE_WIDTH = 1, MIXDOWNSLOPE_WEIGHT_FAR = F, MIXDOWNSLOPE_FRAC_CENTRAL = 0.250000000000000 , DO_BITWISE_EXACT_SUM = F, READ_MIXDOWNSLOPE_MASK = F, MIXDOWNSLOPE_MASK_GFDL = F / NOTE from PE 15: ==>From ocean_mixdownslope_mod: NOT using downslope mixing scheme. &OCEAN_OVERFLOW_NML USE_THIS_MODULE = F, DEBUG_THIS_MODULE = F, OVERFLOW_MU = 1.000000000000000E-004, OVERFLOW_DELTA = 0.333300000000000 , OVERFLOW_UMAX = 1.000000000000000E-002, DO_BITWISE_EXACT_SUM = F, NO_RETURN_FLOW = F, TRANSPORT_UNITS = Sv / NOTE from PE 15: ==>From ocean_overflow_mod: NOT using Campin and Goosse overflow scheme. ==>Warning: ocean_overflow_OFP_init found n_src < 1 for overflow_OFP table. Will NOT use ocean_overflow_OFP. ==>Note from fm_util_mod(fm_util_start_namelist)[ocean_age_tracer_mod(ocean_age _tracer_start)]: Processing namelist ocean_age_tracer/global ==>Note from ocean_age_tracer_mod(ocean_age_tracer_start): Assigning month 1 ==>Note from ocean_age_tracer_mod(set_array): age_global region: 1 ==>Note from ocean_age_tracer_mod(ocean_age_tracer_start): Duplicating month 1 as 2 ==>Note from ocean_age_tracer_mod(ocean_age_tracer_start): Duplicating month 1 as 3 ==>Note from ocean_age_tracer_mod(ocean_age_tracer_start): Duplicating month 1 as 4 ==>Note from ocean_age_tracer_mod(ocean_age_tracer_start): Duplicating month 1 as 5 ==>Note from ocean_age_tracer_mod(ocean_age_tracer_start): Duplicating month 1 as 6 ==>Note from ocean_age_tracer_mod(ocean_age_tracer_start): Duplicating month 1 as 7 ==>Note from ocean_age_tracer_mod(ocean_age_tracer_start): Duplicating month 1 as 8 ==>Note from ocean_age_tracer_mod(ocean_age_tracer_start): Duplicating month 1 as 9 ==>Note from ocean_age_tracer_mod(ocean_age_tracer_start): Duplicating month 1 as 10 ==>Note from ocean_age_tracer_mod(ocean_age_tracer_start): Duplicating month 1 as 11 ==>Note from ocean_age_tracer_mod(ocean_age_tracer_start): Duplicating month 1 as 12 &OCEAN_ADV_VEL_DIAG_NML MAX_CFL_VALUE = 100.000000000000 , LARGE_CFL_VALUE = 10.0000000000000 , VERBOSE_CFL = F, DIAG_STEP = 12 / ==>Maximum T-->U remapping error = 9.558E-06 m/s at (i,j) = (198,175), (lon,lat) = ( -82.00, 65.00) ==>Note: T-->U remapping error will be small (i.e., order 1e-20) only for spherical grids. &OCEAN_TRACER_DIAG_NML TRACER_CONSERVE_DAYS = 1.00000000000000 , DIAG_STEP = 12, PSU2PPT = 1.00486700000000 , DEBUG_DIAGNOSE_MIXINGA = F, DEBUG_DIAGNOSE_MIXINGB = F, DEBUG_DIAGNOSE_MIXINGC = F, DEBUG_DIAGNOSE_MIXINGD = F, SMOOTH_KAPPA_SORT = 0, RHO_GRAD_MIN = 1.000000000000000E-005, RHO_GRAD_MAX = 1.000000000000000E+028, BUOYANCY_CRIT = 3.000000000000000E-004, DO_BITWISE_EXACT_SUM = F, FRAZIL_FACTOR = 1.00000000000000 , SMOOTH_MLD = F, SMOOTH_MLD_FOR_SUBDUCTION = T, DTHETA_CRIT = 2.00000000000000 / Note: tracer and mass/volume conservation tests based on time_tendency==twolevel. Note: Set frazil_factor = 1.00 for computation of heat diagnostics. Be sure this agrees with the value set in nml for ocean_frazil_mod &OCEAN_VELOCITY_DIAG_NML DIAG_STEP = 12, ENERGY_DIAG_STEP = 12, DO_BITWISE_EXACT_SUM = F, DEBUG_THIS_MODULE = F, MAX_CFL_VALUE = 100.000000000000 , LARGE_CFL_VALUE = 10.0000000000000 , VERBOSE_CFL = F, LAND_CELL_NUM_MAX = 100 / &OCEAN_INCREMENT_ETA_NML USE_THIS_MODULE = F, FRACTION_INCREMENT = 1.00000000000000 , DAYS_TO_INCREMENT = 1, SECS_TO_INCREMENT = 0 / &OCEAN_INCREMENT_TRACER_NML USE_THIS_MODULE = F, FRACTION_INCREMENT = 1.00000000000000 , DAYS_TO_INCREMENT = 1, SECS_TO_INCREMENT = 0 / &OCEAN_INCREMENT_VELOCITY_NML USE_THIS_MODULE = F, FRACTION_INCREMENT = 1.00000000000000 , DAYS_TO_INCREMENT = 1, SECS_TO_INCREMENT = 0 / &OCEAN_WAVE_NML WAVEDAMP = -10.0000000000000 , DAMP_WHERE_ICE = T, WRITE_A_RESTART = T, DEBUG_THIS_MODULE = F, USE_TMA = T, FILTER_WAVE_MOM = T, USE_THIS_MODULE = F / ==>Note: Not using the idealized ocean surface wave module. &ODA_NML USE_THIS_MODULE = F, ASSIM_METHOD = NO_ASSIM, ASSIM_START_LAT = -90.0000000000000 , ASSIM_END_LAT = 90.0000000000000 , NK_ASM = 0, ASSIM_INTERVAL = 24, DO_IAU = T, DO_CONVECT_ADJUST = T, MAX_PROFILES = 100000, MAX_SFC_OBS = 100, SAVE_OMF_SNAPS = T, SAVE_OMA_SNAPS = F, ASSIM_LAYOUT = 2*0 / &OCEAN_DRIFTERS_NML USE_THIS_MODULE = F, OUTPUT_INTERVAL = 1 / ======== COMPLETED MOM INITIALIZATION ======== Finished initializing ocean model at 20220718 170005.289 20220718 170005.322: Memuse(MB) at ocean_model_init= 5.128E+02 5.241E+02 3.247E+00 5.187E+02 Starting to initialize data_override at 20220718 170005.322 Finished initializing data_override at 20220718 170005.342 Finished initializing component models at 20220718 170005.342 NOTE from PE 0: vegetation_mod: Reading NetCDF formatted restart file: INPUT/vegetation.res.nc Finished initializing land model at 20220718 170005.412 20220718 170005.413: Memuse(MB) at land_model_init= 2.368E+02 2.446E+02 2.436E+00 2.389E+02 Starting to initialize ice model at 20220718 170005.418 &ICE_MODEL_NML MOM_ROUGH_ICE = 1.000000000000000E-004, HEAT_ROUGH_ICE = 1.000000000000000E-004, P0 = 27500.0000000000 , C0 = 20.0000000000000 , CDW = 3.240000000000000E-003, WD_TURN = 0.000000000000000E+000, KMELT = 240.000000000000 , ALB_SNO = 0.800000000000000 , ALB_ICE = 0.582600000000000 , PEN_ICE = 0.300000000000000 , OPT_DEP_ICE = 0.670000000000000 , NSTEPS_DYN = 72, NSTEPS_ADV = 1, NUM_PART = 6, ATMOS_WINDS = T, SLAB_ICE = F, SPEC_ICE = F, ICE_BULK_SALIN = 5.000000000000000E-003, LAYOUT = 2*0, DO_ICE_RESTORE = F, DO_ICE_LIMIT = F, MAX_ICE_LIMIT = 4.00000000000000 , ICE_RESTORE_TIMESCALE = 5.00000000000000 , SLP2OCEAN = F, CONSERVATION_CHECK = T, T_RANGE_MELT = 10.0000000000000 , CM2_BUGS = T, KS = 0.310000000000000 , H_LO_LIM = 0.000000000000000E+000, VERBOSE = F, DO_ICEBERGS = F, ADD_DIURNAL_SW = F, IO_LAYOUT = 2*0, CHANNEL_VISCOSITY = 0.000000000000000E+000, SMAG_OCN = 0.150000000000000 , SSH_GRAVITY = 9.81000000000000 , CHAN_CFL_LIMIT = 0.250000000000000 , DO_SUN_ANGLE_FOR_ALB = F, MASK_TABLE = INPUT/ice_mask_table , REPRODUCE_SIENA_201303 = T / ==>Note from ice_grid_mod(set_ice_grid): read grid from mosaic version grid NOTE from PE 0: ==>Note from ice_grid_mod: x_boundary_type is cyclic NOTE from PE 0: ==>Note from ice_grid_mod: y_boundary_type is fold_north_edge ice model domain decomposition whalo = 1, ehalo = 1, shalo = 1, nhalo = 1 X-AXIS = 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 Y-AXIS = 200 NOTE from PE 0: mpp_domains_define.inc(mpp_define_io_domain): io domain will not be defined for ice model when one or both entry of io_layout is not positive ice_nohalo domain decomposition whalo = 0, ehalo = 0, shalo = 0, nhalo = 0 X-AXIS = 24 24 24 24 24 24 24 24 24 24 24 24 24 24 24 Y-AXIS = 200 NOTE from PE 0: mpp_domains_define.inc(mpp_define_io_domain): io domain will not be defined for ice_nohalo when one or both entry of io_layout is not positive WARNING from PE 11: ice_model_init: Restart file does not contain flux_sw_* subcomponents! WARNING from PE 4: ice_model_init: Restart file does not contain flux_sw_* subcomponents! WARNING from PE 3: ice_model_init: Restart file does not contain flux_sw_* subcomponents! WARNING from PE 7: ice_model_init: Restart file does not contain flux_sw_* subcomponents! WARNING from PE 9: ice_model_init: Restart file does not contain flux_sw_* subcomponents! WARNING from PE 10: ice_model_init: Restart file does not contain flux_sw_* subcomponents! WARNING from PE 13: ice_model_init: Restart file does not contain flux_sw_* subcomponents! WARNING from PE 12: ice_model_init: Restart file does not contain flux_sw_* subcomponents! WARNING from PE 14: ice_model_init: Restart file does not contain flux_sw_* subcomponents! WARNING from PE 2: ice_model_init: Restart file does not contain flux_sw_* subcomponents! WARNING from PE 6: ice_model_init: Restart file does not contain flux_sw_* subcomponents! WARNING from PE 8: ice_model_init: Restart file does not contain flux_sw_* subcomponents! WARNING from PE 0: ice_model_init: Restart file does not contain flux_sw_* subcomponents! WARNING from PE 0: ice_model_init: Restart file does not contain flux_sw_* subcomponents! WARNING from PE 5: ice_model_init: Restart file does not contain flux_sw_* subcomponents! WARNING from PE 1: ice_model_init: Restart file does not contain flux_sw_* subcomponents! Finished initializing ice model at 20220718 170005.481 20220718 170005.481: Memuse(MB) at ice_model_init= 2.494E+02 2.594E+02 2.635E+00 2.526E+02 Finished initializing component models at 20220718 170005.487 Starting to initialize flux_exchange at 20220718 170005.488 ==>Note from ocean_tpm_util_mod(otpm_set_tracer_package)[ocmip2_cfc_mod(ocmip2_ cfc_flux_init)]: Processing tracer package ocmip2_cfc ==>Note from ocmip2_cfc_mod(ocmip2_cfc_flux_init): No instances ==>Note from ocean_tpm_util_mod(otpm_set_tracer_package)[ocmip2_he_mod(ocmip2_h e_flux_init)]: Processing tracer package ocmip2_he ==>Note from ocmip2_he_mod(ocmip2_he_flux_init): No instances ==>Note from ocean_tpm_util_mod(otpm_set_tracer_package)[ocean_pert_co2_mod(oce an_pert_co2_flux_init)]: Processing tracer package ocean_pert_co2 ==>Note from ocean_pert_co2_mod(ocean_pert_co2_flux_init): No instances ==>Note from ocean_tpm_util_mod(otpm_set_tracer_package)[ocmip2_abiotic_mod(ocm ip2_abiotic_flux_init)]: Processing tracer package ocmip2_abiotic ==>Note from ocmip2_abiotic_mod(ocmip2_abiotic_flux_init): No instances ==>Note from ocean_tpm_util_mod(otpm_set_tracer_package)[ocmip2_biotic_mod(ocmi p2_biotic_flux_init)]: Processing tracer package ocmip2_biotic ==>Note from ocmip2_biotic_mod(ocmip2_biotic_flux_init): No instances ==>Note from ocean_tpm_util_mod(otpm_set_tracer_package)[ocean_bgc_restore_mod( ocean_bgc_restore_flux_init)]: Processing tracer package ocean_bgc_restore ==>Note from ocean_bgc_restore_mod(ocean_bgc_restore_flux_init): No instances ==>Note from ocean_tpm_util_mod(otpm_set_tracer_package)[ocean_ibgc_mod(ocean_i bgc_flux_init)]: Processing tracer package ocean_ibgc ==>Note from ocean_ibgc_mod(ocean_ibgc_flux_init): No instances &GENERIC_TRACER_NML DO_GENERIC_TRACER = T, DO_GENERIC_CFC = T, DO_GENERIC_TOPAZ = F, DO_GENERIC_ERGOM = F, DO_GENERIC_BLING = F, DO_GENERIC_MINIBLING = F, DO_GENERIC_COBALT = F / ==>Note from atmos_ocean_fluxes_mod(aof_set_coupler_flux): Processing coupler fluxes cfc_11_flux ==>Note from atmos_ocean_fluxes_mod(aof_set_coupler_flux): Processing coupler fluxes cfc_12_flux ==>Note from atmos_ocean_fluxes_mod(atmos_ocean_fluxes_init): Processing 2 gas fluxes Dumping fluxes tracer tree fluxes/ cfc_11_flux/ flux_type = 'air_sea_gas_flux_generic' implementation = 'ocmip2' atm_tr_index = 0 mol_wt = 0.000000000000000E+000 ice_restart_file = 'ice_ocmip2_cfc.res.nc' ocean_restart_file = 'ocmip2_cfc_airsea_flux.res.nc' param[1] = 9.360000000000000E-007 param[2] = 9.756100000000000E-006 flag = NULL flux-units = 'mol/m^2/s' flux-long_name = 'Surface flux' deltap-units = 'uatm' deltap-long_name = 'Ocean-air delta pressure' kw-units = 'm/s' kw-long_name = 'Piston velocity' pcair-units = 'mol/mol' pcair-long_name = 'Atmospheric concentration' u10-units = 'm/s' u10-long_name = 'Wind speed at 10 m' psurf-units = 'Pa' psurf-long_name = 'Surface atmospheric pressure' alpha-units = 'mol/m^3/atm' alpha-long_name = 'Solubility w.r.t. atmosphere' csurf-units = 'mol/m^3' csurf-long_name = 'Ocean concentration' sc_no-units = 'dimensionless' sc_no-long_name = 'Schmidt number' cfc_12_flux/ flux_type = 'air_sea_gas_flux_generic' implementation = 'ocmip2' atm_tr_index = 0 mol_wt = 0.000000000000000E+000 ice_restart_file = 'ice_ocmip2_cfc.res.nc' ocean_restart_file = 'ocmip2_cfc_airsea_flux.res.nc' param[1] = 9.360000000000000E-007 param[2] = 9.756100000000000E-006 flag = NULL flux-units = 'mol/m^2/s' flux-long_name = 'Surface flux' deltap-units = 'uatm' deltap-long_name = 'Ocean-air delta pressure' kw-units = 'm/s' kw-long_name = 'Piston velocity' pcair-units = 'mol/mol' pcair-long_name = 'Atmospheric concentration' u10-units = 'm/s' u10-long_name = 'Wind speed at 10 m' psurf-units = 'Pa' psurf-long_name = 'Surface atmospheric pressure' alpha-units = 'mol/m^3/atm' alpha-long_name = 'Solubility w.r.t. atmosphere' csurf-units = 'mol/m^3' csurf-long_name = 'Ocean concentration' sc_no-units = 'dimensionless' sc_no-long_name = 'Schmidt number' ==>Note from flux_exchange_mod(flux_exchange_init): Number of exchanged tracers = 1 Tracer field name :sphum NOTE from PE 0: flux_exchange_mod: tracer "co2" not present in the atmosphere NOTE from xgrid_mod: use second_order conservative exchange NOTE from PE 0: xgrid_mod: reading exchange grid information from mosaic grid file Checked data is array of constant 1 ATM(ATMOCNLND)= 0.999999999965637 0.999999999965637 0.999999999965637 OCN(ATMOCNLND)= 0.703873657789463 0.703873657789463 0.703873657789463 LND(ATMOCNLND)= 0.296126342176174 0.296126342176174 0.296126342176174 Checked data is array of random number between 0 and 1 using second_order ATM(ATMOCNLND)= 0.520351679192166 0.520351679193352 0.520351679193352 OCN(ATMOCNLND)= 0.351937636809641 0.351937636809641 0.351937636812393 LND(ATMOCNLND)= 0.145579763294605 0.145579763294605 0.145579763293996 NOTE from PE 0: xgrid_mod: reading exchange grid information from mosaic grid file Checked data is array of constant 1 LND(LNDOCN)= 0.703873657789463 0.703873657789463 0.703873657789463 OCN(LNDOCN)= 0.703873657789463 0.703873657789463 0.703873657789463 FATAL from PE 15: ==>Error from coupler_types_mod (CT_spawn_1d_3d): Disordered k-dimension index bound list 1 0 FATAL from PE 22: ==>Error from coupler_types_mod (CT_spawn_1d_3d): Disordered k-dimension index bound list 1 0 FATAL from PE 16: ==>Error from coupler_types_mod (CT_spawn_1d_3d): Disordered k-dimension index bound list 1 0 FATAL from PE 18: ==>Error from coupler_types_mod (CT_spawn_1d_3d): Disordered k-dimension index bound list 1 0 FATAL from PE 20: ==>Error from coupler_types_mod (CT_spawn_1d_3d): Disordered k-dimension index bound list 1 0 FATAL from PE 21: ==>Error from coupler_types_mod (CT_spawn_1d_3d): Disordered k-dimension index bound list 1 0 FATAL from PE 17: ==>Error from coupler_types_mod (CT_spawn_1d_3d): Disordered k-dimension index bound list 1 0 FATAL from PE 19: ==>Error from coupler_types_mod (CT_spawn_1d_3d): Disordered k-dimension index bound list 1 0 FATAL from PE 26: ==>Error from coupler_types_mod (CT_spawn_1d_3d): Disordered k-dimension index bound list 1 0 FATAL from PE 28: ==>Error from coupler_types_mod (CT_spawn_1d_3d): Disordered k-dimension index bound list 1 0 FATAL from PE 37: ==>Error from coupler_types_mod (CT_spawn_1d_3d): Disordered k-dimension index bound list 1 0 FATAL from PE 38: ==>Error from coupler_types_mod (CT_spawn_1d_3d): Disordered k-dimension index bound list 1 0 FATAL from PE 39: ==>Error from coupler_types_mod (CT_spawn_1d_3d): Disordered k-dimension index bound list 1 0 FATAL from PE 40: ==>Error from coupler_types_mod (CT_spawn_1d_3d): Disordered k-dimension index bound list 1 0 FATAL from PE 24: ==>Error from coupler_types_mod (CT_spawn_1d_3d): Disordered k-dimension index bound list 1 0 FATAL from PE 25: ==>Error from coupler_types_mod (CT_spawn_1d_3d): Disordered k-dimension index bound list 1 0 FATAL from PE 27: ==>Error from coupler_types_mod (CT_spawn_1d_3d): Disordered k-dimension index bound list 1 0 FATAL from PE 29: ==>Error from coupler_types_mod (CT_spawn_1d_3d): Disordered k-dimension index bound list 1 0 FATAL from PE 30: ==>Error from coupler_types_mod (CT_spawn_1d_3d): Disordered k-dimension index bound list 1 0 FATAL from PE 31: ==>Error from coupler_types_mod (CT_spawn_1d_3d): Disordered k-dimension index bound list 1 0 FATAL from PE 32: ==>Error from coupler_types_mod (CT_spawn_1d_3d): Disordered k-dimension index bound list 1 0 FATAL from PE 33: ==>Error from coupler_types_mod (CT_spawn_1d_3d): Disordered k-dimension index bound list 1 0 FATAL from PE 34: ==>Error from coupler_types_mod (CT_spawn_1d_3d): Disordered k-dimension index bound list 1 0 FATAL from PE 35: ==>Error from coupler_types_mod (CT_spawn_1d_3d): Disordered k-dimension index bound list 1 0 FATAL from PE 41: ==>Error from coupler_types_mod (CT_spawn_1d_3d): Disordered k-dimension index bound list 1 0 FATAL from PE 42: ==>Error from coupler_types_mod (CT_spawn_1d_3d): Disordered k-dimension index bound list 1 0 FATAL from PE 43: ==>Error from coupler_types_mod (CT_spawn_1d_3d): Disordered k-dimension index bound list 1 0 FATAL from PE 44: ==>Error from coupler_types_mod (CT_spawn_1d_3d): Disordered k-dimension index bound list 1 0 FATAL from PE 36: ==>Error from coupler_types_mod (CT_spawn_1d_3d): Disordered k-dimension index bound list 1 0 FATAL from PE 23: ==>Error from coupler_types_mod (CT_spawn_1d_3d): Disordered k-dimension index bound list 1 0 Image PC Routine Line Source fms_CM2M.x 0000000003A66376 Unknown Unknown Unknown fms_CM2M.x 00000000039D94A3 mpp_mod_mp_mpp_er 75 mpp_util_mpi.inc fms_CM2M.x 000000000348C6A9 coupler_types_mod 1344 coupler_types.F90 fms_CM2M.x 000000000348C097 coupler_types_mod 1000 coupler_types.F90 fms_CM2M.x 00000000004791EA flux_exchange_mod 1067 flux_exchange.F90 fms_CM2M.x 0000000000482296 coupler_main_IP_c 1362 coupler_main.F90 fms_CM2M.x 000000000047B0F0 MAIN__ 428 coupler_main.F90 fms_CM2M.x 0000000000414BA2 Unknown Unknown Unknown libc-2.17.so 00002AD9A4F463D5 __libc_start_main Unknown Unknown fms_CM2M.x 0000000000414AA9 Unknown Unknown Unknown Image PC Routine Line Source fms_CM2M.x 0000000003A66376 Unknown Unknown Unknown fms_CM2M.x 00000000039D94A3 mpp_mod_mp_mpp_er 75 mpp_util_mpi.inc fms_CM2M.x 000000000348C6A9 coupler_types_mod 1344 coupler_types.F90 fms_CM2M.x 000000000348C097 coupler_types_mod 1000 coupler_types.F90 fms_CM2M.x 00000000004791EA flux_exchange_mod 1067 flux_exchange.F90 fms_CM2M.x 0000000000482296 coupler_main_IP_c 1362 coupler_main.F90 fms_CM2M.x 000000000047B0F0 MAIN__ 428 coupler_main.F90 fms_CM2M.x 0000000000414BA2 Unknown Unknown Unknown libc-2.17.so 00002B3FCDD5B3D5 __libc_start_main Unknown Unknown fms_CM2M.x 0000000000414AA9 Unknown Unknown Unknown Image PC Routine Line Source fms_CM2M.x 0000000003A66376 Unknown Unknown Unknown fms_CM2M.x 00000000039D94A3 mpp_mod_mp_mpp_er 75 mpp_util_mpi.inc fms_CM2M.x 000000000348C6A9 coupler_types_mod 1344 coupler_types.F90 fms_CM2M.x 000000000348C097 coupler_types_mod 1000 coupler_types.F90 fms_CM2M.x 00000000004791EA flux_exchange_mod 1067 flux_exchange.F90 fms_CM2M.x 0000000000482296 coupler_main_IP_c 1362 coupler_main.F90 fms_CM2M.x 000000000047B0F0 MAIN__ 428 coupler_main.F90 fms_CM2M.x 0000000000414BA2 Unknown Unknown Unknown libc-2.17.so 00002ADCFD9913D5 __libc_start_main Unknown Unknown fms_CM2M.x 0000000000414AA9 Unknown Unknown Unknown Image PC Routine Line Source fms_CM2M.x 0000000003A66376 Unknown Unknown Unknown fms_CM2M.x 00000000039D94A3 mpp_mod_mp_mpp_er 75 mpp_util_mpi.inc fms_CM2M.x 000000000348C6A9 coupler_types_mod 1344 coupler_types.F90 fms_CM2M.x 000000000348C097 coupler_types_mod 1000 coupler_types.F90 fms_CM2M.x 00000000004791EA flux_exchange_mod 1067 flux_exchange.F90 fms_CM2M.x 0000000000482296 coupler_main_IP_c 1362 coupler_main.F90 fms_CM2M.x 000000000047B0F0 MAIN__ 428 coupler_main.F90 fms_CM2M.x 0000000000414BA2 Unknown Unknown Unknown libc-2.17.so 00002AF49631D3D5 __libc_start_main Unknown Unknown fms_CM2M.x 0000000000414AA9 Unknown Unknown Unknown Image PC Routine Line Source fms_CM2M.x 0000000003A66376 Unknown Unknown Unknown fms_CM2M.x 00000000039D94A3 mpp_mod_mp_mpp_er 75 mpp_util_mpi.inc fms_CM2M.x 000000000348C6A9 coupler_types_mod 1344 coupler_types.F90 fms_CM2M.x 000000000348C097 coupler_types_mod 1000 coupler_types.F90 fms_CM2M.x 00000000004791EA flux_exchange_mod 1067 flux_exchange.F90 fms_CM2M.x 0000000000482296 coupler_main_IP_c 1362 coupler_main.F90 fms_CM2M.x 000000000047B0F0 MAIN__ 428 coupler_main.F90 fms_CM2M.x 0000000000414BA2 Unknown Unknown Unknown libc-2.17.so 00002B330E2693D5 __libc_start_main Unknown Unknown fms_CM2M.x 0000000000414AA9 Unknown Unknown Unknown Image PC Routine Line Source fms_CM2M.x 0000000003A66376 Unknown Unknown Unknown fms_CM2M.x 00000000039D94A3 mpp_mod_mp_mpp_er 75 mpp_util_mpi.inc fms_CM2M.x 000000000348C6A9 coupler_types_mod 1344 coupler_types.F90 fms_CM2M.x 000000000348C097 coupler_types_mod 1000 coupler_types.F90 fms_CM2M.x 00000000004791EA flux_exchange_mod 1067 flux_exchange.F90 fms_CM2M.x 0000000000482296 coupler_main_IP_c 1362 coupler_main.F90 fms_CM2M.x 000000000047B0F0 MAIN__ 428 coupler_main.F90 fms_CM2M.x 0000000000414BA2 Unknown Unknown Unknown libc-2.17.so 00002B485B7EE3D5 __libc_start_main Unknown Unknown fms_CM2M.x 0000000000414AA9 Unknown Unknown Unknown Image PC Routine Line Source fms_CM2M.x 0000000003A66376 Unknown Unknown Unknown fms_CM2M.x 00000000039D94A3 mpp_mod_mp_mpp_er 75 mpp_util_mpi.inc fms_CM2M.x 000000000348C6A9 coupler_types_mod 1344 coupler_types.F90 fms_CM2M.x 000000000348C097 coupler_types_mod 1000 coupler_types.F90 fms_CM2M.x 00000000004791EA flux_exchange_mod 1067 flux_exchange.F90 fms_CM2M.x 0000000000482296 coupler_main_IP_c 1362 coupler_main.F90 fms_CM2M.x 000000000047B0F0 MAIN__ 428 coupler_main.F90 fms_CM2M.x 0000000000414BA2 Unknown Unknown Unknown libc-2.17.so 00002B30380F93D5 __libc_start_main Unknown Unknown fms_CM2M.x 0000000000414AA9 Unknown Unknown Unknown Abort(1) on node 16 (rank 16 in comm 0): application called MPI_Abort(MPI_COMM_WORLD, 1) - process 16 Image PC Routine Line Source fms_CM2M.x 0000000003A66376 Unknown Unknown Unknown fms_CM2M.x 00000000039D94A3 mpp_mod_mp_mpp_er 75 mpp_util_mpi.inc fms_CM2M.x 000000000348C6A9 coupler_types_mod 1344 coupler_types.F90 fms_CM2M.x 000000000348C097 coupler_types_mod 1000 coupler_types.F90 fms_CM2M.x 00000000004791EA flux_exchange_mod 1067 flux_exchange.F90 fms_CM2M.x 0000000000482296 coupler_main_IP_c 1362 coupler_main.F90 fms_CM2M.x 000000000047B0F0 MAIN__ 428 coupler_main.F90 fms_CM2M.x 0000000000414BA2 Unknown Unknown Unknown libc-2.17.so 00002AC217E233D5 __libc_start_main Unknown Unknown fms_CM2M.x 0000000000414AA9 Unknown Unknown Unknown Abort(1) on node 15 (rank 15 in comm 0): application called MPI_Abort(MPI_COMM_WORLD, 1) - process 15 Abort(1) on node 18 (rank 18 in comm 0): application called MPI_Abort(MPI_COMM_WORLD, 1) - process 18 Abort(1) on node 20 (rank 20 in comm 0): application called MPI_Abort(MPI_COMM_WORLD, 1) - process 20 Image PC Routine Line Source fms_CM2M.x 0000000003A66376 Unknown Unknown Unknown fms_CM2M.x 00000000039D94A3 mpp_mod_mp_mpp_er 75 mpp_util_mpi.inc fms_CM2M.x 000000000348C6A9 coupler_types_mod 1344 coupler_types.F90 fms_CM2M.x 000000000348C097 coupler_types_mod 1000 coupler_types.F90 fms_CM2M.x 00000000004791EA flux_exchange_mod 1067 flux_exchange.F90 fms_CM2M.x 0000000000482296 coupler_main_IP_c 1362 coupler_main.F90 fms_CM2M.x 000000000047B0F0 MAIN__ 428 coupler_main.F90 fms_CM2M.x 0000000000414BA2 Unknown Unknown Unknown libc-2.17.so 00002B2D87C303D5 __libc_start_main Unknown Unknown fms_CM2M.x 0000000000414AA9 Unknown Unknown Unknown Abort(1) on node 21 (rank 21 in comm 0): application called MPI_Abort(MPI_COMM_WORLD, 1) - process 21 Abort(1) on node 37 (rank 37 in comm 0): application called MPI_Abort(MPI_COMM_WORLD, 1) - process 37 Abort(1) on node 40 (rank 40 in comm 0): application called MPI_Abort(MPI_COMM_WORLD, 1) - process 40 Abort(1) on node 41 (rank 41 in comm 0): application called MPI_Abort(MPI_COMM_WORLD, 1) - process 41 Abort(1) on node 43 (rank 43 in comm 0): application called MPI_Abort(MPI_COMM_WORLD, 1) - process 43 ERROR: Model failed to run to completion