Hi!
I tried and failed to add a 5d array in the registry. It should be xyz and two new non gridded dimensions. It seemed to be fine. It was allocated correctly and I could use it in the code, but i never managed to get the output in the history file. After searching around I now came to the conclusion that it is not possible. I can have at most 4D and I have to package them up, as in moist. Correct?
I have now written a script which will add nxm new fields to the registry, where n and m are my two dimension extents. These are then indirectly accessed by an array of derived type pointers. It works. Its not pretty and a bit slow. I am wondering: do I now do a halo exchange on all of these fields? Can I switch it off? Failed to find that functionality so far.
Many thanks for any hints and comments!
Best wishes,
Markus