Changeset 10622 for branches/UKMO/AMM15_v3_6_STABLE_package_collate_BGC_DA/NEMOGCM/NEMO/TOP_SRC/FABM/trcini_fabm.F90

Timestamp:

2019-02-01T17:27:20+01:00 (5 years ago)

Author:

dford

Message:

Implement biogeochemistry assimilation for FABM-ERSEM.

File:

• branches/UKMO/AMM15_v3_6_STABLE_package_collate_BGC_DA/NEMOGCM/NEMO/TOP_SRC/FABM/trcini_fabm.F90 (modified) (3 diffs)

branches/UKMO/AMM15_v3_6_STABLE_package_collate_BGC_DA/NEMOGCM/NEMO/TOP_SRC/FABM/trcini_fabm.F90

@@ -105 +105 @@
       jp_fabm_o2o  = fabm_state_index( 'O2_o' )
       jp_fabm_o3c  = fabm_state_index( 'O3_c' )
-      jp_fabm_o3a  = fabm_state_index( 'O3_bioalk' )
+      jp_fabm_o3ba = fabm_state_index( 'O3_bioalk' )
+      jp_fabm_r4n  = fabm_state_index( 'R4_n' )
+      jp_fabm_r4c  = fabm_state_index( 'R4_c' )
+      jp_fabm_r4p  = fabm_state_index( 'R4_p' )
+      jp_fabm_r6n  = fabm_state_index( 'R6_n' )
+      jp_fabm_r6c  = fabm_state_index( 'R6_c' )
+      jp_fabm_r6p  = fabm_state_index( 'R6_p' )
+      jp_fabm_r6s  = fabm_state_index( 'R6_s' )
+      jp_fabm_r8n  = fabm_state_index( 'R8_n' )
+      jp_fabm_r8c  = fabm_state_index( 'R8_c' )
+      jp_fabm_r8p  = fabm_state_index( 'R8_p' )
+      jp_fabm_r8s  = fabm_state_index( 'R8_s' )
 
       ! Get indexes for select diagnostic variables
-      jp_fabm_o3ph = fabm_diag_index( 'O3_pH' )
-      jp_fabm_o3pc = fabm_diag_index( 'O3_pCO2' )
+      jp_fabm_o3ta  = fabm_diag_index( 'O3_TA' )
+      jp_fabm_o3ph  = fabm_diag_index( 'O3_pH' )
+      jp_fabm_o3pc  = fabm_diag_index( 'O3_pCO2' )
+      jp_fabm_pgrow = fabm_diag_index( 'p_grow_sum_result' )
+      jp_fabm_ploss = fabm_diag_index( 'p_loss_sum_result' )
+      
+      MLD_MAX(:,:)   = 0.0
+      PGROW_AVG(:,:) = 0.0
+      PLOSS_AVG(:,:) = 0.0
+      PHYT_AVG(:,:)  = 0.0
 
       IF (lwp) THEN
@@ -445 +464 @@
       END DO
       IF (fabm_state_index == -1) THEN
-         CALL ctl_stop( 'Could not find '//TRIM(state_name)//' state variable' )
+         CALL ctl_warn( 'Could not find '//TRIM(state_name)//' state variable' )
       ELSE
          IF (lwp) WRITE(numout,*) 'Index for '//TRIM(state_name)//' is: ', fabm_state_index
@@ -477 +496 @@
       END DO
       IF (fabm_diag_index == -1) THEN
-         CALL ctl_stop( 'Could not find '//TRIM(diag_name)//' diagnostic' )
+         CALL ctl_warn( 'Could not find '//TRIM(diag_name)//' diagnostic' )
       ELSE
          IF (lwp) WRITE(numout,*) 'Index for '//TRIM(diag_name)//' is: ', fabm_diag_index