<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>