Difference between revisions of "SELFE FAQ"
| Line 1: | Line 1: | ||
| − | + | <OL> | |
| + | <LI> <A "#anchorA" NAME="faq1" href="#faq1_2">How do I set up a vertical grid?</A> | ||
| + | <LI> <A "#anchorA" NAME="faq12" href="#faq12_2">My results show platform/compiler/CPU dependency. | ||
| + | <LI> <A "#anchorA" NAME="faq13" href="#faq13_2">Run crashed with a fort.11 error message"QUICKSEARCH: no intersecting edge....". | ||
| + | </A> | ||
| + | <LI> <A "#anchorA" NAME="faq10" href="#faq10_2">How to do a tidal simulation with SELFE?</A> | ||
| + | <LI> <A "#anchorA" NAME="faq5" href="#faq5_2">How do I make GOTM libraries?</A> | ||
| + | <LI> <A "#anchorA" NAME="faq2" href="#faq2_2">How do I choose the flags in lqk.gr3?</A> | ||
| + | <LI> <A "#anchorA" NAME="faq3" href="#faq3_2">My river and plume are too fresh.</A> | ||
| + | <LI> <A "#anchorA" NAME="faq4" href="#faq4_2">I'd like to add some passive tracer transport to SELFE. How can I interface my own code to it?</A> | ||
| + | <LI> <A "#anchorA" NAME="faq6" href="#faq6_2">How to impose river discharge if the depth is negative there?</A> | ||
| + | <LI> <A "#anchorA" NAME="faq8" href="#faq8_2"> How do I set up sflux/ using NARR (North America only) files in the | ||
| + | Utility Library for atmospheric model (nws=2)? </A> | ||
| + | <LI> <A "#anchorA" NAME="faq7" href="#faq7_2"> How do I prepare my own netcdf files for atmospheric model (nws=2)? </A> | ||
| + | <LI> <A "#anchorA" NAME="faq9" href="#faq9_2"> My non-hydrostatic run blows up. </A> | ||
| + | <LI> <A "#anchorA" NAME="faq14" href="#faq14_2"> I have large velocity at open boundary </A> | ||
| + | <LI> <A "#anchorA" NAME="faq15" href="#faq15_2"> My run crashed; how can I find out | ||
| + | why?</A> | ||
| + | </OL> | ||
| + | <HR> | ||
| + | <HR> | ||
| − | < | + | <HR> |
| − | < | + | <p> </p> |
| − | + | ||
| − | + | <p> <b><font color="#FF0000"><a name="faq1_2">How do I set up a vertical grid? </a></font></b></p> | |
| − | + | <p> </p> You should set it up according to your application. The first question to ask is whether a "pure S" model is | |
| − | + | sufficient for your problem. We found it's adequate for most problems (99.9%) except when strong stratification encounters | |
| − | < | + | steep bottom slope (e.g., strongly stratified Columbia River plume flows over a shelf break |
| − | + | with slope in excess of 1/20). In the S part of vgrid.in, adjust the constants h_c, theta_b, theta_f to resolve appropriate | |
| − | + | boundary layers. Note that h_c also determines the boundary between S-zone and sigma-zone (where traditional sigma-coordinates | |
| − | + | are used becasue S-transformation is invalid when h < h_c). | |
| − | + | ||
| − | + | If (in the unlikely case that) a SZ model is needed, h_s may be chosen as the max. depth of the shelf, and make sure there is | |
| − | + | a smooth transition between S and Z coordinates. You can get a preview of some sample z-coordinates at various depths | |
| − | + | by running the code with the pre-processor flag (ipre in param.in) turned on. | |
| − | </ | + | |
| + | <p> <a href="#faq1">go back</a></p> | ||
Revision as of 12:41, 17 September 2012
- <A "#anchorA" NAME="faq1" href="#faq1_2">How do I set up a vertical grid?</A>
- <A "#anchorA" NAME="faq12" href="#faq12_2">My results show platform/compiler/CPU dependency.
- <A "#anchorA" NAME="faq13" href="#faq13_2">Run crashed with a fort.11 error message"QUICKSEARCH: no intersecting edge....". </A>
- <A "#anchorA" NAME="faq10" href="#faq10_2">How to do a tidal simulation with SELFE?</A>
- <A "#anchorA" NAME="faq5" href="#faq5_2">How do I make GOTM libraries?</A>
- <A "#anchorA" NAME="faq2" href="#faq2_2">How do I choose the flags in lqk.gr3?</A>
- <A "#anchorA" NAME="faq3" href="#faq3_2">My river and plume are too fresh.</A>
- <A "#anchorA" NAME="faq4" href="#faq4_2">I'd like to add some passive tracer transport to SELFE. How can I interface my own code to it?</A>
- <A "#anchorA" NAME="faq6" href="#faq6_2">How to impose river discharge if the depth is negative there?</A>
- <A "#anchorA" NAME="faq8" href="#faq8_2"> How do I set up sflux/ using NARR (North America only) files in the Utility Library for atmospheric model (nws=2)? </A>
- <A "#anchorA" NAME="faq7" href="#faq7_2"> How do I prepare my own netcdf files for atmospheric model (nws=2)? </A>
- <A "#anchorA" NAME="faq9" href="#faq9_2"> My non-hydrostatic run blows up. </A>
- <A "#anchorA" NAME="faq14" href="#faq14_2"> I have large velocity at open boundary </A>
- <A "#anchorA" NAME="faq15" href="#faq15_2"> My run crashed; how can I find out why?</A>
<a name="faq1_2">How do I set up a vertical grid? </a>
You should set it up according to your application. The first question to ask is whether a "pure S" model is
sufficient for your problem. We found it's adequate for most problems (99.9%) except when strong stratification encounters
steep bottom slope (e.g., strongly stratified Columbia River plume flows over a shelf break with slope in excess of 1/20). In the S part of vgrid.in, adjust the constants h_c, theta_b, theta_f to resolve appropriate boundary layers. Note that h_c also determines the boundary between S-zone and sigma-zone (where traditional sigma-coordinates are used becasue S-transformation is invalid when h < h_c).
If (in the unlikely case that) a SZ model is needed, h_s may be chosen as the max. depth of the shelf, and make sure there is
a smooth transition between S and Z coordinates. You can get a preview of some sample z-coordinates at various depths by running the code with the pre-processor flag (ipre in param.in) turned on.
<a href="#faq1">go back</a>