Changeset 3668 for palm/trunk/SCRIPTS/palm_csd_files/palm_csd_netcdf_interface.py

Timestamp:

Jan 14, 2019 12:49:24 PM (5 years ago)

Author:

maronga

Message:

removed most_methods circular and lookup. added improved version of palm_csd

File:

• palm/trunk/SCRIPTS/palm_csd_files/palm_csd_netcdf_interface.py (modified) (3 diffs)

palm/trunk/SCRIPTS/palm_csd_files/palm_csd_netcdf_interface.py

@@ -25 +25 @@
 # -----------------
 # $Id$
+# Some improvements and new routines
+# 
+# 3629 2018-12-13 12:18:54Z maronga
 # Some new routines
 # 
@@ -39 +42 @@
 
 from netCDF4 import Dataset
+
+def nc_read_from_file_1d_all(filename, varname):
+
+
+   import numpy as np
+   import sys
+   
+   try:
+      f = open(filename)
+      f.close()
+#      print("Load: " + filename + ".")
+   except FileNotFoundError:
+      print("Error: " + filename + ". No such file. Aborting...")
+      sys.exit(1)
+   
+   nc_file = Dataset(filename, "r+", format="NETCDF4")
+   tmp_array = np.array(nc_file.variables[varname][:] , dtype=type(nc_file.variables[varname]))
+
+   return tmp_array
 
 def nc_read_from_file_2d(filename, varname, x0, x1, y0, y1):
@@ -174 +196 @@
    return 0
 
+def nc_write_global_attribute(filename,attribute,value):
+
+   import datetime
+   import os
+
+   print("Writing attribute " + attribute + " to file...")
+
+   f =  Dataset(filename, "a", format="NETCDF4")
+   
+   f.setncattr(attribute,value)
+
+   f.close()
+
 
 def nc_write_dimension(filename,varname,array,datatype):