Changeset 460 for trunk/NEMO/OPA_SRC/DIA/diaptr.F90
- Timestamp:
- 2006-05-10T19:10:25+02:00 (18 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/OPA_SRC/DIA/diaptr.F90
r406 r460 14 14 !! : flux array; Generic interface: ptr_vj_3d, ptr_vj_2d 15 15 !!---------------------------------------------------------------------- 16 !! History : 17 !! 9.0 ! 03-09 (C. Talandir, G. Madec) Original code 18 !! 9.0 ! 06-01 (A. Biastoch) Allow sub-basins computation 19 !!---------------------------------------------------------------------- 16 20 !! * Modules used 17 21 USE oce ! ocean dynamics and active tracers … … 22 26 USE dianam 23 27 USE phycst 28 USE ioipsl ! NetCDF IPSL library 29 USE daymod 24 30 25 31 IMPLICIT NONE … … 89 95 !! ** Action : - p_fval: i-k-mean poleward flux of pva 90 96 !! 91 !! History :92 !! 9.0 ! 03-09 (G. Madec) Original code93 97 !!---------------------------------------------------------------------- 94 98 !! * arguments … … 98 102 !! * local declarations 99 103 INTEGER :: ji, jj, jk ! dummy loop arguments 100 #if ! defined key_agrif101 INTEGER :: ijpj = jpj ! ???102 #else103 104 INTEGER :: ijpj ! ??? 104 #endif105 105 REAL(wp),DIMENSION(jpj) :: & 106 106 p_fval ! function value 107 107 !!-------------------------------------------------------------------- 108 #if defined key_agrif 108 109 109 ijpj = jpj 110 #endif111 112 110 p_fval(:) = 0.e0 113 111 DO jk = 1, jpkm1 … … 137 135 !! ** Action : - p_fval: i-k-mean poleward flux of pva 138 136 !! 139 !! History :140 !! 9.0 ! 03-09 (G. Madec) Original code141 137 !!---------------------------------------------------------------------- 142 138 !! * arguments … … 146 142 !! * local declarations 147 143 INTEGER :: ji,jj ! dummy loop arguments 148 #if ! defined key_agrif149 INTEGER :: ijpj = jpj ! ???150 #else151 144 INTEGER :: ijpj ! ??? 152 #endif153 145 REAL(wp),DIMENSION(jpj) :: & 154 146 p_fval ! function value 155 147 !!-------------------------------------------------------------------- 156 #if defined key_agrif 148 157 149 ijpj = jpj 158 #endif159 160 150 p_fval(:) = 0.e0 161 151 DO jj = 2, jpjm1 … … 182 172 !! ** Action : - p_fval: i-k-mean poleward flux of pva 183 173 !! 184 !! History :185 !! 9.0 ! 03-09 (G. Madec) Original code186 174 !!---------------------------------------------------------------------- 187 175 !! * arguments … … 204 192 DO jj = 2, jpjm1 205 193 DO ji = fs_2, fs_jpim1 206 p_fval(jj,jk) = p_fval(jj,jk) + pva(ji,jj,jk) * e1v(ji,jj) * fse3v(ji,jj,jk) &207 & * tmask_i(ji,jj+1) * tmask_i(ji,jj)194 p_fval(jj,jk) = p_fval(jj,jk) + pva(ji,jj,jk) * e1v(ji,jj) * fse3v(ji,jj,jk) & 195 & * tmask_i(ji,jj+1) * tmask_i(ji,jj) 208 196 END DO 209 197 END DO 210 198 END DO 211 199 212 IF(lk_mpp) 213 ish(1) = jpj*jpk ; ish2(1)=jpj ; ish2(2)=jpk214 zwork(:)= RESHAPE( p_fval, ish )215 CALL mpp_sum( zwork, jpj*jpk )216 p_fval(:,:)= RESHAPE( zwork,ish2)200 IF(lk_mpp) THEN 201 ish(1) = jpj*jpk ; ish2(1) = jpj ; ish2(2) = jpk 202 zwork(:)= RESHAPE( p_fval, ish ) 203 CALL mpp_sum( zwork, jpj*jpk ) 204 p_fval(:,:)= RESHAPE( zwork, ish2 ) 217 205 END IF 218 206 … … 231 219 !! ** Action : - p_fval: i-k-mean poleward flux of pva 232 220 !! 233 !! History :234 !! 9.0 ! 03-09 (G. Madec) Original code235 !! 9.0 ! 06-01 (A. Biastoch) Allow sub-basins computation236 221 !!---------------------------------------------------------------------- 237 222 !! * arguments … … 260 245 END DO 261 246 p_fval(:,:) = p_fval(:,:) * 0.5 262 IF(lk_mpp) 263 ish(1) = jpj*jpk ; ish2(1)=jpj ; ish2(2)=jpk264 zwork(:)= RESHAPE( p_fval, ish )265 CALL mpp_sum( zwork, jpj*jpk )247 IF(lk_mpp) THEN 248 ish(1) = jpj*jpk ; ish2(1) = jpj ; ish2(2) = jpk 249 zwork(:)= RESHAPE( p_fval, ish ) 250 CALL mpp_sum( zwork, jpj*jpk ) 266 251 p_fval(:,:)= RESHAPE(zwork,ish2) 267 252 END IF … … 491 476 !! ** Purpose : Initialization, namelist read 492 477 !! 493 !! ** Method : 494 !! 495 !! ** input : Namlist namptr 496 !! 497 !! ** Action : 498 !! 499 !! history : 500 !! 9.0 ! 03-08 (Autor Names) Original code 501 !!---------------------------------------------------------------------- 502 !! * local declarations 478 !!---------------------------------------------------------------------- 503 479 REAL(wp), DIMENSION(jpi,jpj,jpk) :: z_1 ! temporary workspace 504 480 … … 529 505 END SUBROUTINE dia_ptr_init 530 506 531 !!---------------------------------------------------------------------532 !! Default option : NetCDF file533 !!---------------------------------------------------------------------534 507 535 508 SUBROUTINE dia_ptr_wri( kt ) … … 541 514 !! ** Method : NetCDF file 542 515 !! 543 !! History : 544 !! 9.0 ! 03-09 (G. Madec) Original code 545 !!---------------------------------------------------------------------- 546 USE ioipsl ! NetCDF IPSL library 547 USE daymod 548 516 !!---------------------------------------------------------------------- 549 517 !! * Arguments 550 518 INTEGER, INTENT(in) :: kt ! ocean time-step index … … 635 603 CALL histbeg(clhstnam, 1, zfoo, jpj, zphi, & 636 604 1, 1, 1, jpj, 0, zjulian, zdt, nhoridz, numptr, domain_id=nidom ) 637 ! Vertical grids : gdept , gdepw605 ! Vertical grids : gdept_0, gdepw_0 638 606 CALL histvert( numptr, "deptht", "Vertical T levels", & 639 "m", jpk, gdept , ndepidzt )607 "m", jpk, gdept_0, ndepidzt ) 640 608 CALL histvert( numptr, "depthw", "Vertical W levels", & 641 "m", jpk, gdepw , ndepidzw )609 "m", jpk, gdepw_0, ndepidzw ) 642 610 643 611 ! Zonal mean T and S
Note: See TracChangeset
for help on using the changeset viewer.