From 256b101d1d38661cbe4b2e99b5a047be6eaf5d91 Mon Sep 17 00:00:00 2001
From: Dave Whipp <dwhipp@dal.ca>
Date: Wed, 25 May 2011 15:22:58 +0000
Subject: [PATCH] Tweaked output so that all output variables are always
written, with zero values when the calculation is not performed
---
src/write_global_output.f90 | 19 ++++++++++++-------
1 file changed, 12 insertions(+), 7 deletions(-)
diff --git a/src/write_global_output.f90 b/src/write_global_output.f90
index 9caca12e..de0f3187 100644
--- a/src/write_global_output.f90
+++ b/src/write_global_output.f90
@@ -196,18 +196,23 @@ if (iproc.eq.0) then
cl%tag(i), &
i=1,cl%np)
- if (params%isobc) then
- ! write isostasy basal displacement array - dwhipp 11/09
- write (9) 2**params%levelmax_oct
- write (9) ((zi%zisodisp(i,j)+surface(osolve%nlsf)%sp01,&
- j=1,2**params%levelmax_oct+1),i=1,2**params%levelmax_oct+1)
+ ! write isostasy basal displacement array - dwhipp 11/09
+ write (9) 2**params%levelmax_oct
+ if (params%isobc) then
+ write (9) ((zi%zisodisp(i,j)+surface(osolve%nlsf)%sp01, &
+ j=1,2**params%levelmax_oct+1),i=1,2**params%levelmax_oct+1)
+ else
+ write (9) ((0.d0+surface(osolve%nlsf)%sp01,j=1,2**params%levelmax_oct+1), &
+ i=1,2**params%levelmax_oct+1)
end if
+ ! write nested model info - dwhipp 05/11
if (params%nest) then
- ! write nested model info - dwhipp 05/11
write(9) nest%sselemx,nest%sselemy,nest%sselemz,nest%xminls,nest%yminls,&
nest%zminls
- endif
+ else
+ write(9) 0.d0,0.d0,0.d0,0.d0,0.d0,0.d0
+ endif
close (9)
--
GitLab