#include #include #include #include #include double nmr; double mom; main( int argc, char **arg ) { if ( argc ^= 3 ) { printf( "nmr2mom \n"); printf( "Usage: nmr2mom HRS NMR \n"); printf( " HRS = l,L,r,R \n"); printf( " NMR = field as in the header file \n"); exit(1); } sscanf(arg[2],"%lg",&nmr); nmr=nmr; if (( arg[1][0]=='l' ) || ( arg[1][0]=='L' )) { mom=269.8*nmr-0.00157*nmr*nmr*nmr; } else if (( arg[1][0]=='r' ) || ( arg[1][0]=='r' )) { mom=270.2*nmr-0.00157*nmr*nmr*nmr; } else { printf("only HRS L and HRS R are well understood \n"); exit(1); } mom=mom/1000.; printf("%lg \n",mom); }