Difference between revisions of "Getdiff"
Jump to navigation
Jump to search
Line 2: | Line 2: | ||
This function returns the diffusivity of ice flow for a point along the glacier. This is determined from | This function returns the diffusivity of ice flow for a point along the glacier. This is determined from | ||
− | <math>D^{'}= -\frac{2fA}{n+2} \left({\rho\,g}\right)^{n} < | + | <math>D^{'}= -\frac{2fA}{n+2} \left({\rho\,g}\right)^{n} </math> |
where standard symbols are defined in Table 1. The diffusivity is scaled and is determined on a staggered grid that has points positioned halfway between adjacent thickness grid points. The function is called for each staggered grid point in turn; this means that it is the responsibility of the calling routine to deal with branching in the system. | where standard symbols are defined in Table 1. The diffusivity is scaled and is determined on a staggered grid that has points positioned halfway between adjacent thickness grid points. The function is called for each staggered grid point in turn; this means that it is the responsibility of the calling routine to deal with branching in the system. | ||
This function is called from calcthck and the returned information is held in array diffus. | This function is called from calcthck and the returned information is held in array diffus. |
Revision as of 17:53, 6 November 2007
function getdiff
This function returns the diffusivity of ice flow for a point along the glacier. This is determined from [math]\displaystyle{ D^{'}= -\frac{2fA}{n+2} \left({\rho\,g}\right)^{n} }[/math]
where standard symbols are defined in Table 1. The diffusivity is scaled and is determined on a staggered grid that has points positioned halfway between adjacent thickness grid points. The function is called for each staggered grid point in turn; this means that it is the responsibility of the calling routine to deal with branching in the system. This function is called from calcthck and the returned information is held in array diffus.