Changeset 10661 for branches/UKMO/AMM15_v3_6_STABLE_package_collate_BGC_DA/NEMOGCM/NEMO/OPA_SRC/ASM/asmbgc.F90

Timestamp:

2019-02-12T18:08:39+01:00 (5 years ago)

Author:

dford

Message:

Add option to apply phosphate increments.

File:

• branches/UKMO/AMM15_v3_6_STABLE_package_collate_BGC_DA/NEMOGCM/NEMO/OPA_SRC/ASM/asmbgc.F90 (modified) (10 diffs)

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

@@ -135 +135 @@
    LOGICAL, PUBLIC :: ln_pno3inc     = .FALSE. !: No profile nitrate                       increment
    LOGICAL, PUBLIC :: ln_psi4inc     = .FALSE. !: No profile silicate                      increment
+   LOGICAL, PUBLIC :: ln_ppo4inc     = .FALSE. !: No profile phosphate                     increment
    LOGICAL, PUBLIC :: ln_pdicinc     = .FALSE. !: No profile dissolved inorganic carbon    increment
    LOGICAL, PUBLIC :: ln_palkinc     = .FALSE. !: No profile alkalinity                    increment
@@ -175 +176 @@
    REAL(wp), DIMENSION(:,:,:),   ALLOCATABLE :: pno3_bkginc     ! pno3 inc
    REAL(wp), DIMENSION(:,:,:),   ALLOCATABLE :: psi4_bkginc     ! psi4 inc
+   REAL(wp), DIMENSION(:,:,:),   ALLOCATABLE :: ppo4_bkginc     ! ppo4 inc
    REAL(wp), DIMENSION(:,:,:),   ALLOCATABLE :: pdic_bkginc     ! pdic inc
    REAL(wp), DIMENSION(:,:,:),   ALLOCATABLE :: palk_bkginc     ! palk inc
@@ -229 +231 @@
 
 #if defined key_hadocc
-      IF ( ln_slchldiainc .OR. ln_slchlnoninc .OR. ln_slphydiainc .OR. &
-         & ln_slphynoninc .OR. ln_psi4inc .OR. ln_pphinc .OR. ln_po2inc ) THEN
-         CALL ctl_stop( ' Cannot assimilate PFTs, Si4, pH or O2 into HadOCC' )
+      IF ( ln_slchldiainc .OR. ln_slchlnoninc .OR. ln_slchlnaninc .OR. &
+         & ln_slchldininc .OR. ln_slchlpicinc .OR. ln_slphydiainc .OR. &
+         & ln_slphynoninc .OR. ln_psi4inc     .OR. ln_pphinc      .OR. &
+         & ln_po2inc      .OR. ln_ppo4inc ) THEN
+         CALL ctl_stop( ' Cannot assimilate PFTs, Si4, PO4, pH or O2 into HadOCC' )
       ENDIF
 #endif
@@ -238 +242 @@
       IF ( .NOT. ln_foam_medusa ) THEN
          CALL ctl_stop( ' MEDUSA data assimilation options not turned on: set ln_foam_medusa' )
+      ENDIF
+
+      IF ( ln_ppo4inc ) THEN
+         CALL ctl_stop( ' Cannot assimilate PO4 into MEDUSA' )
       ENDIF
 #endif
@@ -419 +427 @@
          ALLOCATE( psi4_bkginc(jpi,jpj,jpk) )
          CALL asm_bgc_read_incs_3d( knum, 'bckinpsi4', psi4_bkginc )
+      ENDIF
+      
+      IF ( ln_ppo4inc ) THEN
+         ALLOCATE( ppo4_bkginc(jpi,jpj,jpk) )
+         CALL asm_bgc_read_incs_3d( knum, 'bckinppo4', ppo4_bkginc )
       ENDIF
       
@@ -2441 +2454 @@
          ENDIF
 
