PISCES and psource， read in error

biogeochemistry
1

Hi,
I am working with croco-pisces, and there is only a large river in our model domain.
I am using the psource defined in croco_inter.in, and the model runs very successfully without PISCES.
However, the ‘read in error’ occurs when I turned on the pisces.
Does anybody know how to solve this problem quickly please?
Thanks a lot in advance,
Please see in attachment the error message, and the input files.

croco_inter.in
param.h
cppdefs.h
WechatIMG4712.png

2

Hi

With rivers, biogeochemical tracers are considered as any other tracers. So in the psource section of croco_inter.in, you have to specify the same kind of informations than for temperature and salinity :

psource: Nsrc Isrc Jsrc Dsrc Qbar
[m3/s] Lsrc Tsrc
1
61 148 1 29433. T T 20. 10.82

Here the flags T T 20. 10.82 are for T and S, meaning we want to input values of 20. and 10.82 respectively.
So with Pisces you have 25 biogeochemical tracers. Let’s say you don’t want biogeochemical input (just an example), it would lead to :

psource: Nsrc Isrc Jsrc Dsrc Qbar
[m3/s] Lsrc Tsrc
1
61 148 1 29433. T T 25 * F 20. 10.82 25 * 0.

hope it helps
Rachid