MODFLOW 6  version 6.6.0.dev0
USGS Modular Hydrologic Model
Data Types | Functions/Subroutines
meshmodelmodule Module Reference

This module contains the MeshModelModule. More...

Data Types

type  meshncdimidtype
 type for storing model export dimension ids More...
 
type  meshncvaridtype
 type for storing model export variable ids More...
 
type  meshmodeltype
 base ugrid netcdf export type More...
 
interface  nc_array_export_if
 abstract interfaces for derived ugrid netcd export types More...
 
type  mesh2dmodeltype
 

Functions/Subroutines

subroutine mesh_init (this, modelname, modeltype, modelfname, disenum, nctype, iout)
 initialize More...
 
subroutine mesh_destroy (this)
 initialize More...
 
subroutine add_global_att (this)
 create file (group) attributes More...
 
subroutine export_input_arrays (this, pkgtype, pkgname, mempath, param_dfns)
 write package gridded input data More...
 
subroutine add_pkg_data (this)
 determine packages to write gridded input More...
 
subroutine define_dependent (this)
 create the model layer dependent variables More...
 
subroutine define_gridmap (this)
 create the file grid mapping container variable More...
 
subroutine create_mesh (this)
 create the file mesh container variable More...
 
subroutine, public ncvar_chunk (ncid, varid, chunk_face, nc_fname)
 define variable chunking More...
 
subroutine, public ncvar_deflate (ncid, varid, deflate, shuffle, nc_fname)
 define variable compression More...
 
subroutine, public ncvar_gridmap (ncid, varid, gridmap_name, nc_fname)
 put variable gridmap attributes More...
 
subroutine, public ncvar_mf6attr (ncid, varid, layer, iper, iaux, nc_tag, nc_fname)
 put variable internal attributes More...
 
character(len=linelength) function, public export_varname (varname, layer, iper, iaux)
 build netcdf variable name More...
 

Detailed Description

This module defines a base class for UGRID based (mesh) model netcdf exports. It is dependent on external netcdf libraries.

Function/Subroutine Documentation

◆ add_global_att()

subroutine meshmodelmodule::add_global_att ( class(meshmodeltype), intent(inout)  this)

Definition at line 147 of file MeshNCModel.f90.

148  class(MeshModelType), intent(inout) :: this
149  ! -- file scoped title
150  call nf_verify(nf90_put_att(this%ncid, nf90_global, 'title', &
151  this%annotation%title), this%nc_fname)
152  ! -- source (MODFLOW 6)
153  call nf_verify(nf90_put_att(this%ncid, nf90_global, 'source', &
154  this%annotation%source), this%nc_fname)
155  ! -- export type (MODFLOW 6)
156  call nf_verify(nf90_put_att(this%ncid, nf90_global, 'modflow6_grid', &
157  this%annotation%grid), this%nc_fname)
158  ! -- MODFLOW 6 model type
159  call nf_verify(nf90_put_att(this%ncid, nf90_global, 'modflow6_model', &
160  this%annotation%model), this%nc_fname)
161  ! -- generation datetime
162  call nf_verify(nf90_put_att(this%ncid, nf90_global, 'history', &
163  this%annotation%history), this%nc_fname)
164  ! -- supported conventions
165  call nf_verify(nf90_put_att(this%ncid, nf90_global, 'Conventions', &
166  this%annotation%conventions), &
167  this%nc_fname)
Here is the call graph for this function:

◆ add_pkg_data()

subroutine meshmodelmodule::add_pkg_data ( class(meshmodeltype), intent(inout)  this)

Definition at line 201 of file MeshNCModel.f90.

202  use simvariablesmodule, only: idm_context
203  use memorymanagerextmodule, only: mem_set_value
204  use memoryhelpermodule, only: create_mem_path
205  use idmdfnselectormodule, only: param_definitions
206  use sourcecommonmodule, only: idm_subcomponent_type
207  use idmdfnselectormodule, only: idm_multi_package
208  class(MeshModelType), intent(inout) :: this
209  character(LENCOMPONENTNAME) :: ptype, pname, pkgtype
210  character(len=LENMEMPATH) :: input_mempath
211  type(CharacterStringType), dimension(:), contiguous, &
212  pointer :: pkgtypes => null()
213  type(CharacterStringType), dimension(:), contiguous, &
214  pointer :: pkgnames => null()
215  type(CharacterStringType), dimension(:), contiguous, &
216  pointer :: mempaths => null()
217  type(InputParamDefinitionType), dimension(:), pointer :: param_dfns
218  character(len=LENMEMPATH) :: mempath
219  integer(I4B) :: n
220  integer(I4B), pointer :: export_arrays
221  logical(LGP) :: found
222  !
223  input_mempath = create_mem_path(component=this%modelname, context=idm_context)
224  !
225  ! -- set pointers to model path package info
226  call mem_setptr(pkgtypes, 'PKGTYPES', input_mempath)
227  call mem_setptr(pkgnames, 'PKGNAMES', input_mempath)
228  call mem_setptr(mempaths, 'MEMPATHS', input_mempath)
229  !
230  allocate (export_arrays)
231  !
232  do n = 1, size(mempaths)
233  !
234  ! -- initialize export_arrays
235  export_arrays = 0
236  !
237  ! -- set package attributes
238  mempath = mempaths(n)
239  pname = pkgnames(n)
240  ptype = pkgtypes(n)
241  !
242  ! -- export input arrays
243  if (mempath /= '') then
244  ! -- update export
245  call mem_set_value(export_arrays, 'EXPORT_NC', mempath, found)
246  !
247  if (export_arrays > 0) then
248  pkgtype = idm_subcomponent_type(this%modeltype, ptype)
249  param_dfns => param_definitions(this%modeltype, pkgtype)
250  call this%export_input_arrays(ptype, pname, mempath, param_dfns)
251  end if
252  end if
253  end do
254  !
255  ! -- cleanup
256  deallocate (export_arrays)
memorymanagerextmodule::mem_set_value
Definition: MemoryManagerExt.f90:14
idmdfnselectormodule::param_definitions
type(inputparamdefinitiontype) function, dimension(:), pointer, public param_definitions(component, subcomponent)
Definition: IdmDfnSelector.f90:32
idmdfnselectormodule::idm_multi_package
logical function, public idm_multi_package(component, subcomponent)
Definition: IdmDfnSelector.f90:125
memoryhelpermodule::create_mem_path
character(len=lenmempath) function create_mem_path(component, subcomponent, context)
returns the path to the memory object
Definition: MemoryHelper.f90:22
simvariablesmodule
This module contains simulation variables.
Definition: SimVariables.f90:9
simvariablesmodule::idm_context
character(len=linelength) idm_context
Definition: SimVariables.f90:18
sourcecommonmodule
This module contains the SourceCommonModule.
Definition: SourceCommon.f90:7
sourcecommonmodule::idm_subcomponent_type
character(len=lencomponentname) function, public idm_subcomponent_type(component, subcomponent)
component from package or model type
Definition: SourceCommon.f90:101
Here is the call graph for this function:

◆ create_mesh()

subroutine meshmodelmodule::create_mesh ( class(mesh2dmodeltype), intent(inout)  this)
private

Definition at line 347 of file MeshNCModel.f90.

