Warm start crashes FLEXPART-WRF

[adingwell]

FLEXPART-WRF crashes with the following cryptic error if trying to run with option IPOUT=1:

At line 46 of file skplin.f90 (unit = 1, file = 'flex_nyiragongo_main.input')
Fortran runtime error: End of file

I think the problem is in subroutine readinput in this block of code:

...
        if (iouttype .eq. 0 .or. iouttype .eq.2) then
          read(unitpartin,end=99) itimein,numpart_in, &
            iomode_xycoord_in
        else
          read(unitpartin,*,end=99) itimein,numpart_in, &
            iomode_xycoord_in
        endif
...

After jumping back to 99, readinput will repeat code that has already been run. WRF input files will be read a second time, followed by the gridcheck routines. After that, readinput will try to read OUTGRID setting a second time, but since this has already been done, the initial call to skplin will reach the end of file(?).
This much I'm sure about.

I'm not sure if I got the rest right since this subroutine is a bit hard to follow but hopefully this will help solve the problem. I tried solving the problem by skipping forward instead:

...
        if (iouttype .eq. 0 .or. iouttype .eq.2) then
          read(unitpartin,end=101) itimein,numpart_in, &
            iomode_xycoord_in
        else
          read(unitpartin,*,end=101) itimein,numpart_in, &
            iomode_xycoord_in
        endif

101   nparttot2=maxpart+numpart_in

...

The model will now run without any errors, however, the "old" particles aren't properly loaded. Plotting the particle positions before and after the attempted warm start clearly shows that the model started from scratch.

So, is the warm start feature currently broken or not yet implemented?

[jbrioude]

Hi Adam,
sorry for the late reply.
Can you attach your input file before the warm start and when you try the warm start?

Thanks

[adingwell]

Initial input file

[adingwell]

input for warm start

[jbrioude]

bug fix

[jbrioude]

Hi Adam,
I have added a tar file with different fortran codes to fix the warm start issue you have.
Let me know if that works.

[adingwell]

Hi!
We might have come across something bigger here, I still got a crash with the new code (segmentation fault) so I re-compiled with debug flags and got this:

At line 133 of file interpol_wind.f90
Fortran runtime error: Index '0' of dimension 1 of array 'height' below lower bound of 1

[adingwell]

Updated to next milestone, currently no fix.

[pesei]

Ticket retargeted after milestone closed