Changeset 9435 for branches/UKMO/dev_r5518_GO6_package_asm_3d_bgc/NEMOGCM/NEMO/OPA_SRC/ASM/asmbgc.F90

Timestamp:

2018-03-27T15:20:18+02:00 (6 years ago)

Author:

dford

Message:

Allow assimilation of PFT chlorophyll.

File:

• branches/UKMO/dev_r5518_GO6_package_asm_3d_bgc/NEMOGCM/NEMO/OPA_SRC/ASM/asmbgc.F90 (modified) (2 diffs)

branches/UKMO/dev_r5518_GO6_package_asm_3d_bgc/NEMOGCM/NEMO/OPA_SRC/ASM/asmbgc.F90

@@ -249 +249 @@
       ENDIF
 
+      IF ( ln_phytobal .AND.                                      &
+         & ( ( ln_slchlnoninc .AND.( .NOT. ln_slchldiainc ) ).OR. &
+         &   ( ln_slchldiainc .AND.( .NOT. ln_slchlnoninc ) ) ) ) THEN
+         CALL ctl_stop( ' Cannot calculate phytoplankton balancing increments', &
+            &           ' unless assimilating all model PFTs,')
+      ENDIF
+
       IF ( ln_slphytotinc .AND. ( ln_slphydiainc .OR. ln_slphynoninc ) ) THEN
          CALL ctl_stop( ' Can only assimilate total or PFT surface phytoplankton carbon, not both' )
+      ENDIF
+
+      IF ( ln_phytobal .AND.                                      &
+         & ( ( ln_slphynoninc .AND.( .NOT. ln_slphydiainc ) ).OR. &
+         &   ( ln_slphydiainc .AND.( .NOT. ln_slphynoninc ) ) ) ) THEN
+         CALL ctl_stop( ' Cannot calculate phytoplankton balancing increments', &
+            &           ' unless assimilating all model PFTs,')
       ENDIF
 
@@ -905 +919 @@
 
          ! Un-log any log increments for passing to balancing routines
+         ! Remember that two sets of non-log increments should not be
+         ! expected to be in the same ratio as their log equivalents
+         
          ! Total chlorophyll
          IF ( ln_slchltotinc ) THEN