Page 12 - Demo

Add-on Creator 2%uf09f %uf09f RCP Consult 2021-2025 Page 12!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! PLOT CREATION USING GNUPLOT !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! real*8 function StoneCircularPlot(ier) * bind(C,name='StoneCircularPlot') implicit none!! Include file with internal STRUREL definitions include 'strurel.hf'!GCC$ ATTRIBUTES DLLEXPORT::StoneCircularPlot! Output function for initial and calculated values interface character*200 function outV(a,b,n,m) real*8 a,b integer n,m end function end interface integer, intent(out) :: ier! Local variable: width, heigth, ouput window number, etc. real*8 w, h, t, l, r integer nf character*200 title,state logical Init!! Global flag to manage plot facilities of engine StrurelPlot=.TRUE.!# Uses global names of Variables and Parameters to access their values.!# The sequence is arbitrary. Note: STRUREL is case sensitive but FORTRAN not. if(StrurelPlot) then if(StrurelMode.eq.0.or.StrurelMode.eq.2) then w=8. h=5. title=\ if(StrurelMode.eq.0) then ! 0 - initial deterministic state state=\ nf=1 else ! 2 - final stochastic solution state=\ nf=2 end if! Initialize an Gnuplot environment Init=GnuplotInit() if(Init) then! Create an output window call GnuplotWindowRelative(nf,40,40,.TRUE.)! Set a title call GnuplotWindowTitle(nf, trim(title)//c_) write(buffer,\',a,'''')\ * trim(title),trim(state) call GnuplotPut(nf,trim(buffer)//c_)! Define range for axes write(buffer,\w call GnuplotPut(nf, trim(buffer)//c_) write(buffer,\h call GnuplotPut(nf, trim(buffer)//c_) call GnuplotPut(nf, \_);! Plot a main subtitle of plot write(buffer,\ * \ at 1,4.5 font ''Arial,16''')\ * trim(StrurelName),trim(StrurelModule),trim(StrurelVersion) call GnuplotPut(nf, trim(buffer)//c_) if(StrurelMode.eq.0) then! Input data at state of mean values call GnuplotPut(nf, \ * \ at 1,4.3 font 'Arial,14'\_) else! Final data of stochastic solution write(buffer,\ * \ at 1,4.3 font ''Arial,14''')\