Changeset 9325 for branches/2017/dev_merge_2017/NEMOGCM/TOOLS/MISCELLANEOUS/icb_pp.py

Timestamp:

2018-02-13T17:28:22+01:00 (6 years ago)

Author:

davestorkey

Message:

Small fix to TOOLS/MISCELLANEOUS/icb_pp.py. Ticket #2008.

File:

• branches/2017/dev_merge_2017/NEMOGCM/TOOLS/MISCELLANEOUS/icb_pp.py (modified) (1 diff)

branches/2017/dev_merge_2017/NEMOGCM/TOOLS/MISCELLANEOUS/icb_pp.py

@@ -92 +92 @@
 # Now flatten the lists and reduce to the unique spanning set
 #
-icu = np.concatenate(icu)
+try:
+    icu = np.concatenate(icu)
+except ValueError:
+    # No icebergs: create an empty output file.
+    print 'No icebergs in the model.'
+    fw = Dataset(pathstart+'0000.nc')
+    fo = Dataset(pathout, 'w', format='NETCDF4_CLASSIC')
+    ntrj = fo.createDimension('ntraj', None)
+    icbn = fo.createVariable('iceberg_number', 'i4',('ntraj'))
+    n = 0
+    for key, value in varlist.iteritems() :
+        if key != "iceberg_number" :
+            print 'key is ',key
+            oout = fo.createVariable(key, value.dtype, ('ntraj'),
+                                 zlib=True, complevel=1)
+            oout.long_name = fw.variables[key].getncattr('long_name')
+            oout.units = fw.variables[key].getncattr('units')
+            n = n + 1
+    fw.close()
+    fo.close()    
+    sys.exit()
+
 icu = np.unique(icu)
 times = np.concatenate(times)