Here, the salinity equation is initialised. First, the namelist salt is read from getm.inp. Then, depending on the salt_method, the salinity field is read from a hotstart file (salt_method=0), initialised with a constant value (salt_method=1), initialised and interpolated with horizontally homogeneous salinity from a given salinity profile (salt_method=2), or read in and interpolated from a 3D netCDF field (salt_method=3). Finally, a number of sanity checks are performed for the chosen salinity advection schemes.
Apart from this, there are various options for specific initial conditions which are selected by means of compiler options. USES:
use advection, only: J7 use advection_3d, only: print_adv_settings_3d IMPLICIT NONEINPUT/OUTPUT PARAMETERS:
integer :: i,j,k,n integer :: status NAMELIST /salt/ & salt_method,salt_const,salt_file, & salt_format,salt_name,salt_field_no, & salt_adv_split,salt_adv_hor,salt_adv_ver,salt_AH, & salt_check,min_salt,max_salt