Changeset 239 for trunk/NEMO/OPA_SRC/DIA/diaptr.F90
- Timestamp:
- 2005-03-22T18:55:20+01:00 (19 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/OPA_SRC/DIA/diaptr.F90
r190 r239 345 345 END SUBROUTINE dia_ptr_init 346 346 347 #if defined key_fdir348 !!---------------------------------------------------------------------349 !! 'key_fdir' direct access file350 !!---------------------------------------------------------------------351 352 SUBROUTINE dia_ptr_wri( kt )353 !!---------------------------------------------------------------------354 !! *** ROUTINE dia_ptr_wri ***355 !!356 !! ** Purpose : output of poleward fluxes357 !!358 !! ** Method : NetCDF file359 !!360 !! History :361 !! 9.0 ! 03-09 (G. Madec) Original code362 !!----------------------------------------------------------------------363 !! * Arguments364 INTEGER, INTENT(in) :: kt ! ocean time-step index365 REAL(wp), DIMENSION(jpj) :: zphi, zfoo366 !!----------------------------------------------------------------------367 368 IF( kt == nit000 ) THEN369 370 ! Reference latitude371 ! ------------------372 ! ! =======================373 IF( cp_cfg == "orca" ) THEN ! ORCA configurations374 ! ! =======================375 376 IF( jp_cfg == 05 ) iline = 192 ! i-line that passes near the North Pole377 IF( jp_cfg == 025 ) iline = 384 ! i-line that passes near the North Pole378 IF( jp_cfg == 2 ) iline = 48 ! i-line that passes near the North Pole379 IF( jp_cfg == 4 ) iline = 24 ! i-line that passes near the North Pole380 381 zphi(:) = 0.e0382 DO ji = mi0(iline), mi1(iline)383 zphi(:) = gphiv(ji,:) ! if iline is in the local domain384 END DO385 ! provide the correct zphi to all local domains386 IF( lk_mpp ) CALL mpp_sum( zphi, jpj )387 ! introduce arbitray northernmost grid point to avoid netcdf error388 DO jj=mj0(jpjglo), mj1(jpjglo)389 zphi(jj) = 2*zphi(jj-1)-zphi(jj-2)390 ENDDO391 392 ! ! =======================393 ELSE ! OTHER configurations394 ! ! =======================395 zphi(:) = gphiv(1,:) ! assume lat/lon coordinate, select the first i-line396 !397 ENDIF398 399 ! open the output file400 CALL ctlopn( numptr, 'opaptr.output', 'UNKNOWN', 'UNFORMATTED', 'SEQUENTIAL', 1, numout, lwp, 1 )401 402 ! header of output403 WRITE( numptr ) cexper, no, zdt, nf_ptr, jpj, jpk, zphi404 405 ENDIF406 407 IF( MOD( kt, nf_ptr ) == 0 ) THEN408 IF(lwp) WRITE( numptr ) kt, tn_jk, sn_jk, v_msf, &409 # if defined key_diaeiv410 & pht_eiv, pst_eiv, v_msf_eiv, &411 # endif412 & pht_adv, pht_ldf, pst_adv, pst_ldf413 ENDIF414 415 IF( kt == nitend ) CLOSE( numptr )416 417 END SUBROUTINE dia_ptr_wri418 419 #else420 347 !!--------------------------------------------------------------------- 421 348 !! Default option : NetCDF file … … 598 525 END SUBROUTINE dia_ptr_wri 599 526 600 #endif601 602 527 !!====================================================================== 603 528 END MODULE diaptr
Note: See TracChangeset
for help on using the changeset viewer.