bc_update.f90
dt_adjust.f90
dt_mass
flac
fl_move
fl_node
fl_rheol
78,79d77
< !	   sarc1v=0.978*sarc1v 
< !	   sarc2v=0.8317*sarc1v*poiss/(1-poiss) 				       !G.Ito (Calibrated values)
82d79
< 
86d82
< 	   sarc2p=0.9882*sarc1p*poiss/(1-poiss);
92,93d87
< !	   sarc1v=0.978*sarc1v 
< !	   sarc2v=0.8317*sarc1v*poiss/(1-poiss) 				       !G.Ito (Calibrated values)
96d89
< 
100d92
< 	   sarc2p=0.9882*sarc1p*poiss/(1-poiss);
fl_srate
fl_therm
functions
init_areas
init_bc
init_cord
36c36,37
< dx_init = abs(cord(1,2,1)-cord(1,1,1))
---
> dx_init = abs(cord(1,(nx/2),1)-cord(1,(nx/2) + 1,1))
> !write(*,*) 'd1 = ',iinj,dx_init
init_phase
init_stress
init_temp
12,29d11
< if (mod(nloop,100).eq.0.and.nloop<1000) write(*,*) '>>>>>>>> INIT_TEMP <<<<<<<<<<', nloop  !G.Ito
< if (igeotherm.eq.10) then
<    if (nloop<2) write(*,*) '>> Initial Error Function Temperature Profile << '  !G.Ito
<    uspread=0.0;
<    do ii=1,nofbc
<      if ((nofside(ii).eq.1).or.(nofside(ii).eq.3)) then
< 	uspread=dmax1(dabs(bca(ii)),uspread)
<      endif
<      if (nofside(ii).eq.4) then
<        uspread=dmax1(dabs(bca(ii)),uspread)
<      endif
<    enddo	
<   if (uspread.eq.0.0) then
<     write(*,*) 'igeotherm=10 for error function temperature profile requires a '
<     write(*,*) 'kinematically driven-spreading rate.  Stopping in INIT_TEMP'
<     stop
<   endif
< endif
57c39
< !if (iynts.eq.1) then		G.Ito (commented out)
---
> if (iynts.ge.1) then
84,96d65
< !----------------------------------------------------------------------------------------
< ! Error function temperature profile. G. Ito 7/04
< !----------------------------------------------------------------------------------------
<        if (igeotherm.eq.10) then 
< 	 diff = 1e-6;
< 	 diff = conduct(1)/cp(1)/3300.0
< 	 difft=dabs(xc)/uspread
< 	 difft=dmax1(difft,1e+11) !limit temperatures near ridge axis
<          zz=dabs(yc)/(2*dsqrt(diff*difft))
< 	 temp(j,i)=t_bot-t_bot*erfc(zz)
< !	 if (i.eq.nx/2.or.i.eq.2) write(*,'(I6,5e12.3)') j,yc,xc,difft,uspread,temp(j,i)
<        endif
< 	 	
101c70
< !endif	G.Ito
---
> endif
init_visc
matprops
outflac
161a162
> 
315a317
> 
par
read_more
read_params
rem_cord
33a34
> !cord(1,iinj,2) = (cord(1,iinj-1,2) + cord(1,iinj+1,2)) / 2.
35c36
< write(*,*) cord(1,iinj-3,1),cord(1,iinj-2,1),cord(1,iinj-1,1),cord(1,iinj,1)
---
> !write(*,*) cord(1,iinj-3,1),cord(1,iinj-2,1),cord(1,iinj-1,1),cord(1,iinj,1)
remesh
rem_test
rh_elastic
rh_maxwell
rh_plastic
rh_visc
rmasses
rsflac
saveflac
setflac
user_ab
user_luc
arrays
params
precision