I have some questions about reader_netCDF_CF_unstructured.py code.
If you look at the get_variable function, there is a slight difference between node_variable and face_variable in the part that fetches the data at the closest point to the particle position.
The node_variables says to extract the required data using sigma_ind when dvar.dimensions has 'siglev'.
On the other hand, face_variables use sigma_ind to extract data when 'siglay' is in the variable dimension.
In my opinion, if node_variables, face_variables, or depth dimensions exist, they can all use 'siglay' the same, but what made them set them differently? Is there any special reason?
I think I made this change because siglay was for face_variables and siglev was for node variables, otherwise the dimensions did not match up. But it has been a while since I did this.
I have some questions about reader_netCDF_CF_unstructured.py code.
If you look at the get_variable function, there is a slight difference between node_variable and face_variable in the part that fetches the data at the closest point to the particle position. The node_variables says to extract the required data using sigma_ind when dvar.dimensions has 'siglev'. On the other hand, face_variables use sigma_ind to extract data when 'siglay' is in the variable dimension.
In my opinion, if node_variables, face_variables, or depth dimensions exist, they can all use 'siglay' the same, but what made them set them differently? Is there any special reason?
face_variable interpolating![image](https://github.com/OpenDrift/opendrift/assets/154101128/c56f8a47-ac5a-42b4-a3c4-fdd0d00ce0fa)
node_variable interpolating![image](https://github.com/OpenDrift/opendrift/assets/154101128/d8f84047-547a-4ecb-b931-fb1c941a34ec)