


DEBUG2 - comparison of CIELO and ICE solutions This script is used to debug Cielo and Ice by comparing the results of both solutions


0001 %DEBUG2 - comparison of CIELO and ICE solutions 0002 % 0003 % This script is used to debug Cielo and Ice by comparing 0004 % the results of both solutions 0005 0006 clear @pentaelem/CreateKMatrix 0007 clear @pentaelem/CreatePVector 0008 clear @triaelem/CreateKMatrix 0009 clear @triaelem/CreatePVector 0010 clear @triaelem/CreateKMatrixHoriz 0011 clear @triaelem/CreatePVectorHoriz 0012 clear @icefront/CreatePVector 0013 0014 md.element_debug=1; 0015 md.element_debugid=41; 0016 0017 %first run ice 0018 md.vx=NaN; md.vy=NaN; md.vz=NaN; md.vel=NaN; 0019 md.cluster='none';md.batch=0; 0020 0021 try, 0022 md=solve(md,'diagnostic','ice'); 0023 catch, 0024 end 0025 0026 %now run cielo 0027 md.cluster='wilkes';md.batch=1; 0028 md=solve(md,'diagnostic_horiz','cielo');