[7041] | 1 | MODULE trcwri_c14 |
---|
| 2 | !!====================================================================== |
---|
| 3 | !! *** MODULE trcwri *** |
---|
| 4 | !! MY_SRC : Additional outputs for C14 tracers |
---|
| 5 | !!====================================================================== |
---|
| 6 | !! History : 1.0 ! 2009-05 (C. Ethe) Original code |
---|
| 7 | !! History : 2.0 ! 2015 (A. Mouchet) adapted code for C14 |
---|
| 8 | !!---------------------------------------------------------------------- |
---|
| 9 | #if defined key_top && defined key_iomput |
---|
| 10 | !!---------------------------------------------------------------------- |
---|
| 11 | !! trc_wri_c14 : outputs of ventilation fields |
---|
| 12 | !!---------------------------------------------------------------------- |
---|
| 13 | USE oce_trc ! Ocean variables |
---|
| 14 | USE trc ! passive tracers common variables |
---|
| 15 | USE iom ! I/O manager |
---|
| 16 | USE sms_c14 |
---|
| 17 | |
---|
| 18 | IMPLICIT NONE |
---|
| 19 | PRIVATE |
---|
| 20 | |
---|
| 21 | PUBLIC trc_wri_c14 |
---|
| 22 | ! |
---|
| 23 | ! Standard ratio: 1.176E-12 ; Avogadro's nbr = 6.022E+23 at/mol ; bomb C14 traditionally reported as 1.E+26 atoms |
---|
| 24 | REAL(wp), PARAMETER :: atomc14 = 1.176 * 6.022E-15 ! conversion factor |
---|
| 25 | |
---|
| 26 | |
---|
| 27 | CONTAINS |
---|
| 28 | |
---|
| 29 | SUBROUTINE trc_wri_c14 |
---|
| 30 | !!--------------------------------------------------------------------- |
---|
| 31 | !! *** ROUTINE trc_wri_c14 *** |
---|
| 32 | !! |
---|
| 33 | !! ** Purpose : output additional C14 tracers fields |
---|
| 34 | !!--------------------------------------------------------------------- |
---|
| 35 | CHARACTER (len=20) :: cltra ! short title for tracer |
---|
| 36 | INTEGER :: ji,jj,jk,jn ! dummy loop indexes |
---|
| 37 | REAL(wp) :: zage,zarea,ztemp ! temporary |
---|
[9125] | 38 | REAL(wp), ALLOCATABLE, DIMENSION(:,:) :: zres, z2d ! temporary storage 2D |
---|
| 39 | REAL(wp), ALLOCATABLE, DIMENSION(:,:,:) :: z3d , zz3d ! temporary storage 3D |
---|
[7041] | 40 | !!--------------------------------------------------------------------- |
---|
| 41 | |
---|
| 42 | ! write the tracer concentrations in the file |
---|
| 43 | ! --------------------------------------- |
---|
[7068] | 44 | cltra = TRIM( ctrcnm(jp_c14) ) ! short title for tracer |
---|
| 45 | CALL iom_put( cltra, trn(:,:,:,jp_c14) ) |
---|
[7041] | 46 | |
---|
| 47 | ! compute and write the tracer diagnostic in the file |
---|
| 48 | ! --------------------------------------- |
---|
| 49 | |
---|
| 50 | IF( iom_use("DeltaC14") .OR. iom_use("C14Age") .OR. iom_use("RAge") ) THEN |
---|
| 51 | ! |
---|
[9125] | 52 | ALLOCATE( z2d(jpi,jpj), zres(jpi,jpj) ) |
---|
| 53 | ALLOCATE( z3d(jpi,jpj,jpk), zz3d(jpi,jpj,jpk) ) |
---|
[7041] | 54 | ! |
---|
| 55 | zage = -1._wp / rlam14 / rsiyea ! factor for radioages in year |
---|
| 56 | z3d(:,:,:) = 1._wp |
---|
| 57 | zz3d(:,:,:) = 0._wp |
---|
| 58 | ! |
---|
| 59 | DO jk = 1, jpkm1 |
---|
| 60 | DO jj = 1, jpj |
---|
| 61 | DO ji = 1, jpi |
---|
| 62 | IF( tmask(ji,jj,jk) > 0._wp) THEN |
---|
[7068] | 63 | z3d (ji,jj,jk) = trn(ji,jj,jk,jp_c14) |
---|
[7041] | 64 | zz3d(ji,jj,jk) = LOG( z3d(ji,jj,jk) ) |
---|
| 65 | ENDIF |
---|
| 66 | ENDDO |
---|
| 67 | ENDDO |
---|
| 68 | ENDDO |
---|
| 69 | zres(:,:) = z3d(:,:,1) |
---|
| 70 | |
---|
| 71 | ! Reservoir age [yr] |
---|
| 72 | z2d(:,:) =0._wp |
---|
| 73 | jk = 1 |
---|
| 74 | DO jj = 1, jpj |
---|
| 75 | DO ji = 1, jpi |
---|
| 76 | ztemp = zres(ji,jj) / c14sbc(ji,jj) |
---|
| 77 | IF( ztemp > 0._wp .AND. tmask(ji,jj,jk) > 0._wp ) z2d(ji,jj) = LOG( ztemp ) |
---|
| 78 | ENDDO |
---|
| 79 | ENDDO |
---|
| 80 | ! |
---|
| 81 | z3d(:,:,:) = 1.d03 * ( z3d(:,:,:) - 1._wp ) |
---|
| 82 | CALL iom_put( "DeltaC14" , z3d(:,:,:) ) ! Delta C14 [permil] |
---|
| 83 | CALL iom_put( "C14Age" , zage * zz3d(:,:,:) ) ! Radiocarbon age [yr] |
---|
| 84 | |
---|
| 85 | CALL iom_put( "qtr_c14", rsiyea * qtr_c14(:,:) ) ! Radiocarbon surf flux [./m2/yr] |
---|
| 86 | CALL iom_put( "qint_c14" , qint_c14 ) ! cumulative flux [./m2] |
---|
| 87 | CALL iom_put( "RAge" , zage * z2d(:,:) ) ! Reservoir age [yr] |
---|
| 88 | ! |
---|
[9125] | 89 | DEALLOCATE( z2d, zres, z3d, zz3d ) |
---|
[7041] | 90 | ! |
---|
| 91 | ENDIF |
---|
| 92 | ! |
---|
| 93 | ! 0-D fields |
---|
| 94 | ! |
---|
| 95 | CALL iom_put( "AtmCO2", co2sbc ) ! global atmospheric CO2 [ppm] |
---|
| 96 | |
---|
| 97 | IF( iom_use("AtmC14") ) THEN |
---|
[10425] | 98 | zarea = glob_sum( 'trcwri_c14', e1e2t(:,:) ) ! global ocean surface |
---|
| 99 | ztemp = glob_sum( 'trcwri_c14', c14sbc(:,:) * e1e2t(:,:) ) |
---|
[7041] | 100 | ztemp = ( ztemp / zarea - 1._wp ) * 1000._wp |
---|
| 101 | CALL iom_put( "AtmC14" , ztemp ) ! Global atmospheric DeltaC14 [permil] |
---|
| 102 | ENDIF |
---|
| 103 | IF( iom_use("K_C14") ) THEN |
---|
[10425] | 104 | ztemp = glob_sum ( 'trcwri_c14', exch_c14(:,:) * e1e2t(:,:) ) |
---|
[7041] | 105 | ztemp = rsiyea * ztemp / zarea |
---|
| 106 | CALL iom_put( "K_C14" , ztemp ) ! global mean exchange velocity for C14/C ratio [m/yr] |
---|
| 107 | ENDIF |
---|
| 108 | IF( iom_use("K_CO2") ) THEN |
---|
[10425] | 109 | zarea = glob_sum( 'trcwri_c14', e1e2t(:,:) ) ! global ocean surface |
---|
| 110 | ztemp = glob_sum ( 'trcwri_c14', exch_co2(:,:) * e1e2t(:,:) ) |
---|
[7041] | 111 | ztemp = 360000._wp * ztemp / zarea ! cm/h units: directly comparable with literature |
---|
| 112 | CALL iom_put( "K_CO2", ztemp ) ! global mean CO2 piston velocity [cm/hr] |
---|
| 113 | ENDIF |
---|
| 114 | IF( iom_use("C14Inv") ) THEN |
---|
[10425] | 115 | ztemp = glob_sum( 'trcwri_c14', trn(:,:,:,jp_c14) * cvol(:,:,:) ) |
---|
[7041] | 116 | ztemp = atomc14 * xdicsur * ztemp |
---|
| 117 | CALL iom_put( "C14Inv", ztemp ) ! Radiocarbon ocean inventory [10^26 atoms] |
---|
| 118 | END IF |
---|
| 119 | ! |
---|
| 120 | END SUBROUTINE trc_wri_c14 |
---|
| 121 | |
---|
| 122 | #else |
---|
| 123 | !!---------------------------------------------------------------------- |
---|
| 124 | !! Dummy module : No C14 tracer |
---|
| 125 | !!---------------------------------------------------------------------- |
---|
| 126 | PUBLIC trc_wri_c14 |
---|
| 127 | CONTAINS |
---|
| 128 | SUBROUTINE trc_wri_c14 ! Empty routine |
---|
| 129 | END SUBROUTINE trc_wri_c14 |
---|
| 130 | #endif |
---|
| 131 | |
---|
| 132 | !!---------------------------------------------------------------------- |
---|
[10067] | 133 | !! NEMO/TOP 4.0 , NEMO Consortium (2018) |
---|
[10070] | 134 | !! $Id$ |
---|
[10068] | 135 | !! Software governed by the CeCILL license (see ./LICENSE) |
---|
[7041] | 136 | !!====================================================================== |
---|
| 137 | END MODULE trcwri_c14 |
---|