Changeset 8653


Ignore:
Timestamp:
2017-10-23T18:32:44+02:00 (3 years ago)
Author:
dford
Message:

Updates for getting biogeochemical fields from MEDUSA. See GMED ticket 344.

Location:
branches/UKMO/dev_r5518_obs_oper_update/NEMOGCM/NEMO/OPA_SRC/OBS
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • branches/UKMO/dev_r5518_obs_oper_update/NEMOGCM/NEMO/OPA_SRC/OBS/diaobs.F90

    r8157 r8653  
    646646#elif defined key_medusa && defined key_foam_medusa 
    647647      USE trc, ONLY :  &           ! MEDUSA chlorophyll, fCO2 and pCO2 
    648          & MEDUSA_CHL, & 
    649          & MEDUSA_FCO2, & 
    650          & MEDUSA_PCO2, & 
    651          & MEDUSA_FILL_FLT 
     648         & trn 
     649      USE par_medusa, ONLY: & 
     650         & jpchn, & 
     651         & jpchd 
     652#if defined key_roam 
     653      USE sms_medusa, ONLY: & 
     654         & f2_pco2w, & 
     655         & f2_fco2w 
     656#endif 
    652657#elif defined key_fabm 
    653658      USE fabm 
     
    792797               zsurfvar(:,:) = HADOCC_CHL(:,:,1)    ! (not log) chlorophyll from HadOCC 
    793798#elif defined key_medusa && defined key_foam_medusa 
    794                zsurfvar(:,:) = MEDUSA_CHL(:,:,1)    ! (not log) chlorophyll from HadOCC 
     799               ! Add non-diatom and diatom surface chlorophyll from MEDUSA 
     800               zsurfvar(:,:) = trn(:,:,1,jpchn) + trn(:,:,1,jpchd) 
    795801#elif defined key_fabm 
    796802               chl_3d(:,:,:) = fabm_get_bulk_diagnostic_data(model, jp_fabmdia_chltot) 
     
    829835                     &           ' as HADOCC_FCO2(:,:) == HADOCC_FILL_FLT' ) 
    830836               ENDIF 
    831 #elif defined key_medusa && defined key_foam_medusa 
    832                zsurfmask(:,:) = MEDUSA_FCO2(:,:)    ! fCO2 from MEDUSA 
    833                IF ( ( MINVAL( MEDUSA_FCO2 ) == MEDUSA_FILL_FLT ) .AND. & 
    834                   & ( MAXVAL( MEDUSA_FCO2 ) == MEDUSA_FILL_FLT ) ) THEN 
    835                   zsurfvar(:,:) = obfillflt 
    836                   zsurfmask(:,:) = 0 
    837                   CALL ctl_warn( ' MEDUSA fCO2 values masked out for observation operator', & 
    838                      &           ' as MEDUSA_FCO2(:,:) == MEDUSA_FILL_FLT' ) 
    839                ENDIF 
     837#elif defined key_medusa && defined key_foam_medusa && defined key_roam 
     838               zsurfvar(:,:) = f2_fco2w(:,:) 
    840839#elif defined key_fabm 
    841840               ! First, get pCO2 from FABM 
     
    875874                     &           ' as HADOCC_PCO2(:,:) == HADOCC_FILL_FLT' ) 
    876875               ENDIF 
    877 #elif defined key_medusa && defined key_foam_medusa 
    878                zsurfvar(:,:) = MEDUSA_PCO2(:,:)    ! pCO2 from MEDUSA 
    879                IF ( ( MINVAL( MEDUSA_PCO2 ) == MEDUSA_FILL_FLT ) .AND. & 
    880                   & ( MAXVAL( MEDUSA_PCO2 ) == MEDUSA_FILL_FLT ) ) THEN 
    881                   zsurfvar(:,:) = obfillflt 
    882                   zsurfmask(:,:) = 0 
    883                   CALL ctl_warn( ' MEDUSA pCO2 values masked out for observation operator', & 
    884                      &           ' as MEDUSA_PCO2(:,:) == MEDUSA_FILL_FLT' ) 
    885                ENDIF 
     876#elif defined key_medusa && defined key_foam_medusa && defined key_roam 
     877               zsurfvar(:,:) = f2_pco2w(:,:) 
    886878#elif defined key_fabm 
    887879               pco2_3d(:,:,:) = fabm_get_bulk_diagnostic_data(model, jp_fabm_o3pc) 
  • branches/UKMO/dev_r5518_obs_oper_update/NEMOGCM/NEMO/OPA_SRC/OBS/obs_oper.F90

    r7992 r8653  
    982982            surfdataqc%rmod(jobs,1) = zext(1) 
    983983         ENDIF 
     984          
     985         IF ( zext(1) == obfillflt ) THEN 
     986            ! If the observation value is a fill value, set QC flag to bad 
     987            surfdataqc%nqc(jobs) = 4 
     988         ENDIF 
    984989 
    985990      END DO 
Note: See TracChangeset for help on using the changeset viewer.