/[MITgcm]/MITgcm_contrib/darwin2/pkg/darwin/darwin_init_vari.F
ViewVC logotype

Annotation of /MITgcm_contrib/darwin2/pkg/darwin/darwin_init_vari.F

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


Revision 1.7 - (hide annotations) (download)
Wed Dec 4 21:27:54 2013 UTC (11 years, 7 months ago) by jahn
Branch: MAIN
CVS Tags: ctrb_darwin2_ckpt65w_20160512, ctrb_darwin2_ckpt65j_20150225, ctrb_darwin2_ckpt66g_20170424, ctrb_darwin2_ckpt66k_20171025, ctrb_darwin2_ckpt66n_20180118, ctrb_darwin2_ckpt65v_20160409, ctrb_darwin2_ckpt65s_20160114, ctrb_darwin2_ckpt65_20140718, ctrb_darwin2_ckpt66d_20170214, ctrb_darwin2_ckpt64r_20131210, ctrb_darwin2_ckpt65m_20150615, ctrb_darwin2_ckpt65q_20151118, ctrb_darwin2_ckpt65o_20150914, ctrb_darwin2_ckpt65p_20151023, ctrb_darwin2_ckpt65e_20140929, 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_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_ckpt65f_20141014, ctrb_darwin2_ckpt66b_20161219, ctrb_darwin2_ckpt64u_20140308, 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_ckpt66o_20180209, ctrb_darwin2_ckpt66e_20170314, 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.6: +11 -1 lines
read seaice fraction through exf
- set areamaskfile in data.exf to a file with multiple time records
  and provide timing information
- need to define EXF_SEAICE_FRACTION in EXF_OPTIONS.h for this to work

1 jahn 1.7 C $Header: /u/gcmpack/MITgcm_contrib/darwin2/pkg/darwin/darwin_init_vari.F,v 1.6 2013/12/04 21:21:49 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 jahn 1.7 DO bj = myByLo(myThid), myByHi(myThid)
44     DO bi = myBxLo(myThid), myBxHi(myThid)
45     DO j=1-Oly,sNy+Oly
46     DO i=1-Olx,sNx+Olx
47     FIce(i,j,bi,bj) = 0. _d 0
48     ENDDO
49     ENDDO
50     ENDDO
51     ENDDO
52    
53 jahn 1.1 #ifdef ALLOW_QUOTA
54     call QUOTA_INIT_VARI(myThid )
55     #endif
56    
57     #ifdef ALLOW_MONOD
58     call MONOD_INIT_VARI(myThid )
59     #endif
60    
61 stephd 1.2 #ifdef CALC_RATE_TOTALS
62     c find unused units for darwin_rate output
63 jahn 1.5 IF (MASTER_CPU_THREAD(myThid)) THEN
64     CALL MDSFINDUNIT( DAR_ratetot_unit, mythid )
65     open(DAR_ratetot_unit,file='darwin_ratetot.txt',status='unknown')
66 stephd 1.2 #ifdef ALLOW_CARBON
67 jahn 1.5 CALL MDSFINDUNIT( DAR_fluxtot_unit, mythid )
68     open(DAR_fluxtot_unit,file='darwin_fluxtot.txt',status='unknown')
69 stephd 1.2 #endif
70 jahn 1.5 ENDIF
71 stephd 1.2 #endif
72    
73 jahn 1.1 #endif /* ALLOW_DARWIN */
74    
75     RETURN
76     END

  ViewVC Help
Powered by ViewVC 1.1.22