New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 7803 for branches/UKMO/dev_r5518_v3.6_asm_nemovar_community/NEMOGCM/NEMO/OPA_SRC/ASM/asminc.F90 – NEMO

Ignore:
Timestamp:
2017-03-16T12:18:47+01:00 (7 years ago)
Author:
dford
Message:

Merge in changes to allow a maximum chlorophyll assimilation increment to be specified. See internal UKMO ticket 690.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/UKMO/dev_r5518_v3.6_asm_nemovar_community/NEMOGCM/NEMO/OPA_SRC/ASM/asminc.F90

    r7731 r7803  
    131131   REAL(wp), PUBLIC, DIMENSION(:,:,:,:), ALLOCATABLE :: logchl_balinc  !: Increment to BGC variables from logchl assim 
    132132#endif 
     133   REAL(wp) :: rn_maxchlinc = -999.0  !: maximum absolute non-log chlorophyll increment from logchl assimilation 
     134                                      !: <= 0 implies no maximum applied (switch turned off) 
     135                                      !:  > 0 implies maximum absolute chl increment capped at this value 
    133136 
    134137   INTEGER :: mld_choice        = 4   !: choice of mld criteria to use for physics assimilation 
     
    207210         &                 nitbkg, nitdin, nitiaustr, nitiaufin, niaufn,   & 
    208211         &                 ln_salfix, salfixmin, nn_divdmp, nitavgbkg, mld_choice, & 
    209          &                 mld_choice_bgc 
     212         &                 mld_choice_bgc, rn_maxchlinc 
    210213      !!---------------------------------------------------------------------- 
    211214 
     
    288291         WRITE(numout,*) '      Choice of MLD for physics assimilation                  mld_choice = ', mld_choice 
    289292         WRITE(numout,*) '      Choice of MLD for BGC assimilation                  mld_choice_bgc = ', mld_choice_bgc 
     293         WRITE(numout,*) '      Maximum absolute chlorophyll increment (<=0 = off)    rn_maxchlinc = ', rn_maxchlinc 
    290294      ENDIF 
    291295 
     
    14081412#if defined key_fabm 
    14091413         CALL asm_logchl_bal_ersem( ln_logchlpftinc, nn_asmpfts, mld_choice_bgc, & 
    1410             &                       logchl_bkginc, logchl_balinc ) 
     1414            &                       rn_maxchlinc, logchl_bkginc, logchl_balinc ) 
    14111415#elif defined key_medusa && defined key_foam_medusa 
    14121416         !CALL asm_logchl_bal_medusa() 
Note: See TracChangeset for help on using the changeset viewer.