Rivers in MEDUSA

Rivers in eORCA1

Rivers in eORCA025

Developing the NPD runs and willing to have rivers for the SE-NEMO, we've tested the rivers in MEDUSA-NEMO4.2.2 within the NPD runs (GOSI10)
First, the code itself does what we want. just fixed 3 bugs :

• NPD reads ancillaries in directories (not in the running directory: "./"), and that needed a fix see
• forgot to add the river max depth in the namelist_medusa_ref
• the river nutrient check was in integer instead of REAL

first runs did not work because of something in the newly made river input files.

• tested 1 by 1 seems to come from the N river inputs. but all fields show high values in some specific river mouth, of which some in the Arctic look present very high values (here some example after 5month):
  
  those results are when using the files :
  • river42_total_ALK_flux_LS200km_SP10_OBfix_ORCA025.nc
  • river42_total_DIC_flux_LS200km_SP10_OBfix_ORCA025.nc
  • river42_tot_N_flux_LS200km_SP10_OBfix_ORCA025.nc
  • river42_tot_Si_flux_LS200km_SP10_OBfix_ORCA025.nc
• maybe the we can try with LS 400? here a diff:
  
  Running the LS 400 river nutrients did not solve the problem.
• Andrew added made some special treatment to the Ob fjord in the Russian Arctic. (diff with - without Ob special treatment)
  
  It didn't Solve the problem.
• I've them re-run the last rivers, spreading the nutrients down to 70m deep (max) instead of 10m earlier (see DIN between 2 and 102 mmol/m3 -- max go up to 155mmol/m3 in around the Ob).
  
  Still no Luck.
• Looking at the bathy, spreading the river nutrients down to 70m won't make a difference in the Arctic...
  
  As the sea-floor stands between 8 and 15m deep around there...
• Last hope LS400 with special Ob treatment.
  This test ran longer than the 200km version (7month compared to 2) but still, it eventually failed. same problem around the Ob :
  
  
  And Here LS400 just before it failed
  
  A solution would be to fill (on top of the Ob river mouth) the other problematic fjords, circled in the last plot. The "filling" is only for the river nutrient flux file creation, not at run time. This way, the nutrients would be released in a better mixed location and avoid the accumulation we see in those fjords. Best would be to use the 200km files made this way.
• New Test with more Fjords filled :
  Here is the new modified river input (for N and Si) in the Arctic:
  
  Unfortunately, it failed... again...
  Here are some plots :
  The usual suspects look OK-ish
  
  
  The DIN values are not insane, but this is after 1 month... we might want to fill a tiny bit more in the circled areas
  Further East another beast : This one looks nasty... it's in the middle of probably quite still waters. would be good to fill the suggested part (narrow, and no freshwater).
  
  last bit of the Arctic looks OK
  
• Another tour of improvements,
  Andrew has added an automated way of avoiding narrow/isolated coastal river mouth, to avoid having nutrient accumulating too much
  
  This time the model runs!!! \o/
  But it seems to still accumulate nutrients in some places. Concentration can be quite high (especially DIC get to >7000 mmol/m3 in some places, but the model holds:
  
• Andrew gave another go with the last results in mind to improve the results