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_DARWIN |
8 |
|
|
#ifdef ALLOW_QUOTA |
9 |
|
|
c ==================================================================== |
10 |
|
|
c SUBROUTINE QUOTA_TEMPFUNC |
11 |
|
|
c ==================================================================== |
12 |
|
|
|
13 |
|
|
SUBROUTINE QUOTA_TEMPFUNC( |
14 |
|
|
I Tlocal, |
15 |
|
|
O photo_TempFunction, |
16 |
|
|
O activ_TempFunction, |
17 |
|
|
I myThid) |
18 |
|
|
|
19 |
|
|
implicit none |
20 |
|
|
#include "QUOTA_SIZE.h" |
21 |
|
|
#include "QUOTA.h" |
22 |
|
|
|
23 |
|
|
C !INPUT PARAMETERS: =================================================== |
24 |
|
|
C myThid :: thread number |
25 |
|
|
INTEGER myThid |
26 |
|
|
_RL photo_TempFunction |
27 |
|
|
_RL activ_TempFunction |
28 |
|
|
_RL Tlocal |
29 |
|
|
c local |
30 |
|
|
_RL Temp_P,Temp_A,Tempref |
31 |
|
|
INTEGER jp |
32 |
|
|
c |
33 |
|
|
|
34 |
|
|
! Temp_P = 0.0637 _d 0 |
35 |
|
|
! Temp_A = 0.1 _d 0 |
36 |
|
|
Temp_P = 0.05 _d 0 |
37 |
|
|
Temp_A = 0.05 _d 0 |
38 |
|
|
|
39 |
benw |
1.2 |
Tempref= 20. _d 0 |
40 |
jahn |
1.1 |
c |
41 |
|
|
photo_TempFunction = exp(Temp_P*(Tlocal - Tempref)) |
42 |
|
|
photo_TempFunction = max(photo_TempFunction,1. _d -10) |
43 |
|
|
c |
44 |
|
|
activ_TempFunction = exp(Temp_A*(Tlocal - Tempref)) |
45 |
|
|
activ_TempFunction = max(activ_TempFunction,1. _d -10) |
46 |
|
|
c |
47 |
|
|
RETURN |
48 |
|
|
END |
49 |
|
|
#endif /*ALLOW_QUOTA*/ |
50 |
|
|
#endif /*DARWIN*/ |
51 |
|
|
#endif /*ALLOW_PTRACERS*/ |
52 |
|
|
c ================================================================== |