REMORA
Regional Modeling of Oceans Refined Adaptively
Loading...
Searching...
No Matches
REMORA_ReadFromInitNetcdf.cpp
Go to the documentation of this file.
1#include "REMORA_NCFile.H"
2#include "AMReX_FArrayBox.H"
3#include "REMORA_DataStruct.H"
4
5using namespace amrex;
6
7#ifdef REMORA_USE_NETCDF
8/**
9 * @param lev level of data to read
10 * @param domain simulation domain
11 * @param fname file name to read from
12 * @param NC_temp_fab container for temperature data
13 * @param NC_salt_fab container for salinity data
14 * @param NC_u_fab container for u velocity data
15 * @param NC_v_fab container for v velocity data
16 */
17void
18read_data_from_netcdf (int /*lev*/,
19 const Box& domain,
20 const std::string& fname,
21 FArrayBox& NC_temp_fab, FArrayBox& NC_salt_fab,
22 FArrayBox& NC_xvel_fab, FArrayBox& NC_yvel_fab)
23{
24 amrex::Print() << "Loading initial solution data from NetCDF file " << fname << std::endl;
25
26 Vector<FArrayBox*> NC_fabs;
27 Vector<std::string> NC_names;
28 Vector<enum NC_Data_Dims_Type> NC_dim_types;
29
30 NC_fabs.push_back(&NC_temp_fab); NC_names.push_back("temp"); NC_dim_types.push_back(NC_Data_Dims_Type::Time_BT_SN_WE); // 0
31 NC_fabs.push_back(&NC_salt_fab); NC_names.push_back("salt"); NC_dim_types.push_back(NC_Data_Dims_Type::Time_BT_SN_WE); // 1
32 NC_fabs.push_back(&NC_xvel_fab); NC_names.push_back("u"); NC_dim_types.push_back(NC_Data_Dims_Type::Time_BT_SN_WE); // 2
33 NC_fabs.push_back(&NC_yvel_fab); NC_names.push_back("v"); NC_dim_types.push_back(NC_Data_Dims_Type::Time_BT_SN_WE); // 3
34
35 // Read the netcdf file and fill these FABs
36 BuildFABsFromNetCDFFile<FArrayBox,Real>(domain, fname, NC_names, NC_dim_types, NC_fabs);
37}
38
39/**
40 * @param lev level of data to read
41 * @param domain simulation domain
42 * @param fname file name to read from
43 * @param NC_temp_fab container for temperature data
44 * @param NC_salt_fab container for salinity data
45 * @param NC_u_fab container for u velocity data
46 * @param NC_v_fab container for v velocity data
47 * @param ngrow number of grow cells to read
48 */
49void
50read_data_full_domain_from_netcdf (int /*lev*/,
51 const Box& domain,
52 const std::string& fname,
53 FArrayBox& NC_temp_fab, FArrayBox& NC_salt_fab,
54 FArrayBox& NC_xvel_fab, FArrayBox& NC_yvel_fab,
55 IntVect ngrow)
56{
57 amrex::Print() << "Loading initial solution data from NetCDF file " << fname << std::endl;
58
59 Vector<FArrayBox*> NC_fabs;
60 Vector<std::string> NC_names;
61 Vector<enum NC_Data_Dims_Type> NC_dim_types;
62
63 NC_fabs.push_back(&NC_temp_fab); NC_names.push_back("temp"); NC_dim_types.push_back(NC_Data_Dims_Type::Time_BT_SN_WE); // 0
64 NC_fabs.push_back(&NC_salt_fab); NC_names.push_back("salt"); NC_dim_types.push_back(NC_Data_Dims_Type::Time_BT_SN_WE); // 1
65 NC_fabs.push_back(&NC_xvel_fab); NC_names.push_back("u"); NC_dim_types.push_back(NC_Data_Dims_Type::Time_BT_SN_WE); // 2
66 NC_fabs.push_back(&NC_yvel_fab); NC_names.push_back("v"); NC_dim_types.push_back(NC_Data_Dims_Type::Time_BT_SN_WE); // 3
67
68 // Read the netcdf file and fill these FABs
69 BuildFABsFromNetCDFFile<FArrayBox,Real>(domain, fname, NC_names, NC_dim_types, NC_fabs, false, 0, ngrow);
70}
71
72/**
73 * @param lev level of data to read
74 * @param domain simulation domain
75 * @param fname file name to read from
76 * @param NC_zeta_fab container for sea surface height data
77 */
78void
79read_zeta_from_netcdf (int /*lev*/,
80 const Box& domain,
81 const std::string& fname,
82 FArrayBox& NC_zeta_fab)
83{
84 amrex::Print() << "Loading initial sea surface height from NetCDF file " << fname << std::endl;
85
86 Vector<FArrayBox*> NC_fabs;
87 Vector<std::string> NC_names;
88 Vector<enum NC_Data_Dims_Type> NC_dim_types;
89
90 NC_fabs.push_back(&NC_zeta_fab ) ; NC_names.push_back("zeta") ; NC_dim_types.push_back(NC_Data_Dims_Type::Time_SN_WE); // 0
91
92 // Read the netcdf file and fill these FABs
93 BuildFABsFromNetCDFFile<FArrayBox,Real>(domain, fname, NC_names, NC_dim_types, NC_fabs);
94}
95
96/**
97 * @param lev level of data to read
98 * @param domain simulation domain
99 * @param fname file name to read from
100 * @param NC_zeta_fab container for sea surface height data
101 * @param ngrow number of grow cells to read in, if not default
102 */
103void
104read_zeta_full_domain_from_netcdf (int /*lev*/,
105 const Box& domain,
106 const std::string& fname,
107 FArrayBox& NC_zeta_fab,
108 IntVect ngrow)
109{
110 amrex::Print() << "Loading initial sea surface height from NetCDF file " << fname << std::endl;
111
112 Vector<FArrayBox*> NC_fabs;
113 Vector<std::string> NC_names;
114 Vector<enum NC_Data_Dims_Type> NC_dim_types;
115
116 NC_fabs.push_back(&NC_zeta_fab ) ; NC_names.push_back("zeta") ; NC_dim_types.push_back(NC_Data_Dims_Type::Time_SN_WE); // 0
117
118 // Read the netcdf file and fill these FABs
119 BuildFABsFromNetCDFFile<FArrayBox,Real>(domain, fname, NC_names, NC_dim_types, NC_fabs, false, 0, ngrow);
120}
121
122/**
123 * @param lev level of data to read
124 * @param domain simulation domain
125 * @param fname file name to read from
126 * @param NC_h_fab container for bathymetry data
127 */
128void
129read_bathymetry_from_netcdf (int /*lev*/,
130 const Box& domain,
131 const std::string& fname,
132 FArrayBox& NC_h_fab)
133{
134 amrex::Print() << "Loading initial bathymetry from NetCDF file " << fname << std::endl;
135
136 Vector<FArrayBox*> NC_fabs;
137 Vector<std::string> NC_names;
138 Vector<enum NC_Data_Dims_Type> NC_dim_types;
139
140 NC_fabs.push_back(&NC_h_fab ) ; NC_names.push_back("h") ; NC_dim_types.push_back(NC_Data_Dims_Type::SN_WE); // 0
141
142 // Read the netcdf file and fill these FABs
143 BuildFABsFromNetCDFFile<FArrayBox,Real>(domain, fname, NC_names, NC_dim_types, NC_fabs);
144}
145
146/**
147 * @param lev level of data to read
148 * @param domain simulation domain
149 * @param fname file name to read from
150 * @param NC_pm_fab container for pm data
151 * @param NC_pn_fab container for pn data
152 * @param NC_xr_fab container for x_rho data
153 * @param NC_yr_fab container for y_rho data
154 * @param NC_xu_fab container for x_u data
155 * @param NC_yu_fab container for y_u data
156 * @param NC_xv_fab container for x_v data
157 * @param NC_yv_fab container for y_v data
158 * @param NC_xp_fab container for x_p data
159 * @param NC_yp_fab container for y_p data
160 */
161void
162read_grid_vars_from_netcdf (int /*lev*/,
163 const Box& domain,
164 const std::string& fname,
165 FArrayBox& NC_pm_fab, FArrayBox& NC_pn_fab,
166 FArrayBox& NC_xr_fab, FArrayBox& NC_yr_fab,
167 FArrayBox& NC_xu_fab, FArrayBox& NC_yu_fab,
168 FArrayBox& NC_xv_fab, FArrayBox& NC_yv_fab,
169 FArrayBox& NC_xp_fab, FArrayBox& NC_yp_fab)
170{
171 amrex::Print() << "Loading grid variables from NetCDF file " << fname << std::endl;
172
173 Vector<FArrayBox*> NC_fabs;
174 Vector<std::string> NC_names;
175 Vector<enum NC_Data_Dims_Type> NC_dim_types;
176
177 NC_fabs.push_back(&NC_pm_fab) ; NC_names.push_back("pm") ; NC_dim_types.push_back(NC_Data_Dims_Type::SN_WE); // 1
178 NC_fabs.push_back(&NC_pn_fab) ; NC_names.push_back("pn") ; NC_dim_types.push_back(NC_Data_Dims_Type::SN_WE); // 2
179 NC_fabs.push_back(&NC_xr_fab) ; NC_names.push_back("x_rho") ; NC_dim_types.push_back(NC_Data_Dims_Type::SN_WE); // 3
180 NC_fabs.push_back(&NC_yr_fab) ; NC_names.push_back("y_rho") ; NC_dim_types.push_back(NC_Data_Dims_Type::SN_WE); // 4
181 NC_fabs.push_back(&NC_xu_fab) ; NC_names.push_back("x_u") ; NC_dim_types.push_back(NC_Data_Dims_Type::SN_WE); // 5
182 NC_fabs.push_back(&NC_yu_fab) ; NC_names.push_back("y_u") ; NC_dim_types.push_back(NC_Data_Dims_Type::SN_WE); // 6
183 NC_fabs.push_back(&NC_xv_fab) ; NC_names.push_back("x_v") ; NC_dim_types.push_back(NC_Data_Dims_Type::SN_WE); // 7
184 NC_fabs.push_back(&NC_yv_fab) ; NC_names.push_back("y_v") ; NC_dim_types.push_back(NC_Data_Dims_Type::SN_WE); // 8
185 NC_fabs.push_back(&NC_xp_fab) ; NC_names.push_back("x_psi") ; NC_dim_types.push_back(NC_Data_Dims_Type::SN_WE); // 9
186 NC_fabs.push_back(&NC_yp_fab) ; NC_names.push_back("y_psi") ; NC_dim_types.push_back(NC_Data_Dims_Type::SN_WE); // 10
187
188 // Read the netcdf file and fill these FABs
189 BuildFABsFromNetCDFFile<FArrayBox,Real>(domain, fname, NC_names, NC_dim_types, NC_fabs);
190}
191
192/**
193 * @param domain simulation domain
194 * @param fname file name to read from
195 * @param NC_lonp_fab container for lon_psi data
196 * @param NC_latp_fab container for lat_psi data
197 * @returns whether the file carried spherical psi coordinates and they were read
198 *
199 * ROMS grid files for spherical grids carry lon_psi/lat_psi alongside
200 * x_psi/y_psi, and those degree-valued corner arrays are what a coupled driver
201 * needs when the two components do not share a projected Cartesian frame
202 * (SCRIP/COAWST selects the same way; see Lib/SCRIP_COAWST/read_roms.f).
203 * Idealized grid files carry no spherical coordinates, so this read is optional
204 * and reports whether it happened rather than aborting.
205 */
206bool
207read_spherical_grid_vars_from_netcdf (int /*lev*/,
208 const Box& domain,
209 const std::string& fname,
210 FArrayBox& NC_lonp_fab, FArrayBox& NC_latp_fab)
211{
212 Vector<std::string> NC_names;
213 NC_names.push_back("lon_psi");
214 NC_names.push_back("lat_psi");
215
216 if (!QueryNetCDFHasVars(fname, NC_names)) {
217 amrex::Print() << "Grid file " << fname << " has no lon_psi/lat_psi; "
218 << "spherical psi coordinates unavailable" << std::endl;
219 return false;
220 }
221
222 amrex::Print() << "Loading spherical psi coordinates from NetCDF file " << fname << std::endl;
223
224 Vector<FArrayBox*> NC_fabs;
225 Vector<enum NC_Data_Dims_Type> NC_dim_types;
226
227 NC_fabs.push_back(&NC_lonp_fab) ; NC_dim_types.push_back(NC_Data_Dims_Type::SN_WE);
228 NC_fabs.push_back(&NC_latp_fab) ; NC_dim_types.push_back(NC_Data_Dims_Type::SN_WE);
229
230 BuildFABsFromNetCDFFile<FArrayBox,Real>(domain, fname, NC_names, NC_dim_types, NC_fabs);
231 return true;
232}
233
234/**
235 * @param domain simulation domain at nc_hires_grid_level
236 * @param fname file name to read from
237 * @param NC_h_fab container for bathymetry data
238 * @param ngrow number of grow cells to read in, if not default
239 */
240void
241read_bathymetry_full_domain_from_netcdf (const Box& domain,
242 const std::string& fname,
243 FArrayBox& NC_h_fab, IntVect ngrow)
244{
245 amrex::Print() << "Loading high resolution bathymetry from NetCDF file " << fname << std::endl;
246
247 Vector<FArrayBox*> NC_fabs;
248 Vector<std::string> NC_names;
249 Vector<enum NC_Data_Dims_Type> NC_dim_types;
250
251 NC_fabs.push_back(&NC_h_fab ) ; NC_names.push_back("h") ; NC_dim_types.push_back(NC_Data_Dims_Type::SN_WE); // 0
252
253 // Read the netcdf file and fill these FABs
254 BuildFABsFromNetCDFFile<FArrayBox,Real>(domain, fname, NC_names, NC_dim_types, NC_fabs, false, 0, ngrow);
255}
256
257/**
258 * @param lev level of data to read
259 * @param domain simulation domain
260 * @param fname file name to read from
261 * @param NC_pm_fab container for pm data
262 * @param NC_pn_fab container for pn data
263 * @param ngrow number of grow cells to read in, if not default
264 */
265void
266read_grid_vars_full_domain_from_netcdf (const Box& domain,
267 const std::string& fname,
268 FArrayBox& NC_pm_fab, FArrayBox& NC_pn_fab,
269 IntVect ngrow)
270{
271 amrex::Print() << "Loading high resolution grid variables from NetCDF file " << fname << std::endl;
272
273 Vector<FArrayBox*> NC_fabs;
274 Vector<std::string> NC_names;
275 Vector<enum NC_Data_Dims_Type> NC_dim_types;
276
277 NC_fabs.push_back(&NC_pm_fab) ; NC_names.push_back("pm") ; NC_dim_types.push_back(NC_Data_Dims_Type::SN_WE); // 0
278 NC_fabs.push_back(&NC_pn_fab) ; NC_names.push_back("pn") ; NC_dim_types.push_back(NC_Data_Dims_Type::SN_WE); // 1
279 // Read the netcdf file and fill these FABs
280 BuildFABsFromNetCDFFile<FArrayBox,Real>(domain, fname, NC_names, NC_dim_types, NC_fabs, false, 0, ngrow);
281}
282
283/**
284 * @param lev level of data to read
285 * @param domain simulation domain
286 * @param fname file name to read from
287 * @param NC_fcor_fab container for Coriolis parameter data
288 */
289void
290read_coriolis_from_netcdf (int /*lev*/,
291 const Box& domain,
292 const std::string& fname,
293 FArrayBox& NC_fcor_fab)
294{
295 amrex::Print() << "Loading initial coriolis from NetCDF file " << fname << std::endl;
296
297 Vector<FArrayBox*> NC_fabs;
298 Vector<std::string> NC_names;
299 Vector<enum NC_Data_Dims_Type> NC_dim_types;
300
301 NC_fabs.push_back(&NC_fcor_fab ) ; NC_names.push_back("f") ; NC_dim_types.push_back(NC_Data_Dims_Type::SN_WE); // 0
302
303 // Read the netcdf file and fill these FABs
304 BuildFABsFromNetCDFFile<FArrayBox,Real>(domain, fname, NC_names, NC_dim_types, NC_fabs);
305}
306
307/**
308 * @param lev level of data to read
309 * @param domain simulation domain
310 * @param fname file name to read from
311 * @param NC_mskr_fab container for rho-point land/sea mask data
312 * @param NC_msku_fab container for u-point land/sea mask data
313 * @param NC_mskv_fab container for v-point land/sea mask data
314 */
315void
316read_masks_from_netcdf (int /*lev*/,
317 const Box& domain,
318 const std::string& fname,
319 FArrayBox& NC_mskr_fab,
320 FArrayBox& NC_msku_fab,
321 FArrayBox& NC_mskv_fab)
322{
323 amrex::Print() << "Loading masks from NetCDF file " << fname << std::endl;
324
325 Vector<FArrayBox*> NC_fabs;
326 Vector<std::string> NC_names;
327 Vector<enum NC_Data_Dims_Type> NC_dim_types;
328
329 NC_fabs.push_back(&NC_mskr_fab ) ; NC_names.push_back("mask_rho") ; NC_dim_types.push_back(NC_Data_Dims_Type::SN_WE); // 0
330 NC_fabs.push_back(&NC_msku_fab ) ; NC_names.push_back("mask_u") ; NC_dim_types.push_back(NC_Data_Dims_Type::SN_WE); // 1
331 NC_fabs.push_back(&NC_mskv_fab ) ; NC_names.push_back("mask_v") ; NC_dim_types.push_back(NC_Data_Dims_Type::SN_WE); // 2
332
333 // Read the netcdf file and fill these FABs
334 BuildFABsFromNetCDFFile<FArrayBox,Real>(domain, fname, NC_names, NC_dim_types, NC_fabs);
335}
336
337/**
338 * @param lev level of data to read
339 * @param domain simulation domain
340 * @param fname file name to read from
341 * @param do_m2_clim_nudg whether to do 2d momentum climatology nudging
342 * @param do_m3_clim_nudg whether to do 3d momentum climatology nudging
343 * @param do_temp_clim_nudg whether to do temperature climatology nudging
344 * @param do_salt_clim_nudg whether to do salinity climatology nudging
345 * @param NC_M2NC_fab container for 2d momentum climatology data
346 * @param NC_M3NC_fab container for 3d momentum climatology data
347 * @param NC_TempNC_fab container for temperature climatology data
348 * @param NC_SaltNC_fab container for salinity climatology data
349 */
350void
351read_clim_nudg_coeff_from_netcdf (int /*lev*/,
352 const Box& domain,
353 const std::string& fname,
354 bool do_m2_clim_nudg,
355 bool do_m3_clim_nudg,
356 bool do_temp_clim_nudg,
357 bool do_salt_clim_nudg,
358 FArrayBox& NC_M2NC_fab,
359 FArrayBox& NC_M3NC_fab,
360 FArrayBox& NC_TempNC_fab,
361 FArrayBox& NC_SaltNC_fab)
362{
363 amrex::Print() << "Loading nudging coefficients from NetCDF file " << fname << std::endl;
364
365 Vector<FArrayBox*> NC_fabs;
366 Vector<std::string> NC_names;
367 Vector<enum NC_Data_Dims_Type> NC_dim_types;
368
369 if (do_m3_clim_nudg) {
370 NC_fabs.push_back(&NC_M3NC_fab ); NC_names.push_back("M3_NudgeCoef"); NC_dim_types.push_back(NC_Data_Dims_Type::BT_SN_WE);
371 }
372 if (do_m2_clim_nudg) {
373 NC_fabs.push_back(&NC_M2NC_fab ); NC_names.push_back("M2_NudgeCoef"); NC_dim_types.push_back(NC_Data_Dims_Type::SN_WE);
374 }
375 if (do_temp_clim_nudg) {
376 NC_fabs.push_back(&NC_TempNC_fab ); NC_names.push_back("temp_NudgeCoef"); NC_dim_types.push_back(NC_Data_Dims_Type::BT_SN_WE);
377 }
378 if (do_salt_clim_nudg) {
379 NC_fabs.push_back(&NC_SaltNC_fab ); NC_names.push_back("salt_NudgeCoef"); NC_dim_types.push_back(NC_Data_Dims_Type::BT_SN_WE);
380 }
381 // Read the netcdf file and fill these FABs
382 BuildFABsFromNetCDFFile<FArrayBox,Real>(domain, fname, NC_names, NC_dim_types, NC_fabs);
383}
384
385/**
386 * @param lev level of data to read
387 * @param fnames file name(s) to read from
388 * @param field_name field name to read
389 * @param vec_dat vector to fill data
390 */
391//template <typename DType>
392void read_vec_from_netcdf (int /*lev*/, const amrex::Vector<std::string>& fnames, const std::string& field_name, amrex::Vector<int>& vec_dat)
393{
394 AMREX_ALWAYS_ASSERT(!fnames.empty());
395 const std::string& fname = fnames[0];
396
397 amrex::Print() << "Reading " << field_name << " from NetCDF file" << std::endl;
398
399 // get x-positions and put in array
400 using ARRAY = NDArray<int>;
401 amrex::Vector<ARRAY> array_dat(1);
402 ReadNetCDFFile(fname, {field_name}, array_dat); // filled only on proc 0
403 if (amrex::ParallelDescriptor::IOProcessor())
404 {
405 int n = array_dat[0].get_vshape()[0];
406 for (int i(0); i < n; i++)
407 {
408 vec_dat.push_back((*(array_dat[0].get_data() + i)));
409 }
410 }
411 int nvals = vec_dat.size();
412 int ioproc = amrex::ParallelDescriptor::IOProcessorNumber();
413 amrex::ParallelDescriptor::Bcast(&nvals,1,ioproc);
414 if (!(amrex::ParallelDescriptor::IOProcessor())) {
415 vec_dat.resize(nvals);
416 }
417 amrex::ParallelDescriptor::Bcast(vec_dat.data(), vec_dat.size(), ioproc);
418}
419
420#endif // ROMSX_USE_NETCDF
REMORA_DataStruct.H
AMREX_ALWAYS_ASSERT
AMREX_ALWAYS_ASSERT(!NSPeriodic||!EWPeriodic)
REMORA_NCFile.H
NC_Data_Dims_Type::BT_SN_WE
@ BT_SN_WE
NC_Data_Dims_Type::Time_SN_WE
@ Time_SN_WE
NC_Data_Dims_Type::SN_WE
@ SN_WE
NC_Data_Dims_Type::Time_BT_SN_WE
@ Time_BT_SN_WE
ReadNetCDFFile
void ReadNetCDFFile(const std::string &fname, amrex::Vector< std::string > names, amrex::Vector< NDArray< DType > > &arrays, bool one_time=false, int fill_time=0)
Read in data from netcdf file and save to data arrays.
Definition REMORA_NCFile.H:91
QueryNetCDFHasVars
bool QueryNetCDFHasVars(const std::string &fname, const amrex::Vector< std::string > &var_names)
Helper function for testing whether a file carries every named variable.
Definition REMORA_NCFile.cpp:65
read_clim_nudg_coeff_from_netcdf
void read_clim_nudg_coeff_from_netcdf(int, const Box &domain, const std::string &fname, bool do_m2_clim_nudg, bool do_m3_clim_nudg, bool do_temp_clim_nudg, bool do_salt_clim_nudg, FArrayBox &NC_M2NC_fab, FArrayBox &NC_M3NC_fab, FArrayBox &NC_TempNC_fab, FArrayBox &NC_SaltNC_fab)
helper function to read climatology nudging from netcdf
Definition REMORA_ReadFromInitNetcdf.cpp:351
read_zeta_full_domain_from_netcdf
void read_zeta_full_domain_from_netcdf(int, const Box &domain, const std::string &fname, FArrayBox &NC_zeta_fab, IntVect ngrow)
helper function to read high-resolution full-domain sea surface height from netcdf
Definition REMORA_ReadFromInitNetcdf.cpp:104
read_bathymetry_from_netcdf
void read_bathymetry_from_netcdf(int, const Box &domain, const std::string &fname, FArrayBox &NC_h_fab)
helper function to read bathymetry from netcdf
Definition REMORA_ReadFromInitNetcdf.cpp:129
read_vec_from_netcdf
void read_vec_from_netcdf(int, const amrex::Vector< std::string > &fnames, const std::string &field_name, amrex::Vector< int > &vec_dat)
helper function to read in vector of data from netcdf
Definition REMORA_ReadFromInitNetcdf.cpp:392
read_data_full_domain_from_netcdf
void read_data_full_domain_from_netcdf(int, const Box &domain, const std::string &fname, FArrayBox &NC_temp_fab, FArrayBox &NC_salt_fab, FArrayBox &NC_xvel_fab, FArrayBox &NC_yvel_fab, IntVect ngrow)
helper function for reading in full domain high-resolution initial state data from netcdf
Definition REMORA_ReadFromInitNetcdf.cpp:50
read_grid_vars_from_netcdf
void read_grid_vars_from_netcdf(int, const Box &domain, const std::string &fname, FArrayBox &NC_pm_fab, FArrayBox &NC_pn_fab, FArrayBox &NC_xr_fab, FArrayBox &NC_yr_fab, FArrayBox &NC_xu_fab, FArrayBox &NC_yu_fab, FArrayBox &NC_xv_fab, FArrayBox &NC_yv_fab, FArrayBox &NC_xp_fab, FArrayBox &NC_yp_fab)
helper function to read grid variables from netcdf
Definition REMORA_ReadFromInitNetcdf.cpp:162
read_bathymetry_full_domain_from_netcdf
void read_bathymetry_full_domain_from_netcdf(const Box &domain, const std::string &fname, FArrayBox &NC_h_fab, IntVect ngrow)
helper function to read full-domain high resolution bathymetry from netcdf
Definition REMORA_ReadFromInitNetcdf.cpp:241
read_coriolis_from_netcdf
void read_coriolis_from_netcdf(int, const Box &domain, const std::string &fname, FArrayBox &NC_fcor_fab)
helper function to read coriolis factor from netcdf
Definition REMORA_ReadFromInitNetcdf.cpp:290
read_masks_from_netcdf
void read_masks_from_netcdf(int, const Box &domain, const std::string &fname, FArrayBox &NC_mskr_fab, FArrayBox &NC_msku_fab, FArrayBox &NC_mskv_fab)
helper function for reading in land-sea masks from netcdf
Definition REMORA_ReadFromInitNetcdf.cpp:316
read_zeta_from_netcdf
void read_zeta_from_netcdf(int, const Box &domain, const std::string &fname, FArrayBox &NC_zeta_fab)
helper function to read sea surface height from netcdf
Definition REMORA_ReadFromInitNetcdf.cpp:79
read_spherical_grid_vars_from_netcdf
bool read_spherical_grid_vars_from_netcdf(int, const Box &domain, const std::string &fname, FArrayBox &NC_lonp_fab, FArrayBox &NC_latp_fab)
helper function to read optional spherical psi coordinates from netcdf
Definition REMORA_ReadFromInitNetcdf.cpp:207
read_data_from_netcdf
void read_data_from_netcdf(int, const Box &domain, const std::string &fname, FArrayBox &NC_temp_fab, FArrayBox &NC_salt_fab, FArrayBox &NC_xvel_fab, FArrayBox &NC_yvel_fab)
helper function for reading in initial state data from netcdf
Definition REMORA_ReadFromInitNetcdf.cpp:18
read_grid_vars_full_domain_from_netcdf
void read_grid_vars_full_domain_from_netcdf(const Box &domain, const std::string &fname, FArrayBox &NC_pm_fab, FArrayBox &NC_pn_fab, IntVect ngrow)
helper function to read full-domain high resolution grid variables from netcdf
Definition REMORA_ReadFromInitNetcdf.cpp:266
NDArray
NDArray is the datatype designed to hold any data, including scalars, multidimensional arrays,...
Definition REMORA_NCFile.H:36