# lat variable attributes
nc_lat_var = nc_file.createVariable('lat', 'f4', ('station_nm',))
nc_lat_var.units = 'degree_north'
nc_lat_var.long_name = 'Latitude'
#nc_lat_var.standard_name = 'latitude'
# lon variable attributes
nc_lon_var = nc_file.createVariable('lon', 'f4',('station_nm',))
nc_lon_var.units = 'degree_east'
nc_lon_var.long_name = 'Longitude'
#nc_lon_var.standard_name = 'longitude'
# Create ncdf attributes
nc_file.WML_Conventions = 'CF-1.6'
nc_file.WML_featureType = 'timeSeries'
nc_file.WML_cdm_data_type = 'Station'
nc_file.WML_standard_name_vocabulary = 'CF-1.6'
nc_file.title = nc_title
nc_file.summary = nc_summary
nc_file.id = 'testing_id'
nc_file.naming_authory = 'testing_authority'
nc_file.WML_date_created = nc_date_create
nc_file.WML_creator_name = nc_creator_name
nc_file.creator_email = nc_creator_email
nc_file.project = nc_project
nc_file.processing_level = nc_proc_level
nc_file.WML_profile = 'single variable'
# data
dates = [datetime(2001,3,1)+n*timedelta(hours=12) for n in range(12)]
nc_time[:] = date2num(dates,units=nc_time.units,calendar=nc_time.calendar)