348  class(Mesh2dModelType), intent(inout) :: this
349  !
350  ! -- create mesh container variable
351  call nf_verify(nf90_def_var(this%ncid, this%mesh_name, nf90_int, &
352  this%var_ids%mesh), this%nc_fname)
353  !
354  ! -- assign container variable attributes
355  call nf_verify(nf90_put_att(this%ncid, this%var_ids%mesh, 'cf_role', &
356  'mesh_topology'), this%nc_fname)
357  call nf_verify(nf90_put_att(this%ncid, this%var_ids%mesh, 'long_name', &
358  '2D mesh topology'), this%nc_fname)
359  call nf_verify(nf90_put_att(this%ncid, this%var_ids%mesh, &
360  'topology_dimension', 2), this%nc_fname)
361  call nf_verify(nf90_put_att(this%ncid, this%var_ids%mesh, 'face_dimension', &
362  'nmesh_face'), this%nc_fname)
363  call nf_verify(nf90_put_att(this%ncid, this%var_ids%mesh, &
364  'node_coordinates', 'mesh_node_x mesh_node_y'), &
365  this%nc_fname)
366  call nf_verify(nf90_put_att(this%ncid, this%var_ids%mesh, &
367  'face_coordinates', 'mesh_face_x mesh_face_y'), &
368  this%nc_fname)
369  call nf_verify(nf90_put_att(this%ncid, this%var_ids%mesh, &
370  'face_node_connectivity', 'mesh_face_nodes'), &
371  this%nc_fname)
372 
373  ! -- create mesh x node (mesh vertex) variable
374  call nf_verify(nf90_def_var(this%ncid, 'mesh_node_x', nf90_double, &
375  (/this%dim_ids%nmesh_node/), &
376  this%var_ids%mesh_node_x), this%nc_fname)
377  !
378  ! -- assign mesh x node variable attributes
379  call nf_verify(nf90_put_att(this%ncid, this%var_ids%mesh_node_x, &
380  'units', 'm'), this%nc_fname)
381  call nf_verify(nf90_put_att(this%ncid, this%var_ids%mesh_node_x, &
382  'standard_name', 'projection_x_coordinate'), &
383  this%nc_fname)
384  call nf_verify(nf90_put_att(this%ncid, this%var_ids%mesh_node_x, &
385  'long_name', 'Easting'), this%nc_fname)
386  !
387  if (this%ogc_wkt /= '') then
388  ! -- associate with projection
389  call nf_verify(nf90_put_att(this%ncid, this%var_ids%mesh_node_x, &
390  'grid_mapping', this%gridmap_name), &
391  this%nc_fname)
392  end if
393 
394  ! -- create mesh y node (mesh vertex) variable
395  call nf_verify(nf90_def_var(this%ncid, 'mesh_node_y', nf90_double, &
396  (/this%dim_ids%nmesh_node/), &
397  this%var_ids%mesh_node_y), this%nc_fname)
398  !
399  ! -- assign mesh y variable attributes
400  call nf_verify(nf90_put_att(this%ncid, this%var_ids%mesh_node_y, &
401  'units', 'm'), this%nc_fname)
402  call nf_verify(nf90_put_att(this%ncid, this%var_ids%mesh_node_y, &
403  'standard_name', 'projection_y_coordinate'), &
404  this%nc_fname)
405  call nf_verify(nf90_put_att(this%ncid, this%var_ids%mesh_node_y, &
406  'long_name', 'Northing'), this%nc_fname)
407  !
408  if (this%ogc_wkt /= '') then
409  ! -- associate with projection
410  call nf_verify(nf90_put_att(this%ncid, this%var_ids%mesh_node_y, &
411  'grid_mapping', this%gridmap_name), &
412  this%nc_fname)
413  end if
414 
415  ! -- create mesh x face (cell vertex) variable
416  call nf_verify(nf90_def_var(this%ncid, 'mesh_face_x', nf90_double, &
417  (/this%dim_ids%nmesh_face/), &
418  this%var_ids%mesh_face_x), this%nc_fname)
419  !
420  ! -- assign mesh x face variable attributes
421  call nf_verify(nf90_put_att(this%ncid, this%var_ids%mesh_face_x, &
422  'units', 'm'), this%nc_fname)
423  call nf_verify(nf90_put_att(this%ncid, this%var_ids%mesh_face_x, &
424  'standard_name', 'projection_x_coordinate'), &
425  this%nc_fname)
426  call nf_verify(nf90_put_att(this%ncid, this%var_ids%mesh_face_x, &
427  'long_name', 'Easting'), this%nc_fname)
428  call nf_verify(nf90_put_att(this%ncid, this%var_ids%mesh_face_x, 'bounds', &
429  'mesh_face_xbnds'), this%nc_fname)
430  if (this%ogc_wkt /= '') then
431  ! -- associate with projection
432  call nf_verify(nf90_put_att(this%ncid, this%var_ids%mesh_face_x, &
433  'grid_mapping', this%gridmap_name), &
434  this%nc_fname)
435  end if
436 
437  ! -- create mesh x cell bounds variable
438  call nf_verify(nf90_def_var(this%ncid, 'mesh_face_xbnds', nf90_double, &
439  (/this%dim_ids%max_nmesh_face_nodes, &
440  this%dim_ids%nmesh_face/), &
441  this%var_ids%mesh_face_xbnds), &
442  this%nc_fname)
443  !
444  ! -- create mesh y face (cell vertex) variable
445  call nf_verify(nf90_def_var(this%ncid, 'mesh_face_y', nf90_double, &
446  (/this%dim_ids%nmesh_face/), &
447  this%var_ids%mesh_face_y), this%nc_fname)
448  !
449  ! -- assign mesh y face variable attributes
450  call nf_verify(nf90_put_att(this%ncid, this%var_ids%mesh_face_y, &
451  'units', 'm'), this%nc_fname)
452  call nf_verify(nf90_put_att(this%ncid, this%var_ids%mesh_face_y, &
453  'standard_name', 'projection_y_coordinate'), &
454  this%nc_fname)
455  call nf_verify(nf90_put_att(this%ncid, this%var_ids%mesh_face_y, &
456  'long_name', 'Northing'), this%nc_fname)
457  call nf_verify(nf90_put_att(this%ncid, this%var_ids%mesh_face_y, 'bounds', &
458  'mesh_face_ybnds'), this%nc_fname)
459  !
460  if (this%ogc_wkt /= '') then
461  ! -- associate with projection
462  call nf_verify(nf90_put_att(this%ncid, this%var_ids%mesh_face_y, &
463  'grid_mapping', this%gridmap_name), &
464  this%nc_fname)
465  end if
466 
467  ! -- create mesh y cell bounds variable
468  call nf_verify(nf90_def_var(this%ncid, 'mesh_face_ybnds', nf90_double, &
469  (/this%dim_ids%max_nmesh_face_nodes, &
470  this%dim_ids%nmesh_face/), &
471  this%var_ids%mesh_face_ybnds), &
472  this%nc_fname)
473  !
474  ! -- create mesh face nodes variable
475  call nf_verify(nf90_def_var(this%ncid, 'mesh_face_nodes', nf90_int, &
476  (/this%dim_ids%max_nmesh_face_nodes, &
477  this%dim_ids%nmesh_face/), &
478  this%var_ids%mesh_face_nodes), &
479  this%nc_fname)
480  !
481  ! -- assign variable attributes
482  call nf_verify(nf90_put_att(this%ncid, this%var_ids%mesh_face_nodes, &
483  'cf_role', 'face_node_connectivity'), &
484  this%nc_fname)
485  call nf_verify(nf90_put_att(this%ncid, this%var_ids%mesh_face_nodes, &
486  'long_name', &
487  'Vertices bounding cell (counterclockwise)'), &
488  this%nc_fname)
489  call nf_verify(nf90_put_att(this%ncid, this%var_ids%mesh_face_nodes, &
490  '_FillValue', (/nf90_fill_int/)), &
491  this%nc_fname)
492  call nf_verify(nf90_put_att(this%ncid, this%var_ids%mesh_face_nodes, &
493  'start_index', 1), this%nc_fname)
Here is the call graph for this function:

