source: tags/LIM1D_v3.20/GRAPHICS/open_fields_v7.pro @ 6

Last change on this file since 6 was 6, checked in by vancop, 8 years ago

initial import of v3.20 /Users/ioulianikolskaia/Boulot/CODES/LIM1D/ARCHIVE/TMP/LIM1D_v3.20/
File size: 11.5 KB
Line 
1
2PRO open_fields_v7, indir, exp_id, c_bio_model,            $
3             numt, doy, ts_d, ts_m,                        $
4             h_i, h_s, z_ip, z_ib, s_i, t_i, e_i, PAR, Ra, $
5             dhib, dhisu, dhisi,                           $
6             AoCb, dINb, dIPb, dSib, chla, eoCb,           $
7             Argb, Argbub,                                 $
8             Oxyb, Oxybub,                                 $
9             DICb, Alkb, CO2b, CO2bub,                     $
10             CO2aq, HCO3m, CO32m, pH, pCO2,                $
11             Ikab,                                         $
12             dFeb, aFeb, eFeb,                             $
13             AoCt, dINt, dIPt, dSit, eoCt, Argt,           $
14             DICt, Alkt, CO2t, Ikat,                       $
15             dFet, aFet, eFet,                             $
16             FAoC, FAoCb, FAoCsi,                          $
17             FdIN, FdINb, FdINsi,                          $
18             FCO2_atm, FCO2_bub,                           $
19             syn, rsp, lys, rem,                           $
20             lim_lig, lim_no3, lim_po4, lim_dsi, lim_tem, lim_sal
21
22; working for source_3.05
23; Helsinki, 2014
24; Bruxelles, nov 2014 avec Julie
25; this little nice script reads output from a bio simulation and extract stuff from it
26;
27
28; IN
29; indir, exp_id
30; OUT
31; remainder
32;==================================================================================================
33
34PRINT, ' indir = ', indir
35PRINT, ' exp_id= ', exp_id
36file_no = NCDF_OPEN(indir+"ice.out.nc."+exp_id, /NOWRITE)
37
38numt_no = NCDF_VARID( file_no, 'numt' )
39doy_no  = NCDF_VARID( file_no, 'doy' )
40ts_m_no = NCDF_VARID( file_no, 'ts_m' )
41ts_d_no = NCDF_VARID( file_no, 'ts_d' )
42
43h_i_no  = NCDF_VARID( file_no, 'h_i'  )
44h_s_no  = NCDF_VARID( file_no, 'h_s'  )
45z_ip_no = NCDF_VARID( file_no, 'z_ip' )
46z_ib_no = NCDF_VARID( file_no, 'z_ib' )
47s_i_no  = NCDF_VARID( file_no, 's_i' )
48t_i_no  = NCDF_VARID( file_no, 't_i' )
49e_i_no  = NCDF_VARID( file_no, 'e_i' )
50PAR_no  = NCDF_VARID( file_no, 'PAR' )
51Ra_no   = NCDF_VARID( file_no, 'ra' )
52
53dhib_no = NCDF_VARID( file_no, 'dhib' )
54dhisu_no = NCDF_VARID( file_no, 'dhisu' )
55dhisi_no = NCDF_VARID( file_no, 'dhisi' )
56
57IF ( c_bio_model EQ 'KRILL' ) THEN BEGIN
58
59AoCb_no = NCDF_VARID( file_no, 'AoCb' )
60dINb_no = NCDF_VARID( file_no, 'dINb' )
61dIPb_no = NCDF_VARID( file_no, 'dIPb' )
62dSib_no = NCDF_VARID( file_no, 'dSib' )
63chla_no = NCDF_VARID( file_no, 'Chla_bio' )
64eoCb_no = NCDF_VARID( file_no, 'eoCb' )
65Argb_no   = NCDF_VARID( file_no, 'Argb' )
66Argbub_no = NCDF_VARID( file_no, 'Argbub' )
67Oxyb_no   = NCDF_VARID( file_no, 'Oxyb' )
68Oxybub_no = NCDF_VARID( file_no, 'Oxybub' )
69DICb_no = NCDF_VARID( file_no, 'DICb' )
70Alkb_no = NCDF_VARID( file_no, 'Alkb' )
71CO2b_no = NCDF_VARID( file_no, 'CO2b' )
72CO2bub_no = NCDF_VARID( file_no, 'CO2bub' )
73
74CO2aq_no = NCDF_VARID( file_no, 'CO2aq' )
75HCO3m_no = NCDF_VARID( file_no, 'HCO3m' )
76CO32m_no = NCDF_VARID( file_no, 'CO32m' )
77pH_no    = NCDF_VARID( file_no, 'pH'   )
78pCO2_no  = NCDF_VARID( file_no, 'pCO2' )
79Ikab_no  = NCDF_VARID( file_no, 'Ikab' )
80dFeb_no  = NCDF_VARID( file_no, 'dFeb' )
81aFeb_no  = NCDF_VARID( file_no, 'aFeb' )
82eFeb_no  = NCDF_VARID( file_no, 'eFeb' )
83
84AoCt_no = NCDF_VARID( file_no, 'AoCt' )
85dINt_no = NCDF_VARID( file_no, 'dINt' )
86dIPt_no = NCDF_VARID( file_no, 'dIPt' )
87dSit_no = NCDF_VARID( file_no, 'dSit' )
88eoCt_no = NCDF_VARID( file_no, 'eoCt' )
89Argt_no = NCDF_VARID( file_no, 'Argt' )
90DICt_no = NCDF_VARID( file_no, 'DICt' )
91Alkt_no = NCDF_VARID( file_no, 'Alkt' )
92CO2t_no = NCDF_VARID( file_no, 'CO2t' )
93Ikat_no = NCDF_VARID( file_no, 'Ikat' )
94
95dFet_no = NCDF_VARID( file_no, 'dFet' )
96aFet_no = NCDF_VARID( file_no, 'aFet' )
97eFet_no = NCDF_VARID( file_no, 'eFet' )
98
99FAoC_no = NCDF_VARID( file_no, 'FAoC' )
100FAoCb_no = NCDF_VARID( file_no, 'FAoCb' )
101FAoCsi_no = NCDF_VARID( file_no, 'FAoCsi' )
102
103FdIN_no = NCDF_VARID( file_no, 'FdIN' )
104FdINb_no = NCDF_VARID( file_no, 'FdINb' )
105FdINsi_no = NCDF_VARID( file_no, 'FdINsi' )
106
107FCO2_atm_no = NCDF_VARID( file_no, 'FCO2_atm' )
108FCO2_bub_no = NCDF_VARID( file_no, 'FCO2_bub' )
109
110syn_no  = NCDF_VARID( file_no, 'syn' )
111rsp_no  = NCDF_VARID( file_no, 'rsp' )
112lys_no  = NCDF_VARID( file_no, 'lys' )
113rem_no  = NCDF_VARID( file_no, 'rem' )
114
115lim_lig_no = NCDF_VARID( file_no, 'lim_lig' )
116lim_no3_no = NCDF_VARID( file_no, 'lim_no3' )
117lim_po4_no = NCDF_VARID( file_no, 'lim_po4' )
118lim_dsi_no = NCDF_VARID( file_no, 'lim_dsi' )
119lim_tem_no = NCDF_VARID( file_no, 'lim_tem' )
120lim_sal_no = NCDF_VARID( file_no, 'lim_sal' )
121
122ENDIF
123
124IF ( c_bio_model EQ 'BFMSI' ) THEN BEGIN
125
126S1cb_no = NCDF_VARID( file_no, 'S1cb' )
127S2cb_no = NCDF_VARID( file_no, 'S2cb' )
128dIPb_no = NCDF_VARID( file_no, 'I1pb' )
129dINb_no = NCDF_VARID( file_no, 'I3nb' )
130dSib_no = NCDF_VARID( file_no, 'I5sb' )
131chla1_no = NCDF_VARID( file_no, 'S1lb' )
132chla2_no = NCDF_VARID( file_no, 'S1lb' )
133
134S1ct_no = NCDF_VARID( file_no, 'S1ct' )
135S2ct_no = NCDF_VARID( file_no, 'S2ct' )
136dIPt_no = NCDF_VARID( file_no, 'I1pb' )
137dINt_no = NCDF_VARID( file_no, 'I3nt' )
138dSit_no = NCDF_VARID( file_no, 'I5st' )
139
140FS1c_no = NCDF_VARID( file_no, 'FS1c' )
141FS1cb_no = NCDF_VARID( file_no, 'FS1cb' )
142FS1csi_no = NCDF_VARID( file_no, 'FS1csi' )
143FS1c_no = NCDF_VARID( file_no, 'FS2c' )
144FS1cb_no = NCDF_VARID( file_no, 'FS2cb' )
145FS1csi_no = NCDF_VARID( file_no, 'FS2csi' )
146
147FdIN_no = NCDF_VARID( file_no, 'FI3n' )
148FdINb_no = NCDF_VARID( file_no, 'FI3nb' )
149FdINsi_no = NCDF_VARID( file_no, 'FI3nsi' )
150
151syn_no  = NCDF_VARID( file_no, 'rugC' )
152lys_no  = NCDF_VARID( file_no, 'slc' )
153;rem_no  = NCDF_VARID( file_no, 'rem' )
154
155lim_lig_no = NCDF_VARID( file_no, 'eiSI' )
156lim_no3_no = NCDF_VARID( file_no, 'lim_no3' )
157lim_po4_no = NCDF_VARID( file_no, 'lim_po4' )
158lim_dsi_no = NCDF_VARID( file_no, 'ei5s' )
159lim_tem_no = NCDF_VARID( file_no, 'et' )
160
161ENDIF
162
163IF ( numt_no GT 0 ) THEN NCDF_VARGET, file_no, numt_no, numt
164IF ( doy_no  GT 0 ) THEN NCDF_VARGET, file_no, doy_no, doy
165IF ( ts_m_no GT 0 ) THEN NCDF_VARGET, file_no, ts_m_no, ts_m
166IF ( ts_d_no GT 0 ) THEN NCDF_VARGET, file_no, ts_d_no, ts_d
167
168IF ( h_i_no  GT 0 ) THEN NCDF_VARGET, file_no, h_i_no, h_i
169IF ( h_s_no  GT 0 ) THEN NCDF_VARGET, file_no, h_s_no, h_s
170IF ( z_ip_no GT 0 ) THEN NCDF_VARGET, file_no, z_ip_no, z_ip
171IF ( z_ib_no GT 0 ) THEN NCDF_VARGET, file_no, z_ib_no, z_ib
172IF ( s_i_no  GT 0 ) THEN NCDF_VARGET, file_no, s_i_no, s_i
173IF ( t_i_no  GT 0 ) THEN NCDF_VARGET, file_no, t_i_no, t_i
174IF ( e_i_no  GT 0 ) THEN NCDF_VARGET, file_no, e_i_no, e_i
175IF ( PAR_no  GT 0 ) THEN NCDF_VARGET, file_no, PAR_no, PAR
176IF ( Ra_no   GT 0 ) THEN NCDF_VARGET, file_no, Ra_no, Ra
177
178IF ( dhib_no GT 0 ) THEN NCDF_VARGET, file_no, dhib_no , dhib
179IF ( dhisu_no GT 0) THEN NCDF_VARGET, file_no, dhisu_no, dhisu
180IF ( dhisi_no GT 0) THEN NCDF_VARGET, file_no, dhisi_no, dhisi
181
182IF ( c_bio_model EQ 'KRILL' ) THEN BEGIN
183
184IF ( AoCb_no   GT 0 ) THEN NCDF_VARGET, file_no, AoCb_no, AoCb
185IF ( dINb_no   GT 0 ) THEN NCDF_VARGET, file_no, dINb_no, dINb
186IF ( dIPb_no   GT 0 ) THEN NCDF_VARGET, file_no, dIPb_no, dIPb
187IF ( dsib_no   GT 0 ) THEN NCDF_VARGET, file_no, dsib_no, dsib
188IF ( chla_no   GT 0 ) THEN NCDF_VARGET, file_no, chla_no, chla
189IF ( eoCb_no   GT 0 ) THEN NCDF_VARGET, file_no, eoCb_no, eoCb
190IF ( Argb_no   GT 0 ) THEN NCDF_VARGET, file_no, Argb_no, Argb
191IF ( Argbub_no   GT 0 ) THEN NCDF_VARGET, file_no, Argbub_no, Argbub
192IF ( Oxyb_no   GT 0 ) THEN NCDF_VARGET, file_no, Oxyb_no, Oxyb
193IF ( Oxybub_no   GT 0 ) THEN NCDF_VARGET, file_no, Oxybub_no, Oxybub
194IF ( DICb_no   GT 0 ) THEN NCDF_VARGET, file_no, DICb_no, DICb
195IF ( Alkb_no   GT 0 ) THEN NCDF_VARGET, file_no, Alkb_no, Alkb
196IF ( CO2b_no   GT 0 ) THEN NCDF_VARGET, file_no, CO2b_no, CO2b
197IF ( CO2bub_no   GT 0 ) THEN NCDF_VARGET, file_no, CO2bub_no, CO2bub
198IF ( CO2aq_no  GT 0 ) THEN NCDF_VARGET, file_no, CO2aq_no, CO2aq
199IF ( HCO3m_no  GT 0 ) THEN NCDF_VARGET, file_no, HCO3m_no, HCO3m
200IF ( CO32m_no  GT 0 ) THEN NCDF_VARGET, file_no, CO32m_no, CO32m
201IF ( pCO2_no   GT 0 ) THEN NCDF_VARGET, file_no, pCO2_no, pCO2
202IF ( Ikab_no   GT 0 ) THEN NCDF_VARGET, file_no, Ikab_no, Ikab
203IF ( dFeb_no   GT 0 ) THEN NCDF_VARGET, file_no, dFeb_no, dFeb
204IF ( aFeb_no   GT 0 ) THEN NCDF_VARGET, file_no, aFeb_no, aFeb
205IF ( eFeb_no   GT 0 ) THEN NCDF_VARGET, file_no, eFeb_no, eFeb
206IF ( pH_no     GT 0 ) THEN NCDF_VARGET, file_no, pH_no, pH
207
208IF ( AoCt_no   GT 0 ) THEN NCDF_VARGET, file_no, AoCt_no, AoCt
209IF ( dINt_no   GT 0 ) THEN NCDF_VARGET, file_no, dINt_no, dINt
210IF ( dIPt_no   GT 0 ) THEN NCDF_VARGET, file_no, dIPt_no, dIPt
211IF ( dsit_no   GT 0 ) THEN NCDF_VARGET, file_no, dsit_no, dsit
212IF ( eoCt_no   GT 0 ) THEN NCDF_VARGET, file_no, eoCt_no, eoCt
213IF ( Argt_no   GT 0 ) THEN NCDF_VARGET, file_no, Argt_no, Argt
214IF ( DICt_no   GT 0 ) THEN NCDF_VARGET, file_no, DICt_no, DICt
215IF ( Alkt_no   GT 0 ) THEN NCDF_VARGET, file_no, Alkt_no, Alkt
216IF ( CO2t_no   GT 0 ) THEN NCDF_VARGET, file_no, CO2t_no, CO2t
217IF ( Ikat_no   GT 0 ) THEN NCDF_VARGET, file_no, Ikat_no, Ikat
218IF ( dFet_no   GT 0 ) THEN NCDF_VARGET, file_no, dFet_no, dFet
219IF ( aFet_no   GT 0 ) THEN NCDF_VARGET, file_no, aFet_no, aFet
220IF ( eFet_no   GT 0 ) THEN NCDF_VARGET, file_no, eFet_no, eFet
221
222If ( FAoC_no   GT 0 ) THEN NCDF_VARGET, file_no, FAoC_no, FAoC
223IF ( FAoCb_no  GT 0 ) THEN NCDF_VARGET, file_no, FAoCb_no, FAoCb
224IF ( FAoCsi_no GT 0 ) THEN NCDF_VARGET, file_no, FAoCsi_no, FAoCsi
225
226IF ( FdIN_no   GT 0 ) THEN NCDF_VARGET, file_no, FdIN_no, FdIN
227IF ( FdINb_no  GT 0 ) THEN NCDF_VARGET, file_no, FdINb_no, FdINb
228IF ( FdINsi_no GT 0 ) THEN NCDF_VARGET, file_no, FdINsi_no, FdINsi
229IF ( FCO2_atm_no GT 0 ) THEN NCDF_VARGET, file_no, FCO2_atm_no, FCO2_atm
230IF ( FCO2_bub_no GT 0 ) THEN NCDF_VARGET, file_no, FCO2_bub_no, FCO2_bub
231
232IF ( syn_no    GT 0 ) THEN NCDF_VARGET, file_no, syn_no, syn
233IF ( rsp_no    GT 0 ) THEN NCDF_VARGET, file_no, rsp_no, rsp
234IF ( lys_no    GT 0 ) THEN NCDF_VARGET, file_no, lys_no, lys
235IF ( rem_no    GT 0 ) THEN NCDF_VARGET, file_no, rem_no, rem
236
237IF ( lim_lig_no GT 0 ) THEN NCDF_VARGET, file_no, lim_lig_no, lim_lig
238IF ( lim_no3_no GT 0 ) THEN NCDF_VARGET, file_no, lim_no3_no, lim_no3
239IF ( lim_po4_no GT 0 ) THEN NCDF_VARGET, file_no, lim_po4_no, lim_po4
240IF ( lim_dsi_no GT 0 ) THEN NCDF_VARGET, file_no, lim_dsi_no, lim_dsi
241IF ( lim_tem_no GT 0 ) THEN NCDF_VARGET, file_no, lim_tem_no, lim_tem
242IF ( lim_sal_no GT 0 ) THEN NCDF_VARGET, file_no, lim_sal_no, lim_sal
243
244ENDIF
245
246IF ( c_bio_model EQ 'BFMSI' ) THEN BEGIN
247
248NCDF_VARGET, file_no, S1cb_no, S1cb
249NCDF_VARGET, file_no, S2cb_no, S2cb
250AoCb = S1cb + S2cb
251NCDF_VARGET, file_no, dINb_no, dINb
252NCDF_VARGET, file_no, dIPb_no, dIPb
253NCDF_VARGET, file_no, dsib_no, dsib
254NCDF_VARGET, file_no, chla1_no, chla1
255NCDF_VARGET, file_no, chla2_no, chla2
256chla = chla1 + chla2
257
258NCDF_VARGET, file_no, S1ct_no, S1ct
259NCDF_VARGET, file_no, S2ct_no, S2ct
260AoCT = S1ct + S2ct
261NCDF_VARGET, file_no, dINt_no, dINt
262NCDF_VARGET, file_no, dIPt_no, dIPt
263NCDF_VARGET, file_no, dsit_no, dsit
264
265NCDF_VARGET, file_no, FS1c_no, FS1c
266NCDF_VARGET, file_no, FS1cb_no, FS1cb
267NCDF_VARGET, file_no, FS1csi_no, FS1csi
268NCDF_VARGET, file_no, FS1c_no, FS2c
269NCDF_VARGET, file_no, FS1cb_no, FS2cb
270NCDF_VARGET, file_no, FS1csi_no, FS2csi
271
272FAoC   = FS1c + FS2c
273FAoCb  = FS1cb + FS2cb
274FAoCsi = FS1csi + FS2csi
275
276NCDF_VARGET, file_no, FdIN_no, FdIN
277NCDF_VARGET, file_no, FdINb_no, FdINb
278NCDF_VARGET, file_no, FdINsi_no, FdINsi
279
280NCDF_VARGET, file_no, syn_no, syn
281NCDF_VARGET, file_no, rsp_no, rsp
282NCDF_VARGET, file_no, lys_no, lys
283NCDF_VARGET, file_no, rem_no, rem
284
285NCDF_VARGET, file_no, lim_lig_no, lim_lig
286NCDF_VARGET, file_no, lim_no3_no, lim_no3
287NCDF_VARGET, file_no, lim_po4_no, lim_po4
288NCDF_VARGET, file_no, lim_dsi_no, lim_dsi
289NCDF_VARGET, file_no, lim_tem_no, lim_tem
290;NCDF_VARGET, file_no, lim_sal_no, lim_sal
291lim_sal = lim_lig & lim_sal[*] = 1
292
293ENDIF
294
295NCDF_CLOSE, file_no
296
297;==================================================================================================
298
299END
Note: See TracBrowser for help on using the repository browser.