/[MITgcm]/MITgcm_contrib/darwin2/pkg/darwin/DARWIN_IO.h
ViewVC logotype

Annotation of /MITgcm_contrib/darwin2/pkg/darwin/DARWIN_IO.h

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


Revision 1.4 - (hide annotations) (download)
Thu Sep 22 19:23:31 2011 UTC (13 years, 9 months ago) by jahn
Branch: MAIN
CVS Tags: ctrb_darwin2_ckpt63l_20120405, ctrb_darwin2_ckpt63f_20111201, ctrb_darwin2_ckpt63o_20120629, ctrb_darwin2_ckpt63c_20111011, ctrb_darwin2_ckpt63i_20120124, ctrb_darwin2_ckpt63m_20120506, ctrb_darwin2_ckpt63e_20111107, ctrb_darwin2_ckpt63j_20120217, ctrb_darwin2_ckpt63g_20111220, ctrb_darwin2_ckpt63h_20111230, ctrb_darwin2_ckpt63d_20111107, ctrb_darwin2_ckpt63n_20120604, ctrb_darwin2_ckpt63k_20120317
Changes since 1.3: +9 -2 lines
File MIME type: text/plain
add namelist parameters darwin_ironFileUnits darwin_dustFrac

1 jahn 1.4 C $Header: /u/gcmpack/MITgcm_contrib/darwin2/pkg/darwin/DARWIN_IO.h,v 1.3 2011/09/22 19:18:37 jahn Exp $
2 stephd 1.2 C $Name: $
3 jahn 1.1
4     #include "DARWIN_OPTIONS.h"
5    
6     c DARWIN_IO.h
7     c --------------------------------------------------------------------
8     c Description: input and diagnostic output for DARWIN pkg
9     c Stephanie Dutkiewicz: Spring 2006
10     c NOTE: additional diagnostics still needed
11     c --------------------------------------------------------------------
12    
13     C-- COMMON /DARWIN_FILENAMES/
14     C darwin_iceFile :: file name of seaice fraction
15     C darwin_ironFile :: file name of aeolian iron flux
16     C darwin_PARFile :: file name of Photosynthetically Active Radiation at surface
17     C darwin_nutWVelFile :: file name of nutrient wvel
18     C darwin_waterabsorbFile :: file name of water absorption/scattering spectra
19     C darwin_phytoabsorbFile :: file name of phyto absorption/scattering spectra
20     C darwin_particleabsorbFile :: file name of particle absorption/scattering spectra
21     C darwin_surfacespecFile :: file name of spectrum of incident light
22     C darwin_acdomFile :: file name of CDOM absorption spectrum
23 jahn 1.3 C darwin_PARunits :: units of data in PAR file; one of Ein/m2/d, uEin/m2/s, W/m2
24     C darwin_PARFileConv :: conversion factor from PAR file units to uEin/m2/s
25     C darwin_W_to_uEins :: conversion factor from W to uEin/s
26 jahn 1.4 C darwin_ironUnits :: units of data in iron file; one of mol Fe/m2/s, kg dust/m2/s
27     C darwin_dustFrac :: fraction of dust by weight that is iron
28 jahn 1.1 C darwin_forcingPeriod :: perioDARWIN forcing parameter specific for DARWIN (seconds)
29     C darwin_forcingCycle :: perioDARWIN forcing parameter specific for DARWIN (seconds)
30     C darwin_seed :: seed for the random number generator
31    
32     COMMON /DARWIN_FILENAMES/
33     & darwin_iceFile,
34     & darwin_ironFile,
35     & darwin_PARFile,
36     & darwin_nutWVelFile,
37     & darwin_PO4_relaxFile, darwin_NO3_relaxFile,
38     & darwin_FeT_relaxFile, darwin_Si_relaxFile,
39     & darwin_relaxscale,
40     & darwin_PO4_fluxFile, darwin_NO3_FluxFile,
41     & darwin_FeT_fluxFile, darwin_Si_fluxFile,
42     & darwin_oasim_edFile, darwin_oasim_esFile,
43     & darwin_waterabsorbFile, darwin_phytoabsorbFile,
44     & darwin_surfacespecFile, darwin_acdomFile,
45     & darwin_particleabsorbFile,
46 jahn 1.4 & darwin_PARunits, darwin_ironUnits,
47     & darwin_PARFileConv, darwin_ironFileConv,
48     & darwin_dustFrac,
49 jahn 1.3 & darwin_W_to_uEins,
50 jahn 1.1 & darwin_forcingPeriod, darwin_forcingCycle,
51     & darwin_seed
52    
53     CHARACTER*(MAX_LEN_FNAM) darwin_iceFile
54     CHARACTER*(MAX_LEN_FNAM) darwin_ironFile
55     CHARACTER*(MAX_LEN_FNAM) darwin_PARFile
56     CHARACTER*(MAX_LEN_FNAM) darwin_NutWVelFile
57     CHARACTER*(MAX_LEN_FNAM) darwin_PO4_relaxFile
58     CHARACTER*(MAX_LEN_FNAM) darwin_NO3_relaxFile
59     CHARACTER*(MAX_LEN_FNAM) darwin_FeT_relaxFile
60     CHARACTER*(MAX_LEN_FNAM) darwin_Si_relaxFile
61     CHARACTER*(MAX_LEN_FNAM) darwin_PO4_fluxFile
62     CHARACTER*(MAX_LEN_FNAM) darwin_NO3_fluxFile
63     CHARACTER*(MAX_LEN_FNAM) darwin_FeT_fluxFile
64     CHARACTER*(MAX_LEN_FNAM) darwin_Si_fluxFile
65     CHARACTER*(MAX_LEN_FNAM) darwin_oasim_edFile
66     CHARACTER*(MAX_LEN_FNAM) darwin_oasim_esFile
67     CHARACTER*(MAX_LEN_FNAM) darwin_waterabsorbFile
68     CHARACTER*(MAX_LEN_FNAM) darwin_phytoabsorbFile
69     CHARACTER*(MAX_LEN_FNAM) darwin_particleabsorbFile
70     CHARACTER*(MAX_LEN_FNAM) darwin_surfacespecFile
71     CHARACTER*(MAX_LEN_FNAM) darwin_acdomFile
72 jahn 1.3 CHARACTER*(MAX_LEN_FNAM) darwin_PARunits
73 jahn 1.4 CHARACTER*(MAX_LEN_FNAM) darwin_ironUnits
74 jahn 1.1
75     _RL darwin_relaxscale
76     _RL darwin_forcingPeriod
77     _RL darwin_forcingCycle
78 jahn 1.3 _RL darwin_PARFileConv
79     _RL darwin_W_to_uEins
80 jahn 1.4 _RL darwin_ironFileConv
81     _RL darwin_dustFrac
82 jahn 1.1 INTEGER darwin_seed
83    
84     c INPUT:
85     c fice - ice fraction
86     c inputFe - aeolian input of iron
87     c sur_par - surface PAR
88     COMMON /DARWIN_INPUT/
89     & fice, inputFe, sur_par
90     #ifdef NUT_SUPPLY
91     & ,nut_wvel
92     #endif
93     #ifdef RELAX_NUTS
94     & ,po4_obs, no3_obs, fet_obs, si_obs
95     #endif
96     #ifdef FLUX_NUTS
97     & ,po4_flx, no3_flx, fet_flx, si_flx
98     #endif
99     c
100     _RL fice(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx, nSy)
101     _RL inputFe(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx, nSy)
102     _RL sur_par(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx, nSy)
103     #ifdef NUT_SUPPLY
104     _RL nut_wvel(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy)
105     #endif
106     #ifdef RELAX_NUTS
107     _RL po4_obs(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy)
108     _RL no3_obs(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy)
109     _RL fet_obs(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy)
110     _RL si_obs(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy)
111     #endif
112     #ifdef FLUX_NUTS
113     _RL po4_flx(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy)
114     _RL no3_flx(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy)
115     _RL fet_flx(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy)
116     _RL si_flx(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy)
117     #endif
118     c
119     c OUPUT DIAGNOSTICS
120     c PPave - average primary production
121     c Nfixave - average N fixation
122     c Zoograzave - average zooplankton consumption
123     c Parave - average PAR
124     c Chlave - average chlorophyll
125     c npzd_timeave - time for averaging
126    
127     COMMON /DARWIN_OUPUT/
128     & PPave, Nfixave, Zoograzave,
129     & PARave, Chlave, Denitave,
130     c ANNA_TAVE
131     #ifdef WAVES_DIAG_PCHL
132     & Pchlave,
133     #endif
134     #ifdef DAR_DIAG_ACDOM
135     & aCDOMave,
136     #endif
137     c ANNA end TAVE
138     #ifdef DAR_DIAG_RSTAR
139     & Rstarave, RNstarave,
140     #endif
141     #ifdef DAR_DIAG_DIVER
142     & Diver1ave, Diver2ave, Diver3ave, Diver4ave,
143     #endif
144     #ifdef DAR_DIAG_GROW
145     & Growave, Growsqave,
146     #endif
147     #ifdef ALLOW_DIAZ
148     #ifdef DAR_DIAG_NFIXP
149     & NfixPave,
150     #endif
151     #endif
152     & dar_timeave
153     c
154     _RL PPave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy)
155     _RL Nfixave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy)
156     _RL Zoograzave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy)
157     _RL PARave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy)
158     _RL Chlave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy)
159     _RL Denitave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy)
160     c ANNA_TAVE
161     #ifdef WAVES_DIAG_PCHL
162     _RL Pchlave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy,npmax)
163     #endif
164     #ifdef DAR_DIAG_ACDOM
165     _RL aCDOMave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx,nSy)
166     #endif
167     c ANNA end TAVE
168     #ifdef DAR_DIAG_RSTAR
169     _RL Rstarave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy,npmax)
170     _RL RNstarave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy,npmax)
171     #endif
172     #ifdef DAR_DIAG_DIVER
173     _RL Diver1ave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy)
174     _RL Diver2ave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy)
175     _RL Diver3ave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy)
176     _RL Diver4ave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy)
177     #endif
178     #ifdef DAR_DIAG_GROW
179     _RL Growave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy,npmax)
180     _RL Growsqave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy,npmax)
181     #endif
182     #ifdef ALLOW_DIAZ
183     #ifdef DAR_DIAG_NFIXP
184     _RL NfixPave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy,npmax)
185     #endif
186     #endif
187     _RL DAR_timeave(nSx,nSy,nR)
188    
189     #ifdef CHECK_CONS
190     COMMON /DARWIN_CHECK_CONS/
191     & DAR_cons_unit1, DAR_cons_unit2,
192     & DAR_cons_unit3, DAR_cons_unit4,
193     & DAR_cons_unit5, DAR_cons_unit6,
194     & DAR_cons_unit7
195     INTEGER DAR_cons_unit1
196     INTEGER DAR_cons_unit2
197     INTEGER DAR_cons_unit3
198     INTEGER DAR_cons_unit4
199     INTEGER DAR_cons_unit5
200     INTEGER DAR_cons_unit6
201     INTEGER DAR_cons_unit7
202     #endif
203    
204 stephd 1.2 #ifdef CALC_RATE_TOTALS
205     COMMON /DARWIN_RATE_TOT/
206     & DAR_ratetot_unit, DAR_fluxtot_unit
207     INTEGER DAR_ratetot_unit
208     INTEGER DAR_fluxtot_unit
209     #endif
210    
211 jahn 1.1 #ifdef ALLOW_PAR_DAY
212     C PARday :: array for accumulating/storing daily-averaged PAR
213     COMMON /DARWIN_PAR_DAY/ PARday,
214     & darwin_PARavPeriod, darwin_PARnav
215     _RL PARday(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx,nSy,2)
216     _RL darwin_PARavPeriod
217     INTEGER darwin_PARnav
218     #endif
219    
220     #ifdef ALLOW_CARBON
221     C-- COMMON /DIC_FILENAMES/
222     C DIC_windFile :: file name of wind speeds
223     C DIC_atmospFile :: file name of atmospheric pressure
224     C dic_pCO2 :: Atmospheric pCO2 to be rad in data.dic
225     C dic_int* :: place holder to read in a integer number, set at run time
226    
227     COMMON /DIC_FILENAMES/
228     & DIC_windFile, DIC_atmospFile,
229     & dic_pCO2, dic_int1, dic_int2, dic_int3, dic_int4
230     CHARACTER*(MAX_LEN_FNAM) DIC_windFile
231     CHARACTER*(MAX_LEN_FNAM) DIC_atmospFile
232     _RL dic_pCO2
233     INTEGER dic_int1
234     INTEGER dic_int2
235     INTEGER dic_int3
236     INTEGER dic_int4
237     #endif
238    

  ViewVC Help
Powered by ViewVC 1.1.22