◆ define_dependent()

subroutine meshmodelmodule::define_dependent ( class(meshmodeltype), intent(inout)  this)

Definition at line 261 of file MeshNCModel.f90.

262  class(MeshModelType), intent(inout) :: this
263  character(len=LINELENGTH) :: varname, longname
264  integer(I4B) :: k
265  !
266  ! -- create a dependent variable for each layer
267  do k = 1, this%nlay
268  !
269  ! -- initialize names
270  varname = ''
271  longname = ''
272  !
273  ! -- set layer variable and longnames
274  write (varname, '(a,i0)') trim(this%xname)//'_l', k
275  write (longname, '(a,i0,a)') trim(this%annotation%longname)// &
276  ' (layer ', k, ')'
277  !
278  ! -- create the netcdf dependent layer variable
279  call nf_verify(nf90_def_var(this%ncid, varname, nf90_double, &
280  (/this%dim_ids%nmesh_face, &
281  this%dim_ids%time/), &
282  this%var_ids%dependent(k)), &
283  this%nc_fname)
284  !
285  ! -- apply chunking parameters
286  if (this%chunking_active) then
287  call nf_verify(nf90_def_var_chunking(this%ncid, &
288  this%var_ids%dependent(k), &
289  nf90_chunked, &
290  (/this%chunk_face, &
291  this%chunk_time/)), &
292  this%nc_fname)
293  end if
294  ! -- deflate and shuffle
295  call ncvar_deflate(this%ncid, this%var_ids%dependent(k), this%deflate, &
296  this%shuffle, this%nc_fname)
297  !
298  ! -- assign variable attributes
299  call nf_verify(nf90_put_att(this%ncid, this%var_ids%dependent(k), &
300  'units', 'm'), this%nc_fname)
301  call nf_verify(nf90_put_att(this%ncid, this%var_ids%dependent(k), &
302  'standard_name', this%annotation%stdname), &
303  this%nc_fname)
304  call nf_verify(nf90_put_att(this%ncid, this%var_ids%dependent(k), &
305  'long_name', longname), this%nc_fname)
306  call nf_verify(nf90_put_att(this%ncid, this%var_ids%dependent(k), &
307  '_FillValue', (/dhnoflo/)), &
308  this%nc_fname)
309  call nf_verify(nf90_put_att(this%ncid, this%var_ids%dependent(k), &
310  'mesh', this%mesh_name), this%nc_fname)
311  call nf_verify(nf90_put_att(this%ncid, this%var_ids%dependent(k), &
312  'location', 'face'), this%nc_fname)
313  !
314  ! -- add grid mapping
315  call ncvar_gridmap(this%ncid, this%var_ids%dependent(k), &
316  this%gridmap_name, this%nc_fname)
317  end do
Here is the call graph for this function:

