Commit 0f66ee57 authored by ldemaine's avatar ldemaine
Browse files

ajout de sortie pour vériefier get_slopes

parent 7a99b3b4
......@@ -7426,12 +7426,24 @@ c$$$ double precision double_hA, double_hB, double_hC
parameter (pi=3.14159265358979d0)
parameter (degtorad=pi/180.0d0)
write(*,*) '|-----------------------------------------|'
write(*,*) '| Input get_slopes'
write(*,*) '|-----------------------------------------|'
write(*,*) '|'
write(*,*) '| latitude', latitude
write(*,*) '| longitude', longitude
write(*,*) '| distance', distance
write(*,*) '| zradius', zradius
write(*,*) '|-----------------------------------------|'
write(*,*) '|'
! 1. Computation of the coordinates of A, B and C
distance = distance*1000
write(*,*) '| distance', distance
c$$$ distance_local = distance
c$$$ write(*,*) 'distance:', distance
c$$$ write(*,*) 'distance_local:', distance_local
lonA = longitude
latA = latitude + distance/(Zradius*degtorad)
......@@ -7443,12 +7455,14 @@ c$$$ write(*,*) 'distance_local:', distance_local
lonC = longitude-distance*sqrt(3.)/
& (2*Zradius*cos(latBC*degtorad)*degtorad)
c$$$ write(*,*) 'lonA, lonB, lonC'
c$$$ write(*,*) lonA, lonB, lonC
c$$$
c$$$ write(*,*) 'latA, latB, latC'
c$$$ write(*,*) latA, latBC, latBC
c$$$
write(*,*) '|'
write(*,*) '| lonA', lonA
write(*,*) '| lonB', lonB
write(*,*) '| lonC', lonC
write(*,*) '| latA', latA
write(*,*) '| latB', latBC
write(*,*) '| latC', latBC
c$$$ write(*,*) 'xlon-lonA, xlon-lonB, xlon-lonC'
c$$$ write(*,*) longitude-lonA, longitude-lonB, longitude-lonC
c$$$
......@@ -7460,6 +7474,10 @@ c$$$ write(*,*) latitude-latA, latitude-latBC, latitude-latBC
call mola(dset,latBC,lonB,hB,ier)
call mola(dset,latBC,lonC,hC,ier)
write(*,*) '|'
write(*,*) '| hA', hA
write(*,*) '| hB', hB
write(*,*) '| hC', hC
c$$$ double_hA, double_hB, double_hC = hA, hB, hC
......@@ -7474,13 +7492,18 @@ c$$$ double_hA, double_hB, double_hC = hA, hB, hC
& hB*hC)), 3*distance)
endif
theta_s = theta_s/degtorad
write(*,*) '|'
write(*,*) '| theta_s', theta_s
psy_s = atan2( hC-hB,
& sqrt(3.)*(hB + hC - 2*hA))
psy_s = psy_s/degtorad
write(*,*) '| psy_s', psy_s
return
......@@ -7755,7 +7778,7 @@ c$$$ real,intent(in) :: pratio
parameter (degtorad=pi/180.0d0)
open(1, file='Test_slope/TEST_5.txt', status = 'new')
open(1, file='Test_slope/TEST_5.txt', status = 'old')
write(1,*) ''
write(1,*) '|-----------------------------------------|'
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment