1 |
|
C $Header$ |
2 |
|
C $Name$ |
3 |
|
|
4 |
#include "CPP_OPTIONS.h" |
#include "CPP_OPTIONS.h" |
5 |
|
|
60 |
jmax = sny |
jmax = sny |
61 |
imin = 1 |
imin = 1 |
62 |
imax = snx |
imax = snx |
63 |
doglobalread = .false. |
doglobalread = .false. |
64 |
ladinit = .false. |
ladinit = .false. |
65 |
equal = .true. |
equal = .true. |
66 |
if (equal) then |
if (equal) then |
67 |
fac = 1.d0 |
fac = 1.d0 |
68 |
else |
else |
74 |
DO j = jmin, jmax |
DO j = jmin, jmax |
75 |
DO i = imin, imax |
DO i = imin, imax |
76 |
DO k=1,Nr |
DO k=1,Nr |
77 |
adtheta(i,j,k,bi,bj) = |
adtheta(i,j,k,bi,bj) = |
78 |
& g_objf_state_final(i,j,bi,bj,k) |
& g_objf_state_final(i,j,bi,bj,k) |
79 |
adsalt(i,j,k,bi,bj) = |
adsalt(i,j,k,bi,bj) = |
80 |
& g_objf_state_final(i,j,bi,bj,1*Nr+k) |
& g_objf_state_final(i,j,bi,bj,1*Nr+k) |
81 |
aduvel(i,j,k,bi,bj) = |
aduvel(i,j,k,bi,bj) = |
82 |
& g_objf_state_final(i,j,bi,bj,2*Nr+k) |
& g_objf_state_final(i,j,bi,bj,2*Nr+k) |
83 |
advvel(i,j,k,bi,bj) = |
advvel(i,j,k,bi,bj) = |
84 |
& g_objf_state_final(i,j,bi,bj,3*Nr+k) |
& g_objf_state_final(i,j,bi,bj,3*Nr+k) |
85 |
END DO |
END DO |
86 |
adetan(i,j,bi,bj) = |
adetan(i,j,bi,bj) = |
87 |
& g_objf_state_final(i,j,bi,bj,4*Nr+1) |
& g_objf_state_final(i,j,bi,bj,4*Nr+1) |
88 |
END DO |
END DO |
89 |
END DO |
END DO |
103 |
end do |
end do |
104 |
end do |
end do |
105 |
il = ilnblnk(xx_etan_file) |
il = ilnblnk(xx_etan_file) |
106 |
write(unit=fnamegeneric(1:80),fmt='(2a,i10.10)') |
write(unit=fnamegeneric(1:80),fmt='(2a,i10.10)') |
107 |
& xx_etan_file(1:il),'.',optimcycle |
& xx_etan_file(1:il),'.',optimcycle |
108 |
call adactive_read_xy_loc( fnamegeneric,1,doglobalread,ladinit, |
call adactive_read_xy_loc( fnamegeneric,1,doglobalread,ladinit, |
109 |
& optimcycle,mythid,tmpfld2d ) |
& optimcycle,mythid,tmpfld2d ) |
121 |
end do |
end do |
122 |
end do |
end do |
123 |
il = ilnblnk(xx_vvel_file) |
il = ilnblnk(xx_vvel_file) |
124 |
write(unit=fnamegeneric(1:80),fmt='(2a,i10.10)') |
write(unit=fnamegeneric(1:80),fmt='(2a,i10.10)') |
125 |
& xx_vvel_file(1:il),'.',optimcycle |
& xx_vvel_file(1:il),'.',optimcycle |
126 |
call adactive_read_xyz( fnamegeneric,1,doglobalread,ladinit, |
call adactive_read_xyz( fnamegeneric,1,doglobalread,ladinit, |
127 |
& optimcycle,mythid,tmpfld3d ) |
& optimcycle,mythid,tmpfld3d ) |
139 |
end do |
end do |
140 |
end do |
end do |
141 |
il = ilnblnk(xx_uvel_file) |
il = ilnblnk(xx_uvel_file) |
142 |
write(unit=fnamegeneric(1:80),fmt='(2a,i10.10)') |
write(unit=fnamegeneric(1:80),fmt='(2a,i10.10)') |
143 |
& xx_uvel_file(1:il),'.',optimcycle |
& xx_uvel_file(1:il),'.',optimcycle |
144 |
call adactive_read_xyz( fnamegeneric,1,doglobalread,ladinit, |
call adactive_read_xyz( fnamegeneric,1,doglobalread,ladinit, |
145 |
& optimcycle,mythid,tmpfld3d ) |
& optimcycle,mythid,tmpfld3d ) |
157 |
end do |
end do |
158 |
end do |
end do |
159 |
il = ilnblnk(xx_salt_file) |
il = ilnblnk(xx_salt_file) |
160 |
write(unit=fnamegeneric(1:80),fmt='(2a,i10.10)') |
write(unit=fnamegeneric(1:80),fmt='(2a,i10.10)') |
161 |
& xx_salt_file(1:il),'.',optimcycle |
& xx_salt_file(1:il),'.',optimcycle |
162 |
call adactive_read_xyz_loc( fnamegeneric,1,doglobalread,ladinit, |
call adactive_read_xyz_loc( fnamegeneric,1,doglobalread,ladinit, |
163 |
& optimcycle,mythid,tmpfld3d ) |
& optimcycle,mythid,tmpfld3d ) |
175 |
end do |
end do |
176 |
end do |
end do |
177 |
il = ilnblnk(xx_theta_file) |
il = ilnblnk(xx_theta_file) |
178 |
write(unit=fnamegeneric(1:80),fmt='(2a,i10.10)') |
write(unit=fnamegeneric(1:80),fmt='(2a,i10.10)') |
179 |
& xx_theta_file(1:il),'.',optimcycle |
& xx_theta_file(1:il),'.',optimcycle |
180 |
call adactive_read_xyz_loc( fnamegeneric,1,doglobalread,ladinit, |
call adactive_read_xyz_loc( fnamegeneric,1,doglobalread,ladinit, |
181 |
& optimcycle,mythid,tmpfld3d ) |
& optimcycle,mythid,tmpfld3d ) |
182 |
|
|
183 |
#endif |
#endif |
184 |
|
|
185 |
end |
end |