Changeset 8455


Ignore:
Timestamp:
2017-08-23T12:57:45+02:00 (3 years ago)
Author:
dford
Message:

Update getting chlorophyll, pCO2 and fCO2 from MEDUSA.

File:
1 edited

Legend:

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

    r8157 r8455  
    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               zsurfvar(:,:) = trn(:,:,1,jpchn) + trn(:,:,1,jpchd) 
    795800#elif defined key_fabm 
    796801               chl_3d(:,:,:) = fabm_get_bulk_diagnostic_data(model, jp_fabmdia_chltot) 
     
    829834                     &           ' as HADOCC_FCO2(:,:) == HADOCC_FILL_FLT' ) 
    830835               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 
     836#elif defined key_medusa && defined key_foam_medusa && defined key_roam 
     837               zsurfvar(:,:) = f2_fco2w(:,:) 
    840838#elif defined key_fabm 
    841839               ! First, get pCO2 from FABM 
     
    875873                     &           ' as HADOCC_PCO2(:,:) == HADOCC_FILL_FLT' ) 
    876874               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 
     875#elif defined key_medusa && defined key_foam_medusa && defined key_roam 
     876               zsurfvar(:,:) = f2_pco2w(:,:) 
    886877#elif defined key_fabm 
    887878               pco2_3d(:,:,:) = fabm_get_bulk_diagnostic_data(model, jp_fabm_o3pc) 
Note: See TracChangeset for help on using the changeset viewer.