<p><b>ringler@lanl.gov</b> 2010-03-18 21:43:32 -0600 (Thu, 18 Mar 2010)</p><p><br>
file changed: module_scvt.F<br>
<br>
added hyperbolic function to create variable resolution meshes<br>
</p><hr noshade><pre><font color="gray">Modified: trunk/grid_gen/global_scvt/src/module_scvt.F
===================================================================
--- trunk/grid_gen/global_scvt/src/module_scvt.F        2010-03-18 20:34:08 UTC (rev 147)
+++ trunk/grid_gen/global_scvt/src/module_scvt.F        2010-03-19 03:43:32 UTC (rev 148)
@@ -205,7 +205,7 @@
       real :: r, norm, t_cent
       real :: r1
       real :: pi
-!     real :: width, trans_center, min_val
+      real :: width, trans_center, min_val
 
       pi = 4.0*atan(1.0)
 
@@ -214,20 +214,26 @@
       ! Uniform Density Function
        density_for_point = 1.0
 
-      !Target Density Function
-!     p_local%lat = c_lat_rad
-!     p_local%lon = c_lon_rad
-!     call convert_lx(rx, ry, rz, 1.0, p)
-!     call convert_lx(prx, pry, prz, 1.0, p_local)
-!     r = acos(rx*prx + ry*pry + rz*prz)
+      !Target Density Function based on hyperbolic tangent
+      !                p_local%lat = latitude (radians) center of high-resolution region
+      !         p_local%lon = longitude (radians) center of high-resolution region
+      !         width = width of transition zone 
+      !         trans_center = width (radians) of high resolution zone
+      !         minval = minimum density value. to have grid spacing vary by a factor of 8
+      !                  set minval = (1.0 / 8.0)**4
 
-      ! hyperbolic tangent density function
-!     width = 0.15
-!     t_cent = trans_center*pi
-!     min_val = 0.0001
-!     norm = 1.0/(1.0-min_val)
-!     density_for_point = ((tanh((t_cent-r)*(1.0/width))+1.0)/2)/norm + min_val
+  !   p_local%lat = pii/4.0
+  !   p_local%lon = 1.25*pii
+  !   call convert_lx(rx, ry, rz, 1.0, p)
+  !   call convert_lx(prx, pry, prz, 1.0, p_local)
+  !   r = acos(rx*prx + ry*pry + rz*prz)
 
+  !   width = 0.15
+  !   trans_center = pi/6.0
+  !   min_val = (1.0/8.0)**4
+  !   norm = 1.0/(1.0-min_val)
+  !   density_for_point = ((tanh((trans_center-r)*(1.0/width))+1.0)/2)/norm + min_val
+
    end function density_for_point
 
 end module scvt

</font>
</pre>