simp:true;
fortran_float_precision:double;
load("%INCLUDE_/share/vecfem3/fortra.lisp");
NK:%NK_;
ARRAY(DB,NK);
ARRAY(B,NK);
ARRAY(U0,NK);
ARRAY(U0ZZ,NK);
ARRAY(IM,3);

ARRAY(F0,NK,3);
ARRAY(F1,NK,3,3);
ARRAY(F0U,NK,NK,3);
ARRAY(F0UT,NK,NK,3);
ARRAY(F0UX,NK,NK,3,3);
ARRAY(F0UTX,NK,NK,3,3);
ARRAY(F1U,NK,NK,3,3);
ARRAY(F1UT,NK,NK,3,3);
ARRAY(F1UX,NK,NK,3,3,3);
ARRAY(F1UTX,NK,NK,3,3,3);

ARRAY(HF0,NK,3);
ARRAY(HF1,NK,3,3);
ARRAY(HF0U,NK,NK,3);
ARRAY(HF0UT,NK,NK,3);
ARRAY(HF0UX,NK,NK,3,3);
ARRAY(HF0UTX,NK,NK,3,3);
ARRAY(HF1U,NK,NK,3,3);
ARRAY(HF1UT,NK,NK,3,3);
ARRAY(HF1UX,NK,NK,3,3,3);
ARRAY(HF1UTX,NK,NK,3,3,3);

ARRAY(dudxzz,NK,3);
ARRAY(dutdxzz,NK,3);
ARRAY(uzz,NK);
ARRAY(utzz,NK);
ARRAY(dvdxzz,NK,3);
ARRAY(vzz,NK);

FOR k:1 THRU NK DO
   (DB[k]:0,
   U0[k]:0);

FOR class:0 THRU 3 DO
   IM[class]:0;
