Ignore:
Timestamp:
2018-11-15T17:27:18+01:00 (23 months ago)
Author:
smasson
Message:

dev_r10164_HPC09_ESIWACE_PREP_MERGE: action 2: add generic glob_min/max/sum and locmin/max, complete timing and report (including bdy and icb), see #2133

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2018/dev_r10164_HPC09_ESIWACE_PREP_MERGE/src/TOP/C14/trcwri_c14.F90

    r10070 r10314  
    9696     
    9797      IF( iom_use("AtmC14") ) THEN 
    98          zarea = glob_sum( e1e2t(:,:) )           ! global ocean surface 
    99          ztemp = glob_sum( c14sbc(:,:) * e1e2t(:,:) ) 
     98         zarea = glob_sum( 'trcwri_c14', e1e2t(:,:) )           ! global ocean surface 
     99         ztemp = glob_sum( 'trcwri_c14', c14sbc(:,:) * e1e2t(:,:) ) 
    100100         ztemp = ( ztemp / zarea - 1._wp ) * 1000._wp 
    101101         CALL iom_put( "AtmC14" , ztemp )   ! Global atmospheric DeltaC14 [permil] 
    102102      ENDIF 
    103103      IF( iom_use("K_C14") ) THEN 
    104          ztemp = glob_sum ( exch_c14(:,:) * e1e2t(:,:) ) 
     104         ztemp = glob_sum ( 'trcwri_c14', exch_c14(:,:) * e1e2t(:,:) ) 
    105105         ztemp = rsiyea * ztemp / zarea 
    106106         CALL iom_put( "K_C14" , ztemp )   ! global mean exchange velocity for C14/C ratio [m/yr] 
    107107      ENDIF 
    108108      IF( iom_use("K_CO2") ) THEN 
    109          zarea = glob_sum( e1e2t(:,:) )           ! global ocean surface 
    110          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(:,:) ) 
    111111         ztemp = 360000._wp * ztemp / zarea       ! cm/h units: directly comparable with literature 
    112112         CALL iom_put( "K_CO2", ztemp )  !  global mean CO2 piston velocity [cm/hr] 
    113113      ENDIF 
    114114      IF( iom_use("C14Inv") ) THEN 
    115          ztemp = glob_sum( trn(:,:,:,jp_c14) * cvol(:,:,:) ) 
     115         ztemp = glob_sum( 'trcwri_c14', trn(:,:,:,jp_c14) * cvol(:,:,:) ) 
    116116         ztemp = atomc14 * xdicsur * ztemp 
    117117         CALL iom_put( "C14Inv", ztemp )  !  Radiocarbon ocean inventory [10^26 atoms] 
Note: See TracChangeset for help on using the changeset viewer.