I\'m making a function call to a fortran subroutine:
CALL EGFpar(T_e,X_e,Y_e,cpi(1:nx,1:Nspec),WEG,IWEG)
Here IWEG is initialized with intege