Changeset 7791 for branches/UKMO/dev_r5518_v3.6_asm_nemovar_community_maxchlinc/NEMOGCM/NEMO/OPA_SRC/ASM/asminc.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/asminc.F90 (modified) (4 diffs)

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

@@ -131 +131 @@
    REAL(wp), PUBLIC, DIMENSION(:,:,:,:), ALLOCATABLE :: logchl_balinc  !: Increment to BGC variables from logchl assim
 #endif
+   REAL(wp) :: rn_maxchlinc = -999.0  !: maximum absolute non-log chlorophyll increment from logchl assimilation
+                                      !: <= 0 implies no maximum applied (switch turned off)
+                                      !:  > 0 implies maximum absolute chl increment capped at this value
 
    INTEGER :: mld_choice        = 4   !: choice of mld criteria to use for physics assimilation
@@ -207 +210 @@
          &                 nitbkg, nitdin, nitiaustr, nitiaufin, niaufn,   &
          &                 ln_salfix, salfixmin, nn_divdmp, nitavgbkg, mld_choice, &
-         &                 mld_choice_bgc
+         &                 mld_choice_bgc, rn_maxchlinc
       !!----------------------------------------------------------------------
 
@@ -288 +291 @@
          WRITE(numout,*) '      Choice of MLD for physics assimilation                  mld_choice = ', mld_choice
          WRITE(numout,*) '      Choice of MLD for BGC assimilation                  mld_choice_bgc = ', mld_choice_bgc
+         WRITE(numout,*) '      Maximum absolute chlorophyll increment (<=0 = off)    rn_maxchlinc = ', rn_maxchlinc
       ENDIF
 
@@ -1408 +1412 @@
 #if defined key_fabm
          CALL asm_logchl_bal_ersem( ln_logchlpftinc, nn_asmpfts, mld_choice_bgc, &
-            &                       logchl_bkginc, logchl_balinc )
+            &                       rn_maxchlinc, logchl_bkginc, logchl_balinc )
 #elif defined key_medusa && defined key_foam_medusa
          !CALL asm_logchl_bal_medusa()