CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > General Forums > Main CFD Forum

Can someone help to check this code on Roe solver

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   June 14, 2019, 21:20
Default Can someone help to check this code on Roe solver
  #1
Member
 
Join Date: Feb 2019
Posts: 59
Rep Power: 5
cfdnewb123 is on a distinguished road
The code is found here,
https://ir.cwi.nl/pub/10964/10964D.pdf

I am trying to write my own code for Roe solver on Matlab and I have been using this paper as a guide. However when I run the code from this paper on Fortran 95, it ran into some floating point error. As I am not familiar with Fortran, I am wondering whether can someone help me check the validity of the code over here.

I am working on the Sod shock tube in particular.
cfdnewb123 is offline   Reply With Quote

Old   June 15, 2019, 03:34
Default
  #2
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,062
Rep Power: 64
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Quote:
Originally Posted by cfdnewb123 View Post
The code is found here,
https://ir.cwi.nl/pub/10964/10964D.pdf

I am trying to write my own code for Roe solver on Matlab and I have been using this paper as a guide. However when I run the code from this paper on Fortran 95, it ran into some floating point error. As I am not familiar with Fortran, I am wondering whether can someone help me check the validity of the code over here.

I am working on the Sod shock tube in particular.

First, depending on the compiler, use all the options to debug the code. Disable all optimizations and use the checkbound for arrays and other setting for the floating points errors in the link procedure
FMDenaro is offline   Reply With Quote

Old   June 15, 2019, 04:20
Default
  #3
Member
 
Join Date: Feb 2019
Posts: 59
Rep Power: 5
cfdnewb123 is on a distinguished road
Quote:
Originally Posted by FMDenaro View Post
First, depending on the compiler, use all the options to debug the code. Disable all optimizations and use the checkbound for arrays and other setting for the floating points errors in the link procedure
I am sorry I don't really get this as I am still new to Fortran. I am using FTN95 (https://www.silverfrost.com/). I have try to debug the code but as I am using this paper as I guide, it is really hard for me to know whether the code is correct or not.
cfdnewb123 is offline   Reply With Quote

Old   June 15, 2019, 04:23
Default
  #4
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,062
Rep Power: 64
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Quote:
Originally Posted by cfdnewb123 View Post
I am sorry I don't really get this as I am still new to Fortran. I am using FTN95 (https://www.silverfrost.com/). I have try to debug the code but as I am using this paper as I guide, it is really hard for me to know whether the code is correct or not.



If you compile correctly the source code, there is no conflict in the FTN95 language. You need to check the error in the run-time procedure.
Use the specific debug options in compiling the source, see here https://www.silverfrost.com/ftn95-he...lassified.aspx
FMDenaro is offline   Reply With Quote

Reply

Tags
fortran, riemann solvers, roe solver, shock, shock-tube

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
The ANSYS CFX solver exited with return code 1 kola77 CFX 24 April 11, 2022 07:32
Compilation error for OpenFoam FSI solver with third party FEA code danny261083 OpenFOAM Running, Solving & CFD 1 October 6, 2017 09:50
Changing Temperature of patch within code of solver SalvusApfel OpenFOAM Pre-Processing 1 February 24, 2017 07:29
modify the source code of simpleSRFFoam solver Pirlo OpenFOAM Programming & Development 0 October 19, 2009 08:49
OpenFOAM15 paraFoam bug koen OpenFOAM Bugs 19 June 30, 2009 10:46


All times are GMT -4. The time now is 03:19.