Changeset 7479
- Timestamp:
- 2016-12-08T14:25:38+01:00 (8 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
r6857 r7479 1446 1446 & MEDUSA_FILL_FLT 1447 1447 #elif defined key_fabm 1448 !USE ??? ! ERSEM chlorophyll, fCO2 and pCO2 1448 USE fabm 1449 USE trcsms_fabm, ONLY: model ! ERSEM chlorophyll, fCO2 and pCO2 1449 1450 #endif 1450 1451 #if defined key_spm … … 1487 1488 maskpco2 ! array for special pco2 mask 1488 1489 INTEGER :: jn ! loop index 1490 #if defined key_fabm 1491 INTEGER :: chl_index 1492 REAL(wp), DIMENSION(jpi,jpj,jpk) :: logchl_3d 1493 #endif 1489 1494 CHARACTER(LEN=20) :: datestr=" ",timestr=" " 1490 1495 … … 1599 1604 1600 1605 #if defined key_hadocc 1601 logchl(:,:) 1606 logchl(:,:) = HADOCC_CHL(:,:,1) ! (not log) chlorophyll from HadOCC 1602 1607 #elif defined key_medusa && defined key_foam_medusa 1603 logchl(:,:) 1608 logchl(:,:) = MEDUSA_CHL(:,:,1) ! (not log) chlorophyll from HadOCC 1604 1609 #elif defined key_fabm 1605 !logchl(:,:) = ??? ! (not log) chlorophyll from ERSEM 1606 CALL ctl_stop( ' Trying to run logchl observation operator', & 1607 & ' but not properly implemented for FABM-ERSEM yet' ) 1610 DO jn = 1, SIZE(model%diagnostic_variables) 1611 IF (TRIM(model%diagnostic_variables(jn)%standard_variable%name) == 'total_chlorophyll') THEN 1612 chl_index = jn 1613 EXIT 1614 ENDIF 1615 END DO 1616 logchl_3d(:,:,:) = fabm_get_bulk_diagnostic_data(model, chl_index) 1617 logchl(:,:) = logchl_3d(:,:,1) 1608 1618 #else 1609 1619 CALL ctl_stop( ' Trying to run logchl observation operator', & … … 1758 1768 INTEGER :: jfbini 1759 1769 CHARACTER(LEN=20) :: datestr=" ",timestr=" " 1760 CHARACTER(LEN= 10) :: cdtmp1770 CHARACTER(LEN=20) :: cdtmp 1761 1771 !----------------------------------------------------------------------- 1762 1772 ! Depending on switches call various observation output routines
Note: See TracChangeset
for help on using the changeset viewer.