Actually I had this script beacuse an accident, when I was trying to make an equipotential line in order to modeling a geoelectrical method. but I found that my fail in this script can make it looks like a multilevel helix 😀

phi = 3.14;
k=1;
for m = n:n:r;
    i(k)=m;
    k=k+1;    
end
V = ((I*rho)/(4*phi))*(1./i);
a = [0:1:360];
for p=1:length(i);
    x(:,p) = i(p)*cos(a.*pi/180);
    y(:,p) = i(p)*sin(a.*pi/180);
    z(:,p) = linspace(0,100,length(a));
end
figure 
for nl=1:length(i);
    plot3(x(:,:,nl),y(:,:,nl),z(:,:,nl))
    %text(x(1,nl),y(1,nl),[num2str(V(nl),2)])
    %title(‘Nilai equipotential pada tiap lintasan’)
    hold on
    grid on
end

phi = 3.14;
k=1;for m = n:n:r;    i(k)=m;    k=k+1;    end
V = ((I*rho)/(4*phi))*(1./i);
a = [0:1:360];
for p=1:length(i);    x(:,p) = i(p)*cos(a.*pi/180);    y(:,p) = i(p)*sin(a.*pi/180);    z(:,p) = linspace(0,100,length(a));end
figure for nl=1:length(i);    plot3(x(:,:,nl),y(:,:,nl),z(:,:,nl))  

 %text(x(1,nl),y(1,nl),[num2str(V(nl),2)])    

%title(‘Nilai equipotential pada tiap lintasan’)    hold on    grid on

end


 

I Don’t know what can I do when Iam stuck installing some programs in my windows 7. I think that previous version of windows (XP) is the most compatible.
Actually I almost desperate and going to install dual OS in order to install some programs in my laptop.

like ArcGis 9.2, I found that in windows 7 it can’t run. finally, It is no problem when my friend accomodates me by bought ArcGis 9.3 that running well in windows 7. problem fixed.

I got a problem again,  when I was installing google sketch up. I found that an error occured before the installation finished. However I must find the solution , and I start googling. In hours, I found that the error message that say’

” An error occurred during the installation of assembly ‘Microsoft.VC80.ATL,type=”win32″,version=”8.0.50727.4053″,publicKeyToken=”1fc8b3b9a1e18e3b”,processorArchitecture=”x86″‘. Please refer to Help and Support for more information. HRESULT: 0x800736FD.

can be fixed with registry editor.Actually I don’t know what the message means, but in my searching i found that error caused by “running out registry space”.

so the step for solution is :

1. go to registry editor

2. find the HKEY_LOCAL_MACHINE\System\CurrentControlSet

3. Increase the registry space by changing value 5d00000(in hex) which is 97517568 to value fffffffff (in hex value).

I restart my machine, reinstall google sketch up and finally… it works. Thaks to my reference for it’ 

http://geo.geek.nz/development/running-out-of-registry-space-with-microsoft-windows-7-and-arcgis-9-3-1-or-other-programs/

prototipe untuk program downhole seismik //prototipe for downhole seismic program

% u/ downhole seismic

% look for velocity in vertical position

load tz.dat

load vp2a.dat

load vs2.dat

tp=tz(:,2);

ts=tz(:,3);

z=tz(:,1);

figure(1)

plot(z,tp,’*b’),title(‘waktu tempuh gel P terhadap kedalaman’),xlabel(‘depth (m)’),ylabel(‘tp (ms)’)

figure(2)

plot(z,ts,’*g’),title(‘waktu tempuh gel s terhadap kedalaman’),xlabel(‘depth (m)’),ylabel(‘ts (ms)’)

% mencari kecepatan gel P dan gel S

n=length(z)-1;

for i=1:1:n

vp(i)=(z(i+1)-z(i))/(tp(i+1)-tp(i));

vs(i)=(z(i+1)-z(i))/(ts(i+1)-ts(i));

r(i,:)=z(i):z(i+1)

end

