Changeset 8124


Ignore:
Timestamp:
2017-06-02T08:06:17+02:00 (3 years ago)
Author:
alistairsellar
Message:

Introduce a ramp to the minimum term so that it scales with fq1. It cannot drop
below 0.57.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/UKMO/test_r7519_MEDUSA_anderson_DMSminramp/NEMOGCM/NEMO/TOP_SRC/MEDUSA/trcdms_medusa.F90

    r6719 r8124  
    110110        !! of diatoms and non-diatoms) 
    111111        Qterm = xdin / (xdin + 0.5) 
    112         fq1 = log10(CHL * Jterm * Qterm) 
     112        fq1 = max(0.0, log10(CHL * Jterm * Qterm)) 
    113113        if (fq1 > 1.72) then 
    114114           dms_andr = (8.24 * (fq1 - 1.72)) + 2.29 
    115115        else 
    116            dms_andr = 2.29 
     116           dms_andr = fq1 + 0.57 
    117117        endif 
    118118! 
Note: See TracChangeset for help on using the changeset viewer.