Query about need to change in Cgrid_transformation_tools.py

croco tools (pre/post-processing)
1

Hi,
In mostly python base croco_pytools have a one file called Cgrid_transformation_tools.py inside the prepro/Modules . During reading the code found the Quadratic output shows very suspicious. Its very confused me.
can you please check:
line (519-524)

    var_rho[iz, :, :] = 0.5625 * (
        var_w[iz + 1, :, :] + var_w[iz, :, :]
    ) - 0.0625 * (var_w[iz + 2 :, :] + var_w[iz - 1, :, :])
var_rho[0, :, :] = (
    -0.125 * var_w[1, :, :] + 0.75 * var_w[1, :, :] + 0.375 * var_w[0, :, :]
)

maybe need to change:

    var_rho[iz, :, :] = 0.5625 * (
        var_w[iz + 1, :, :] + var_w[iz, :, :]
    ) - 0.0625 * (var_w[iz + 2, :, :] + var_w[iz - 1, :, :])
var_rho[0, :, :] = (
    -0.125 * var_w[2, :, :] + 0.75 * var_w[1, :, :] + 0.375 * var_w[0, :, :]
)

I added only one comma (,) and change one 1 to 2 two. you can check it above. Is it correct ?
can you please check it. then I also change (very confusing for me).
please let me know if that changes needed or not.

Thank you for your kind understanding and help.
Best regards.

2

Thanks for the report and for pointing this out.

This issue has already been raised on the CROCO forum here: Bug in w2rho() loop body / croco_pytools v2.0.4

It is currently an open work item in our GitLab repository. It has not been addressed yet, but it is already identified and will be investigated.

Regards,

3

Thank you for your kind responses
Best regards