_RL a,g,agrav,akwnmb

      g      = 1.0

        gstar = g*frsf*frsf/(frsf*frsf+a*a)

      _RL  getcon,g,r,ak

      integer i,j, g

      DO g = 1,MNC_MAX_ID

        mnc_cw_gname(g)(1:MNC_MAX_CHAR) = blank(1:MNC_MAX_CHAR)

        mnc_cw_vname(g)(1:MNC_MAX_CHAR) = blank(1:MNC_MAX_CHAR)

        mnc_cw_vgind(g) = 0

        mnc_cw_vfmv(g) = 0

        mnc_cw_ndim(g) = 0

        mnc_cw_fgnm(g)(1:MNC_MAX_CHAR) = blank(1:MNC_MAX_CHAR)

        mnc_cw_fgud(g) = 0

        mnc_cw_fgis(g) = 0

        mnc_cw_fgci(g) = 1

          mnc_cw_dn(i,g)(1:MNC_MAX_CHAR) = blank(1:MNC_MAX_CHAR)

          mnc_cw_vtnm(i,g)(1:MNC_MAX_CHAR) = blank(1:MNC_MAX_CHAR)

          mnc_cw_vinm(i,g)(1:MNC_MAX_CHAR) = blank(1:MNC_MAX_CHAR)

          mnc_cw_vdnm(i,g)(1:MNC_MAX_CHAR) = blank(1:MNC_MAX_CHAR)

          mnc_cw_vtat(i,g)(1:MNC_MAX_CHAR) = blank(1:MNC_MAX_CHAR)

          mnc_cw_dims(i,g) = 0

          mnc_cw_is(i,g) = 0

          mnc_cw_ie(i,g) = 0

          mnc_cw_viat(i,g) = 0

          mnc_cw_vdat(i,g) = 0.0D0

          mnc_cw_vnat(i,g) = 0

           mnc_cw_vbij(i,g) = 0

           mnc_cw_vmvi(i,g) = 0

           mnc_cw_vmvr(i,g) = 0.0

           mnc_cw_vmvd(i,g) = 0.0D0

        PtrISt(g)som_P(a,b,c,d,e,f)

    (e.g.: shared_code) which is both in the include path of the

    (e.g.: shared_code) which is both in the include path of the

    (e.g.: shared_code) which is both in the include path of the

  model (e.g., where your SIZE.h is), and comment out these lines.