Good morning everyone,
I am facing some problems setting up some transport sections using the diadct.f90 module. My section_ijglobal.diadct file reads
-68. -54.5 -60. -64.7 00 okstrpond noice ACC_Drake_Passage
27. -73. 27. -30. 00 nostrpond noice ACC_LeCap
116. -70. 116. -33.8 00 nostrpond noice ACC_WAUS
147.3 -70. 147.3 -41. 00 okstrpond noice ACC_TAS
167.4 -73. 167.4 -45.5 00 nostrpond noice ACC_NZ
and it is located in the EXP00 folder. This is nothing else than the first five lines of the file /tools/SECTIONS_DIADCT/run/list_sections.ascii_global, since I am doing this in an ORCA1 configuration (this configuration, that is actually based on de Lavergne et al. 2024) I thought this would suffice.
Running NEMO I got
At line 387 of file .../BLD/ppsrc/nemo/diadct.f90 (unit = 20, file = 'section_ijglobal.diadct')
Fortran runtime error: End of file
and the code stops. Inside ocean.output the diadct section reports
diadct_init: compute transports through sections
~~~~~~~~~~~~~~~~~~~~~
Calculate transport thru sections: ln_diadct = T
Frequency of computation: nn_dct = 15
Frequency of write: nn_dctwri = 15
Debug all sections
debuging for section number: 1
===>>> : E R R O R
===========
diadct: read sections : Problem of section number: isec= **** and jsec= 0001
isec 538976302
and I have a bunch of ocean.output_XXXX with the same error.
Does anyone know what I’m doing wrong?
Thanks in advance,
Francesco
PS: I’m doing this in NEMO 5.0.1, this version is apparently not present in the possible tags for this discourse