/*
 * DO NOT MANUALLY UPDATE THIS FILE.
 * If needed: make cpp_headers
 *
 * fortran2c.py - automatic conversion
 *       nucstf.inc -> nucstf.h
 */

#ifndef __NUCSTF__
#define __NUCSTF__

#include "types.h"

#ifdef __cplusplus
	extern "C" {
#endif

//
//=== nucstf ===========================================================*
//
//----------------------------------------------------------------------*
//                                                                      *
//     Copyright (C) 2003-2019:  CERN & INFN                            *
//     All Rights Reserved.                                             *
//                                                                      *
//     NUClear STatement Functions:                                     *
//                                                                      *
//     Created on  14 January 1992  by    Alfredo Ferrari & Paola Sala  *
//                                                   Infn - Milan       *
//                                                                      *
//     Included in the following routines:                              *
//                                                                      *
//                        FPFRNC                                        *
//                        FRADNC                                        *
//                        FRHINC                                        *
//                        FRHONC                                        *
//                        SBCOMP                                        *
//                                                                      *
//----------------------------------------------------------------------*
//
//  NUClear STatement Functions: nucleus with a costant core and
//  linearly decreasing skin and halo
//  rho (r):
//     Statement functions:
//  other
//     RHCORE (R) = RHOCEN
//  Skin with a linear and a cubic term:
//     RHSKIN (R) = RHOCEN * ( ( RADSKN - R ) / SKNEFF + SK3PAR * ( R
//    &           - CPARWS )**3 )
//  Skin with a linear, a cubic and a quartic term:
//     RHSKIN (R) = RHOCEN * ( ( RADSKN - R ) / SKNEFF
//    &           + ( R - CPARWS )**3 * ( SK3PAR + SK4PAR
//    &           * ( R - CPARWS ) ) )
//     RHHALO (R) = RHOCEN * HABPAR * ( RADTOT - R )**1.5D+00
// / r      2
// | 4 pi r' rho (r') dr':
// / 0
//     Statement functions:
//  Other:
//     RICORE (R) = 4.D+00 * PI / 3.D+00 * RHOCEN * R**3
//  Skin with a linear and a cubic term:
//     RISKIN (R) = RICORE (RADIU0) + 4.D+00 * PI / 3.D+00 * RHOCEN
//    &           * ( ( RADSKN / SKNEFF - CPARWS**3 * SK3PAR )
//    &           * ( R**3 - RADIU0**3 ) - 0.75D+00 * ( 1.D+00 / SKNEFF
//    &           - 3.D+00 * SK3PAR * CPARWS**2 ) * ( R**4 - RADIU0**4 )
//    &           - 1.8D+00 * SK3PAR * CPARWS * ( R**5 - RADIU0**5 )
//    &           + 0.5D+00 * SK3PAR * ( R**6 - RADIU0**6 ) )
//  Skin with a linear, a cubic and a quartic term:
//     RISKIN (R) = RICORE (RADIU0) + 4.D+00 * PI / 3.D+00 * RHOCEN
//    &           * ( ( RADSKN / SKNEFF - CPARWS**3 * ( SK3PAR - CPARWS
//    &           * SK4PAR ) ) * ( R**3 - RADIU0**3 ) - 0.75D+00
//    &           * ( 1.D+00 / SKNEFF - CPARWS**2 * ( 3.D+00 * SK3PAR
//    &           - 4.D+00 * SK4PAR * CPARWS ) ) * ( R**4 - RADIU0**4 )
//    &           - 1.8D+00 * CPARWS * ( SK3PAR - 2.D+00 * SK4PAR
//    &           * CPARWS ) * ( R**5 - RADIU0**5 )
//    &           - ( 2.D+00 * SK4PAR * CPARWS - 0.5D+00 * SK3PAR )
//    &           * ( R**6 - RADIU0**6 ) + 3.D+00 / 7.D+00 * SK4PAR
//    &           * ( R**7 - RADIU0**7 ) )
//     RIHALO (R) = RISKIN (RADIU1) + 8.D+00 * PI * RHOCEN
//    &           * HABPAR * ( ( 0.2D+00 * RADTOT**2 - 2.D+00 / 7.D+00
//    &           * RADTOT * HALODP + HALODP**2 / 9.D+00 )
//    &           * ( SQRT (HALODP) )**5 - ( 0.2D+00 * RADTOT**2
//    &           - 2.D+00 / 7.D+00 * RADTOT * ( RADTOT - R )
//    &           + ( RADTOT - R )**2 / 9.D+00 )
//    &           * ( RADTOT - R )**2.5D+00 )
//  r (rho):
//     Statement functions:
//  p_f (rho,i), i=1 proton, 2 neutron:
//     Statement functions:
//  end NUClear STatement Functions

#ifdef __cplusplus
}
#endif
#endif