Ask Your Question


asked 2022-02-03 12:45:54 +0200

Hi there,

I'm having difficulty getting point sources to work for my model. I have followed the tutorial instructions to the letter, but I get the following runtime issue, whether using constant flow specified via, or time varying flow via the runoff netcdf file.

forrtl: severe (408): fort: (2): Subscript #1 of the array ISRC has value 1 which is greater than the upper bound of 0

I can't be the first person to encounter such a problem, but I see no reference to the problem on the CROCO forum or in the mailing list communications.

Your advice would be greatly appreciated.

Best regards, Joe

edit retag flag offensive close merge delete

4 answers

Sort by ยป oldest newest most voted

answered 2022-03-30 14:06:02 +0200

GildasC gravatar image

Hi all,

The problem was due to a bug in the v1.2 It has been corrected into the v1.2.1 newly released.It is available on Here the change log for the v1.2.1

  • ONLINE bulk: wind speed was not properly computed, with a consequent underestimation of fluxes when using bulk ONLINE. This is now corrected

  • Runoff : Msrc parameter in param.h was wrongly set to 0, it is now corrected and set to 100

  • Plurimonth_scripts: DT was wrongly intialized with nesting. This is now corrected.

  • In 2D cases, the size of wrthis array was wrong. This is now corrected. A typo in def_diags has been corrected

  • CVTK: some continuous integration tests have been updated

Regards, Gildas

edit flag offensive delete link more

answered 2022-02-08 16:21:53 +0200

this post is marked as community wiki

This post is a wiki. Anyone with karma >5 is welcome to improve it.

Just to close this off, I found that this issue was linked to Msrc in param.h file which is referred to in the Ifremer specific section of the tutorials document for v1.2, whereas I was only looking at section 15 of the document and matlab scripts.

I sent out the same query and the corresponding solution via last week, apologies for the delay updating the forum.

Best, Joe

edit flag offensive delete link more

answered 2022-02-08 16:16:48 +0200

MartinV gravatar image

updated 2022-02-08 16:17:15 +0200

Regarding the ISRC problem, my first guess would be to check the number of sources in your file and n_qbar variable in the runoff netcdf file. Best, Martin

edit flag offensive delete link more

answered 2022-02-05 21:05:44 +0200

Ansong gravatar image

I've also been unsuccessful trying to get PSOURCE/ANA_PSOURCE as well as PSOURCE/PSOURCE_NCFILE to work for the BENGUELA test case. The code runs well but I don't get any river flow when I plot the results. But the RIVER test case works fine. I've checked to make sure the location of the point sources are on land points, and also used the make_runoff.m code to create a netcdf file. But to no avail.

Has anyone been successful with the PSOURCE setup for the BENGUELA test case? Please share your setup with us. Thanks.

Joseph Ansong

edit flag offensive delete link more
Login/Signup to Answer

Question Tools

1 follower


Asked: 2022-02-03 12:45:54 +0200

Seen: 334 times

Last updated: Mar 30 '22