◆ define_gridmap()

subroutine meshmodelmodule::define_gridmap ( class(meshmodeltype), intent(inout)  this)
private

Definition at line 322 of file MeshNCModel.f90.

323  class(MeshModelType), intent(inout) :: this
324  integer(I4B) :: var_id
325  !
326  ! -- was projection info provided
327  if (this%ogc_wkt /= '') then
328  !
329  ! -- create projection variable
330  call nf_verify(nf90_redef(this%ncid), this%nc_fname)
331  call nf_verify(nf90_def_var(this%ncid, this%gridmap_name, nf90_int, &
332  var_id), this%nc_fname)
333  ! -- cf-conventions prefers 'crs_wkt'
334  !call nf_verify(nf90_put_att(this%ncid, var_id, 'crs_wkt', this%ogc_wkt), &
335  ! this%nc_fname)
336  ! -- QGIS recognizes 'wkt'
337  call nf_verify(nf90_put_att(this%ncid, var_id, 'wkt', this%ogc_wkt), &
338  this%nc_fname)
339  call nf_verify(nf90_enddef(this%ncid), this%nc_fname)
340  call nf_verify(nf90_put_var(this%ncid, var_id, 1), &
341  this%nc_fname)
342  end if
Here is the call graph for this function:

◆ export_input_arrays()

subroutine meshmodelmodule::export_input_arrays ( class(meshmodeltype), intent(inout)  this,
character(len=*), intent(in)  pkgtype,
character(len=*), intent(in)  pkgname,
character(len=*), intent(in)  mempath,
type(inputparamdefinitiontype), dimension(:), intent(in), pointer  param_dfns 
)
private

Definition at line 172 of file MeshNCModel.f90.

173  use memorymanagermodule, only: get_isize
174  class(MeshModelType), intent(inout) :: this
175  character(len=*), intent(in) :: pkgtype
176  character(len=*), intent(in) :: pkgname
177  character(len=*), intent(in) :: mempath
178  type(InputParamDefinitionType), dimension(:), pointer, &
179  intent(in) :: param_dfns
180  type(InputParamDefinitionType), pointer :: idt
181  integer(I4B) :: iparam, isize
182  !
183  ! -- export griddata block parameters
184  do iparam = 1, size(param_dfns)
185  ! -- assign param definition pointer
186  idt => param_dfns(iparam)
187  ! -- for now
188  if (idt%blockname == 'GRIDDATA') then
189  ! -- veriy variable is allocated
190  call get_isize(idt%mf6varname, mempath, isize)
191  !
192  if (isize > 0) then
193  call this%export_input_array(pkgtype, pkgname, mempath, idt)
194  end if
195  end if
196  end do
memorymanagermodule::get_isize
subroutine, public get_isize(name, mem_path, isize)
@ brief Get the number of elements for this variable
Definition: MemoryManager.f90:265
Here is the call graph for this function:

