Issue connected to Tnudg_SSS

[c.drinkorn]

I am trying to set up Kate's branch for an Arctic application including ice and sediments. When I first used the code as I received it from github, there was a compilation error about the dimension of Rval in inp_par.F Obviously, there is a major difference between the SVN version and Kate's branch regarding the handling of the keyword value loading from the .in files. I found a mismatch of dimension between decode_line in inp_par.F and read_sedpar.F (sediment.inp.h) I changed the dimension value to 100 in the decode_line function and the compilation went smooth. However, since then I am receiving a segmentation fault regarding the handling of Tnudg_SSS. I had to add this value to my .in file first since it was missing there. (The .in file is from my previous application using the SVN code) However, this didn't solve the problem. Here is an excerpt of the error backtrace output:

Code: Select all

==== backtrace ====
 2 0x00000000000575cc mxm_handle_error()  /scrap/jenkins/workspace/hpc-power-pack/label/r-vmb-rhel6-u7-x86-64-MOFED-CHECKER/hpcx_root/src/hpcx-v1.8.0-gcc-OFED-3.18-redhat6.7-x86_64/mxm-v3.6/src/mxm/util/debug/debug.c:641
 3 0x000000000005773c mxm_error_signal_handler()  /scrap/jenkins/workspace/hpc-power-pack/label/r-vmb-rhel6-u7-x86-64-MOFED-CHECKER/hpcx_root/src/hpcx-v1.8.0-gcc-OFED-3.18-redhat6.7-x86_64/mxm-v3.6/src/mxm/util/debug/debug.c:616
 4 0x0000000000032510 killpg()  ??:0
 5 0x00000000024d48c2 load_r()  /mnt/lustre01/scratch/b/b380636/Arctic20km/obcatmforceriversinisedbulkice/Build/inp_par.f90:1547
 6 0x0000000002711e80 read_phypar()  /mnt/lustre01/scratch/b/b380636/Arctic20km/obcatmforceriversinisedbulkice/Build/read_phypar.f90:674
 7 0x00000000024f0b87 inp_par()  /mnt/lustre01/scratch/b/b380636/Arctic20km/obcatmforceriversinisedbulkice/Build/inp_par.f90:94
 8 0x000000000040cffe roms_initialize()  /mnt/lustre01/scratch/b/b380636/Arctic20km/obcatmforceriversinisedbulkice/Build/ocean_control.f90:86
 9 0x000000000040ca86 ocean()  /mnt/lustre01/scratch/b/b380636/Arctic20km/obcatmforceriversinisedbulkice/Build/master.f90:95
10 0x000000000040c85e main()  ??:0
11 0x000000000001ed1d __libc_start_main()  ??:0
12 0x000000000040c769 _start()  ??:0

In load_r the line referred here is the one where Vout(i) is set to Vinp(i). Maybe there is still some dimension mismatch?
I could provide a lot of output and configuration here but since I don't know exactly where to look, I'd rather provide these when asked for them in particular.

Also, I am not so familiar with the ice code and really want to keep it simple. I haven't exactly figured, yet what the most essential ice simulation setting could be. I also tried to implement the ice module into the SVN version of ROMS but this didn't really mainly also due to the above reasons (lack of experience with the ice code - no clue what the essential ice code is).

Any help, hints or advice will be highly appreciated!

[kate]

If you are using my code, you really, really need to start from one of my ocean.in files. In particular, you need not only NFFILES for the number of forcing files, but also NBCFILES and NCLMFILES for the number of boundary and climatology files as well. Maybe even NSSFFILES for the rivers if the code is recent enough.

I'm afraid I've gotten behind on the Hernan updates. The latest COAWST code has both sea ice and the Hernan updates, not to mention active sediment users.

[c.drinkorn]

If you are using my code, you really, really need to start from one of my ocean.in files. In particular, you need not only NFFILES for the number of forcing files, but also NBCFILES and NCLMFILES for the number of boundary and climatology files as well. Maybe even NSSFFILES for the rivers if the code is recent enough.

I'm afraid I've gotten behind on the Hernan updates. The latest COAWST code has both sea ice and the Hernan updates, not to mention active sediment users.

Hi Kate, thank you for the hint about the COAWST code!
I will have a look into my .in file and compare it with the ocean.in provided in the User/External section of your branch. However, I did already add the number of files variables. I think, I remembered this from another post and added these right away. So this can't be the issue. I will try to find the problem as long as this day lasts but then check out COAWST.
For a better understanding: Why is ROMS calling the nudging time scale for sea surface salinity at all? For now I have closed all ice boundaries and use an analytical ice initial. Where and why does it want to nudge? I am not providing any climatologies....

[c.drinkorn]

Unfortunately, I am reproducing the error with one of your ocean.in files. I took my time to look into the referred code lines but there is nothing suspicious or inconsistent. I guess, I'll switch to the COAWST code, then, unless anyone has an idea what produces this problem.