From a877014f79685437298b1690c3b2592f8ae57e41 Mon Sep 17 00:00:00 2001 From: Lorenzo Figini Date: Thu, 15 Nov 2012 14:52:31 +0000 Subject: [PATCH] added index_rt to fort.33 (ray trajectories) output --- src/gray.f | 10 +++++----- 1 file changed, 5 insertions(+), 5 deletions(-) diff --git a/src/gray.f b/src/gray.f index 8322d92..a62bdc3 100644 --- a/src/gray.f +++ b/src/gray.f @@ -479,7 +479,7 @@ c print dIds in MA/m, Jcd in MA/m^2, ppabs and pdjki in MW/m^3 if(istpl.eq.istpl0) then if(j.eq.nray.and.tauv(j,k,i).le.taucr) then write(33,111) i,j,k,stm,xxm,yym,rrm,zzm, - . psjki(j,k,i),taujk,anpl,alphav(j,k,i) + . psjki(j,k,i),taujk,anpl,alphav(j,k,i),dble(index_rt) end if c if(k.eq.ktx) write(33,*) ' ' end if @@ -503,7 +503,7 @@ c write(8,*) ' #istep j k xt yt zt rt psin' write(9,*) ' #istep j k xt yt zt rt psin' write(17,*) ' #sst Dr_11 Dr_Nr1 Di_Nr1' - write(33,*) ' #i j k sst x y R z psi tauv Npl alpha' + write(33,*) ' #i j k sst x y R z psi tauv Npl alpha index_rt' write(12,*) ' #i sst psi w1 w2' write(7,*)'#Icd Pa Jphimx dPdVmx '// .'rhotj rhotjava rhotp rhotpav '// @@ -3921,7 +3921,7 @@ c r0m=r0/1.0d2 z0m=z0/1.0d2 write(33,111) izero,j,k,zero,x0m,y0m,r0m,z0m, - . -1.0d0,zero,zero,zero + . -1.0d0,zero,zero,zero,one end if if(j.eq.1.and.k.eq.1) then write(17,99) zero,zero,zero,zero @@ -4076,7 +4076,7 @@ c r0m=r0/1.0d2 z0m=z0/1.0d2 write(33,111) izero,j,k,zero,x0m,y0m,r0m,z0m, - . -1.0d0,zero,zero,zero + . -1.0d0,zero,zero,zero,one end if if(j.eq.1.and.k.eq.1) then write(17,99) zero,zero,zero,zero @@ -4175,7 +4175,7 @@ c r0m=r0/1.0d2 z0m=z0/1.0d2 write(33,111) izero,j,k,zero,x0m,y0m,r0m,z0m, - . -1.0d0,zero,zero,zero + . -1.0d0,zero,zero,zero,one end if if(j.eq.1.and.k.eq.1) then write(17,99) zero,zero,zero,zero