Difference between revisions of "SELFE FAQ"

From ccrmwiki
Jump to: navigation, search
Line 1: Line 1:
The list of Q/A below applies to both serial and parallel versions.
+
<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>
  
<UL>
+
<HR>
     <LI>How do I set up a vertical grid?
+
     <p> </p>
    <LI>My results show platform/compiler/CPU dependency.
+
 
    <LI>Run crashed with a fort.11 error message"QUICKSEARCH: no intersecting edge....".
+
<p> <b><font color="#FF0000"><a name="faq1_2">How do I set up a vertical grid? </a></font></b></p>
    <LI>How to do a tidal simulation with SELFE?
+
     <p> </p> You should set it up according to your application. The first question to ask is whether a "pure S" model is
    <LI>How do I make GOTM libraries? 
+
sufficient for your problem. We found it's adequate for most problems (99.9%) except when strong stratification encounters
     <LI>My river and plume are too fresh
+
steep bottom slope (e.g., strongly stratified Columbia River plume flows over a shelf break
    <LI>I'd like to add some passive tracer transport to SELFE. How can I interface my own code to it?
+
  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
    <LI>How to impose river discharge if the depth is negative there?
+
boundary layers. Note that h_c also determines the boundary between S-zone and sigma-zone (where traditional sigma-coordinates
    <LI> How do I set up sflux/ using NARR (North America only) files in the Utility Library for atmospheric model (nws=2)?
+
are used becasue S-transformation is invalid when h &lt h_c).
    <LI>How do I prepare my own netcdf files for atmospheric model (nws=2)?
+
 
    <LI>My non-hydrostatic run blows up.
+
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
    <LI> I have large velocity at open boundary
+
a smooth transition between S and Z coordinates. You can get a preview of some sample z-coordinates at various depths
    <!--LI>My run crashed; how can I find out why? -->
+
by running the code with the pre-processor flag (ipre in param.in) turned on.
</UL>
+
 
 +
<p> <a href="#faq1">go back</a></p>

Revision as of 12:41, 17 September 2012

  1. <A "#anchorA" NAME="faq1" href="#faq1_2">How do I set up a vertical grid?</A>
  2. <A "#anchorA" NAME="faq12" href="#faq12_2">My results show platform/compiler/CPU dependency.
  3. <A "#anchorA" NAME="faq13" href="#faq13_2">Run crashed with a fort.11 error message"QUICKSEARCH: no intersecting edge....". </A>
  4. <A "#anchorA" NAME="faq10" href="#faq10_2">How to do a tidal simulation with SELFE?</A>
  5. <A "#anchorA" NAME="faq5" href="#faq5_2">How do I make GOTM libraries?</A>
  6. <A "#anchorA" NAME="faq2" href="#faq2_2">How do I choose the flags in lqk.gr3?</A>
  7. <A "#anchorA" NAME="faq3" href="#faq3_2">My river and plume are too fresh.</A>
  8. <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>
  9. <A "#anchorA" NAME="faq6" href="#faq6_2">How to impose river discharge if the depth is negative there?</A>
  10. <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>
  11. <A "#anchorA" NAME="faq7" href="#faq7_2"> How do I prepare my own netcdf files for atmospheric model (nws=2)? </A>
  12. <A "#anchorA" NAME="faq9" href="#faq9_2"> My non-hydrostatic run blows up. </A>
  13. <A "#anchorA" NAME="faq14" href="#faq14_2"> I have large velocity at open boundary </A>
  14. <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 &lt 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>