1 |
C $Header: /u/gcmpack/MITgcm_contrib/darwin2/pkg/darwin/darwin_init_fixed.F,v 1.2 2013/12/04 21:21:32 jahn Exp $ |
2 |
C $Name: $ |
3 |
|
4 |
#include "DARWIN_OPTIONS.h" |
5 |
|
6 |
C !INTERFACE: ========================================================== |
7 |
SUBROUTINE DARWIN_INIT_FIXED(myThid ) |
8 |
|
9 |
C !DESCRIPTION: |
10 |
C calls subroutines that initialized fixed variables for any tracer |
11 |
C experiment |
12 |
|
13 |
C !USES: =============================================================== |
14 |
IMPLICIT NONE |
15 |
#include "SIZE.h" |
16 |
#include "EEPARAMS.h" |
17 |
#include "PARAMS.h" |
18 |
#ifdef ALLOW_SPECTRAL |
19 |
#include "SPECTRAL_SIZE.h" |
20 |
#include "SPECTRAL_PARAMS.h" |
21 |
#endif |
22 |
|
23 |
C !INPUT PARAMETERS: =================================================== |
24 |
C myThid :: thread number |
25 |
INTEGER myThid |
26 |
CEOP |
27 |
|
28 |
#ifdef ALLOW_DARWIN |
29 |
|
30 |
C !LOCAL VARIABLES: |
31 |
#ifdef ALLOW_SPECTRAL |
32 |
INTEGER ilam |
33 |
_RL planck, c, hc, oavo, hcoavo, rlamm |
34 |
#endif |
35 |
|
36 |
_BEGIN_MASTER(myThid) |
37 |
|
38 |
#ifdef ALLOW_SPECTRAL |
39 |
c Quanta conversion |
40 |
planck = 6.6256 _d -34 !Plancks constant J sec |
41 |
c = 2.998 _d 8 !speed of light m/sec |
42 |
hc = 1.0/(planck*c) |
43 |
oavo = 1.0/6.023 _d 23 ! 1/Avogadros number |
44 |
hcoavo = hc*oavo |
45 |
do ilam = 1,tlam |
46 |
rlamm = darwin_waves(ilam)*1 _d -9 !lambda in m |
47 |
WtouEins(ilam) = 1 _d 6*rlamm*hcoavo !Watts to uEin/s conversion |
48 |
enddo |
49 |
#endif |
50 |
|
51 |
_END_MASTER(myThid) |
52 |
|
53 |
#ifdef ALLOW_QUOTA |
54 |
call QUOTA_INIT_FIXED(myThid) |
55 |
#endif |
56 |
#ifdef ALLOW_MONOD |
57 |
call MONOD_INIT_FIXED(myThid) |
58 |
#endif |
59 |
|
60 |
_BARRIER |
61 |
|
62 |
#endif /* ALLOW_DARWIN */ |
63 |
|
64 |
RETURN |
65 |
END |