#include?
#include?“diom_globals.h“
#include?“diom_procedures.h“

void?diatom_point_test2?（?REAL?*x?REAL?*y?REAL?*z?REAL?*dr?REAL?*value
??INTEGER?*ival?）
{
/*?
??DIATOM_POINT_TEST2

??Discussion:

????This?is?a?modified?version?of?DIATOM_POINT_TEST?as?supplied?by?David
????Crawford.??In?this?version?the?arguments?X?Y?Z?and?DR?are?pointers.
????This?facilitates?calling?the?routine?from?FORTRAN.

????This?routine?functions?as?an?interface?between?higher?level?FORTRAN
????routines?and?the?DIATOM?library?which?is?written?primarily?in?C.

??Modified:

????12?April?2001

??Parameters:

????Input?REAL?*X?*Y?*Z?pointers?to?the?（XYZ）?coordinates?of?a?point.

????Input?REAL?*DR?pointer?to?a?tolerance?used?when?determining?if?a
????point?is?in?on?or?outside?the?region.

????Input?REAL?*VALUE?pointer?to?an?output?quantity?a?material?density?
????to?be?computed?by?this?routine.

????Input?INTEGER?*IVAL?pointer?to?an?output?quantity?an?index?to?be?
????computed?by?this?routine?which?is?1?if?the?point?（XYZ）?is?in?the?
????region?-1?if?it?is?on?the?boundary?of?the?region?and?0?if?it?is?
????outside?the?region.
*/
?static?struct?diom?*curr_diom;
?struct?geom_object?*curr_geom;
?REAL?dt?=?0;
?INTEGER?i;
?int?init?=?1;
?REAL?io;
?INTEGER?kbad;
?INTEGER?mat_id;
?REAL?*mmtmp;
?REAL?time?=?0;
?REAL?tmp;
?REAL?vf;
?REAL?vf_all;
?REAL?x1;
?REAL?x2;
?REAL?y1;
?REAL?y2;
?REAL?z1;
?REAL?z2;

/*?printf?（?“DIATOM_POINT_TEST2:?xyz?=?%f?%f?%f.\n“?*x?*y?*z?）;?*/
/*??printf?（?“DIATOM_POINT_TEST2:?dr?=?%f.\n“?*dr?）;?*/

?x1?=?*x?-?*dr;
?y1?=?*y?-?*dr;
?z1?=?*z?-?*dr;

?x2?=?*x?+?*dr;
?y2?=?*y?+?*dr;
?z2?=?*z?+?*dr;

?curr_diom?=?head_diom;

?if?（?curr_diom?!=?NULL?）?{
??curr_geom?=?curr_diom->geom_list;
?}?else?{
??*ival?=?0;
??printf?（?“\n“?）;
??printf?（?“DIATOM_POINT_TEST2:?Bad?News!\n“?）;
??printf?（?“??The?data?structure?is?not?set?up.\n“?）;
??return;
?}

?vf_all=0;
?*value=0;

?while?（curr_diom?!=?NULL）?{
??vf=0;
??if?（（（init?==?1）?&&?（curr_diom->off?==?time））?||?
??????（（init?==?1）?&&?（curr_diom->off?speed?==?0））?||
??????（（init?==?1）?&&?（curr_diom->anywhere?==?2））?||
??????（（init?==?0）?&&?（time+dt?>=?curr_diom->on）?&&?（time?off）））?{
???if?（diom_brick_VF?!=?NULL）?{
?????if?（（curr_diom->anywhere?==?1）?||?（（curr_diom->anywhere?==?2）?&&?（init?==?1））?||
?????????（in_ureg（number_dimensionsx1?y1?z1?x2?y2?z2）?!=?0））?{
???????curr_geom?=?curr_diom->geom_list;
???????while?（curr_geom?!=?NULL）?{
????????mmtmp?=?curr_geom->minmax;
????????if?（（mmtmp[0]x1）?&&
??????????（（（mmtmp[1]y1））?||?（number_dimensions<2））?&&
??????????（（（mmtmp[2]z1））?||?（number_dimensions<3）））?{
?????????/*curr_geom->GSUB（curr_geom->G&curr_geom->NGxyzdrdrdr&io）;
?????????vf?+=?curr_geom->VFADD*fabs（io）;*/
?????????vf?+=?curr_geom->VFADD*diom_brick_VF（curr_geom?x1?y1?z1?x2?y2?z2?0）;
?????????if?（curr_diom->IP[1]???????????if?（vf???????????

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件????????6544??2010-12-16?07:10??yiweipenguan.f90
?????文件??????124585??2010-12-16?07:36??cvt_basis\cvt_basis.f90
?????文件?????????838??2010-12-16?07:36??cvt_basis\cvt_run_01.inp
?????文件???????10727??2010-12-16?07:36??cvt_basis\cvt_run_01.out
?????文件?????????280??2010-12-16?07:36??cvt_basis\raw.txt
?????文件?????????152??2010-12-16?07:36??cvt_basis\gnuplot_raw.inp
?????文件???????12506??2010-12-16?07:36??cvt_basis\raw.eps
?????文件?????????280??2010-12-16?07:36??cvt_basis\raw2.txt
?????文件?????????155??2010-12-16?07:36??cvt_basis\gnuplot_raw2.inp
?????文件???????12746??2010-12-16?07:36??cvt_basis\raw2.eps
?????文件?????????800??2010-12-16?07:37??cvt_basis\cvt_run_02.inp
?????文件????????5944??2010-12-16?07:37??cvt_basis\cvt_run_02.out
?????文件???????48749??2010-12-16?07:37??cvt_basis\gen_001.txt
?????文件???????48749??2010-12-16?07:37??cvt_basis\gen_002.txt
?????文件???????48749??2010-12-16?07:37??cvt_basis\gen_003.txt
?????文件???????48749??2010-12-16?07:37??cvt_basis\gen_004.txt
?????文件???????48749??2010-12-16?07:37??cvt_basis\gen_005.txt
?????文件???????48749??2010-12-16?07:37??cvt_basis\gen_006.txt
?????文件???????48749??2010-12-16?07:37??cvt_basis\gen_007.txt
?????文件???????48749??2010-12-16?07:37??cvt_basis\gen_008.txt
?????文件????????1075??2010-12-16?07:37??cvt_basis\cvt_run_03.inp
?????文件????????5107??2010-12-16?07:37??cvt_basis\cvt_run_03.out
?????文件???????82345??2010-12-16?07:34??cvt\cvt.f90
?????文件???????14799??2010-12-16?07:34??cvt\cvt_main.f90
?????文件????????4857??2010-12-16?07:34??cvt\diatom_setup.c
?????文件?????????279??2010-12-16?07:34??cvt\diatom_test.in
?????文件????????3711??2010-12-16?07:34??cvt\diatom_point_test2.c
?????文件????????5080??2010-12-16?07:34??cvt\diom_globals.h
?????文件????????7667??2010-12-16?07:34??cvt\diom_procedures.h
?????文件????????2057??2010-12-16?07:34??cvt\cvt_01.out
?????文件????????2127??2010-12-16?07:34??cvt\cvt_02.out
............此處省略73個文件信息