add warning when wkt provided for rotated grid with structured export

Author: mjreno

doc/mf6io/mf6ivar/dfn/utl-ncf.dfn

@@ -96,7 +96,7 @@ shape (ncpl)
 optional true
 reader readarray
 longname cell center latitude
-description cell center latitude. Only supported for NETCDF\_STRUCTURED export types.
+description cell center latitude. Only supported for NETCDF\_STRUCTURED export type.
 
 block griddata
 name longitude
@@ -105,4 +105,4 @@ shape (ncpl)
 optional true
 reader readarray
 longname cell center longitude
-description cell center longitude. Only supported for NETCDF\_STRUCTURED export types.
+description cell center longitude. Only supported for NETCDF\_STRUCTURED export type.

src/Utilities/Export/DisNCStructured.f90

@@ -9,7 +9,7 @@ module DisNCStructuredModule
 
   use KindModule, only: DP, I4B, LGP
   use ConstantsModule, only: LINELENGTH, LENBIGLINE, LENCOMPONENTNAME, &
-                             LENMEMPATH, LENVARNAME
+                             LENMEMPATH, LENVARNAME, DNODATA, DZERO
   use SimVariablesModule, only: errmsg, warnmsg
   use SimModule, only: store_error, store_warning, store_error_filename
   use MemoryManagerModule, only: mem_setptr
@@ -127,6 +127,7 @@ subroutine dis_export_init(this, modelname, modeltype, modelfname, nc_fname, &
     ! initialize base class
     call this%NCModelExportType%init(modelname, modeltype, modelfname, nc_fname, &
                                      disenum, nctype, iout)
+
     ! update values from input context
     if (this%ncf_mempath /= '') then
       call mem_set_value(this%chunk_z, 'CHUNK_Z', this%ncf_mempath, found)
@@ -163,6 +164,14 @@ subroutine dis_export_init(this, modelname, modeltype, modelfname, nc_fname, &
         call mem_setptr(this%latitude, 'LATITUDE', this%ncf_mempath)
         call mem_setptr(this%longitude, 'LONGITUDE', this%ncf_mempath)
       end if
+
+      if (this%wkt /= '') then
+        if (this%dis%angrot /= DZERO) then
+          write (warnmsg, '(a)') 'WKT parameter set with structured rotated &
+            &grid. Projected coordinates will have grid local values.'
+          call store_warning(warnmsg)
+        end if
+      end if
     end if
 
     if (this%dis%lenuni == 1) then
@@ -380,7 +389,6 @@ end subroutine export_input_arrays
   !> @brief netcdf export package dynamic input with ilayer index variable
   !<
   subroutine package_step_ilayer(this, export_pkg, ilayer_varname, ilayer)
-    use ConstantsModule, only: DNODATA, DZERO
     use TdisModule, only: kper
     use NCModelExportModule, only: ExportPackageType
     use DefinitionSelectModule, only: get_param_definition_type
@@ -528,7 +536,6 @@ end subroutine package_step
   !<
   subroutine export_layer_3d(this, export_pkg, idt, ilayer_read, ialayer, &
                              dbl1d, nc_varname, input_attr, iaux)
-    use ConstantsModule, only: DNODATA, DZERO
     use TdisModule, only: kper
     use NCModelExportModule, only: ExportPackageType
     class(DisNCStructuredType), intent(inout) :: this
@@ -878,7 +885,6 @@ end subroutine add_proj_data
   !> @brief add grid coordinates
   !<
   subroutine add_grid_data(this)
-    use ConstantsModule, only: DZERO
     class(DisNCStructuredType), intent(inout) :: this
     integer(I4B) :: ibnd, n !, k, i, j
     real(DP), dimension(:, :), pointer, contiguous :: dbl2d