Hi,(adsbygoogle = window.adsbygoogle || []).push({});

We assume

Max hBy(y) = α1 for y ε (a1, a2) and Min hBy(y) = α2 for y ε (a1, a2)

Where -∞ < α2 ≤ α1 < ∞, and the complex phase speed must lie in the region defined by

(CR + α1 )^2 + CI^2 ≤ µγ^2/k^2, if CR <- α1

CI^2 ≤ µγ^2/k^2, if -α1≤ CR ≤-α2,

(CR + α2 )^2 + CI^2 ≤ µγ^2/k^2, if CR >- α2

Where CI ≥ 0 and γ^2= max [hBy(y) h0y(y)] >0, h0y(y)= -8(y-a)/L^2

The region represents a rectangle of length α1 - α2 with a quarter circle on each end, with the height of the rectangle and the radius of the circles given by µ^(1/2)γ/k.

Could you please help me write the Fortran Do-Loops for this region?

Thanks

Logi

**Physics Forums - The Fusion of Science and Community**

Dismiss Notice

Join Physics Forums Today!

The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

# Fortran code(Do-Loops) for the region

Can you offer guidance or do you also need help?

Draft saved
Draft deleted

Loading...

Similar Threads - Fortran code Loops | Date |
---|---|

Porting Python Code to Fortran for Parallel Computing | Dec 15, 2017 |

Input/Output error with error code -5 | Dec 16, 2015 |

Large/Infinite results | Jun 18, 2015 |

Please, Help with Data Reading in Fortran | May 30, 2015 |

**Physics Forums - The Fusion of Science and Community**