Changeset 3879
- Timestamp:
- 2013-04-22T11:00:51+02:00 (11 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2012/dev_v3_4_STABLE_2012/NEMOGCM/NEMO/TOP_SRC/PISCES/p4zflx.F90
r3294 r3879 112 112 ! then the first atmospheric CO2 record read is at years(1) 113 113 zyr_dec = REAL( nyear + nn_offset, wp ) + REAL( nday_year, wp ) / REAL( nyear_len(1), wp ) 114 jm = 2115 DO WHILE( jm <= nmaxrec .AND. years(jm -1) < zyr_dec .AND. years(jm) >=zyr_dec ) ; jm = jm + 1 ; END DO114 jm = 1 115 DO WHILE( jm <= nmaxrec .AND. years(jm) < zyr_dec ) ; jm = jm + 1 ; END DO 116 116 iind = jm ; iindm1 = jm - 1 117 117 zdco2dt = ( atcco2h(iind) - atcco2h(iindm1) ) / ( years(iind) - years(iindm1) + rtrn ) … … 201 201 t_oce_co2_flx = t_oce_co2_flx + glob_sum( oce_co2(:,:) ) ! Cumulative Total Flux of Carbon 202 202 IF( kt == nitend ) THEN 203 t_atm_co2_flx = glob_sum( satmco2(:,:) * patm(:,:) *e1e2t(:,:) ) ! Total atmospheric pCO2203 t_atm_co2_flx = glob_sum( satmco2(:,:) * e1e2t(:,:) ) ! Total atmospheric pCO2 204 204 ! 205 205 t_oce_co2_flx = (-1.) * t_oce_co2_flx * 12. / 1.e15 ! Conversion in PgC ; negative for out of the ocean
Note: See TracChangeset
for help on using the changeset viewer.