Might not need to interpolate initial conditions

[andrew-c-ross]

Initial conditions are currently interpolated and flooded during the setup process:

CEFI-regional-MOM6/tools/initial/write_bgc_initial.py

Line 22 in 2934bf2

def interpolate_flood(target_grid, ztarget, ds, xdim=None, ydim=None, periodic=True):

This can be time-consuming, especially as the size and resolution of the model grids increases. MOM has the ability to read z/lat/lon files and interpolate on the fly. Using this for the physics and BGC initial conditions should involve setting the MOM_input option INIT_LAYERS_FROM_Z_FILE = True and setting the input data source and some other properties in the field_table, for example:

no3_src_file = INPUT/woa18_all_n00_01.nc
no3_src_var_name = n_an
no3_src_var_unit = micromoles_per_kg
no3_src_var_record = 1
no3_src_var_gridspec = none
no3_valid_min = 0.0

This is what I do for the NWA25 with BGC because of the huge size of the model grid. It should also work for the other grids.

As discussed today with @amoebaliz and @uwagura .