| 35 |
|
#define G_BOXZ 14 |
| 36 |
|
#define G_SAMPLETIME 15 |
| 37 |
|
#define G_STATUSTIME 16 |
| 38 |
< |
#define G_RRF 17 |
| 38 |
> |
#define G_ECR 17 |
| 39 |
|
#define G_DIELECTRIC 18 |
| 40 |
|
#define G_TEMPSET 19 |
| 41 |
|
#define G_THERMALTIME 20 |
| 42 |
< |
#define G_USEPBC 21 |
| 43 |
< |
#define G_MIXINGRULE 22 |
| 42 |
> |
#define G_USEPBC 21 |
| 43 |
> |
#define G_MIXINGRULE 22 |
| 44 |
> |
#define G_EST 23 |
| 45 |
> |
#define G_USERF 24 |
| 46 |
|
|
| 47 |
|
|
| 48 |
|
Globals::Globals(){ |
| 72 |
|
addHash( "boxZ", G_BOXZ ); |
| 73 |
|
addHash( "sampleTime", G_SAMPLETIME ); |
| 74 |
|
addHash( "statusTime", G_STATUSTIME ); |
| 75 |
< |
addHash( "rRF", G_RRF ); |
| 75 |
> |
addHash( "electrostaticCutoffRadius", G_ECR ); |
| 76 |
|
addHash( "dielectric", G_DIELECTRIC ); |
| 77 |
|
addHash( "tempSet", G_TEMPSET ); |
| 78 |
|
addHash( "thermalTime", G_THERMALTIME ); |
| 79 |
|
addHash( "mixingRule", G_MIXINGRULE); |
| 80 |
< |
addHash( "periodicBox", G_USEPBC); |
| 80 |
> |
addHash( "usePeriodicBoundaryConditions", G_USEPBC); |
| 81 |
> |
addHash( "electrostaticSkinThickness", G_EST ); |
| 82 |
> |
addHash( "useReactionField", G_USERF ); |
| 83 |
|
|
| 84 |
|
// define some default values |
| 85 |
|
|
| 86 |
|
strcpy( mixingRule,"standard"); //default mixing rules to standard. |
| 87 |
|
usePBC = 1; //default periodic boundry conditions to on |
| 88 |
+ |
useRF = 0; |
| 89 |
|
|
| 85 |
– |
|
| 90 |
|
have_force_field = 0; |
| 91 |
|
have_n_components = 0; |
| 92 |
|
have_target_temp = 0; |
| 106 |
|
have_sample_time = 0; |
| 107 |
|
have_status_time = 0; |
| 108 |
|
have_thermal_time = 0; |
| 109 |
< |
have_rrf = 0; |
| 109 |
> |
have_ecr = 0; |
| 110 |
|
have_dielectric = 0; |
| 111 |
|
have_tempSet = 0; |
| 112 |
+ |
have_est = 0; |
| 113 |
|
} |
| 114 |
|
|
| 115 |
|
Globals::~Globals(){ |
| 644 |
|
} |
| 645 |
|
break; |
| 646 |
|
|
| 647 |
< |
case G_RRF: |
| 647 |
> |
case G_ECR: |
| 648 |
|
switch( the_type ){ |
| 649 |
|
|
| 650 |
|
case STRING: |
| 651 |
|
the_event->err_msg = |
| 652 |
< |
strdup( "Global error. rRF is not a double or int.\n" ); |
| 652 |
> |
strdup( "Global error. electrostaticCutoffRadius is not a double or int.\n" ); |
| 653 |
|
return 0; |
| 654 |
|
break; |
| 655 |
|
|
| 656 |
|
case DOUBLE: |
| 657 |
< |
rRF = the_event->evt.asmt.rhs.dval; |
| 658 |
< |
have_rrf = 1; |
| 657 |
> |
ecr = the_event->evt.asmt.rhs.dval; |
| 658 |
> |
have_ecr = 1; |
| 659 |
|
return 1; |
| 660 |
|
break; |
| 661 |
|
|
| 662 |
|
case INT: |
| 663 |
< |
rRF = (double)the_event->evt.asmt.rhs.ival; |
| 664 |
< |
have_rrf = 1; |
| 663 |
> |
ecr = (double)the_event->evt.asmt.rhs.ival; |
| 664 |
> |
have_ecr = 1; |
| 665 |
|
return 1; |
| 666 |
|
break; |
| 667 |
|
|
| 668 |
|
default: |
| 669 |
|
the_event->err_msg = |
| 670 |
< |
strdup( "Global error. rRF unrecognized.\n" ); |
| 670 |
> |
strdup( "Global error. electrostaticCutoffRadius unrecognized.\n" ); |
| 671 |
> |
return 0; |
| 672 |
> |
break; |
| 673 |
> |
} |
| 674 |
> |
break; |
| 675 |
> |
|
| 676 |
> |
case G_EST: |
| 677 |
> |
switch( the_type ){ |
| 678 |
> |
|
| 679 |
> |
case STRING: |
| 680 |
> |
the_event->err_msg = |
| 681 |
> |
strdup( "Global error. electrostaticSkinThickness is not a double or int.\n" ); |
| 682 |
|
return 0; |
| 683 |
|
break; |
| 684 |
+ |
|
| 685 |
+ |
case DOUBLE: |
| 686 |
+ |
est = the_event->evt.asmt.rhs.dval; |
| 687 |
+ |
have_est = 1; |
| 688 |
+ |
return 1; |
| 689 |
+ |
break; |
| 690 |
+ |
|
| 691 |
+ |
case INT: |
| 692 |
+ |
est = (double)the_event->evt.asmt.rhs.ival; |
| 693 |
+ |
have_est = 1; |
| 694 |
+ |
return 1; |
| 695 |
+ |
break; |
| 696 |
+ |
|
| 697 |
+ |
default: |
| 698 |
+ |
the_event->err_msg = |
| 699 |
+ |
strdup( "Global error. electrostaticSkinThickness unrecognized.\n" ); |
| 700 |
+ |
return 0; |
| 701 |
+ |
break; |
| 702 |
|
} |
| 703 |
|
break; |
| 704 |
|
|
| 757 |
|
else if( !strcmp( "false", the_event->evt.asmt.rhs.sval )) usePBC = 0; |
| 758 |
|
else{ |
| 759 |
|
the_event->err_msg = |
| 760 |
< |
strdup( "Global error. periodicBoundryConditions was not \"true\" or \"false\".\n" ); |
| 760 |
> |
strdup( "Global error. usePeriodicBoundaryConditions was not \"true\" or \"false\".\n" ); |
| 761 |
|
return 0; |
| 762 |
|
} |
| 763 |
|
return 1; |
| 764 |
|
} |
| 765 |
|
|
| 766 |
|
the_event->err_msg = |
| 767 |
< |
strdup( "Global error. tempSet was not \"true\" or \"false\".\n" ); |
| 767 |
> |
strdup( "Global error. usePeriodicBoundaryConditions was not \"true\" or \"false\".\n" ); |
| 768 |
|
return 0; |
| 769 |
|
break; |
| 770 |
+ |
|
| 771 |
+ |
case G_USERF: |
| 772 |
+ |
if( the_type == STRING ){ |
| 773 |
+ |
|
| 774 |
+ |
if( !strcmp( "true", the_event->evt.asmt.rhs.sval )) useRF = 1; |
| 775 |
+ |
else if( !strcmp( "false", the_event->evt.asmt.rhs.sval )) useRF = 0; |
| 776 |
+ |
else{ |
| 777 |
+ |
the_event->err_msg = |
| 778 |
+ |
strdup( "Global error. useReactionField was not \"true\" or \"false\".\n" ); |
| 779 |
+ |
return 0; |
| 780 |
+ |
} |
| 781 |
+ |
return 1; |
| 782 |
+ |
} |
| 783 |
+ |
|
| 784 |
+ |
the_event->err_msg = |
| 785 |
+ |
strdup( "Global error. useReactionField was not \"true\" or \"false\".\n" ); |
| 786 |
+ |
return 0; |
| 787 |
+ |
break; |
| 788 |
+ |
|
| 789 |
+ |
|
| 790 |
|
} |
| 791 |
|
} |
| 792 |
|
|