/[MITgcm]/MITgcm_contrib/darwin2/pkg/monod/monod_tempfunc.F
ViewVC logotype

Annotation of /MITgcm_contrib/darwin2/pkg/monod/monod_tempfunc.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph

Revision 1.3 - (hide annotations) (download)
Thu Mar 28 20:14:42 2013 UTC (12 years, 3 months ago) by stephd
Branch: MAIN
CVS Tags: ctrb_darwin2_ckpt64k_20130723, ctrb_darwin2_ckpt65w_20160512, ctrb_darwin2_ckpt65j_20150225, ctrb_darwin2_ckpt66g_20170424, ctrb_darwin2_ckpt64h_20130528, ctrb_darwin2_ckpt66k_20171025, ctrb_darwin2_ckpt66n_20180118, ctrb_darwin2_ckpt65v_20160409, ctrb_darwin2_ckpt65s_20160114, ctrb_darwin2_ckpt65_20140718, ctrb_darwin2_ckpt64m_20130820, ctrb_darwin2_ckpt66d_20170214, ctrb_darwin2_ckpt64r_20131210, ctrb_darwin2_ckpt65m_20150615, ctrb_darwin2_ckpt65q_20151118, ctrb_darwin2_ckpt65o_20150914, ctrb_darwin2_ckpt64f_20130405, ctrb_darwin2_ckpt65p_20151023, ctrb_darwin2_ckpt64n_20130826, ctrb_darwin2_ckpt65e_20140929, ctrb_darwin2_ckpt64o_20131024, ctrb_darwin2_ckpt64v_20140411, ctrb_darwin2_ckpt64z_20140711, ctrb_darwin2_ckpt65l_20150504, ctrb_darwin2_ckpt65z_20160929, ctrb_darwin2_ckpt65n_20150729, ctrb_darwin2_ckpt64y_20140622, ctrb_darwin2_ckpt65d_20140915, ctrb_darwin2_ckpt64t_20140202, ctrb_darwin2_ckpt66h_20170602, ctrb_darwin2_ckpt64i_20130622, ctrb_darwin2_ckpt64s_20140105, ctrb_darwin2_ckpt64x_20140524, ctrb_darwin2_ckpt65x_20160612, ctrb_darwin2_ckpt66f_20170407, ctrb_darwin2_ckpt65g_20141120, ctrb_darwin2_ckpt65k_20150402, ctrb_darwin2_ckpt64w_20140502, ctrb_darwin2_ckpt66a_20161020, ctrb_darwin2_ckpt64g_20130503, ctrb_darwin2_ckpt64l_20130806, ctrb_darwin2_ckpt65f_20141014, ctrb_darwin2_ckpt66b_20161219, ctrb_darwin2_ckpt64u_20140308, ctrb_darwin2_ckpt64j_20130704, ctrb_darwin2_ckpt65i_20150123, ctrb_darwin2_ckpt66j_20170815, ctrb_darwin2_ckpt65y_20160801, ctrb_darwin2_ckpt66c_20170121, ctrb_darwin2_ckpt65a_20140728, ctrb_darwin2_ckpt65b_20140812, ctrb_darwin2_ckpt65t_20160221, ctrb_darwin2_ckpt64p_20131118, ctrb_darwin2_ckpt66o_20180209, ctrb_darwin2_ckpt66e_20170314, ctrb_darwin2_ckpt64q_20131118, ctrb_darwin2_ckpt64p_20131024, ctrb_darwin2_ckpt65u_20160315, ctrb_darwin2_ckpt65r_20151221, ctrb_darwin2_ckpt66i_20170718, ctrb_darwin2_ckpt65c_20140830, ctrb_darwin2_ckpt66l_20171025, ctrb_darwin2_ckpt65h_20141217, ctrb_darwin2_ckpt66m_20171213, HEAD
Changes since 1.2: +7 -7 lines 
o comment out so grazing tempfunction is 1 always - to keep backwards compatible
1 jahn 1.1
2     #include "CPP_OPTIONS.h"
3     #include "PTRACERS_OPTIONS.h"
4     #include "DARWIN_OPTIONS.h"
5    
6     #ifdef ALLOW_PTRACERS
7     #ifdef ALLOW_MONOD
8    
9     c ====================================================================
10     c SUBROUTINE MONOD_TEMPFUNC
11     c ====================================================================
12    
13     SUBROUTINE MONOD_TEMPFUNC(
14     I Tlocal,
15     O phytoTempFunction,
16     #ifndef QUOTA
17     O zooTempFunction,
18     #endif
19     O reminTempFunction,
20 stephd 1.2 O mortPTempFunction,
21     O mortZTempFunction,
22     O mortZ2TempFunction,
23 jahn 1.1 I myThid)
24    
25     implicit none
26     #ifdef QUOTA
27     #include "QUOTA_SIZE.h"
28     #include "QUOTA.h"
29     #else
30     #include "MONOD_SIZE.h"
31     #include "MONOD.h"
32     #endif
33     #include "DARWIN_PARAMS.h"
34    
35     C !INPUT PARAMETERS: ===================================================
36     C myThid :: thread number
37     INTEGER myThid
38     _RL phytoTempFunction(npmax)
39     #ifndef QUOTA
40     _RL zooTempFunction(nzmax)
41     #endif
42     _RL reminTempFunction
43 stephd 1.2 _RL mortPTempFunction
44     _RL mortZTempFunction
45     _RL mortZ2TempFunction
46 jahn 1.1 _RL Tlocal
47    
48     c local
49     _RL Tkel
50     _RL TempAe, Tempref, TempCoeff
51     INTEGER np, nz
52    
53     c TEMP_VERSION 1 used in Follows et al (2007) - without max of 1
54     c TEMP_VERSION 1 used in Dutkiewicz et al (2009), Hickman et al,
55     c Monteiro etal, Barton et al
56     c TEMP_RANGE gives all phyto a specific temperature range
57     c if undefined, then full Eppley/Arrenhius curve used
58     c
59     #if TEMP_VERSION == 1
60     c ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
61     c +++++++++++++++++++ VERSION 1 +++++++++++++++++++++++++++++++++++++++
62     c steph's version 1 (pseudo-Eppley)
63     c plankton growth function (unitless)
64     do np=1, npmax
65     cswd -- this gives Eppley curve only
66     phytoTempFunction(np) =
67     & (phytoTempExp1(np)**Tlocal)
68     #ifdef TEMP_RANGE
69     cswd -- temperature range
70     phytoTempFunction(np) = phytoTempFunction(np) *
71     & (exp(- phytoTempExp2(np)*
72     & (abs(Tlocal - phytoTempOptimum(np)))**
73     & phytoDecayPower(np)))
74     #endif
75     phytoTempFunction(np) = phytoTempFunction(np) - tempnorm
76     phytoTempFunction(np) = phytoTempCoeff(np)*
77     & max(phytoTempFunction(np), 1. _d -10)
78     phytoTempFunction(np) = min(phytoTempFunction(np),1. _d 0)
79     enddo
80     #ifndef QUOTA
81     do nz = 1,nzmax
82     c zooTempFunction(nz) = zooTempCoeff(nz)*EXP(
83     c & zooTempExp(nz)*(Tlocal - zooTempOptimum(nz)))
84     zooTempFunction(nz) = 1.0 _d 0
85     end do
86     #endif
87     reminTempFunction = 1.0 _d 0
88 stephd 1.2 mortPTempFunction = 1.0 _d 0
89     mortZTempFunction = 1.0 _d 0
90     mortZ2TempFunction = 1.0 _d 0
91 jahn 1.1 c ++++++++++++++ END VERSION 1 +++++++++++++++++++++++++++++++++++++++
92     c++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
93     #elif TEMP_VERSION == 2
94     c ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
95     c +++++++++++++++++++ VERSION 2 +++++++++++++++++++++++++++++++++++++++
96     c steph's version 2 (pseudo-Arrenhius)
97     Tkel=273.15 _d 0
98     TempAe=-4000. _d 0
99     Tempref=293.15 _d 0
100     tempnorm= 0. _d 0
101     TempCoeff=0.5882 _d 0
102     do np=1, npmax
103     phytoTempCoeff(np) = TempCoeff
104     cswd -- this gives Arrenhius curve only
105     phytoTempFunction(np) =
106     & exp(TempAe*(1. _d 0/(Tlocal+Tkel) -
107     & 1. _d 0/(Tempref) ) )
108     #ifdef TEMP_RANGE
109     cswd -- temperature range
110     phytoTempFunction(np) = phytoTempFunction(np) *
111     & (exp(- phytoTempExp2(np)*
112     & (abs(Tlocal - phytoTempOptimum(np)))**
113     & phytoDecayPower(np)))
114     #endif
115     phytoTempFunction(np) = phytoTempFunction(np) - tempnorm
116     phytoTempFunction(np) = phytoTempCoeff(np)*
117     & max(phytoTempFunction(np), 1. _d -10)
118     enddo
119     #ifndef QUOTA
120     do nz = 1,nzmax
121 stephd 1.3 c zooTempFunction(nz) =
122     c & exp(TempAe*(1/(Tlocal+Tkel) -
123     c & 1/(Tempref) ) )
124     c zooTempFunction(nz) = zooTempFunction(nz) - tempnorm
125     c zooTempFunction(nz) = TempCoeff*
126     c & max(zooTempFunction(nz), 1. _d -10)
127     zooTempFunction(nz) = 1. _d 0
128 jahn 1.1 end do
129     #endif
130     reminTempFunction = exp(TempAe*(1/(Tlocal+Tkel) -
131     & 1/(Tempref) ) )
132     reminTempFunction = reminTempFunction - tempnorm
133     reminTempFunction = TempCoeff*
134     & max(reminTempFunction, 1. _d -10)
135     c reminTempFunction = 1. _d 0
136 stephd 1.2 c mortPTempFunction = reminTempFunction
137     c mortZTempFunction = reminTempFunction
138     c mortZ2TempFunction = reminTempFunction
139     mortPTempFunction = 1.0 _d 0
140     mortZTempFunction = 1.0 _d 0
141     mortZ2TempFunction = 1.0 _d 0
142 jahn 1.1 c ++++++++++++++ END VERSION 2 +++++++++++++++++++++++++++++++++++++++
143     c++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
144     #else
145     #error "TEMP_VERSION must be 1 or 2. Define in DARWIN_OPTIONS.h"
146     #endif
147     c
148     c ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
149     c +++++++++++++++++++ NO TEMP LIMITATION +++++++++++++++++++++++++++++
150     #ifdef NOTEMP
151     do np=1, npmax
152     phytoTempFunction(np) = 1.0 _d 0
153     enddo
154     #ifndef QUOTA
155     do nz = 1,nzmax
156     zooTempFunction(nz) = 1.0 _d 0
157     end do
158     #endif
159     reminTempFunction = 1.0 _d 0
160 stephd 1.2 mortTempFunction = 1.0 _d 0
161     mort2TempFunction = 1.0 _d 0
162 jahn 1.1 #endif
163     c +++++++++++++ END NO TEMP LIMITATION +++++++++++++++++++++++++++++
164     c ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
165     c
166    
167     RETURN
168     END
169     #endif /*MONOD*/
170     #endif /*ALLOW_PTRACERS*/
171     c ==================================================================
  ViewVC Help
Powered by ViewVC 1.1.22