% V as depth function with range

k=[r(:,1)

r(:,2)]

g=sort(k)

% to be continue……………………

vp2=vp2a;

figure(3)

plot(g-23,vp2),xlabel(‘kedalaman (m)’),ylabel(‘vs’),title(‘Vs terhadap kedalaman’)

vpt=vp’;

vst=vs’;

vs2=vs2;

figure(4)

plot(g,vs2),xlabel(‘kedalaman (m)’),ylabel(‘vs’),title(‘Vs terhadap kedalaman’)


(program ini sebagian hanya dummy dan tidak terpakai)

Program ini berisi menampilkan lapisan 2 dimensi yang mempunyai nilai densitas atau kecepatan // this program contain showing a 2-D layer that have a value density or velocity

Program ini prototipe untuk tomografi // This is Prototipe for tomography

Pada penjalanan program kita diminta untuk menuliskan koordinat yang akan di chek atau dicari kecepatannya, hasilnya ditampilkan di layar

script :

% PROGRAM LAYER TOMOGRAFI

clear;

clc;

a=0:100;b=0:100;

[x,y]=meshgrid(a,b);

plot(x,y,’.white’),hold on % horee..coba coba kata “white” eh bisa

plot(0,10,’*r’);plot(0,30,’*r’);plot(0,60,’*r’);plot(0,90,’*r’),hold on;

% gambar template

a=0:100;b=80:100;

[x,y]=meshgrid(a,b);

plot(x,y,’.b’),hold on

a=50:100;b=60:80;

[x,y]=meshgrid(a,b);

plot(x,y,’.b’),hold on

a=0:50;b=60:80;

[x,y]=meshgrid(a,b);

plot(x,y,’.g’)

a=0:100;b=40:60;

[x,y]=meshgrid(a,b);

plot(x,y,’.g’)

a=0:100;b=20:40;

[x,y]=meshgrid(a,b);

plot(x,y,’.c’)

a=0:100;b=0:20;

[x,y]=meshgrid(a,b);

plot(x,y,’.m’)

a=input(‘masukkan nilai a ‘)

b=input(‘masukkan nilai b ‘)

% inisiasi kecepatan tiap lapisan

if (b<=20)

v=3500;

elseif b<=40

v=3000;

elseif ((b<=80)&(a<=50))||((b<=60)&(a>50))

v=2500;

else

v=2000;

end

disp([‘maka kecepatan adalah =’ num2str(v)])

Image1

Image2

Mencari nilai koordinat

Ada file”ref.dat” yang berisi nomor titik yang belum ada koordinatnya. Ada file kedua “val.dat”yang berupa file koordinat. Link dari dua data tersebut adalah ada nomer point( point number), yang jika nomernya sama, maka koordinatnya juga akan sama.

// There is file “ref.dat” contain number of points. And there is 2nd file “val.dat” contain coordinate. Link of both of them is the point number. If the number is same, so the coordinate will copied

Jumlah data pada “ref.dat” jauh lebih besar dari “val.dat” karena merupakan titik tiap ketinggian dengan interval tertentu

// Number of points of ” ref.dat” larger than “val.dat” because this is points at each level altitude with certain interval

% program untuk mencari koordinat

clear

clc

load ref.dat % file ini sebagai input nomer titik yang dicari koordinatnya

load val.dat % file ini adalah input nilai koordinat tiap point

a=ref(:,4);

b=val(:,3);c=val(:,1);d=val(:,2);

for j=1:1:length(a)

for i=1:1:length(b)

if a(j)==b(i)

x(j)=c(i);y(j)=d(i);

else

continue

end

end

end

%h=[a x’ y’];

%save line.dat h -ascii

n=[ref(:,1) ref(:,2) ref(:,3) a x’ y’];

save hasil.dat n -ascii


hasil disimpan dalam file “hasil.dat”

most click

  • None

Calendar

June 2022
M T W T F S S
 12345
6789101112
13141516171819
20212223242526
27282930  

category