NetCDF-3 Classic trim-2x3z.nc {

  dimensions:
    time = 1 ;
    m = 4 ;
    n = 5 ;
    Layer = 4 ;
    LayerInterf = 5 ;
    bounds2 = 2 ;
    bounds4 = 4 ;
    grid_m = 5 ;
    grid_n = 6 ;

  variables:
    // Preference 'PRESERVE_FVD':  false,
    // dimensions consistent with ncBrowse, not with native MATLAB netcdf package.
    double time(time), shape = [1]
      :standard_name = "time" 
      :long_name = "time" 
      :units = "days since 1970-01-01 00:00:00 +00:00" 
      :axis = "T" 
      :delft3d_name = "map-info-series:ITMAPC map-const:ITDATE map-const:DT map-const:TUNIT" 
      :actual_range = 0 0 
    single m(m), shape = [4]
      :standard_name = "projection_x_coordinate" 
      :long_name = "grid cell centres, x-coordinate" 
      :units = "m" 
      :axis = "X" 
      :_FillValue = NaN f
      :actual_range = 200 400 
      :grid_mapping = "CRS" 
      :delft3d_name = "map-const:XZ map-const:KCS TEMPOUT:XWAT TEMPOUT:CODW map-const:COORDINATES" 
    single n(n), shape = [5]
      :standard_name = "projection_y_coordinate" 
      :long_name = "grid cell centres, y-coordinate" 
      :units = "m" 
      :axis = "Y" 
      :_FillValue = NaN f
      :actual_range = 150 350 
      :grid_mapping = "CRS" 
      :delft3d_name = "map-const:YZ map-const:KCS TEMPOUT:YWAT TEMPOUT:CODW map-const:COORDINATES" 
    int32 CRS([]), shape = [1]
      :name = "Amersfoort / RD New" 
      :epsg = 28992 
      :epsg_name = "Oblique Stereographic" 
      :grid_mapping_name = "" 
      :semi_major_axis = 6.3774e+006 
      :semi_minor_axis = 6.35608e+006 
      :inverse_flattening = 299.153 
      :latitude_of_projection_origin = 52.0922 
      :longitude_of_projection_origin = 5.23155 
      :false_easting = 155000 
      :false_northing = 463000 
      :scale_factor_at_projection_origin = 0.999908 
      :proj4_params = "+proj=sterea +lat_0=52.15616055555555 +lon_0=5.38763888888889 +k=0.999908 +x_0=155000 +y_0=463000 +ellps=bessel +units=m +towgs84=565.4174,50.3319,465.5542,-0.398957388243134,0.343987817378283,-1.87740163998045,4.0725 +no_defs" 
      :EPSG_code = "EPSG:28992" 
      :projection_name = "Dutch rijksdriekhoek system" 
      :wkt = "" 
      :comment = "value is equal to EPSG code" 
    single x(m,n), shape = [4 5]
      :standard_name = "projection_x_coordinate" 
      :long_name = "grid cell centres, x-coordinate" 
      :units = "m" 
      :axis = "X" 
      :_FillValue = NaN f
      :actual_range = 200 400 
      :grid_mapping = "CRS" 
      :delft3d_name = "map-const:XZ map-const:KCS TEMPOUT:XWAT TEMPOUT:CODW map-const:COORDINATES" 
      :bounds = "grid_x" 
    single y(m,n), shape = [4 5]
      :standard_name = "projection_y_coordinate" 
      :long_name = "grid cell centres, y-coordinate" 
      :units = "m" 
      :axis = "Y" 
      :_FillValue = NaN f
      :actual_range = 150 350 
      :grid_mapping = "CRS" 
      :delft3d_name = "map-const:YZ map-const:KCS TEMPOUT:YWAT TEMPOUT:CODW map-const:COORDINATES" 
      :bounds = "grid_y" 
    single grid_x(m,n,bounds4), shape = [4 5 4]
      :standard_name = "projection_x_coordinate" 
      :long_name = "grid cell corners, x-coordinate" 
      :units = "m" 
      :axis = "X" 
      :_FillValue = NaN f
      :actual_range = 100 500 
      :grid_mapping = "CRS" 
      :delft3d_name = "map-const:XCOR TEMPOUT:CODB map-const:COORDINATES" 
      :comment = "OpenEarth Matlab function nc_cf_bounds2cor.m reshapes it to a regular 2D array" 
    single grid_y(m,n,bounds4), shape = [4 5 4]
      :standard_name = "projection_y_coordinate" 
      :long_name = "grid cell corners, y-coordinate" 
      :axis = "Y" 
      :_FillValue = NaN f
      :actual_range = 100 400 
      :grid_mapping = "CRS" 
      :delft3d_name = "map-const:YCOR TEMPOUT:CODB map-const:COORDINATES" 
      :comment = "OpenEarth Matlab function nc_cf_bounds2cor.m reshapes it to a regular 2D array" 
    single longitude(m,n), shape = [4 5]
      :standard_name = "longitude" 
      :long_name = "grid cell centers, longitude-coordinate" 
      :units = "degrees_east" 
      :axis = "X" 
      :_FillValue = NaN f
      :actual_range = 3.3161 3.31881 
      :delft3d_name = "map-const:XZ map-const:KCS TEMPOUT:XWAT map-const:CODW map-const:COORDINATES" 
      :bounds = "grid_longitude" 
    single latitude(m,n), shape = [4 5]
      :standard_name = "latitude" 
      :long_name = "grid cell centers, latitude-coordinate" 
      :units = "degrees_north" 
      :axis = "Y" 
      :_FillValue = NaN f
      :actual_range = 47.9762 47.978 
      :delft3d_name = "map-const:YZ map-const:KCS map-const:YWAT map-const:CODW map-const:COORDINATES" 
      :bounds = "grid_latitude" 
    single grid_longitude(m,n,bounds4), shape = [4 5 4]
      :standard_name = "longitude" 
      :long_name = "grid cell corners, longitude-coordinate" 
      :units = "degrees_east" 
      :axis = "X" 
      :_FillValue = NaN f
      :actual_range = 3.31475 3.32017 
      :delft3d_name = "map-const:XCOR map-const:CODB map-const:COORDINATES" 
    single grid_latitude(m,n,bounds4), shape = [4 5 4]
      :standard_name = "latitude" 
      :long_name = "grid cell corners, latitude-coordinate" 
      :units = "degrees_north" 
      :axis = "Y" 
      :_FillValue = NaN f
      :actual_range = 47.9757 47.9785 
      :delft3d_name = "map-const:YCOR map-const:CODB map-const:COORDINATES" 
    single k(Layer), shape = [4]
      :long_name = "layer index" 
      :units = "1" 
      :axis = "Z" 
      :positive = "up" 
      :comment = "The surface layer has index k=1, the bottom layer has index kmax." 
      :delft3d_name = "map-const:KMAX map-const:LAYER_MODEL" 
    single Layer(Layer), shape = [4]
      :long_name = "z at layer midpoints" 
      :standard_name = "altitude" 
      :units = "m" 
      :positive = "up" 
      :actual_range = -114 -24 
      :comment = "The bottom layer has index k=1 and is the bottom depth, the surface layer has index kmax and is z=free water surface." 
      :delft3d_name = "map-const:KMAX map-const:LAYER_MODEL map-const:ZK" 
    single LayerInterf(LayerInterf), shape = [5]
      :long_name = "z at layer interfaces" 
      :units = "m" 
      :positive = "up" 
      :actual_range = -120 0 
      :comment = "The bottom layer has index k=1 and is the bottom depth, the surface layer has index kmax and is z=free water surface." 
      :delft3d_name = "map-const:KMAX map-const:LAYER_MODEL map-const:ZK" 
    single grid_depth(grid_m,grid_n), shape = [5 6]
      :standard_name = "sea_floor_depth_below_sea_level" 
      :long_name = "depth of cell corners" 
      :units = "m" 
      :positive = "down" 
      :_FillValue = NaN f
      :actual_range = NaN NaN 
      :delft3d_name = "map-const:DP map-const:DP0 map-const:DPS map-const:DRYFLP" 
      :comment = "" 
    single depth(m,n), shape = [4 5]
      :standard_name = "sea_floor_depth_below_sea_level" 
      :long_name = "depth of cell centers" 
      :units = "m" 
      :positive = "down" 
      :coordinates = "x y" 
      :_FillValue = NaN f
      :actual_range = -120 -120 
      :delft3d_name = "map-const:DPS0 map-const:KCS" 
      :comment = "" 
    single zactive(m,n), shape = [4 5]
      :long_name = "active water-level point" 
      :units = "1" 
      :coordinates = "x y" 
      :_FillValue = NaN f
      :actual_range = 0 1 
      :delft3d_name = "map-const:KCS" 
    single area(m,n), shape = [4 5]
      :standard_name = "" 
      :long_name = "area of grid cells" 
      :units = "m2" 
      :coordinates = "x y" 
      :_FillValue = NaN f
      :actual_range = 20000 20000 
      :delft3d_name = "map-const:GSQS map-const:KCS" 
      :comment = "This is the exact area spanned geometrically by the 4 corner points. This is not identical to the area GSQS used internally in Delft3D for mass-conservation!" 
    single waterlevel(time,m,n), shape = [1 4 5]
      :standard_name = "sea_surface_height" 
      :long_name = "water level" 
      :units = "m" 
      :positive = "up" 
      :coordinates = "x y" 
      :_FillValue = NaN f
      :actual_range = 0 0 
      :delft3d_name = "map-series:S1 map-const:KCS" 
    single velocity_x(time,Layer,m,n), shape = [1 4 4 5]
      :standard_name = "sea_water_x_velocity" 
      :long_name = "velocity, x-component" 
      :units = "m/s" 
      :coordinates = "x y" 
      :_FillValue = NaN f
      :actual_range = 0.03 0.06 
      :delft3d_name = "map-series:U1 map-series:V1 map-const:ALFAS map-const:KCU map-const:KCV map-const:KFU map-const:KFV map-const:KCS" 
    single velocity_y(time,Layer,m,n), shape = [1 4 4 5]
      :standard_name = "sea_water_y_velocity" 
      :long_name = "velocity, y-component" 
      :units = "m/s" 
      :coordinates = "x y" 
      :_FillValue = NaN f
      :actual_range = 0.04 0.08 
      :delft3d_name = "map-series:U1 map-series:V1 map-const:ALFAS map-const:KCU map-const:KCV map-const:KFU map-const:KFV map-const:KCS" 
    single velocity_omega(time,LayerInterf,m,n), shape = [1 5 4 5]
      :standard_name = "upward_sea_water_velocity" 
      :long_name = "relative layer velocity, z-component" 
      :units = "m/s" 
      :positive = "up" 
      :coordinates = "x y" 
      :_FillValue = NaN f
      :actual_range = Inf -Inf 
      :delft3d_name = "map-series:W map-const:KCS" 
    single velocity_z(time,Layer,m,n), shape = [1 4 4 5]
      :standard_name = "upward_sea_water_velocity" 
      :long_name = "velocity, z-component" 
      :units = "m/s" 
      :positive = "up" 
      :coordinates = "x y" 
      :_FillValue = NaN f
      :actual_range = 0 0 
      :delft3d_name = "map-series:WPHY map-const:KCS" 
    single tau_x(time,m,n), shape = [1 4 5]
      :standard_name = "surface_downward_x_stress" 
      :long_name = "bed shear stress, x-component" 
      :units = "N/m2" 
      :coordinates = "x y" 
      :_FillValue = NaN f
      :actual_range = 0.00895443 0.0174341 
      :delft3d_name = "map-series:TAUKSI map-const:TAUETA map-const:ALFAS map-const:KCS" 
    single tau_y(time,m,n), shape = [1 4 5]
      :standard_name = "surface_downward_y_stress" 
      :long_name = "bed shear stress, y-component" 
      :units = "N/m2" 
      :coordinates = "x y" 
      :_FillValue = NaN f
      :actual_range = 0.0149559 0.0205304 
      :delft3d_name = "map-series:TAUKSI map-const:TAUETA map-const:ALFAS map-const:KCS" 
    single density(time,Layer,m,n), shape = [1 4 4 5]
      :standard_name = "sea_water_density" 
      :long_name = "density" 
      :units = "kg/m3" 
      :coordinates = "x y" 
      :_FillValue = NaN f
      :actual_range = 1023.66 1025.42 
      :delft3d_name = "map-series:RHO map-const:KCS" 
    single salinity(time,Layer,m,n), shape = [1 4 4 5]
      :standard_name = "sea_water_salinity" 
      :long_name = "salinity" 
      :units = "ppt" 
      :coordinates = "x y" 
      :_FillValue = NaN f
      :actual_range = 31 34 
      :delft3d_name = "map-series:R1 map-const:KCS map-const:NAMCON map-const:LSTCI" 
    single temperature(time,Layer,m,n), shape = [1 4 4 5]
      :standard_name = "sea_water_temperature" 
      :long_name = "temperature" 
      :units = "degree_Celsius" 
      :coordinates = "x y" 
      :_FillValue = NaN f
      :actual_range = 11 14 
      :delft3d_name = "map-series:R1 map-const:KCS map-const:NAMCON map-const:LSTCI" 
    single tke(time,LayerInterf,m,n), shape = [1 5 4 5]
      :standard_name = "" 
      :long_name = "turbulent kinetic energy" 
      :units = "m2/s2" 
      :coordinates = "x y" 
      :_FillValue = NaN f
      :actual_range = 1e-007 1e-007 
      :delft3d_name = "map-series:RTUR1 map-const:KCS map-const:NAMCON map-const:LSTCI map-const:LTUR" 
    single eps(time,LayerInterf,m,n), shape = [1 5 4 5]
      :standard_name = "" 
      :long_name = "turbulent energy dissipation" 
      :units = "m2/s3" 
      :coordinates = "x y" 
      :_FillValue = NaN f
      :actual_range = 0 2.22597e+008 
      :delft3d_name = "map-series:RTUR1 map-const:KCS map-const:NAMCON map-const:LSTCI map-const:LTUR" 
    single viscosity_z(time,LayerInterf,m,n), shape = [1 5 4 5]
      :standard_name = "" 
      :long_name = "Vertical eddy viscosity-3D" 
      :units = "m^2/s" 
      :coordinates = "x y" 
      :_FillValue = NaN f
      :actual_range = 0 0 
      :delft3d_name = "map-series:VICWW map-const:KCS" 
    single diffusivity_z(time,LayerInterf,m,n), shape = [1 5 4 5]
      :standard_name = "" 
      :long_name = "Vertical eddy diffusivity-3D" 
      :units = "m^2/s" 
      :coordinates = "x y" 
      :_FillValue = NaN f
      :actual_range = 0 0 
      :delft3d_name = "map-series:DICWW map-const:KCS" 
    single Ri(time,LayerInterf,m,n), shape = [1 5 4 5]
      :standard_name = "" 
      :long_name = "Richardson number" 
      :units = "-" 
      :coordinates = "x y" 
      :_FillValue = NaN f
      :actual_range = 0 0 
      :delft3d_name = "map-series:RICH map-const:KCS" 

  //global Attributes:
      :title = "" 
      :institution = "" 
      :source = "Delft3D trim file" 
      :history = "Original filename: trim-2x3z.dat, Deltares, Delft3D-FLOW Version 4.00.04.757, Aug 23 2011, 18:06:51, file version: 3.54.28, file date:2013-01-07 15:51:53, transformation to netCDF: $HeadURL: https://svn.oss.deltares.nl/repos/openearthtools/trunk/matlab/applications/delft3d/vs_trim2nc.m $ $Id: vs_trim2nc.m 7814 2012-12-10 16:59:37Z boer_g $" 
      :references = "http://svn.oss.deltares.nl" 
      :email = "" 
      :comment = "" 
      :version = "Deltares, Delft3D-FLOW Version 4.00.04.757, Aug 23 2011, 18:06:51, file version: 3.54.28" 
      :Conventions = "CF-1.6" 
      :terms_for_use = "These data can be used freely for research purposes provided that the following source is acknowledged: " 
      :disclaimer = "This data is made available in the hope that it will be useful, but WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE." 
      :delft3d_description = "2x3 cells, 4 layers, 5 times  
test for netCDF conversion    
waterlevel initially 0        
|velocity|=1m/s at 60deg at t0
salinity: last digit=k        
temperature: last digit=k     
<line 7>                      
<line 8>                      
<line 9>                      
<line 10>                     
" 


}