My aim is to simulate CROCO coupled with PISCES over the northern Indian Ocean. I defined dust input as True and PSOURCE also activated. In Climatology RUN, after the time step of 127 month onwards my DIC values becoming NaN. Here attached the image of croco_avg.00128.nc month map and namelist_pisces_ref also. Any suggestions or comments please.
What version of CROCO are you using? I had a similar issue with the atmospheric dust input in an older version. Perhaps to identify the problem, you can turn off the rivers first and just check the effect of turning the dust input on or off.
The river input can also be a source of DIC values becoming NaN when salinity values near river mouths become very low. For a related post, please see the link below:
I implemented the fix in the code suggested here and it seemed to help.
Did you check that temperature and salinity are ok ?
Maybe you can do a test without rivers ?
You should activate in namelist_pisces_cfg,
ln_ironsed = .true. ! boolean for Fe input from sediments
to get the correct input of Fe.