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

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

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

Revision 1.4 - (hide annotations) (download)
Wed Dec 4 21:27:54 2013 UTC (11 years, 11 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.3: +55 -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.4 C $Header: /u/gcmpack/MITgcm_contrib/darwin2/pkg/darwin/darwin_init_fixed.F,v 1.3 2013/12/04 21:21:49 jahn Exp $
2 jahn 1.2 C $Name: $
3 jahn 1.1
4     #include "DARWIN_OPTIONS.h"
5 jahn 1.4 #ifdef ALLOW_EXF
6     #include "EXF_OPTIONS.h"
7     #endif
8 jahn 1.1
9     C !INTERFACE: ==========================================================
10     SUBROUTINE DARWIN_INIT_FIXED(myThid )
11    
12     C !DESCRIPTION:
13     C calls subroutines that initialized fixed variables for any tracer
14     C experiment
15    
16     C !USES: ===============================================================
17     IMPLICIT NONE
18     #include "SIZE.h"
19     #include "EEPARAMS.h"
20     #include "PARAMS.h"
21 jahn 1.4 #ifdef ALLOW_EXF
22     #include "EXF_PARAM.h"
23     #endif
24     #include "DARWIN_SIZE.h"
25     #include "DARWIN_IO.h"
26 jahn 1.1 #ifdef ALLOW_SPECTRAL
27     #include "SPECTRAL_SIZE.h"
28     #include "SPECTRAL_PARAMS.h"
29     #endif
30    
31     C !INPUT PARAMETERS: ===================================================
32     C myThid :: thread number
33     INTEGER myThid
34     CEOP
35    
36     #ifdef ALLOW_DARWIN
37    
38     C !LOCAL VARIABLES:
39 jahn 1.4 CHARACTER*(MAX_LEN_MBUF) msgBuf
40 jahn 1.1 #ifdef ALLOW_SPECTRAL
41     INTEGER ilam
42     _RL planck, c, hc, oavo, hcoavo, rlamm
43     #endif
44    
45 jahn 1.2 _BEGIN_MASTER(myThid)
46    
47 jahn 1.1 #ifdef ALLOW_SPECTRAL
48     c Quanta conversion
49     planck = 6.6256 _d -34 !Plancks constant J sec
50     c = 2.998 _d 8 !speed of light m/sec
51     hc = 1.0/(planck*c)
52     oavo = 1.0/6.023 _d 23 ! 1/Avogadros number
53     hcoavo = hc*oavo
54     do ilam = 1,tlam
55     rlamm = darwin_waves(ilam)*1 _d -9 !lambda in m
56     WtouEins(ilam) = 1 _d 6*rlamm*hcoavo !Watts to uEin/s conversion
57     enddo
58     #endif
59    
60 jahn 1.4 C which ice area fraction to use (by preference):
61     C 1. read from darwin_iceFile if given
62     C 2. read from areamaskfile through exf if given
63     C 3. use iceMask from THSICE if compiled and used
64     C 4. use AREA from seaice if compiled and used
65     darwin_useiceFile = .FALSE.
66     darwin_useareamask = .FALSE.
67     darwin_useiceMask = .FALSE.
68     darwin_useAREA = .FALSE.
69     darwin_haveIce = .FALSE.
70     IF ( darwin_iceFile .NE. ' ' ) THEN
71     darwin_useiceFile = .TRUE.
72     darwin_haveIce = .TRUE.
73     ENDIF
74     #ifdef ALLOW_EXF
75     IF (.NOT.darwin_haveIce.AND.useEXF.AND.areamaskFile.NE.' ') THEN
76     darwin_useareamask = .TRUE.
77     darwin_haveIce = .TRUE.
78     #ifndef EXF_SEAICE_FRACTION
79     WRITE(msgBuf,'(2A)') 'S/R DARWIN_INIT_FIXED:',
80     & ' you specified an areamaskfile but EXF_SEAICE_FRACTION'
81     CALL PRINT_ERROR( msgBuf , 1)
82     WRITE(msgBuf,'(2A)') 'S/R DARWIN_INIT_FIXED:',
83     & ' is not defined. You need to define it in EXF_OPTIONS.h'
84     CALL PRINT_ERROR( msgBuf , 1)
85     WRITE(msgBuf,'(2A)') 'S/R DARWIN_INIT_FIXED:',
86     & ' in order to use areamaskfile.'
87     CALL PRINT_ERROR( msgBuf , 1)
88     STOP 'DARWIN_INIT_FIXED'
89     #endif
90     ENDIF
91     #endif
92     #ifdef ALLOW_THSICE
93     IF ( .NOT.darwin_haveIce .AND. useTHSICE ) THEN
94     darwin_useiceMask = .TRUE.
95     darwin_haveIce = .TRUE.
96     ENDIF
97     #endif
98     #ifdef ALLOW_SEAICE
99     IF ( .NOT.darwin_haveIce .AND. useSEAICE ) THEN
100     darwin_useAREA = .TRUE.
101     darwin_haveIce = .TRUE.
102     ENDIF
103     #endif
104    
105 jahn 1.2 _END_MASTER(myThid)
106    
107 jahn 1.1 #ifdef ALLOW_QUOTA
108     call QUOTA_INIT_FIXED(myThid)
109     #endif
110     #ifdef ALLOW_MONOD
111     call MONOD_INIT_FIXED(myThid)
112     #endif
113    
114 jahn 1.3 _BARRIER
115    
116 jahn 1.1 #endif /* ALLOW_DARWIN */
117    
118     RETURN
119     END
  ViewVC Help
Powered by ViewVC 1.1.22