Skip to content
Snippets Groups Projects
Commit d04c0e10 authored by Douglas Guptill's avatar Douglas Guptill
Browse files

split long lines

parent 6c565651
No related branches found
No related tags found
No related merge requests found
Hide whitespace changes
Inline Side-by-side
Showing
with 54 additions and 27 deletions
define_surface.f90
+ 54
− 27
View file @ d04c0e10
...@@ -71,15 +71,24 @@ if (params%irestart.eq.0) then ...@@ -71,15 +71,24 @@ if (params%irestart.eq.0) then
surface(is)%itype=surface(1)%itype surface(is)%itype=surface(1)%itype
surface(is)%surface_type=surface(1)%surface_type surface(is)%surface_type=surface(1)%surface_type
surface(is)%rand=surface(1)%rand surface(is)%rand=surface(1)%rand
allocate (surface(is)%x(surface(is)%nsurface),stat=threadinfo%err) ; call heap (threadinfo,'surface(is)%x','define_surface',size(surface(is)%x),'dp',+1) allocate (surface(is)%x(surface(is)%nsurface),stat=threadinfo%err)
allocate (surface(is)%y(surface(is)%nsurface),stat=threadinfo%err) ; call heap (threadinfo,'surface(is)%y','define_surface',size(surface(is)%y),'dp',+1) call heap (threadinfo,'surface(is)%x','define_surface',size(surface(is)%x),'dp',+1)
allocate (surface(is)%z(surface(is)%nsurface),stat=threadinfo%err) ; call heap (threadinfo,'surface(is)%z','define_surface',size(surface(is)%z),'dp',+1) allocate (surface(is)%y(surface(is)%nsurface),stat=threadinfo%err)
allocate (surface(is)%xn(surface(is)%nsurface),stat=threadinfo%err) ; call heap (threadinfo,'surface(is)%xn','define_surface',size(surface(is)%xn),'dp',+1) call heap (threadinfo,'surface(is)%y','define_surface',size(surface(is)%y),'dp',+1)
allocate (surface(is)%yn(surface(is)%nsurface),stat=threadinfo%err) ; call heap (threadinfo,'surface(is)%yn','define_surface',size(surface(is)%yn),'dp',+1) allocate (surface(is)%z(surface(is)%nsurface),stat=threadinfo%err)
allocate (surface(is)%zn(surface(is)%nsurface),stat=threadinfo%err) ; call heap (threadinfo,'surface(is)%zn','define_surface',size(surface(is)%zn),'dp',+1) call heap (threadinfo,'surface(is)%z','define_surface',size(surface(is)%z),'dp',+1)
allocate (surface(is)%r(surface(is)%nsurface),stat=threadinfo%err) ; call heap (threadinfo,'surface(is)%r','define_surface',size(surface(is)%r),'dp',+1) allocate (surface(is)%xn(surface(is)%nsurface),stat=threadinfo%err)
allocate (surface(is)%s(surface(is)%nsurface),stat=threadinfo%err) ; call heap (threadinfo,'surface(is)%s','define_surface',size(surface(is)%s),'dp',+1) call heap (threadinfo,'surface(is)%xn','define_surface',size(surface(is)%xn),'dp',+1)
allocate (surface(is)%icon(3,surface(is)%nt),stat=threadinfo%err) ; call heap (threadinfo,'surface(is)%icon','define_surface',size(surface(is)%icon),'int',+1) allocate (surface(is)%yn(surface(is)%nsurface),stat=threadinfo%err)
call heap (threadinfo,'surface(is)%yn','define_surface',size(surface(is)%yn),'dp',+1)
allocate (surface(is)%zn(surface(is)%nsurface),stat=threadinfo%err)
call heap (threadinfo,'surface(is)%zn','define_surface',size(surface(is)%zn),'dp',+1)
allocate (surface(is)%r(surface(is)%nsurface),stat=threadinfo%err)
call heap (threadinfo,'surface(is)%r','define_surface',size(surface(is)%r),'dp',+1)
allocate (surface(is)%s(surface(is)%nsurface),stat=threadinfo%err)
call heap (threadinfo,'surface(is)%s','define_surface',size(surface(is)%s),'dp',+1)
allocate (surface(is)%icon(3,surface(is)%nt),stat=threadinfo%err)
call heap (threadinfo,'surface(is)%icon','define_surface',size(surface(is)%icon),'int',+1)
surface(is)%x=surface(1)%x surface(is)%x=surface(1)%x
surface(is)%y=surface(1)%y surface(is)%y=surface(1)%y
surface(is)%z=surface(1)%z surface(is)%z=surface(1)%z
...@@ -95,15 +104,24 @@ if (params%irestart.eq.0) then ...@@ -95,15 +104,24 @@ if (params%irestart.eq.0) then
surface(is)%nsurface=(2**surface(is)%levelt+2)*(2**surface(is)%levelt+2) !opla surface(is)%nsurface=(2**surface(is)%levelt+2)*(2**surface(is)%levelt+2) !opla
! surface(is)%nt=(2**surface(is)%levelt)*(2**surface(is)%levelt)*2 ! surface(is)%nt=(2**surface(is)%levelt)*(2**surface(is)%levelt)*2
surface(is)%nt=(2**surface(is)%levelt +1)*(2**surface(is)%levelt +1)*2 ! opla surface(is)%nt=(2**surface(is)%levelt +1)*(2**surface(is)%levelt +1)*2 ! opla
allocate (surface(is)%x(surface(is)%nsurface),stat=threadinfo%err) ; call heap (threadinfo,'surface(is)%x','define_surface',size(surface(is)%x),'dp',+1) allocate (surface(is)%x(surface(is)%nsurface),stat=threadinfo%err)
allocate (surface(is)%y(surface(is)%nsurface),stat=threadinfo%err) ; call heap (threadinfo,'surface(is)%y','define_surface',size(surface(is)%y),'dp',+1) call heap (threadinfo,'surface(is)%x','define_surface',size(surface(is)%x),'dp',+1)
allocate (surface(is)%z(surface(is)%nsurface),stat=threadinfo%err) ; call heap (threadinfo,'surface(is)%z','define_surface',size(surface(is)%z),'dp',+1) allocate (surface(is)%y(surface(is)%nsurface),stat=threadinfo%err)
allocate (surface(is)%xn(surface(is)%nsurface),stat=threadinfo%err) ; call heap (threadinfo,'surface(is)%xn','define_surface',size(surface(is)%xn),'dp',+1) call heap (threadinfo,'surface(is)%y','define_surface',size(surface(is)%y),'dp',+1)
allocate (surface(is)%yn(surface(is)%nsurface),stat=threadinfo%err) ; call heap (threadinfo,'surface(is)%yn','define_surface',size(surface(is)%yn),'dp',+1) allocate (surface(is)%z(surface(is)%nsurface),stat=threadinfo%err)
allocate (surface(is)%zn(surface(is)%nsurface),stat=threadinfo%err) ; call heap (threadinfo,'surface(is)%zn','define_surface',size(surface(is)%zn),'dp',+1) call heap (threadinfo,'surface(is)%z','define_surface',size(surface(is)%z),'dp',+1)
allocate (surface(is)%r(surface(is)%nsurface),stat=threadinfo%err) ; call heap (threadinfo,'surface(is)%r','define_surface',size(surface(is)%r),'dp',+1) allocate (surface(is)%xn(surface(is)%nsurface),stat=threadinfo%err)
allocate (surface(is)%s(surface(is)%nsurface),stat=threadinfo%err) ; call heap (threadinfo,'surface(is)%s','define_surface',size(surface(is)%s),'dp',+1) call heap (threadinfo,'surface(is)%xn','define_surface',size(surface(is)%xn),'dp',+1)
allocate (surface(is)%icon(3,surface(is)%nt),stat=threadinfo%err) ; call heap (threadinfo,'surface(is)%icon','define_surface',size(surface(is)%icon),'int',+1) allocate (surface(is)%yn(surface(is)%nsurface),stat=threadinfo%err)
call heap (threadinfo,'surface(is)%yn','define_surface',size(surface(is)%yn),'dp',+1)
allocate (surface(is)%zn(surface(is)%nsurface),stat=threadinfo%err)
call heap (threadinfo,'surface(is)%zn','define_surface',size(surface(is)%zn),'dp',+1)
allocate (surface(is)%r(surface(is)%nsurface),stat=threadinfo%err)
call heap (threadinfo,'surface(is)%r','define_surface',size(surface(is)%r),'dp',+1)
allocate (surface(is)%s(surface(is)%nsurface),stat=threadinfo%err)
call heap (threadinfo,'surface(is)%s','define_surface',size(surface(is)%s),'dp',+1)
allocate (surface(is)%icon(3,surface(is)%nt),stat=threadinfo%err)
call heap (threadinfo,'surface(is)%icon','define_surface',size(surface(is)%icon),'int',+1)
call create_surf(surface(is),is,params%debug) call create_surf(surface(is),is,params%debug)
endif endif
current_time=0.d0 current_time=0.d0
...@@ -147,15 +165,24 @@ else ...@@ -147,15 +165,24 @@ else
read (9) surface(is)%nsurface, & read (9) surface(is)%nsurface, &
activation_time, & activation_time, &
surface(is)%nt surface(is)%nt
allocate (surface(is)%x(surface(is)%nsurface),stat=threadinfo%err) ; call heap (threadinfo,'surface(is)%x','define_surface',size(surface(is)%x),'dp',+1) allocate (surface(is)%x(surface(is)%nsurface),stat=threadinfo%err)
allocate (surface(is)%y(surface(is)%nsurface),stat=threadinfo%err) ; call heap (threadinfo,'surface(is)%y','define_surface',size(surface(is)%y),'dp',+1) call heap (threadinfo,'surface(is)%x','define_surface',size(surface(is)%x),'dp',+1)
allocate (surface(is)%z(surface(is)%nsurface),stat=threadinfo%err) ; call heap (threadinfo,'surface(is)%z','define_surface',size(surface(is)%z),'dp',+1) allocate (surface(is)%y(surface(is)%nsurface),stat=threadinfo%err)
allocate (surface(is)%xn(surface(is)%nsurface),stat=threadinfo%err) ; call heap (threadinfo,'surface(is)%xn','define_surface',size(surface(is)%xn),'dp',+1) call heap (threadinfo,'surface(is)%y','define_surface',size(surface(is)%y),'dp',+1)
allocate (surface(is)%yn(surface(is)%nsurface),stat=threadinfo%err) ; call heap (threadinfo,'surface(is)%yn','define_surface',size(surface(is)%yn),'dp',+1) allocate (surface(is)%z(surface(is)%nsurface),stat=threadinfo%err)
allocate (surface(is)%zn(surface(is)%nsurface),stat=threadinfo%err) ; call heap (threadinfo,'surface(is)%zn','define_surface',size(surface(is)%zn),'dp',+1) call heap (threadinfo,'surface(is)%z','define_surface',size(surface(is)%z),'dp',+1)
allocate (surface(is)%r(surface(is)%nsurface),stat=threadinfo%err) ; call heap (threadinfo,'surface(is)%r','define_surface',size(surface(is)%x),'dp',+1) allocate (surface(is)%xn(surface(is)%nsurface),stat=threadinfo%err)
allocate (surface(is)%s(surface(is)%nsurface),stat=threadinfo%err) ; call heap (threadinfo,'surface(is)%s','define_surface',size(surface(is)%x),'dp',+1) call heap (threadinfo,'surface(is)%xn','define_surface',size(surface(is)%xn),'dp',+1)
allocate (surface(is)%icon(3,surface(is)%nt),stat=threadinfo%err) ; call heap (threadinfo,'surface(is)%icon','define_surface',size(surface(is)%icon),'int',+1) allocate (surface(is)%yn(surface(is)%nsurface),stat=threadinfo%err)
call heap (threadinfo,'surface(is)%yn','define_surface',size(surface(is)%yn),'dp',+1)
allocate (surface(is)%zn(surface(is)%nsurface),stat=threadinfo%err)
call heap (threadinfo,'surface(is)%zn','define_surface',size(surface(is)%zn),'dp',+1)
allocate (surface(is)%r(surface(is)%nsurface),stat=threadinfo%err)
call heap (threadinfo,'surface(is)%r','define_surface',size(surface(is)%x),'dp',+1)
allocate (surface(is)%s(surface(is)%nsurface),stat=threadinfo%err)
call heap (threadinfo,'surface(is)%s','define_surface',size(surface(is)%x),'dp',+1)
allocate (surface(is)%icon(3,surface(is)%nt),stat=threadinfo%err)
call heap (threadinfo,'surface(is)%icon','define_surface',size(surface(is)%icon),'int',+1)
read (9) (surface(is)%r(i), & read (9) (surface(is)%r(i), &
surface(is)%s(i), & surface(is)%s(i), &
......
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