Partel issue when nplan=1

I encountered a crash in partel due to a non-allocated variable IKLE3D_P at line 221 of mod_write_solution.f Looking at the code we could saw at line 213 a test "IF(NPLAN.EQ.0)" then line 221 is only reached if NPLAN>0 The problem is at line 151 where there is the test "IF(NPLAN.GT.1)" then IKLE3D_P is allocated So in case NPLAN=1, IKLE3D_P is not allocated but accessed !!!

According to Selafin description, NPLAN value is IPARAM(7) and the description mention : if IPARAM(7) is not 0: the value corresponds to the number of planes on the vertical (in prisms.)

So there is clearly an ambiguous definition of NPLAN and the associated tests and I suspect some border effects using python scripts of not. IMO a mesh file generated by BK as NPLAN=0 but when using python scripts to alter the mesh the generated file has NPLAN=1...

I'm just surprised this kind of problem was not identified earlier ...