pkg/darwin/darwin_init_vari.F

C $Header: /u/gcmpack/MITgcm_contrib/darwin2/pkg/darwin/darwin_init_vari.F,v 1.5 2013/12/04 21:19:12 jahn Exp $
C $Name:  $

#include "DARWIN_OPTIONS.h"

C !INTERFACE: ==========================================================
      SUBROUTINE DARWIN_INIT_VARI(myThid )

C !DESCRIPTION:
C calls subroutines that initialize any time dependent variables
C for any tracer experiment

C !USES: ===============================================================
      IMPLICIT NONE

#include "SIZE.h"
#include "EEPARAMS.h"
#include "DARWIN_SIZE.h"
#include "DARWIN_IO.h"

C !INPUT PARAMETERS: ===================================================
C  myThid               :: thread number
      INTEGER myThid
CEOP

#ifdef ALLOW_DARWIN

C !FUNCTIONS: ==========================================================
      LOGICAL MASTER_CPU_THREAD
      EXTERNAL MASTER_CPU_THREAD

C !LOCAL VARIABLES: ====================================================
      INTEGER i,j,bi,bj

c set up ecosystem coefficients
c
c seed randomization
      _BEGIN_MASTER( myThid )
      CALL DARWIN_RANDOM_INIT(darwin_seed, myThid)
      _END_MASTER( myThid )
      _BARRIER

#ifdef ALLOW_QUOTA
      call QUOTA_INIT_VARI(myThid )
#endif

#ifdef ALLOW_MONOD
      call MONOD_INIT_VARI(myThid )
#endif

#ifdef CALC_RATE_TOTALS
c find unused units for darwin_rate output
      IF (MASTER_CPU_THREAD(myThid)) THEN
       CALL MDSFINDUNIT( DAR_ratetot_unit, mythid )
       open(DAR_ratetot_unit,file='darwin_ratetot.txt',status='unknown')
#ifdef ALLOW_CARBON
       CALL MDSFINDUNIT( DAR_fluxtot_unit, mythid )
       open(DAR_fluxtot_unit,file='darwin_fluxtot.txt',status='unknown')
#endif
      ENDIF
#endif

#endif /* ALLOW_DARWIN */

      RETURN
      END
1	jahn	1.6	C $Header: /u/gcmpack/MITgcm_contrib/darwin2/pkg/darwin/darwin_init_vari.F,v 1.5 2013/12/04 21:19:12 jahn Exp $
2	stephd	1.2	C $Name: $
3	jahn	1.1
4			#include "DARWIN_OPTIONS.h"
5
6			C !INTERFACE: ==========================================================
7			SUBROUTINE DARWIN_INIT_VARI(myThid )
8
9			C !DESCRIPTION:
10			C calls subroutines that initialize any time dependent variables
11			C for any tracer experiment
12
13			C !USES: ===============================================================
14			IMPLICIT NONE
15
16			#include "SIZE.h"
17			#include "EEPARAMS.h"
18			#include "DARWIN_SIZE.h"
19			#include "DARWIN_IO.h"
20
21			C !INPUT PARAMETERS: ===================================================
22			C myThid :: thread number
23			INTEGER myThid
24			CEOP
25
26			#ifdef ALLOW_DARWIN
27
28	jahn	1.5	C !FUNCTIONS: ==========================================================
29			LOGICAL MASTER_CPU_THREAD
30			EXTERNAL MASTER_CPU_THREAD
31
32			C !LOCAL VARIABLES: ====================================================
33			INTEGER i,j,bi,bj
34	jahn	1.1
35			c set up ecosystem coefficients
36			c
37			c seed randomization
38	jahn	1.4	_BEGIN_MASTER( myThid )
39	jahn	1.1	CALL DARWIN_RANDOM_INIT(darwin_seed, myThid)
40	jahn	1.4	_END_MASTER( myThid )
41	jahn	1.6	_BARRIER
42	jahn	1.1
43			#ifdef ALLOW_QUOTA
44			call QUOTA_INIT_VARI(myThid )
45			#endif
46
47			#ifdef ALLOW_MONOD
48			call MONOD_INIT_VARI(myThid )
49			#endif
50
51	stephd	1.2	#ifdef CALC_RATE_TOTALS
52			c find unused units for darwin_rate output
53	jahn	1.5	IF (MASTER_CPU_THREAD(myThid)) THEN
54			CALL MDSFINDUNIT( DAR_ratetot_unit, mythid )
55			open(DAR_ratetot_unit,file='darwin_ratetot.txt',status='unknown')
56	stephd	1.2	#ifdef ALLOW_CARBON
57	jahn	1.5	CALL MDSFINDUNIT( DAR_fluxtot_unit, mythid )
58			open(DAR_fluxtot_unit,file='darwin_fluxtot.txt',status='unknown')
59	stephd	1.2	#endif
60	jahn	1.5	ENDIF
61	stephd	1.2	#endif
62
63	jahn	1.1	#endif /* ALLOW_DARWIN */
64
65			RETURN
66			END