◆ export_varname()

character(len=linelength) function, public meshmodelmodule::export_varname ( character(len=*), intent(in)  varname,
integer(i4b), intent(in), optional  layer,
integer(i4b), intent(in), optional  iper,
integer(i4b), intent(in), optional  iaux 
)

Definition at line 576 of file MeshNCModel.f90.

577  use inputoutputmodule, only: lowcase
578  character(len=*), intent(in) :: varname
579  integer(I4B), optional, intent(in) :: layer
580  integer(I4B), optional, intent(in) :: iper
581  integer(I4B), optional, intent(in) :: iaux
582  character(len=LINELENGTH) :: vname
583  !
584  vname = ''
585  !
586  if (varname /= '') then
587  vname = varname
588  call lowcase(vname)
589  if (present(layer)) then
590  if (layer > 0) then
591  write (vname, '(a,i0)') trim(vname)//'_l', layer
592  end if
593  end if
594  if (present(iper)) then
595  if (iper > 0) then
596  write (vname, '(a,i0)') trim(vname)//'_p', iper
597  end if
598  end if
599  if (present(iaux)) then
600  if (iaux > 0) then
601  write (vname, '(a,i0)') trim(vname)//'a', iaux
602  end if
603  end if
604  end if
inputoutputmodule::lowcase
subroutine, public lowcase(word)
Convert to lower case.
Definition: InputOutput.f90:211
Here is the call graph for this function:
Here is the caller graph for this function:

◆ mesh_destroy()

subroutine meshmodelmodule::mesh_destroy ( class(meshmodeltype), intent(inout)  this)

Definition at line 135 of file MeshNCModel.f90.

136  use memorymanagerextmodule, only: mem_set_value
137  class(MeshModelType), intent(inout) :: this
138  !
139  call nf_verify(nf90_close(this%ncid), this%nc_fname)
140  !
141  deallocate (this%chunk_face)
142  nullify (this%chunk_face)
Here is the call graph for this function:

◆ mesh_init()

subroutine meshmodelmodule::mesh_init ( class(meshmodeltype), intent(inout)  this,
character(len=*), intent(in)  modelname,
character(len=*), intent(in)  modeltype,
character(len=*), intent(in)  modelfname,
integer(i4b), intent(in)  disenum,
integer(i4b), intent(in)  nctype,
integer(i4b), intent(in)  iout 
)
private

Definition at line 98 of file MeshNCModel.f90.

100  use memorymanagerextmodule, only: mem_set_value
101  class(MeshModelType), intent(inout) :: this
102  character(len=*), intent(in) :: modelname
103  character(len=*), intent(in) :: modeltype
104  character(len=*), intent(in) :: modelfname
105  integer(I4B), intent(in) :: disenum
106  integer(I4B), intent(in) :: nctype
107  integer(I4B), intent(in) :: iout
108  logical(LGP) :: found
109  !
110  ! -- initialize base class
111  call this%NCModelExportType%init(modelname, modeltype, modelfname, disenum, &
112  nctype, iout)
113  !
114  ! -- allocate and initialize
115  allocate (this%chunk_face)
116  this%chunk_face = -1
117  !
118  ! -- update values from input context
119  if (this%ncf_mempath /= '') then
120  call mem_set_value(this%chunk_face, 'CHUNK_FACE', this%ncf_mempath, found)
121  end if
122  !
123  if (this%chunk_time > 0 .and. this%chunk_face > 0) then
124  this%chunking_active = .true.
125  end if
126  !
127  ! -- create the netcdf file
128  call nf_verify(nf90_create(this%nc_fname, &
129  iand(nf90_clobber, nf90_netcdf4), this%ncid), &
130  this%nc_fname)
Here is the call graph for this function:

◆ ncvar_chunk()

subroutine, public meshmodelmodule::ncvar_chunk ( integer(i4b), intent(in)  ncid,
integer(i4b), intent(in)  varid,
integer(i4b), intent(in)  chunk_face,
character(len=*), intent(in)  nc_fname 
)

