anchor.diff:                  if (any(shape(var_stor) /= correct_dims)) &
                       call error("Mismatched dimensions for variable storage!")
               else
  -                allocate(var_stor(correct_dims(2), correct_dims(1)))
  +                allocate(var_stor(correct_dims(1), correct_dims(2)))
               end if
           end subroutine nc_diag_read_assert_var_dims_2d_byte
           
  
positive.code:                 if (any(shape(var_stor) /= correct_dims)) &
                      call error("Mismatched dimensions for variable storage!")
              else
                  allocate(var_stor(correct_dims(2), correct_dims(1)))
              end if
          end subroutine nc_diag_read_assert_var_dims_2d_long
          
  
negative.code:                 if (any(shape(var_stor) /= correct_dims)) &
                      call error("Mismatched dimensions for variable storage!")
              else
                  allocate(var_stor(correct_dims(1), correct_dims(2)))
              end if
          end subroutine nc_diag_read_assert_var_dims_2d_float
          
  