Issue with PSOURCE/PSOURCE_NCFILE/PSOURCE_NCFILE_TS

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 croco.in, 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

add a comment

0

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

GildasC
26

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 croco-ocean.org 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

add a comment

0

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

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 croco-users@inria.fr last week, apologies for the delay updating the forum.

Best, Joe

edit flag offensive delete link

add a comment

0

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

MartinV
41

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 roms.in file and n_qbar variable in the runoff netcdf file. Best, Martin

edit flag offensive delete link

add a comment

0

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

Ansong
1

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