/[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.9 - (hide annotations) (download)
Fri Dec 27 17:29:00 2013 UTC (11 years, 6 months ago) by jahn
Branch: MAIN
CVS Tags: ctrb_darwin2_ckpt65j_20150225, ctrb_darwin2_ckpt65_20140718, ctrb_darwin2_ckpt65e_20140929, ctrb_darwin2_ckpt64v_20140411, ctrb_darwin2_ckpt64z_20140711, ctrb_darwin2_ckpt65l_20150504, ctrb_darwin2_ckpt64y_20140622, ctrb_darwin2_ckpt65d_20140915, ctrb_darwin2_ckpt64t_20140202, ctrb_darwin2_ckpt64s_20140105, ctrb_darwin2_ckpt64x_20140524, ctrb_darwin2_ckpt65g_20141120, ctrb_darwin2_ckpt65k_20150402, ctrb_darwin2_ckpt64w_20140502, ctrb_darwin2_ckpt65f_20141014, ctrb_darwin2_ckpt64u_20140308, ctrb_darwin2_ckpt65i_20150123, ctrb_darwin2_ckpt65a_20140728, ctrb_darwin2_ckpt65b_20140812, ctrb_darwin2_ckpt65c_20140830, ctrb_darwin2_ckpt65h_20141217
Changes since 1.8: +2 -2 lines
File MIME type: text/plain
use simpler (no level index) cumulative-time counter for timeave

1 jahn 1.9 C $Header: /u/gcmpack/MITgcm_contrib/darwin2/pkg/darwin/DARWIN_IO.h,v 1.8 2013/12/04 21:27:54 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 jahn 1.8 COMMON/DARWIN_IO_L/
85     & darwin_useiceFile,
86     & darwin_useareamask,
87     & darwin_useiceMask,
88     & darwin_useAREA,
89     & darwin_haveIce
90     LOGICAL darwin_useiceFile
91     LOGICAL darwin_useareamask
92     LOGICAL darwin_useiceMask
93     LOGICAL darwin_useAREA
94     LOGICAL darwin_haveIce
95    
96 jahn 1.1 c INPUT:
97     c fice - ice fraction
98     c inputFe - aeolian input of iron
99     c sur_par - surface PAR
100     COMMON /DARWIN_INPUT/
101     & fice, inputFe, sur_par
102     #ifdef NUT_SUPPLY
103     & ,nut_wvel
104     #endif
105     #ifdef RELAX_NUTS
106     & ,po4_obs, no3_obs, fet_obs, si_obs
107     #endif
108     #ifdef FLUX_NUTS
109     & ,po4_flx, no3_flx, fet_flx, si_flx
110     #endif
111     c
112     _RL fice(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx, nSy)
113     _RL inputFe(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx, nSy)
114     _RL sur_par(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx, nSy)
115     #ifdef NUT_SUPPLY
116     _RL nut_wvel(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy)
117     #endif
118     #ifdef RELAX_NUTS
119     _RL po4_obs(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy)
120     _RL no3_obs(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy)
121     _RL fet_obs(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy)
122     _RL si_obs(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy)
123     #endif
124     #ifdef FLUX_NUTS
125     _RL po4_flx(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy)
126     _RL no3_flx(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy)
127     _RL fet_flx(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy)
128     _RL si_flx(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy)
129     #endif
130     c
131     c OUPUT DIAGNOSTICS
132     c PPave - average primary production
133     c Nfixave - average N fixation
134     c Zoograzave - average zooplankton consumption
135     c Parave - average PAR
136     c Chlave - average chlorophyll
137 benw 1.5 c npzd_timeave - time for averaging
138     c Nlim,Flim,etc- Quota model average limitation factors
139 jahn 1.1
140     COMMON /DARWIN_OUPUT/
141     & PPave, Nfixave, Zoograzave,
142     & PARave, Chlave, Denitave,
143 benw 1.5 #ifdef QUOTA_DIAG_LIMIT
144     & Nlimave,Flimave,Ilimave,Tlimave,
145     #endif
146 jahn 1.1 c ANNA_TAVE
147     #ifdef WAVES_DIAG_PCHL
148     & Pchlave,
149     #endif
150     #ifdef DAR_DIAG_ACDOM
151     & aCDOMave,
152     #endif
153 stephd 1.6 #ifdef DAR_DIAG_PARW
154     & chl2cave,
155     #endif
156     #ifdef DAR_DIAG_EK
157     & Ekave, EkoverEave, acclimave,
158     #endif
159 jahn 1.1 c ANNA end TAVE
160     #ifdef DAR_DIAG_RSTAR
161     & Rstarave, RNstarave,
162     #endif
163     #ifdef DAR_DIAG_DIVER
164     & Diver1ave, Diver2ave, Diver3ave, Diver4ave,
165     #endif
166     #ifdef DAR_DIAG_GROW
167     & Growave, Growsqave,
168     #endif
169     #ifdef ALLOW_DIAZ
170     #ifdef DAR_DIAG_NFIXP
171     & NfixPave,
172     #endif
173     #endif
174     & dar_timeave
175     c
176     _RL PPave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy)
177     _RL Nfixave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy)
178     _RL Zoograzave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy)
179     _RL PARave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy)
180     _RL Chlave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy)
181     _RL Denitave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy)
182 benw 1.5 #ifdef QUOTA_DIAG_LIMIT
183     _RL Nlimave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy, npmax)
184     _RL Flimave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy, npmax)
185     _RL Ilimave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy, npmax)
186     _RL Tlimave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy)
187     #endif
188 jahn 1.1 c ANNA_TAVE
189     #ifdef WAVES_DIAG_PCHL
190     _RL Pchlave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy,npmax)
191     #endif
192     #ifdef DAR_DIAG_ACDOM
193     _RL aCDOMave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx,nSy)
194     #endif
195 stephd 1.6 #ifdef DAR_DIAG_PARW
196     _RL chl2cave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy,npmax)
197     #endif
198     #ifdef DAR_DIAG_EK
199     _RL Ekave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy,npmax)
200     _RL EkoverEave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy,npmax)
201     _RL acclimave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy,npmax)
202     #endif
203 jahn 1.1 c ANNA end TAVE
204     #ifdef DAR_DIAG_RSTAR
205     _RL Rstarave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy,npmax)
206     _RL RNstarave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy,npmax)
207     #endif
208     #ifdef DAR_DIAG_DIVER
209     _RL Diver1ave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy)
210     _RL Diver2ave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy)
211     _RL Diver3ave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy)
212     _RL Diver4ave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy)
213     #endif
214     #ifdef DAR_DIAG_GROW
215     _RL Growave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy,npmax)
216     _RL Growsqave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy,npmax)
217     #endif
218     #ifdef ALLOW_DIAZ
219     #ifdef DAR_DIAG_NFIXP
220     _RL NfixPave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx, nSy,npmax)
221     #endif
222     #endif
223 jahn 1.9 _RL DAR_timeave(nSx,nSy)
224 jahn 1.1
225     #ifdef CHECK_CONS
226     COMMON /DARWIN_CHECK_CONS/
227 jahn 1.7 & DAR_cons_unitC, DAR_cons_unitN,
228     & DAR_cons_unitP, DAR_cons_unitF,
229     & DAR_cons_unitS, DAR_cons_unitA,
230     & DAR_cons_unitO
231     INTEGER DAR_cons_unitC
232     INTEGER DAR_cons_unitN
233     INTEGER DAR_cons_unitP
234     INTEGER DAR_cons_unitF
235     INTEGER DAR_cons_unitS
236     INTEGER DAR_cons_unitA
237     INTEGER DAR_cons_unitO
238 jahn 1.1 #endif
239    
240 stephd 1.2 #ifdef CALC_RATE_TOTALS
241     COMMON /DARWIN_RATE_TOT/
242     & DAR_ratetot_unit, DAR_fluxtot_unit
243     INTEGER DAR_ratetot_unit
244     INTEGER DAR_fluxtot_unit
245     #endif
246    
247 jahn 1.1 #ifdef ALLOW_PAR_DAY
248     C PARday :: array for accumulating/storing daily-averaged PAR
249     COMMON /DARWIN_PAR_DAY/ PARday,
250     & darwin_PARavPeriod, darwin_PARnav
251     _RL PARday(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx,nSy,2)
252     _RL darwin_PARavPeriod
253     INTEGER darwin_PARnav
254     #endif
255    
256     #ifdef ALLOW_CARBON
257     C-- COMMON /DIC_FILENAMES/
258     C DIC_windFile :: file name of wind speeds
259     C DIC_atmospFile :: file name of atmospheric pressure
260     C dic_pCO2 :: Atmospheric pCO2 to be rad in data.dic
261     C dic_int* :: place holder to read in a integer number, set at run time
262    
263     COMMON /DIC_FILENAMES/
264     & DIC_windFile, DIC_atmospFile,
265     & dic_pCO2, dic_int1, dic_int2, dic_int3, dic_int4
266     CHARACTER*(MAX_LEN_FNAM) DIC_windFile
267     CHARACTER*(MAX_LEN_FNAM) DIC_atmospFile
268     _RL dic_pCO2
269     INTEGER dic_int1
270     INTEGER dic_int2
271     INTEGER dic_int3
272     INTEGER dic_int4
273     #endif
274    

  ViewVC Help
Powered by ViewVC 1.1.22