+         IF ( ln_ppo4inc ) THEN
+#if defined key_fabm
+            it = jp_fabm_n1p
+
+            IF ( ln_phytobal ) THEN
+               ppo4_bkginc(:,:,:) = ppo4_bkginc(:,:,:) - phyto2d_balinc(:,:,:,it)
+            ENDIF
+            IF ( ln_plchltotinc .OR. ln_pchltotinc ) THEN
+               ppo4_bkginc(:,:,:) = ppo4_bkginc(:,:,:) - phyto3d_balinc(:,:,:,it)
+            ENDIF
+            IF ( ln_sfco2inc .OR. ln_spco2inc ) THEN
+               ppo4_bkginc(:,:,:) = ppo4_bkginc(:,:,:) - pco2_balinc(:,:,:,it)
+            ENDIF
+            IF ( ln_pphinc ) THEN
+               ppo4_bkginc(:,:,:) = ppo4_bkginc(:,:,:) - ph_balinc(:,:,:,it)
+            ENDIF
+#else
+            CALL ctl_stop ( ' bgc3d_asm_inc: no compatible BGC model defined' )
+#endif
+         ENDIF
+
          IF ( ln_pdicinc ) THEN
 #if defined key_hadocc || defined key_medusa || defined key_fabm
@@ -2585 +2619 @@
             ENDIF
 
+            IF ( ln_ppo4inc ) THEN
+#if defined key_fabm
+               WHERE( ppo4_bkginc(:,:,:) > 0.0_wp .OR. &
+                  &   trn(:,:,:,jp_fabm_n1p) + ppo4_bkginc(:,:,:) * zincwgt > 0.0_wp )
+                  trn(:,:,:,jp_fabm_n1p) = trn(:,:,:,jp_fabm_n1p) + ppo4_bkginc(:,:,:) * zincwgt
+                  trb(:,:,:,jp_fabm_n1p) = trb(:,:,:,jp_fabm_n1p) + ppo4_bkginc(:,:,:) * zincwgt
+               END WHERE
+#else
+               CALL ctl_stop ( ' bgc3d_asm_inc: no compatible BGC model defined' )
+#endif
+            ENDIF
+
             IF ( ln_pdicinc ) THEN
 #if defined key_hadocc
@@ -2654 +2700 @@
                IF ( ln_pno3inc ) DEALLOCATE( pno3_bkginc )
                IF ( ln_psi4inc ) DEALLOCATE( psi4_bkginc )
+               IF ( ln_ppo4inc ) DEALLOCATE( ppo4_bkginc )
                IF ( ln_pdicinc ) DEALLOCATE( pdic_bkginc )
                IF ( ln_palkinc ) DEALLOCATE( palk_bkginc )
@@ -2718 +2765 @@
             ENDIF
 
+            IF ( ln_ppo4inc ) THEN
+#if defined key_fabm
+               WHERE( ppo4_bkginc(:,:,:) > 0.0_wp .OR. &
+                  &   trn(:,:,:,jp_fabm_n1p) + ppo4_bkginc(:,:,:) > 0.0_wp )
+                  trn(:,:,:,jp_fabm_n1p) = trn(:,:,:,jp_fabm_n1p) + ppo4_bkginc(:,:,:)
+                  trb(:,:,:,jp_fabm_n1p) = trn(:,:,:,jp_fabm_n1p)
+               END WHERE
+#else
+               CALL ctl_stop ( ' bgc3d_asm_inc: no compatible BGC model defined' )
+#endif
+            ENDIF
+
             IF ( ln_pdicinc ) THEN
 #if defined key_hadocc
@@ -2786 +2845 @@
             IF ( ln_pno3inc ) DEALLOCATE( pno3_bkginc )
             IF ( ln_psi4inc ) DEALLOCATE( psi4_bkginc )
+            IF ( ln_ppo4inc ) DEALLOCATE( ppo4_bkginc )
             IF ( ln_pdicinc ) DEALLOCATE( pdic_bkginc )
             IF ( ln_palkinc ) DEALLOCATE( palk_bkginc )