Changeset 7791 for branches/UKMO/dev_r5518_v3.6_asm_nemovar_community_maxchlinc/NEMOGCM/NEMO/OPA_SRC/ASM/asmlogchlbal_ersem.F90

Timestamp:

2017-03-14T16:13:59+01:00 (7 years ago)

Author:

dford

Message:

Code to implement maximum chlorophyll increment option.

File:

• branches/UKMO/dev_r5518_v3.6_asm_nemovar_community_maxchlinc/NEMOGCM/NEMO/OPA_SRC/ASM/asmlogchlbal_ersem.F90 (modified) (5 diffs)

branches/UKMO/dev_r5518_v3.6_asm_nemovar_community_maxchlinc/NEMOGCM/NEMO/OPA_SRC/ASM/asmlogchlbal_ersem.F90

@@ -30 +30 @@
 
    SUBROUTINE asm_logchl_bal_ersem( ln_logchlpftinc, npfts, mld_choice_bgc, &
-      &                             logchl_bkginc, logchl_balinc )
+      &                             rn_maxchlinc, logchl_bkginc, logchl_balinc )
       !!---------------------------------------------------------------------------
       !!                    ***  ROUTINE asm_logchl_bal_ersem  ***
@@ -49 +49 @@
       INTEGER,  INTENT(in   )                               :: npfts
       INTEGER,  INTENT(in   )                               :: mld_choice_bgc
+      REAL(wp), INTENT(in   )                               :: rn_maxchlinc
       REAL(wp), INTENT(in   ), DIMENSION(jpi,jpj,npfts)     :: logchl_bkginc
       REAL(wp), INTENT(  out), DIMENSION(jpi,jpj,jpk,jptra) :: logchl_balinc
@@ -65 +66 @@
       ! 2) Take 10^log10(analysis) to get analysis
       ! 3) Subtract background from analysis to get chl incs
+      ! If rn_maxchlinc > 0 then cap total absolute chlorophyll increment at that value
       !
       ! Only apply increments if all of Chl1-4 background values are > 0
@@ -103 +105 @@
                         &                                             trn(ji,jj,1,jp_fabm_m1+jp_fabm_chl4)
                   ENDIF
+                  IF (rn_maxchlinc > 0.0) THEN
+                     chl_inc = logchl_balinc(ji,jj,1,jp_fabm_m1+jp_fabm_chl1) + &
+                               logchl_balinc(ji,jj,1,jp_fabm_m1+jp_fabm_chl2) + &
+                               logchl_balinc(ji,jj,1,jp_fabm_m1+jp_fabm_chl3) + &
+                               logchl_balinc(ji,jj,1,jp_fabm_m1+jp_fabm_chl4)
+                     IF ( ABS(chl_inc) > rn_maxchlinc ) THEN
+                        chl_tot = ABS(chl_inc) / rn_maxchlinc
+                        logchl_balinc(ji,jj,1,jp_fabm_m1+jp_fabm_chl1) = logchl_balinc(ji,jj,1,jp_fabm_m1+jp_fabm_chl1) / chl_tot
+                        logchl_balinc(ji,jj,1,jp_fabm_m1+jp_fabm_chl2) = logchl_balinc(ji,jj,1,jp_fabm_m1+jp_fabm_chl2) / chl_tot
+                        logchl_balinc(ji,jj,1,jp_fabm_m1+jp_fabm_chl3) = logchl_balinc(ji,jj,1,jp_fabm_m1+jp_fabm_chl3) / chl_tot
+                        logchl_balinc(ji,jj,1,jp_fabm_m1+jp_fabm_chl4) = logchl_balinc(ji,jj,1,jp_fabm_m1+jp_fabm_chl4) / chl_tot
+                     ENDIF
+                  ENDIF
                ENDIF
             END DO
@@ -124 +139 @@
                      &      trn(ji,jj,1,jp_fabm_m1+jp_fabm_chl3) + trn(ji,jj,1,jp_fabm_m1+jp_fabm_chl4)
                   chl_inc = 10**( LOG10( chl_tot ) + logchl_bkginc(ji,jj,1) ) - chl_tot
+                  IF (rn_maxchlinc > 0.0) THEN
+                     chl_inc = MAX( -1.0 * rn_maxchlinc, MIN( chl_inc, rn_maxchlinc ) )
+                  ENDIF
                   logchl_balinc(ji,jj,1,jp_fabm_m1+jp_fabm_chl1) = chl_inc * trn(ji,jj,1,jp_fabm_m1+jp_fabm_chl1) / chl_tot
                   logchl_balinc(ji,jj,1,jp_fabm_m1+jp_fabm_chl2) = chl_inc * trn(ji,jj,1,jp_fabm_m1+jp_fabm_chl2) / chl_tot