Changeset 1329 for trunk/NEMO/TOP_SRC/PISCES/p4zsink.F90
- Timestamp:
- 2009-02-20T11:32:52+01:00 (15 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/TOP_SRC/PISCES/p4zsink.F90
r1264 r1329 279 279 END DO 280 280 281 # if defined key_trc_dia3d281 #if defined key_trc_diaadd 282 282 zrfact2 = 1.e3 * rfact2r 283 trc2d(:,:, jp_pcs0_2d + 4) = sinking (:,:,iksed+1) * zrfact2 284 trc2d(:,:, jp_pcs0_2d + 5) = sinking2(:,:,iksed+1) * zrfact2 285 trc2d(:,:, jp_pcs0_2d + 6) = sinkfer (:,:,iksed+1) * zrfact2 286 trc2d(:,:, jp_pcs0_2d + 8) = sinksil (:,:,iksed+1) * zrfact2 287 trc2d(:,:, jp_pcs0_2d + 9) = sinkcal (:,:,iksed+1) * zrfact2 288 trc3d(:,:,:,jp_pcs0_3d + 11) = sinking (:,:,:) * zrfact2 289 trc3d(:,:,:,jp_pcs0_3d + 12) = sinking2(:,:,:) * zrfact2 290 trc3d(:,:,:,jp_pcs0_3d + 13) = sinksil (:,:,:) * zrfact2 291 trc3d(:,:,:,jp_pcs0_3d + 14) = sinkcal (:,:,:) * zrfact2 292 trc3d(:,:,:,jp_pcs0_3d + 15) = znum3d (:,:,:) 293 trc3d(:,:,:,jp_pcs0_3d + 16) = wsbio3 (:,:,:) 294 trc3d(:,:,:,jp_pcs0_3d + 17) = wsbio4 (:,:,:) 295 # endif 283 DO jj = 1, jpj 284 DO ji = 1, jpi 285 trc2d(ji,jj, jp_pcs0_2d + 4) = sinking (ji,jj,iksed+1) * zrfact2 * tmask(ji,jj,1) 286 trc2d(ji,jj, jp_pcs0_2d + 5) = sinking2(ji,jj,iksed+1) * zrfact2 * tmask(ji,jj,1) 287 trc2d(ji,jj, jp_pcs0_2d + 6) = sinkfer (ji,jj,iksed+1) * zrfact2 * tmask(ji,jj,1) 288 trc2d(ji,jj, jp_pcs0_2d + 7) = sinksil (ji,jj,iksed+1) * zrfact2 * tmask(ji,jj,1) 289 trc2d(ji,jj, jp_pcs0_2d + 8) = sinkcal (ji,jj,iksed+1) * zrfact2 * tmask(ji,jj,1) 290 ENDDO 291 ENDDO 292 # if defined key_trc_dia3d 293 DO jk = 1, jpk 294 DO jj = 1, jpj 295 DO ji = 1, jpi 296 trc3d(ji,jj,jk,jp_pcs0_3d + 11) = sinking (ji,jj,jk) * zrfact2 * tmask(ji,jj,jk) 297 trc3d(ji,jj,jk,jp_pcs0_3d + 12) = sinking2(ji,jj,jk) * zrfact2 * tmask(ji,jj,jk) 298 trc3d(ji,jj,jk,jp_pcs0_3d + 13) = sinksil (ji,jj,jk) * zrfact2 * tmask(ji,jj,jk) 299 trc3d(ji,jj,jk,jp_pcs0_3d + 14) = sinkcal (ji,jj,jk) * zrfact2 * tmask(ji,jj,jk) 300 trc3d(ji,jj,jk,jp_pcs0_3d + 15) = znum3d (ji,jj,jk) * tmask(ji,jj,jk) 301 trc3d(ji,jj,jk,jp_pcs0_3d + 16) = wsbio3 (ji,jj,jk) * tmask(ji,jj,jk) 302 trc3d(ji,jj,jk,jp_pcs0_3d + 17) = wsbio4 (ji,jj,jk) * tmask(ji,jj,jk) 303 ENDDO 304 ENDDO 305 ENDDO 306 # endif 307 308 #endif 296 309 ! 297 310 IF(ln_ctl) THEN ! print mean trends (used for debugging) … … 598 611 # if defined key_trc_diaadd 599 612 zrfact2 = 1.e3 * rfact2r 600 trc2d(:,:, jp_pcs0_2d + 4) = sinking (:,:,iksed+1) * zrfact2 601 trc2d(:,:, jp_pcs0_2d + 5) = sinking2(:,:,iksed+1) * zrfact2 602 trc2d(:,:, jp_pcs0_2d + 6) = sinkfer (:,:,iksed+1) * zrfact2 603 trc2d(:,:, jp_pcs0_2d + 7) = sinkfer2(:,:,iksed+1) * zrfact2 604 trc2d(:,:, jp_pcs0_2d + 8) = sinksil (:,:,iksed+1) * zrfact2 605 trc2d(:,:, jp_pcs0_2d + 9) = sinkcal (:,:,iksed+1) * zrfact2 613 DO jj = 1, jpj 614 DO ji = 1, jpi 615 trc2d(ji,jj, jp_pcs0_2d + 4) = sinking (ji,jj,iksed+1) * zrfact2 * tmask(ji,jj,1) 616 trc2d(ji,jj, jp_pcs0_2d + 5) = sinking2(ji,jj,iksed+1) * zrfact2 * tmask(ji,jj,1) 617 trc2d(ji,jj, jp_pcs0_2d + 6) = sinkfer (ji,jj,iksed+1) * zrfact2 * tmask(ji,jj,1) 618 trc2d(ji,jj, jp_pcs0_2d + 7) = sinkfer2(ji,jj,iksed+1) * zrfact2 * tmask(ji,jj,1) 619 trc2d(ji,jj, jp_pcs0_2d + 8) = sinksil (ji,jj,iksed+1) * zrfact2 * tmask(ji,jj,1) 620 trc2d(ji,jj, jp_pcs0_2d + 9) = sinkcal (ji,jj,iksed+1) * zrfact2 * tmask(ji,jj,1) 621 ENDDO 622 ENDDO 606 623 # endif 607 624 !
Note: See TracChangeset
for help on using the changeset viewer.