72 |
integer startrec |
integer startrec |
73 |
integer endrec |
integer endrec |
74 |
|
|
75 |
|
_RL nwetc3d |
76 |
|
|
77 |
character*(max_len_prec) record |
character*(max_len_prec) record |
78 |
character*(max_len_mbuf) msgbuf |
character*(max_len_mbuf) msgbuf |
79 |
|
|
88 |
namelist /ctrl_nml/ |
namelist /ctrl_nml/ |
89 |
& xx_theta_file, |
& xx_theta_file, |
90 |
& xx_salt_file, |
& xx_salt_file, |
91 |
|
& xx_tr1_file, |
92 |
& xx_tauu_file, |
& xx_tauu_file, |
93 |
& xx_tauv_file, |
& xx_tauv_file, |
94 |
& xx_sflux_file, |
& xx_sflux_file, |
116 |
c-- Set default values. |
c-- Set default values. |
117 |
xx_theta_file = ' ' |
xx_theta_file = ' ' |
118 |
xx_salt_file = ' ' |
xx_salt_file = ' ' |
119 |
|
xx_tr1_file = ' ' |
120 |
xx_tauu_file = ' ' |
xx_tauu_file = ' ' |
121 |
xx_tauv_file = ' ' |
xx_tauv_file = ' ' |
122 |
xx_sflux_file = ' ' |
xx_sflux_file = ' ' |
261 |
_END_MASTER( mythid ) |
_END_MASTER( mythid ) |
262 |
#endif /* ALLOW_TAUV0_CONTROL */ |
#endif /* ALLOW_TAUV0_CONTROL */ |
263 |
|
|
264 |
|
cph( |
265 |
|
cph index 7-10 reserved for atmos. state, |
266 |
|
cph index 11-14 reserved for open boundaries, |
267 |
|
cph index 15-16 reserved for mixing coeff. |
268 |
|
cph index 17 reserved for passive tracer TR1 |
269 |
|
cph index 18,19 reserved for sst, sss |
270 |
|
cph) |
271 |
|
|
272 |
|
#ifdef ALLOW_TR10_CONTROL |
273 |
|
_BEGIN_MASTER( mythid ) |
274 |
|
ncvarindex(17) = 117 |
275 |
|
ncvarrecs(17) = 1 |
276 |
|
ncvarxmax(17) = snx |
277 |
|
ncvarymax(17) = sny |
278 |
|
ncvarnrmax(17) = nr |
279 |
|
ncvargrd(17) = 'c' |
280 |
|
_END_MASTER( mythid ) |
281 |
|
#endif /* ALLOW_TR10_CONTROL */ |
282 |
|
|
283 |
#ifdef ALLOW_SST0_CONTROL |
#ifdef ALLOW_SST0_CONTROL |
284 |
_BEGIN_MASTER( mythid ) |
_BEGIN_MASTER( mythid ) |
285 |
ncvarindex(7) = 107 |
ncvarindex(18) = 118 |
286 |
ncvarrecs(7) = 1 |
ncvarrecs(18) = 1 |
287 |
ncvarxmax(7) = snx |
ncvarxmax(18) = snx |
288 |
ncvarymax(7) = sny |
ncvarymax(18) = sny |
289 |
ncvarnrmax(7) = 1 |
ncvarnrmax(18) = 1 |
290 |
ncvargrd(7) = 'c' |
ncvargrd(18) = 'c' |
291 |
_END_MASTER( mythid ) |
_END_MASTER( mythid ) |
292 |
#endif /* ALLOW_SST0_CONTROL */ |
#endif /* ALLOW_SST0_CONTROL */ |
293 |
|
|
294 |
#ifdef ALLOW_SSS0_CONTROL |
#ifdef ALLOW_SSS0_CONTROL |
295 |
_BEGIN_MASTER( mythid ) |
_BEGIN_MASTER( mythid ) |
296 |
ncvarindex(8) = 108 |
ncvarindex(19) = 119 |
297 |
ncvarrecs(8) = 1 |
ncvarrecs(19) = 1 |
298 |
ncvarxmax(8) = snx |
ncvarxmax(19) = snx |
299 |
ncvarymax(8) = sny |
ncvarymax(19) = sny |
300 |
ncvarnrmax(8) = 1 |
ncvarnrmax(19) = 1 |
301 |
ncvargrd(8) = 'c' |
ncvargrd(19) = 'c' |
302 |
_END_MASTER( mythid ) |
_END_MASTER( mythid ) |
303 |
#endif /* ALLOW_SSS0_CONTROL */ |
#endif /* ALLOW_SSS0_CONTROL */ |
304 |
|
|
382 |
endif |
endif |
383 |
enddo |
enddo |
384 |
|
|
385 |
|
cph( |
386 |
|
print *, 'ph-wet 1: nvarlength = ', nvarlength |
387 |
|
print *, 'ph-wet 2: surface wet C = ', nwetctile(1,1,1) |
388 |
|
print *, 'ph-wet 3: surface wet W = ', nwetwtile(1,1,1) |
389 |
|
print *, 'ph-wet 4: surface wet S = ', nwetstile(1,1,1) |
390 |
|
nwetc3d = 0 |
391 |
|
do k = 1, Nr |
392 |
|
nwetc3d = nwetc3d + nwetctile(1,1,k) |
393 |
|
end do |
394 |
|
print *, 'ph-wet 5: 3D center wet points = ', nwetc3d |
395 |
|
do i = 1, 6 |
396 |
|
print *, 'ph-wet 6: no recs for i = ', i, ncvarrecs(i) |
397 |
|
end do |
398 |
|
print *, 'ph-wet 7: ', |
399 |
|
& 2*nwetc3d + |
400 |
|
& ncvarrecs(3)*nwetctile(1,1,1) + |
401 |
|
& ncvarrecs(4)*nwetctile(1,1,1) + |
402 |
|
& ncvarrecs(5)*nwetwtile(1,1,1) + |
403 |
|
& ncvarrecs(6)*nwetstile(1,1,1) |
404 |
|
cph) |
405 |
|
|
406 |
c |
c |
407 |
c Summation of wet point counters |
c Summation of wet point counters |
408 |
c |
c |