Definition at line 498 of file MeshNCModel.f90.

499  integer(I4B), intent(in) :: ncid
500  integer(I4B), intent(in) :: varid
501  integer(I4B), intent(in) :: chunk_face
502  character(len=*), intent(in) :: nc_fname
503  !
504  ! -- apply chunking parameters
505  if (chunk_face > 0) then
506  call nf_verify(nf90_def_var_chunking(ncid, varid, nf90_chunked, &
507  (/chunk_face/)), nc_fname)
508  end if
Here is the call graph for this function:
Here is the caller graph for this function:

◆ ncvar_deflate()

subroutine, public meshmodelmodule::ncvar_deflate ( integer(i4b), intent(in)  ncid,
integer(i4b), intent(in)  varid,
integer(i4b), intent(in)  deflate,
integer(i4b), intent(in)  shuffle,
character(len=*), intent(in)  nc_fname 
)

Definition at line 513 of file MeshNCModel.f90.

514  integer(I4B), intent(in) :: ncid
515  integer(I4B), intent(in) :: varid
516  integer(I4B), intent(in) :: deflate
517  integer(I4B), intent(in) :: shuffle
518  character(len=*), intent(in) :: nc_fname
519  ! -- deflate and shuffle
520  if (deflate >= 0) then
521  call nf_verify(nf90_def_var_deflate(ncid, varid, shuffle=shuffle, &
522  deflate=1, deflate_level=deflate), &
523  nc_fname)
524  end if
Here is the call graph for this function:
Here is the caller graph for this function:

◆ ncvar_gridmap()

subroutine, public meshmodelmodule::ncvar_gridmap ( integer(i4b), intent(in)  ncid,
integer(i4b), intent(in)  varid,
character(len=*), intent(in)  gridmap_name,
character(len=*), intent(in)  nc_fname 
)

Definition at line 529 of file MeshNCModel.f90.

530  integer(I4B), intent(in) :: ncid
531  integer(I4B), intent(in) :: varid
532  character(len=*), intent(in) :: gridmap_name
533  character(len=*), intent(in) :: nc_fname
534  !
535  if (gridmap_name /= '') then
536  call nf_verify(nf90_put_att(ncid, varid, 'coordinates', &
537  'mesh_face_x mesh_face_y'), nc_fname)
538  call nf_verify(nf90_put_att(ncid, varid, 'grid_mapping', &
539  gridmap_name), nc_fname)
540  end if
Here is the call graph for this function:
Here is the caller graph for this function:

◆ ncvar_mf6attr()

subroutine, public meshmodelmodule::ncvar_mf6attr ( integer(i4b), intent(in)  ncid,
integer(i4b), intent(in)  varid,
integer(i4b), intent(in)  layer,
integer(i4b), intent(in)  iper,
integer(i4b), intent(in)  iaux,
character(len=*), intent(in)  nc_tag,
character(len=*), intent(in)  nc_fname 
)

Definition at line 545 of file MeshNCModel.f90.

546  integer(I4B), intent(in) :: ncid
547  integer(I4B), intent(in) :: varid
548  integer(I4B), intent(in) :: layer
549  integer(I4B), intent(in) :: iper
550  integer(I4B), intent(in) :: iaux
551  character(len=*), intent(in) :: nc_tag
552  character(len=*), intent(in) :: nc_fname
553  !
554  if (nc_tag /= '') then
555  call nf_verify(nf90_put_att(ncid, varid, 'modflow6_input', &
556  nc_tag), nc_fname)
557  if (layer > 0) then
558  call nf_verify(nf90_put_att(ncid, varid, 'modflow6_layer', &
559  layer), nc_fname)
560  end if
561  !
562  if (iper > 0) then
563  call nf_verify(nf90_put_att(ncid, varid, 'modflow6_iper', &
564  iper), nc_fname)
565  end if
566  !
567  if (iaux > 0) then
568  call nf_verify(nf90_put_att(ncid, varid, 'modflow6_iaux', &
569  iaux), nc_fname)
570  end if
571  end if
Here is the call graph for this function:
Here is the caller graph for this function: