Changeset 5208 for branches/2014/dev_r4650_UKMO11_restart_functionality/NEMOGCM/NEMO/TOP_SRC/PISCES/trcwri_pisces.F90
- Timestamp:
- 2015-04-13T15:08:59+02:00 (9 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2014/dev_r4650_UKMO11_restart_functionality/NEMOGCM/NEMO/TOP_SRC/PISCES/trcwri_pisces.F90
r4292 r5208 30 30 !! ** Purpose : output passive tracers fields 31 31 !!--------------------------------------------------------------------- 32 CHARACTER (len=20) :: cltra 33 REAL(wp) :: zrfact 34 INTEGER :: jn 32 CHARACTER (len=20) :: cltra 33 REAL(wp) :: zfact 34 INTEGER :: ji, jj, jk, jn 35 REAL(wp), DIMENSION(jpi,jpj) :: zdic, zo2min, zdepo2min 35 36 !!--------------------------------------------------------------------- 36 37 … … 40 41 DO jn = jp_pcs0, jp_pcs1 41 42 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 ) 43 CALL iom_put( cltra, trn(:,:,:,jn) ) 48 44 END DO 49 45 #else 50 46 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+647 zfact = 1.0e+6 48 IF( jn == jpno3 .OR. jn == jpnh4 ) zfact = rno3 * 1.0e+6 49 IF( jn == jppo4 ) zfact = po4r * 1.0e+6 54 50 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 51 IF( iom_use( cltra ) ) CALL iom_put( cltra, trn(:,:,:,jn) * zfact ) 60 52 END DO 53 54 IF( iom_use( "INTDIC" ) ) THEN ! DIC content in kg/m2 55 zdic(:,:) = 0. 56 DO jk = 1, jpkm1 57 zdic(:,:) = zdic(:,:) + trn(:,:,jk,jpdic) * fse3t(:,:,jk) * tmask(:,:,jk) * 12. 58 ENDDO 59 CALL iom_put( 'INTDIC', zdic ) 60 ENDIF 61 ! 62 IF( iom_use( "O2MIN" ) .OR. iom_use ( "ZO2MIN" ) ) THEN ! Oxygen minimum concentration and depth 63 zo2min (:,:) = trn(:,:,1,jpoxy) * tmask(:,:,1) 64 zdepo2min(:,:) = fsdepw(:,:,1) * tmask(:,:,1) 65 DO jk = 2, jpkm1 66 DO jj = 1, jpj 67 DO ji = 1, jpi 68 IF( tmask(ji,jj,jk) == 1 ) then 69 IF( trn(ji,jj,jk,jpoxy) < zo2min(ji,jj) ) then 70 zo2min (ji,jj) = trn(ji,jj,jk,jpoxy) 71 zdepo2min(ji,jj) = fsdepw(ji,jj,jk) 72 ENDIF 73 ENDIF 74 END DO 75 END DO 76 END DO 77 ! 78 CALL iom_put('O2MIN' , zo2min ) ! oxygen minimum concentration 79 CALL iom_put('ZO2MIN', zdepo2min ) ! depth of oxygen minimum concentration 80 ! 81 ENDIF 61 82 #endif 62 83 !
Note: See TracChangeset
for help on using the changeset viewer.