c ion-param.inc
       PARAMETER (  NI = 3, NJ = 24, NK = 32 )
#include "ion-rest-of-param.inc"