Ignore:
Timestamp:
2019-02-01T17:27:20+01:00 (18 months ago)
Author:
dford
Message:

Implement biogeochemistry assimilation for FABM-ERSEM.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/UKMO/AMM15_v3_6_STABLE_package_collate_BGC_DA/NEMOGCM/NEMO/TOP_SRC/FABM/trcini_fabm.F90

    r10390 r10622  
    105105      jp_fabm_o2o  = fabm_state_index( 'O2_o' ) 
    106106      jp_fabm_o3c  = fabm_state_index( 'O3_c' ) 
    107       jp_fabm_o3a  = fabm_state_index( 'O3_bioalk' ) 
     107      jp_fabm_o3ba = fabm_state_index( 'O3_bioalk' ) 
     108      jp_fabm_r4n  = fabm_state_index( 'R4_n' ) 
     109      jp_fabm_r4c  = fabm_state_index( 'R4_c' ) 
     110      jp_fabm_r4p  = fabm_state_index( 'R4_p' ) 
     111      jp_fabm_r6n  = fabm_state_index( 'R6_n' ) 
     112      jp_fabm_r6c  = fabm_state_index( 'R6_c' ) 
     113      jp_fabm_r6p  = fabm_state_index( 'R6_p' ) 
     114      jp_fabm_r6s  = fabm_state_index( 'R6_s' ) 
     115      jp_fabm_r8n  = fabm_state_index( 'R8_n' ) 
     116      jp_fabm_r8c  = fabm_state_index( 'R8_c' ) 
     117      jp_fabm_r8p  = fabm_state_index( 'R8_p' ) 
     118      jp_fabm_r8s  = fabm_state_index( 'R8_s' ) 
    108119 
    109120      ! Get indexes for select diagnostic variables 
    110       jp_fabm_o3ph = fabm_diag_index( 'O3_pH' ) 
    111       jp_fabm_o3pc = fabm_diag_index( 'O3_pCO2' ) 
     121      jp_fabm_o3ta  = fabm_diag_index( 'O3_TA' ) 
     122      jp_fabm_o3ph  = fabm_diag_index( 'O3_pH' ) 
     123      jp_fabm_o3pc  = fabm_diag_index( 'O3_pCO2' ) 
     124      jp_fabm_pgrow = fabm_diag_index( 'p_grow_sum_result' ) 
     125      jp_fabm_ploss = fabm_diag_index( 'p_loss_sum_result' ) 
     126       
     127      MLD_MAX(:,:)   = 0.0 
     128      PGROW_AVG(:,:) = 0.0 
     129      PLOSS_AVG(:,:) = 0.0 
     130      PHYT_AVG(:,:)  = 0.0 
    112131 
    113132      IF (lwp) THEN 
     
    445464      END DO 
    446465      IF (fabm_state_index == -1) THEN 
    447          CALL ctl_stop( 'Could not find '//TRIM(state_name)//' state variable' ) 
     466         CALL ctl_warn( 'Could not find '//TRIM(state_name)//' state variable' ) 
    448467      ELSE 
    449468         IF (lwp) WRITE(numout,*) 'Index for '//TRIM(state_name)//' is: ', fabm_state_index 
     
    477496      END DO 
    478497      IF (fabm_diag_index == -1) THEN 
    479          CALL ctl_stop( 'Could not find '//TRIM(diag_name)//' diagnostic' ) 
     498         CALL ctl_warn( 'Could not find '//TRIM(diag_name)//' diagnostic' ) 
    480499      ELSE 
    481500         IF (lwp) WRITE(numout,*) 'Index for '//TRIM(diag_name)//' is: ', fabm_diag_index 
Note: See TracChangeset for help on using the changeset viewer.