| ... | @@ -69,3 +69,35 @@ first runs did not work because of something in the newly made river input files |
... | @@ -69,3 +69,35 @@ first runs did not work because of something in the newly made river input files |
|
|
<img src="uploads/96e46a2d5ce64abc58705d6772fa9d65/image.png" width=49%>
|
|
<img src="uploads/96e46a2d5ce64abc58705d6772fa9d65/image.png" width=49%>
|
|
|
<img src="uploads/76151e2d9390d4ae42d5de21a2fe9f07/image.png" width=49%>
|
|
<img src="uploads/76151e2d9390d4ae42d5de21a2fe9f07/image.png" width=49%>
|
|
|
* Andrew gave another go with the last results in mind to improve the results
|
|
* Andrew gave another go with the last results in mind to improve the results
|
|
|
|
This Time the automation has been run on both DIN/SIL and DIC/ALK riverine input sets.
|
|
|
|
It seems to runs better, the model runs for 10years before failing. The run use the `fct` advection scheme for the passive tracers (same all physics) which is the one having the most chance of failing - just to avoid having to retest everything again. and this set of files :
|
|
|
|
```fortran
|
|
|
|
&nammsa_river
|
|
|
|
ln_read_rivn = .true.
|
|
|
|
ln_read_rivsi = .true.
|
|
|
|
ln_read_rivc = .true.
|
|
|
|
ln_read_rivalk = .true.
|
|
|
|
rivn_kind = 2
|
|
|
|
rivsi_kind = 2
|
|
|
|
rivc_kind = 2
|
|
|
|
rivalk_kind = 2
|
|
|
|
rivdepmax = 70.0
|
|
|
|
cn_dir = './INPUT/INIT_MEDUSA/' ! root directory for the dust data location
|
|
|
|
!___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________!
|
|
|
|
! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights filename ! rotation ! land/sea mask !
|
|
|
|
! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! ! pairing ! filename !
|
|
|
|
sn_riv_n = 'river42_tot_N_flux_LS200km_SP10_5ITER_ORCA025.nc' , -12 , 'RIV_N' , .true. , .true. , 'yearly' , '' , '' , ''
|
|
|
|
sn_riv_s = 'river42_tot_Si_flux_LS200km_SP10_5ITER_ORCA025.nc' , -12 , 'RIV_Si' , .true. , .true. , 'yearly' , '' , '' , ''
|
|
|
|
sn_riv_c = 'river42_total_DIC_flux_LS200km_SP10_5ITER_ORCA025.nc', -12 , 'RIV_C' , .true. , .true. , 'yearly' , '' , '' , ''
|
|
|
|
sn_riv_a = 'river42_total_ALK_flux_LS200km_SP10_5ITER_ORCA025.nc', -12 , 'RIV_ALK' , .true. , .true. , 'yearly' , '' , '' , ''
|
|
|
|
/
|
|
|
|
```
|
|
|
|
The error is no more in the Arctic rivers.
|
|
|
|
It still is a river problem, but around china/korea. 1 pixel switch to crazy DIC :
|
|
|
|
<img src="uploads/7ecaf7642c7d63e863e79882d8166cca/image.png" width=49%>
|
|
|
|
<img src="uploads/758795535d9eb6a77df22fd2a3d586db/image.png" width=49%>
|
|
|
|
Because of the new C limitation, we shouldn't have anything quite like that happening...
|
|
|
|
* A quick test could be to add the nutrient down to a shallower depth of 10m instead of 70.
|
|
|
|
* we could also strengthen the C security we've added in the C uptake...
|
|
|
|
I quickly check by changing the river depth. we'll see.
|
|
|
|
With river, the total production in 1986 is of 31 Pg-C per y still too low for our liking. this will need some tweaking. |
