1 |
heimbach |
1.1 |
*** ad_taf_output.f 2005-09-29 15:17:45.000000000 -0400 |
2 |
|
|
--- ad_taf_output.fmod 2005-09-29 15:40:24.000000000 -0400 |
3 |
|
|
*************** |
4 |
|
|
*** 78470,78492 **** |
5 |
|
|
close(unit=76) |
6 |
|
|
else |
7 |
|
|
idivbeg = nchklev_3 |
8 |
|
|
! idivend = nchklev_3+(-1) |
9 |
|
|
endif |
10 |
|
|
endif |
11 |
|
|
call mpi_bcast( idivbeg,1,mpi_integer,0,mpi_comm_world,iers ) |
12 |
|
|
call mpi_bcast( idivend,1,mpi_integer,0,mpi_comm_world,iers ) |
13 |
|
|
|
14 |
|
|
C---------------------------------------------- |
15 |
|
|
- C write divided adjoint control file |
16 |
|
|
- C---------------------------------------------- |
17 |
|
|
- if (iproc .eq. 0) then |
18 |
|
|
- open(unit=76,file='divided.ctrl',form='formatted') |
19 |
|
|
- idivaux = 2*idivend-idivbeg |
20 |
|
|
- write(unit=76,fmt=*) idivend,idivaux |
21 |
|
|
- close(unit=76) |
22 |
|
|
- endif |
23 |
|
|
- |
24 |
|
|
- C---------------------------------------------- |
25 |
|
|
C OPEN TAPE onetape |
26 |
|
|
C---------------------------------------------- |
27 |
|
|
call adopen ( mythid, |
28 |
|
|
--- 78470,78482 ---- |
29 |
|
|
close(unit=76) |
30 |
|
|
else |
31 |
|
|
idivbeg = nchklev_3 |
32 |
|
|
! idivend = MAX(nchklev_3-1,0) |
33 |
|
|
endif |
34 |
|
|
endif |
35 |
|
|
call mpi_bcast( idivbeg,1,mpi_integer,0,mpi_comm_world,iers ) |
36 |
|
|
call mpi_bcast( idivend,1,mpi_integer,0,mpi_comm_world,iers ) |
37 |
|
|
|
38 |
|
|
C---------------------------------------------- |
39 |
|
|
C OPEN TAPE onetape |
40 |
|
|
C---------------------------------------------- |
41 |
|
|
call adopen ( mythid, |
42 |
|
|
*************** |
43 |
|
|
*** 82624,82632 **** |
44 |
|
|
--- 82614,82642 ---- |
45 |
|
|
$1368 ) |
46 |
|
|
|
47 |
|
|
endif |
48 |
|
|
+ |
49 |
|
|
+ C---------------------------------------------- |
50 |
|
|
+ C update divided adjoint control file |
51 |
|
|
+ C---------------------------------------------- |
52 |
|
|
+ if (iproc .eq. 0) then |
53 |
|
|
+ open(unit=76,file='divided.ctrl',form='formatted') |
54 |
|
|
+ write(unit=76,fmt=*) ilev_3-1,idivend |
55 |
|
|
+ close(unit=76) |
56 |
|
|
+ endif |
57 |
|
|
+ |
58 |
|
|
end do |
59 |
|
|
|
60 |
|
|
C---------------------------------------------- |
61 |
|
|
+ C write divaided adjoint control file |
62 |
|
|
+ C---------------------------------------------- |
63 |
|
|
+ if (iproc .eq. 0) then |
64 |
|
|
+ open(unit=76,file='divided.ctrl',form='formatted') |
65 |
|
|
+ idivaux = MAX(2*idivend-idivbeg,0) |
66 |
|
|
+ write(unit=76,fmt=*) idivend,idivaux |
67 |
|
|
+ close(unit=76) |
68 |
|
|
+ endif |
69 |
|
|
+ |
70 |
|
|
+ C---------------------------------------------- |
71 |
|
|
C write snapshot |
72 |
|
|
C---------------------------------------------- |
73 |
|
|
if (idivend .ge. 1) then |