This post describes the changes I made to increase AttSW in the Salish Sea. Note that these changes are in LO_roms_source_alt/npzd_banas/fennel.h. The previous version of fennel.h without my edits is called “fennel_LO_ORIG.h.” You can find both of these files in perigee in /data1/auroral/LO_roms_source_alt/npzd_banas.

Figure 1 shows a mp of the intended changes. Outside of the Salish Sea region, AttSW = 0.05 m-1. Inside of the Salish Sea (ish) region, AttSW = 0.15 m-1.


Fig 1. Map of AttSW values.


The following section enumerates all of the changes I made to fennel.h.

Changes to fennel.h

(1.) Add the grid lon and lat as inputs to the function call of fennel_tile


(2.) Add lonr and latr as expected inputs to the fennel_tile subroutine


(3.) Declare lonr and latr as an array of real numbers that spans the i and j indices of the tile


(4.) Declare a new local variable called ‘AttSW_region’


(5.) If lonr and latr are within the Salish Sea, AttSW_region = AttSW * 3. Otherwise, AttSW_region = AttSW. I also added some print statements for testing, but have since commented them out. Note that AttSW is still defined by the user in bio_Fennel_BLANK.in. The initial value is still 0.05 m-1, and simply gets multiplied by a factor of 3 inside of the Salish Sea.


(6.) Update calculations of light attenuation to use the new regional AttSW_region rather than AttSW


Results

The resulting print statement in the log file demonstrated success:


Fig 2. Print statement output in the log file listing the different AttSW_region values.


Figure 3 shows the corresponding location of these test points on a map.


Fig 3. Location of AttSW test points.