- Timestamp:
- 2017-01-31T13:06:14+01:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/UKMO/dev_r4650_general_vert_coord_obsoper_surf_bgc/NEMOGCM/NEMO/OPA_SRC/OBS/diaobs.F90
r7489 r7625 1423 1423 & tmask, umask, vmask 1424 1424 USE phycst, ONLY : & ! Physical constants 1425 & rday 1425 & rday, & 1426 & rt0 1426 1427 USE oce, ONLY : & ! Ocean dynamics and tracers variables 1427 1428 & tsn, & … … 1455 1456 USE par_spm, ONLY: & ! ERSEM/SPM sediments 1456 1457 & jp_spm 1458 USE trc, ONLY : & 1459 & trn 1457 1460 #endif 1458 1461 IMPLICIT NONE … … 1493 1496 #if defined key_fabm 1494 1497 REAL(wp), DIMENSION(jpi,jpj,jpk) :: logchl_3d 1498 REAL(wp), DIMENSION(jpi,jpj,jpk) :: pco2_3d 1495 1499 #endif 1496 1500 CHARACTER(LEN=20) :: datestr=" ",timestr=" " … … 1674 1678 ENDIF 1675 1679 #elif defined key_fabm 1676 !fco2(:,:) = ??? ! fCO2 from ERSEM 1677 CALL ctl_stop( ' Trying to run fco2 observation operator', & 1678 & ' but not properly implemented for FABM-ERSEM yet' ) 1680 ! First, get pCO2 from FABM 1681 pco2_3d(:,:,:) = fabm_get_bulk_diagnostic_data(model, jp_fabm_o3pc) 1682 pco2(:,:) = pco2_3d(:,:,1) 1683 ! Now, convert pCO2 to fCO2, based on SST in K 1684 ! This follows the method used by SOCAT to convert pCO2 to fCO2, see: 1685 ! http://www.socat.info/Uniform_format_surface_fCO2_database.pdf 1686 fco2(:,:) = pco2(:,:) * EXP((-1636.75 + & 1687 & 12.0408 * (tsn(:,:,1,jp_tem)+rt0) - & 1688 & 0.0327957 * (tsn(:,:,1,jp_tem)+rt0)*(tsn(:,:,1,jp_tem)+rt0) + & 1689 & 0.0000316528 * (tsn(:,:,1,jp_tem)+rt0)*(tsn(:,:,1,jp_tem)+rt0)*(tsn(:,:,1,jp_tem)+rt0) + & 1690 & 2.0 * (57.7 - 0.118 * (tsn(:,:,1,jp_tem)+rt0))) / & 1691 & (82.0578 * (tsn(:,:,1,jp_tem)+rt0))) 1679 1692 #else 1680 1693 CALL ctl_stop( ' Trying to run fco2 observation operator', & … … 1710 1723 ENDIF 1711 1724 #elif defined key_fabm 1712 !pco2(:,:) = ??? ! pCO2 from ERSEM 1713 CALL ctl_stop( ' Trying to run pCO2 observation operator', & 1714 & ' but not properly implemented for FABM-ERSEM yet' ) 1725 pco2_3d(:,:,:) = fabm_get_bulk_diagnostic_data(model, jp_fabm_o3pc) 1726 pco2(:,:) = pco2_3d(:,:,1) 1715 1727 #else 1716 1728 CALL ctl_stop( ' Trying to run pCO2 observation operator', &
Note: See TracChangeset
for help on using the changeset viewer.