Changeset 7156 for branches/NERC
- Timestamp:
- 2016-10-28T13:04:10+02:00 (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/NERC/dev_r5518_NOC_MEDUSA_Stable/NEMOGCM/NEMO/TOP_SRC/MEDUSA/trcbio_medusa.F90
r7017 r7156 305 305 REAL(wp) :: f_fco2w, f_BetaD, f_rhosw, f_opres, f_insitut, f_pco2atm, f_fco2atm 306 306 REAL(wp) :: f_schmidtco2, f_kwco2, f_K0, f_co2starair, f_dpco2, f_kwo2 307 !! jpalm 14-07-2016: convert CO2flux diag in kg/m2/s 308 REAL, PARAMETER :: weight_CO2_mol = 0.0440095 309 REAL, PARAMETER :: secs_in_day = 86400.0 310 REAL, PARAMETER :: CO2flux_conv = weight_CO2_mol / secs_in_day 311 307 !! jpalm 14-07-2016: convert CO2flux diag from mmol-C/m2/s to kg-CO2/m2/s 308 REAL, PARAMETER :: CO2flux_conv = 44.0095 * 1.e-6 !! CO2 weight per mol (kg/mmol) 312 309 !! 313 310 INTEGER :: iters … … 1499 1496 # endif 1500 1497 !! 1498 !! JPALM -- 28-10-2016 -- save co2flux in mmol/m2/s in 2D field 1499 !! before it is converted to mmol/m3/d 1500 IF ( lk_iomput ) THEN 1501 IF( med_diag%CO2FLUX%dgsave ) THEN 1502 f_co2flux2d(ji,jj) = f_co2flux 1503 ENDIF 1504 ENDIF 1505 !! 1501 1506 !! mmol/m2/s -> mmol/m3/d; correct for sea-ice; divide through by layer thickness 1502 1507 f_co2flux = (1. - fr_i(ji,jj)) * f_co2flux * 86400. / fthk … … 1578 1583 !! 1579 1584 !! store 2D outputs 1585 !! JPALM -- 28-10-2016 -- removed co2flux from the list below 1586 !! it is now a little bit above, because a 1587 !! conversion is done in between. 1580 1588 IF ( lk_iomput ) THEN 1581 1589 IF( med_diag%ATM_PCO2%dgsave ) THEN … … 1584 1592 IF( med_diag%OCN_PCO2%dgsave ) THEN 1585 1593 f_pco2w2d(ji,jj) = f_pco2w 1586 ENDIF1587 IF( med_diag%CO2FLUX%dgsave ) THEN1588 f_co2flux2d(ji,jj) = f_co2flux1589 1594 ENDIF 1590 1595 IF( med_diag%TCO2%dgsave ) THEN … … 4371 4376 ENDIF 4372 4377 IF( med_diag%CO2FLUX%dgsave ) THEN 4373 !! convert from m ol/m2/day to kg/m2/s4378 !! convert from mmol-C/m2/s to kg-CO2/m2/s 4374 4379 f_co2flux2d(:,:) = CO2flux_conv * f_co2flux2d(:,:) 4375 4380 CALL lbc_lnk(f_co2flux2d(:,:),'T',1. )
Note: See TracChangeset
for help on using the changeset viewer.