Changeset 6225 for branches/2014/dev_r4704_NOC5_MPP_BDY_UPDATE/NEMOGCM/NEMO/TOP_SRC/PISCES/trcwri_pisces.F90
- Timestamp:
- 2016-01-08T10:35:19+01:00 (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2014/dev_r4704_NOC5_MPP_BDY_UPDATE/NEMOGCM/NEMO/TOP_SRC/PISCES/trcwri_pisces.F90
r4292 r6225 21 21 PUBLIC trc_wri_pisces 22 22 23 # include "top_substitute.h90" 23 !!---------------------------------------------------------------------- 24 !! NEMO/TOP 3.3 , NEMO Consortium (2010) 25 !! $Id: trcnam.F90 5836 2015-10-26 14:49:40Z cetlod $ 26 !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 27 !!---------------------------------------------------------------------- 24 28 CONTAINS 25 29 … … 30 34 !! ** Purpose : output passive tracers fields 31 35 !!--------------------------------------------------------------------- 32 CHARACTER (len=20) :: cltra 33 REAL(wp) :: zrfact 34 INTEGER :: jn 36 CHARACTER (len=20) :: cltra 37 REAL(wp) :: zfact 38 INTEGER :: ji, jj, jk, jn 39 REAL(wp), DIMENSION(jpi,jpj) :: zdic, zo2min, zdepo2min 35 40 !!--------------------------------------------------------------------- 36 41 … … 40 45 DO jn = jp_pcs0, jp_pcs1 41 46 cltra = TRIM( ctrcnm(jn) ) ! short title for tracer 42 IF( lk_vvl ) THEN 43 CALL iom_put( cltra, trn(:,:,:,jn) * fse3t_n(:,:,:) ) 44 ELSE 45 CALL iom_put( cltra, trn(:,:,:,jn) ) 46 ENDIF 47 CALL iom_put( cltra, trn(:,:,:,jn) * zrfact ) 47 CALL iom_put( cltra, trn(:,:,:,jn) ) 48 48 END DO 49 49 #else 50 50 DO jn = jp_pcs0, jp_pcs1 51 z rfact = 1.0e+652 IF( jn == jpno3 .OR. jn == jpnh4 ) z rfact = rno3 * 1.0e+653 IF( jn == jppo4 ) z rfact = po4r * 1.0e+651 zfact = 1.0e+6 52 IF( jn == jpno3 .OR. jn == jpnh4 ) zfact = rno3 * 1.0e+6 53 IF( jn == jppo4 ) zfact = po4r * 1.0e+6 54 54 cltra = TRIM( ctrcnm(jn) ) ! short title for tracer 55 IF( lk_vvl ) THEN 56 CALL iom_put( cltra, trn(:,:,:,jn) * fse3t_n(:,:,:) * zrfact ) 57 ELSE 58 CALL iom_put( cltra, trn(:,:,:,jn) * zrfact ) 59 ENDIF 55 IF( iom_use( cltra ) ) CALL iom_put( cltra, trn(:,:,:,jn) * zfact ) 60 56 END DO 57 58 IF( iom_use( "INTDIC" ) ) THEN ! DIC content in kg/m2 59 zdic(:,:) = 0. 60 DO jk = 1, jpkm1 61 zdic(:,:) = zdic(:,:) + trn(:,:,jk,jpdic) * e3t_n(:,:,jk) * tmask(:,:,jk) * 12. 62 ENDDO 63 CALL iom_put( 'INTDIC', zdic ) 64 ENDIF 65 ! 66 IF( iom_use( "O2MIN" ) .OR. iom_use ( "ZO2MIN" ) ) THEN ! Oxygen minimum concentration and depth 67 zo2min (:,:) = trn(:,:,1,jpoxy) * tmask(:,:,1) 68 zdepo2min(:,:) = gdepw_n(:,:,1) * tmask(:,:,1) 69 DO jk = 2, jpkm1 70 DO jj = 1, jpj 71 DO ji = 1, jpi 72 IF( tmask(ji,jj,jk) == 1 ) then 73 IF( trn(ji,jj,jk,jpoxy) < zo2min(ji,jj) ) then 74 zo2min (ji,jj) = trn(ji,jj,jk,jpoxy) 75 zdepo2min(ji,jj) = gdepw_n(ji,jj,jk) 76 ENDIF 77 ENDIF 78 END DO 79 END DO 80 END DO 81 ! 82 CALL iom_put('O2MIN' , zo2min ) ! oxygen minimum concentration 83 CALL iom_put('ZO2MIN', zdepo2min ) ! depth of oxygen minimum concentration 84 ! 85 ENDIF 61 86 #endif 62 87 !
Note: See TracChangeset
for help on using the changeset viewer.