MARBL ERROR: abs(Jint_Fetot) exceeds Jint_Fetot_thres

[Zhu Zhu]

Hi all,
I am running experiment using fully coupled CESM2.1.5. The model runs fine with the default B1850G configuration. I want to conduct an experiment with prescribed Antarctic sea ice (2000-2100). I changed the Antarctic sea ice and associated fluxes passed to coupler and run with B1850G. After running for about 30 years, the model crashed with an error:

(Task 191, block 1) Message from (lon, lat) ( 310.438, 1.199), which is global (i,j) (312, 192). Level: 3
(Task 191, block 1) MARBL WARNING (marbl_interior_tendency_mod:compute_large_detritus_prod): dz*DOP_loss_P_bal= 0.425E-005 exceeds Jint_Ptot_thres= 0.271E-013
(Task 191, block 1) Message from (lon, lat) ( 310.438, 1.199), which is global (i,j) (312, 192). Level: 3
(Task 191, block 1) MARBL WARNING (marbl_interior_tendency_mod:compute_large_detritus_prod): dz*DOP_loss_P_bal= 0.162E-005 exceeds Jint_Ptot_thres= 0.271E-013
max rss=2302.8 MB
max rss=2302.8 MB
max rss=2302.8 MB
(Task 31, block 1) Message from (lon, lat) ( 301.438, -63.194), which is global (i,j) (304, 31). Level: 1
(Task 31, block 1) MARBL ERROR (marbl_diagnostics_mod:store_diagnostics_iron_fluxes): abs(Jint_Fetot)= 0.954E-017 exceeds Jint_Fetot_thres= 0.951E-017
(Task 31, block 1) MARBL ERROR (marbl_diagnostics_mod:marbl_diagnostics_interior_tendency_compute): Error reported from store_diagnostics_iron_fluxes
(Task 31, block 1) MARBL ERROR (marbl_interior_tendency_mod:marbl_interior_tendency_compute): Error reported from marbl_diagnostics_interior_tendency_compute()
(Task 31, block 1) MARBL ERROR (marbl_interface:interior_tendency_compute): Error reported from marbl_interior_tendency_compute()
(Task 31, block 1) MARBL ERROR (ecosys_driver:ecosys_driver_set_interior): Error reported from marbl_instances(1)%set_interior_forcing()
ERROR reported from MARBL library

Then I changed the Jint_Fetot_thres to 0.951E-015 and continue run the simulation, but it exceeds the new threshhold after about 20 years:

(Task 201, block 1) Message from (lon, lat) ( 181.063, 3.895), which is global (i,j) (197, 202). Level: 15
(Task 201, block 1) MARBL WARNING (marbl_interior_tendency_mod:compute_large_detritus_prod): dz*DOP_loss_P_bal= 0.406E-013 exceeds Jint_Ptot_thres= 0.271E-013
(Task 238, block 1) Message from (lon, lat) ( 292.448, 12.231), which is global (i,j) (296, 234). Level: 14
(Task 238, block 1) MARBL WARNING (marbl_interior_tendency_mod:compute_large_detritus_prod): dz*DOP_loss_P_bal= 0.239E-011 exceeds Jint_Ptot_thres= 0.271E-013
(Task 27, block 1) Message from (lon, lat) ( 219.313, -66.399), which is global (i,j) (231, 25). Level: 1
(Task 27, block 1) MARBL ERROR (marbl_diagnostics_mod:store_diagnostics_iron_fluxes): abs(Jint_Fetot)= 0.711E-014 exceeds Jint_Fetot_thres= 0.951E-015
(Task 27, block 1) MARBL ERROR (marbl_diagnostics_mod:marbl_diagnostics_interior_tendency_compute): Error reported from store_diagnostics_iron_fluxes
(Task 27, block 1) MARBL ERROR (marbl_interior_tendency_mod:marbl_interior_tendency_compute): Error reported from marbl_diagnostics_interior_tendency_compute()
(Task 27, block 1) MARBL ERROR (marbl_interface:interior_tendency_compute): Error reported from marbl_interior_tendency_compute()
(Task 27, block 1) MARBL ERROR (ecosys_driver:ecosys_driver_set_interior): Error reported from marbl_instances(1)%set_interior_forcing()
ERROR reported from MARBL library

I am not sure why this error occurs, since I did not change any variables related to iron. Would it be reasonable to relax Jint_Fetot_thres further (e.g., to 0.951E-013) in order to continue the run? Any suggestions would be very helpful. Thanks a lot!

Best regards,
Zhu Zhu