- Timestamp:
- 2016-10-21T17:38:13+02:00 (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2016/dev_r7012_ROBUST5_CNRS/NEMOGCM/NEMO/TOP_SRC/C14/trcsms_c14.F90
r7041 r7068 46 46 ! freshwater fluxes which should not impact the C14/C ratio 47 47 ! 48 ! => Delta-C14= ( trn(...jp_c14 0) -1)*1000.48 ! => Delta-C14= ( trn(...jp_c14) -1)*1000. 49 49 !! 50 50 !!---------------------------------------------------------------------- … … 121 121 ! Flux of C-14 from air-to-sea; units: (C14/C ratio) x m/s 122 122 ! already masked 123 qtr_c14(:,:) = exch_c14(:,:) * ( c14sbc(:,:) - trb(:,:,1,jp_c14 0) )123 qtr_c14(:,:) = exch_c14(:,:) * ( c14sbc(:,:) - trb(:,:,1,jp_c14) ) 124 124 125 125 ! cumulation of air-to-sea flux at each time step 126 126 qint_c14(:,:) = qint_c14(:,:) + qtr_c14(:,:) * rdttrc 127 127 ! 128 ! Add the surface flux to the trend of jp_c14 0128 ! Add the surface flux to the trend of jp_c14 129 129 DO jj = 1, jpj 130 130 DO ji = 1, jpi 131 tra(ji,jj,1,jp_c14 0) = tra(ji,jj,1,jp_c140) + qtr_c14(ji,jj) / e3t_n(ji,jj,1)131 tra(ji,jj,1,jp_c14) = tra(ji,jj,1,jp_c14) + qtr_c14(ji,jj) / e3t_n(ji,jj,1) 132 132 END DO 133 133 END DO 134 134 ! 135 ! Computation of decay effects on jp_c14 0135 ! Computation of decay effects on jp_c14 136 136 DO jk = 1, jpk 137 137 DO jj = 1, jpj 138 138 DO ji = 1, jpi 139 139 ! 140 tra(ji,jj,jk,jp_c14 0) = tra(ji,jj,jk,jp_c140) - rlam14 * trb(ji,jj,jk,jp_c140) * tmask(ji,jj,jk)140 tra(ji,jj,jk,jp_c14) = tra(ji,jj,jk,jp_c14) - rlam14 * trb(ji,jj,jk,jp_c14) * tmask(ji,jj,jk) 141 141 ! 142 142 END DO … … 159 159 ENDIF 160 160 161 IF( l_trdtrc ) CALL trd_trc( tra(:,:,:,jp_c14 0), jp_c14, jptra_sms, kt ) ! save trends161 IF( l_trdtrc ) CALL trd_trc( tra(:,:,:,jp_c14), 1, jptra_sms, kt ) ! save trends 162 162 ! 163 163 IF( nn_timing == 1 ) CALL timing_stop('trc_sms_c14')
Note: See TracChangeset
for help on using the changeset viewer.