Perturbation output

[lola-cmks]

Hello everyone, I am trying to read LIS_DAPERT binary output files generated by my data assimilation run in Fortran but have been encountering issues. Specifically, I am unsure about the correct way to interpret the binary format, as my attempts have resulted in incorrect or empty data.

What I Have Tried

1. I attempted to read the files using unformatted and stream access with convert='big_endian'.
2. I structured my Fortran code assuming the following data dimensions: (variable, lat, lon, ensemble)
   where: variable: Includes perturbed forcing, model states and observations

Despite this, the data does not appear correctly structured when I inspect the output.

Questions

1. Can someone confirm if the assumed structure of the binary file is correct?
2. How many variables are actually included in the LIS DAPERT output?
3. Is there an existing Fortran script available to correctly read these files?
4. What is the best way to verify the correct byte structure of these files (e.g., precision, ordering)?

Any guidance or example code would be greatly appreciated!

Thanks in advance.