Changeset 10314 for NEMO/branches/2018/dev_r10164_HPC09_ESIWACE_PREP_MERGE/src/TOP/C14/trcwri_c14.F90
- Timestamp:
- 2018-11-15T17:27:18+01:00 (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2018/dev_r10164_HPC09_ESIWACE_PREP_MERGE/src/TOP/C14/trcwri_c14.F90
r10070 r10314 96 96 97 97 IF( iom_use("AtmC14") ) THEN 98 zarea = glob_sum( e1e2t(:,:) ) ! global ocean surface99 ztemp = glob_sum( c14sbc(:,:) * e1e2t(:,:) )98 zarea = glob_sum( 'trcwri_c14', e1e2t(:,:) ) ! global ocean surface 99 ztemp = glob_sum( 'trcwri_c14', c14sbc(:,:) * e1e2t(:,:) ) 100 100 ztemp = ( ztemp / zarea - 1._wp ) * 1000._wp 101 101 CALL iom_put( "AtmC14" , ztemp ) ! Global atmospheric DeltaC14 [permil] 102 102 ENDIF 103 103 IF( iom_use("K_C14") ) THEN 104 ztemp = glob_sum ( exch_c14(:,:) * e1e2t(:,:) )104 ztemp = glob_sum ( 'trcwri_c14', exch_c14(:,:) * e1e2t(:,:) ) 105 105 ztemp = rsiyea * ztemp / zarea 106 106 CALL iom_put( "K_C14" , ztemp ) ! global mean exchange velocity for C14/C ratio [m/yr] 107 107 ENDIF 108 108 IF( iom_use("K_CO2") ) THEN 109 zarea = glob_sum( e1e2t(:,:) ) ! global ocean surface110 ztemp = glob_sum ( exch_co2(:,:) * e1e2t(:,:) )109 zarea = glob_sum( 'trcwri_c14', e1e2t(:,:) ) ! global ocean surface 110 ztemp = glob_sum ( 'trcwri_c14', exch_co2(:,:) * e1e2t(:,:) ) 111 111 ztemp = 360000._wp * ztemp / zarea ! cm/h units: directly comparable with literature 112 112 CALL iom_put( "K_CO2", ztemp ) ! global mean CO2 piston velocity [cm/hr] 113 113 ENDIF 114 114 IF( iom_use("C14Inv") ) THEN 115 ztemp = glob_sum( trn(:,:,:,jp_c14) * cvol(:,:,:) )115 ztemp = glob_sum( 'trcwri_c14', trn(:,:,:,jp_c14) * cvol(:,:,:) ) 116 116 ztemp = atomc14 * xdicsur * ztemp 117 117 CALL iom_put( "C14Inv", ztemp ) ! Radiocarbon ocean inventory [10^26 atoms]
Note: See TracChangeset
for help on using the changeset viewer.