41 |
c == local variables == |
c == local variables == |
42 |
|
|
43 |
_RL fcref |
_RL fcref |
44 |
_RL fcpert |
_RL fcpertplus, fcpertminus |
45 |
_RL xxmemo_ref |
_RL xxmemo_ref |
46 |
_RL xxmemo_pert |
_RL xxmemo_pert |
47 |
_RL gfd |
_RL gfd |
94 |
call print_message( msgbuf, standardmessageunit, |
call print_message( msgbuf, standardmessageunit, |
95 |
& SQUEEZE_RIGHT , mythid) |
& SQUEEZE_RIGHT , mythid) |
96 |
|
|
97 |
write(msgbuf,'(a,a)') |
write(msgbuf,'(a,6(1x,a15))') |
98 |
&' I X(I) FC FC1-FC2/EPS ', |
& 'grdchk output: ', 'I', 'ITILEPOS', 'JTILEPOS', |
99 |
&' GRAD(FC) DIFFERENCE/RATIO ' |
& 'LAYER', 'X(I)', 'X(I)+/-EPS' |
100 |
|
call print_message( msgbuf, standardmessageunit, |
101 |
|
& SQUEEZE_RIGHT , mythid) |
102 |
|
write(msgbuf,'(a,6(1x,a15))') |
103 |
|
& 'grdchk output: ', 'FC', 'FC1', 'FC2', 'FC1-FC2/(2*EPS)', |
104 |
|
& 'GRAD(FC)', '1-FDGRD/ADGRD' |
105 |
call print_message( msgbuf, standardmessageunit, |
call print_message( msgbuf, standardmessageunit, |
106 |
& SQUEEZE_RIGHT , mythid) |
& SQUEEZE_RIGHT , mythid) |
107 |
|
|
116 |
xxmemo_ref = xxmemref (i) |
xxmemo_ref = xxmemref (i) |
117 |
xxmemo_pert = xxmempert (i) |
xxmemo_pert = xxmempert (i) |
118 |
adxxmemo = adxxmem (i) |
adxxmemo = adxxmem (i) |
119 |
fcref = fcrmem (i) |
fcref = fcrmem (i) |
120 |
fcpert = fcpmem (i) |
fcpertplus = fcppmem (i) |
121 |
|
fcpertminus = fcppmem (i) |
122 |
gfd = gfdmem (i) |
gfd = gfdmem (i) |
123 |
ratio = ratiomem(i) |
ratio = ratiomem(i) |
124 |
itile = bimem (i) |
itile = bimem (i) |
129 |
icomp = icompmem(i) |
icomp = icompmem(i) |
130 |
ierr = ierrmem (i) |
ierr = ierrmem (i) |
131 |
|
|
132 |
write(msgbuf,'(A,4(I8),2(1x,D15.9))') |
write(msgbuf,'(A,4(I16),2(1x,D15.9))') |
133 |
& 'grdchk output: ', i, itilepos, jtilepos, layer, |
& 'grdchk output: ', i, itilepos, jtilepos, layer, |
134 |
& xxmemo_ref, xxmemo_pert |
& xxmemo_ref, xxmemo_pert |
135 |
call print_message( msgbuf, standardmessageunit, |
call print_message( msgbuf, standardmessageunit, |
136 |
& SQUEEZE_RIGHT , mythid) |
& SQUEEZE_RIGHT , mythid) |
137 |
if ( ierr .eq. 0 ) then |
if ( ierr .eq. 0 ) then |
138 |
write(msgbuf,'(A,5(1x,D15.9))') |
write(msgbuf,'(A,6(1x,D15.9))') |
139 |
& 'grdchk output: ', |
& 'grdchk output: ', |
140 |
& fcref, fcpert, gfd, adxxmemo, ratio |
& fcref, fcpertplus, fcpertminus, gfd, adxxmemo, ratio |
141 |
call print_message( msgbuf, standardmessageunit, |
call print_message( msgbuf, standardmessageunit, |
142 |
& SQUEEZE_RIGHT , mythid) |
& SQUEEZE_RIGHT , mythid) |
143 |
else |
else |