diff --git a/CHANGELOG.md b/CHANGELOG.md
index c79027aa14aae72985ec10f62c7124cbe2ca493f..8561210f1d6cebc0da7853955423016879ff46aa 100644
--- a/CHANGELOG.md
+++ b/CHANGELOG.md
@@ -3,6 +3,7 @@
 ### 1.1.1
 * Release date: ???
 * Changes and new features:
+  * Write 2D string data to save variables from shapefiles after doing a spatial join ([#49](https://earth.bsc.es/gitlab/es/NES/-/issues/49))
   * Bugs fixing:
     * Bug on `cell_methods` serial write ([#53](https://earth.bsc.es/gitlab/es/NES/-/issues/53))
     * Horizontal Interpolation Conservative: Improvement on memory usage when calculating the weight matrix ([#54](https://earth.bsc.es/gitlab/es/NES/-/issues/54))
diff --git a/nes/nc_projections/default_nes.py b/nes/nc_projections/default_nes.py
index da83b2b923929004314b5a81417bb33b94b2d4fb..7b0759831cd17644ac85f83e64e03354738acea3 100644
--- a/nes/nc_projections/default_nes.py
+++ b/nes/nc_projections/default_nes.py
@@ -88,7 +88,7 @@ class Nes(object):
     """
     def __init__(self, comm=None, path=None, info=False, dataset=None, xarray=False, parallel_method='Y',
                  avoid_first_hours=0, avoid_last_hours=0, first_level=0, last_level=None, create_nes=False,
-                 balanced=False, times=None, **kwargs):
+                 balanced=False, times=None, strlen=75, **kwargs):
         """
         Initialize the Nes class
 
@@ -246,6 +246,9 @@ class Nes(object):
         # Writing options
         self.zip_lvl = 0
 
+        # Get string length
+        self.strlen = strlen
+
         # Dimensions information
         self._var_dim = None
         self._lat_dim = None
@@ -318,6 +321,7 @@ class Nes(object):
             del self.lat_bnds
             del self._lon_bnds
             del self.lon_bnds
+            del self.strlen
             del self.shapefile
             for cell_measure in self.cell_measures.keys():
                 if self.cell_measures[cell_measure]['data'] is not None:
@@ -1425,6 +1429,7 @@ class Nes(object):
         """
 
         units = self.__parse_time_unit(time.units)
+
         if not hasattr(time, 'calendar'):
             calendar = 'standard'
         else:
@@ -1479,7 +1484,6 @@ class Nes(object):
             if self.master:
                 nc_var = self.netcdf.variables['time']
                 time_data, units, calendar = self.__parse_time(nc_var)
-
                 time = num2date(time_data, units, calendar=calendar)
                 time = [aux.replace(second=0, microsecond=0) for aux in time]
             else:
@@ -1797,10 +1801,21 @@ class Nes(object):
                           self.read_axis_limits['x_min']:self.read_axis_limits['x_max']]
             data = data.reshape(1, 1, data.shape[-2], data.shape[-1])
         elif len(var_dims) == 3:
-            data = nc_var[self.read_axis_limits['t_min']:self.read_axis_limits['t_max'],
-                          self.read_axis_limits['y_min']:self.read_axis_limits['y_max'],
-                          self.read_axis_limits['x_min']:self.read_axis_limits['x_max']]
-            data = data.reshape(data.shape[-3], 1, data.shape[-2], data.shape[-1])
+            if 'strlen' in var_dims:
+                data = nc_var[self.read_axis_limits['y_min']:self.read_axis_limits['y_max'],
+                              self.read_axis_limits['x_min']:self.read_axis_limits['x_max'],
+                              :]
+                data_aux = np.empty(shape=(data.shape[0], data.shape[1]), dtype=np.object)         
+                for lat_n in range(data.shape[0]):
+                    for lon_n in range(data.shape[1]):
+                        data_aux[lat_n, lon_n] = ''.join(
+                            data[lat_n, lon_n].tostring().decode('ascii').replace('\x00', ''))
+                data = data_aux.reshape(1, 1, data_aux.shape[-2], data_aux.shape[-1])
+            else:
+                data = nc_var[self.read_axis_limits['t_min']:self.read_axis_limits['t_max'],
+                              self.read_axis_limits['y_min']:self.read_axis_limits['y_max'],
+                              self.read_axis_limits['x_min']:self.read_axis_limits['x_max']]
+                data = data.reshape(data.shape[-3], 1, data.shape[-2], data.shape[-1])
         elif len(var_dims) == 4:
             data = nc_var[self.read_axis_limits['t_min']:self.read_axis_limits['t_max'],
                           self.read_axis_limits['z_min']:self.read_axis_limits['z_max'],
@@ -2077,6 +2092,9 @@ class Nes(object):
         netcdf.createDimension('lon', len(self._lon['data']))
         netcdf.createDimension('lat', len(self._lat['data']))
 
+        # Create string length dimension
+        netcdf.createDimension('strlen', self.strlen)
+
         return None
 
     def _create_dimension_variables(self, netcdf):
@@ -2197,12 +2215,63 @@ class Nes(object):
 
         for i, (var_name, var_dict) in enumerate(self.variables.items()):
             if var_dict['data'] is not None:
+                
+                # Get dimensions
+                if var_dict['data'].shape == 4:
+                    var_dims = ('time', 'lev',) + self._var_dim
+                else:
+                    var_dims = self._var_dim
+
+                # Get data type
+                if 'dtype' in var_dict.keys():
+                    var_dtype = var_dict['dtype']
+                    if var_dtype != var_dict['data'].dtype:
+                        msg = "WARNING!!! "
+                        msg += "Different data types for variable {0}. ".format(var_name)
+                        msg += "Input dtype={0}. Data dtype={1}.".format(var_dtype, 
+                                                                        var_dict['data'].dtype)
+                        warnings.warn(msg)
+                        try:
+                            var_dict['data'] = var_dict['data'].astype(var_dtype)
+                        except Exception as e:  # TODO: Detect exception
+                            raise e("It was not possible to cast the data to the input dtype.")
+                else:
+                    var_dtype = var_dict['data'].dtype
+
+                # Transform objects into strings
+                if var_dtype == np.dtype(object):
+                    var_dict['data'] = var_dict['data'].astype(str)
+                    var_dtype = var_dict['data'].dtype
+
+                # Convert list of strings to chars for parallelization
+                try:
+                    unicode_type = len(max(var_dict['data'].flatten(), key=len))
+                    if ((var_dict['data'].dtype == np.dtype('<U' + str(unicode_type))) or 
+                        (var_dict['data'].dtype == np.str)):
+                        # Add strlen as a dimension if needed
+                        if 'strlen' not in var_dims:
+                            var_dims += ('strlen',)
+                        data_aux = np.chararray(shape=(var_dict['data'].shape[0], 
+                                                       var_dict['data'].shape[1], 
+                                                       self.strlen))
+                        for lat_n in range(var_dict['data'].shape[0]):
+                            for lon_n in range(var_dict['data'].shape[1]):
+                                chars = [v for v in var_dict['data'][lat_n][lon_n]]
+                                data = chars + ['']*(self.strlen-len(chars))
+                                for char_n, char in enumerate(data):
+                                    data_aux[lat_n, lon_n, char_n] = char.encode('ascii', 'ignore')
+                        
+                        var_dict['data'] = data_aux
+                        var_dtype = 'S1'
+                except:
+                    pass
+
                 if self.info:
                     print("Rank {0:03d}: Writing {1} var ({2}/{3})".format(
                         self.rank, var_name, i + 1, len(self.variables)))
                 try:
                     if not chunking:
-                        var = netcdf.createVariable(var_name, var_dict['data'].dtype, ('time', 'lev',) + self._var_dim,
+                        var = netcdf.createVariable(var_name, var_dtype, var_dims,
                                                     zlib=self.zip_lvl > 0, complevel=self.zip_lvl)
                     else:
                         if self.balanced:
@@ -2212,7 +2281,7 @@ class Nes(object):
                         else:
                             chunk_size = None
                         chunk_size = self.comm.bcast(chunk_size, root=0)
-                        var = netcdf.createVariable(var_name, var_dict['data'].dtype, ('time', 'lev',) + self._var_dim,
+                        var = netcdf.createVariable(var_name, var_dtype, var_dims,
                                                     zlib=self.zip_lvl > 0, complevel=self.zip_lvl,
                                                     chunksizes=chunk_size)
                     if self.info:
@@ -2228,31 +2297,47 @@ class Nes(object):
                         if att_name == 'data':
                             if self.info:
                                 print("Rank {0:03d}: Filling {1})".format(self.rank, var_name))
-                            try:
-                                var[self.write_axis_limits['t_min']:self.write_axis_limits['t_max'],
-                                    self.write_axis_limits['z_min']:self.write_axis_limits['z_max'],
-                                    self.write_axis_limits['y_min']:self.write_axis_limits['y_max'],
-                                    self.write_axis_limits['x_min']:self.write_axis_limits['x_max']] = att_value
-                            except ValueError:
-                                var[self.write_axis_limits['t_min']:self.write_axis_limits['t_max'],
-                                    0,
-                                    self.write_axis_limits['y_min']:self.write_axis_limits['y_max'],
-                                    self.write_axis_limits['x_min']:self.write_axis_limits['x_max']] = att_value
-                                # msg = "*WARNING* '{0}' variable is a 3D field. Setting it on first (0) layer.".format(
-                                #     var_name)
-                                # warn(msg)
-                            except IndexError:
-                                raise IndexError("Different shapes. out_shape={0}, data_shp={1}".format(
+                            if len(att_value.shape) == 4:
+                                try:
                                     var[self.write_axis_limits['t_min']:self.write_axis_limits['t_max'],
                                         self.write_axis_limits['z_min']:self.write_axis_limits['z_max'],
                                         self.write_axis_limits['y_min']:self.write_axis_limits['y_max'],
-                                        self.write_axis_limits['x_min']:self.write_axis_limits['x_max']].shape,
-                                    att_value.shape))
+                                        self.write_axis_limits['x_min']:self.write_axis_limits['x_max']] = att_value
+                                except ValueError:
+                                    var[self.write_axis_limits['t_min']:self.write_axis_limits['t_max'],
+                                        0,
+                                        self.write_axis_limits['y_min']:self.write_axis_limits['y_max'],
+                                        self.write_axis_limits['x_min']:self.write_axis_limits['x_max']] = att_value
+                                    # msg = "*WARNING* '{0}' variable is a 3D field. Setting it on first (0) layer.".format(
+                                    #     var_name)
+                                    # warn(msg)
+                                except IndexError:
+                                    raise IndexError("Different shapes. out_shape={0}, data_shp={1}".format(
+                                        var[self.write_axis_limits['t_min']:self.write_axis_limits['t_max'],
+                                            self.write_axis_limits['z_min']:self.write_axis_limits['z_max'],
+                                            self.write_axis_limits['y_min']:self.write_axis_limits['y_max'],
+                                            self.write_axis_limits['x_min']:self.write_axis_limits['x_max']].shape,
+                                        att_value.shape))
+                            elif len(att_value.shape) == 3:
+                                if 'strlen' in var_dims:
+                                    try:
+                                        var[self.write_axis_limits['y_min']:self.write_axis_limits['y_max'],
+                                            self.write_axis_limits['x_min']:self.write_axis_limits['x_max'],
+                                            :] = att_value
+                                    except IndexError:
+                                        raise IndexError("Different shapes. out_shape={0}, data_shp={1}".format(
+                                            var[self.write_axis_limits['y_min']:self.write_axis_limits['y_max'],
+                                                self.write_axis_limits['x_min']:self.write_axis_limits['x_max'], 
+                                                :].shape,
+                                            att_value.shape))
+                                else:
+                                    raise NotImplementedError('It is not possible to write 3D variables.')
                             if self.info:
                                 print("Rank {0:03d}: Var {1} data ({2}/{3})".format(self.rank, var_name, i + 1,
                                                                                     len(self.variables)))
                         elif att_name not in ['chunk_size', 'var_dims', 'dimensions']:
                             var.setncattr(att_name, att_value)
+
                     self._set_var_crs(var)
                     if self.info:
                         print("Rank {0:03d}: Var {1} completed ({2}/{3})".format(self.rank, var_name, i + 1,
diff --git a/nes/nc_projections/points_nes.py b/nes/nc_projections/points_nes.py
index 9c3d859b48ad1c157095d42b570e9977ddd1f0e0..5a2edaa418dae585cdfaaffc6ab6a38f8cb848f3 100644
--- a/nes/nc_projections/points_nes.py
+++ b/nes/nc_projections/points_nes.py
@@ -369,7 +369,7 @@ class PointsNes(Nes):
 
                     # Convert list of strings to chars for parallelization
                     try:
-                        unicode_type = len(max(var_dict['data'], key=len))
+                        unicode_type = len(max(var_dict['data'].flatten(), key=len))
                         if ((var_dict['data'].dtype == np.dtype('<U' + str(unicode_type))) or 
                             (var_dict['data'].dtype == np.str)):
                             # Add strlen as a dimension if needed
diff --git a/nes/nc_projections/points_nes_ghost.py b/nes/nc_projections/points_nes_ghost.py
index 1d95fa18baa925d5b930d411cea97d1737e7c5dd..d7cf59891e6355341d3c7e8506dd25e60bf8ba28 100644
--- a/nes/nc_projections/points_nes_ghost.py
+++ b/nes/nc_projections/points_nes_ghost.py
@@ -351,7 +351,7 @@ class PointsNesGHOST(PointsNes):
 
                     # Convert list of strings to chars for parallelization
                     try:
-                        unicode_type = len(max(var_dict['data'], key=len))
+                        unicode_type = len(max(var_dict['data'].flatten(), key=len))
                         if ((var_dict['data'].dtype == np.dtype('<U' + str(unicode_type))) or 
                             (var_dict['data'].dtype == np.str)):
                             # Add strlen as a dimension if needed
diff --git a/nes/nc_projections/points_nes_providentia.py b/nes/nc_projections/points_nes_providentia.py
index 737392f31da97abc6ca00aa3075127dc36305b51..9d576ad2923dad041c7f448513913855003f2fc7 100644
--- a/nes/nc_projections/points_nes_providentia.py
+++ b/nes/nc_projections/points_nes_providentia.py
@@ -387,7 +387,7 @@ class PointsNesProvidentia(PointsNes):
 
                     # Convert list of strings to chars for parallelization
                     try:
-                        unicode_type = len(max(var_dict['data'], key=len))
+                        unicode_type = len(max(var_dict['data'].flatten(), key=len))
                         if ((var_dict['data'].dtype == np.dtype('<U' + str(unicode_type))) or 
                             (var_dict['data'].dtype == np.str)):
                             # Add strlen as a dimension if needed
diff --git a/tests/2.1-test_spatial_join.py b/tests/2.1-test_spatial_join.py
index 2daab1e66d0803011c917727a3e684de3433cacf..ede7851c737d6b591f9187d393cbccdb820c36e6 100644
--- a/tests/2.1-test_spatial_join.py
+++ b/tests/2.1-test_spatial_join.py
@@ -4,6 +4,7 @@ import sys
 from mpi4py import MPI
 import pandas as pd
 import timeit
+import numpy as np
 from nes import *
 
 comm = MPI.COMM_WORLD
@@ -47,6 +48,21 @@ comm.Barrier()
 result.loc['calculate', test_name] = timeit.default_timer() - st_time
 print('SPATIAL JOIN - Rank {0:03d} - Shapefile: \n{1}'.format(rank, nessy.shapefile))
 
+# ADD TIMEZONES
+timezones = nessy.shapefile['tzid'].values.reshape([nessy.lat['data'].shape[0], 
+                                                    nessy.lon['data'].shape[-1]])
+nessy.variables['tz'] = {'data': timezones,}
+
+# WRITE
+st_time = timeit.default_timer()
+nessy.to_netcdf(test_name.replace(' ', '_') + "_{0:03d}.nc".format(size), info=True)
+comm.Barrier()
+result.loc['write', test_name] = timeit.default_timer() - st_time
+
+# REOPEN
+nessy = open_netcdf(test_name.replace(' ', '_') + "_{0:03d}.nc".format(size), parallel_method=parallel_method)
+nessy.load()
+
 comm.Barrier()
 if rank == 0:
     print(result.loc[:, test_name])
@@ -78,9 +94,23 @@ nessy = from_shapefile(shapefile_path, method='centroid', projection=projection,
                        n_lat=n_lat, n_lon=n_lon)
 comm.Barrier()
 result.loc['calculate', test_name] = timeit.default_timer() - st_time
-
 print('FROM SHAPEFILE - Rank {0:03d} - Shapefile: \n{1}'.format(rank, nessy.shapefile))
 
+# ADD TIMEZONES
+timezones = nessy.shapefile['tzid'].values.reshape([nessy.lat['data'].shape[0], 
+                                                    nessy.lon['data'].shape[-1]])
+nessy.variables['tz'] = {'data': timezones,}
+
+# WRITE
+st_time = timeit.default_timer()
+nessy.to_netcdf(test_name.replace(' ', '_') + "_{0:03d}.nc".format(size), info=True)
+comm.Barrier()
+result.loc['write', test_name] = timeit.default_timer() - st_time
+
+# REOPEN
+nessy = open_netcdf(test_name.replace(' ', '_') + "_{0:03d}.nc".format(size), parallel_method=parallel_method)
+nessy.load()
+
 comm.Barrier()
 if rank == 0:
     print(result.loc[:, test_name])
@@ -108,6 +138,21 @@ comm.Barrier()
 result.loc['calculate', test_name] = timeit.default_timer() - st_time
 print('SPATIAL JOIN - Rank {0:03d} - Shapefile: \n{1}'.format(rank, nessy.shapefile))
 
+# ADD TIMEZONES
+timezones = nessy.shapefile['tzid'].values.reshape([nessy.lat['data'].shape[0], 
+                                                    nessy.lon['data'].shape[-1]])
+nessy.variables['tz'] = {'data': timezones,}
+
+# WRITE
+st_time = timeit.default_timer()
+nessy.to_netcdf(test_name.replace(' ', '_') + "_{0:03d}.nc".format(size), info=True)
+comm.Barrier()
+result.loc['write', test_name] = timeit.default_timer() - st_time
+
+# REOPEN
+nessy = open_netcdf(test_name.replace(' ', '_') + "_{0:03d}.nc".format(size), parallel_method=parallel_method)
+nessy.load()
+
 comm.Barrier()
 if rank == 0:
     print(result.loc[:, test_name])
@@ -139,9 +184,23 @@ nessy = from_shapefile(shapefile_path, method='nearest', projection=projection,
                        n_lat=n_lat, n_lon=n_lon)
 comm.Barrier()
 result.loc['calculate', test_name] = timeit.default_timer() - st_time
-
 print('FROM SHAPEFILE - Rank {0:03d} - Shapefile: \n{1}'.format(rank, nessy.shapefile))
 
+# ADD TIMEZONES
+timezones = nessy.shapefile['tzid'].values.reshape([nessy.lat['data'].shape[0], 
+                                                    nessy.lon['data'].shape[-1]])
+nessy.variables['tz'] = {'data': timezones,}
+
+# WRITE
+st_time = timeit.default_timer()
+nessy.to_netcdf(test_name.replace(' ', '_') + "_{0:03d}.nc".format(size), info=True)
+comm.Barrier()
+result.loc['write', test_name] = timeit.default_timer() - st_time
+
+# REOPEN
+nessy = open_netcdf(test_name.replace(' ', '_') + "_{0:03d}.nc".format(size), parallel_method=parallel_method)
+nessy.load()
+
 comm.Barrier()
 if rank == 0:
     print(result.loc[:, test_name])
@@ -170,6 +229,21 @@ comm.Barrier()
 result.loc['calculate', test_name] = timeit.default_timer() - st_time
 print('SPATIAL JOIN - Rank {0:03d} - Shapefile: \n{1}'.format(rank, nessy.shapefile))
 
+# ADD TIMEZONES
+timezones = nessy.shapefile['tzid'].values.reshape([nessy.lat['data'].shape[0], 
+                                                    nessy.lon['data'].shape[-1]])
+nessy.variables['tz'] = {'data': timezones,}
+
+# WRITE
+st_time = timeit.default_timer()
+nessy.to_netcdf(test_name.replace(' ', '_') + "_{0:03d}.nc".format(size), info=True)
+comm.Barrier()
+result.loc['write', test_name] = timeit.default_timer() - st_time
+
+# REOPEN
+nessy = open_netcdf(test_name.replace(' ', '_') + "_{0:03d}.nc".format(size), parallel_method=parallel_method)
+nessy.load()
+
 comm.Barrier()
 if rank == 0:
     print(result.loc[:, test_name])
@@ -202,9 +276,23 @@ nessy = from_shapefile(shapefile_path, method='intersection', projection=project
                        n_lat=n_lat, n_lon=n_lon)
 comm.Barrier()
 result.loc['calculate', test_name] = timeit.default_timer() - st_time
-
 print('FROM SHAPEFILE - Rank {0:03d} - Shapefile: \n{1}'.format(rank, nessy.shapefile))
 
+# ADD TIMEZONES
+timezones = nessy.shapefile['tzid'].values.reshape([nessy.lat['data'].shape[0], 
+                                                    nessy.lon['data'].shape[-1]])
+nessy.variables['tz'] = {'data': timezones,}
+
+# WRITE
+st_time = timeit.default_timer()
+nessy.to_netcdf(test_name.replace(' ', '_') + "_{0:03d}.nc".format(size), info=True)
+comm.Barrier()
+result.loc['write', test_name] = timeit.default_timer() - st_time
+
+# REOPEN
+nessy = open_netcdf(test_name.replace(' ', '_') + "_{0:03d}.nc".format(size), parallel_method=parallel_method)
+nessy.load()
+
 comm.Barrier()
 if rank == 0:
     print(result.loc[:, test_name])
diff --git a/tests/2.3-test_bounds.py b/tests/2.3-test_bounds.py
index 3e1c85a3ae0e9eec4e0d40c8c2678f3008dcd037..3896a6f49c65abcd350232d5a4642f0714ff2c87 100644
--- a/tests/2.3-test_bounds.py
+++ b/tests/2.3-test_bounds.py
@@ -142,9 +142,7 @@ nessy_5.to_netcdf(test_name.replace(' ', '_') + "_{0:03d}.nc".format(size), info
 comm.Barrier()
 result.loc['write', test_name] = timeit.default_timer() - st_time
 
-# REOPENcomm.Barrier()
-result.loc['write', test_name] = timeit.default_timer() - st_time
-
+# REOPEN
 nessy_6 = open_netcdf(test_name.replace(' ', '_') + "_{0:03d}.nc".format(size), info=True)
 
 # LOAD DATA AND EXPLORE BOUNDS
diff --git a/tutorials/5.Geospatial/5.2.Spatial_Join.ipynb b/tutorials/5.Geospatial/5.2.Spatial_Join.ipynb
index 8ebefeef0f2536b9f35f135f3a3cdc3085aa81d7..98c4665cccc3e7bca1d5badb56e54ce5828c98e2 100644
--- a/tutorials/5.Geospatial/5.2.Spatial_Join.ipynb
+++ b/tutorials/5.Geospatial/5.2.Spatial_Join.ipynb
@@ -16,6 +16,7 @@
     "from nes import *\n",
     "import geopandas as gpd\n",
     "import pandas as pd\n",
+    "import numpy as np\n",
     "import matplotlib.pyplot as plt\n",
     "\n",
     "pd.options.mode.chained_assignment = None"
@@ -82,7 +83,7 @@
      "name": "stderr",
      "output_type": "stream",
      "text": [
-      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2666: UserWarning: Shapefile does not exist. It will be created now.\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2814: UserWarning: Shapefile does not exist. It will be created now.\n",
       "  warnings.warn(msg)\n"
      ]
     }
@@ -219,7 +220,7 @@
    "cell_type": "markdown",
    "metadata": {},
    "source": [
-    "### Add data of last timestep"
+    "### Save timezones data as variable"
    ]
   },
   {
@@ -228,7 +229,7 @@
    "metadata": {},
    "outputs": [],
    "source": [
-    "grid.keep_vars('O3')"
+    "timezones = grid.shapefile['tzid'].values.reshape([grid.lat['data'].shape[0], grid.lon['data'].shape[-1]])"
    ]
   },
   {
@@ -236,13 +237,565 @@
    "execution_count": 9,
    "metadata": {},
    "outputs": [],
+   "source": [
+    "grid.variables['tz'] = {'data': timezones,}"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 11,
+   "metadata": {},
+   "outputs": [
+    {
+     "name": "stderr",
+     "output_type": "stream",
+     "text": [
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable lmp was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable IM was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable JM was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable LM was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable IHRST was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable I_PAR_STA was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable J_PAR_STA was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable NPHS was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable NCLOD was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable NHEAT was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable NPREC was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable NRDLW was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable NRDSW was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable NSRFC was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable AVGMAXLEN was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable MDRMINout was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable MDRMAXout was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable MDIMINout was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable MDIMAXout was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable IDAT was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable DXH was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable SG1 was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable SG2 was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable DSG1 was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable DSG2 was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable SGML1 was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable SGML2 was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable SLDPTH was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable ISLTYP was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable IVGTYP was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable NCFRCV was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable NCFRST was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable FIS was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable GLAT was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable GLON was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable PD was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable VLAT was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable VLON was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable ACPREC was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable CUPREC was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable MIXHT was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable PBLH was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable RLWTOA was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable RSWIN was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable U10 was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable USTAR was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable V10 was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable RMOL was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable T2 was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable relative_humidity_2m was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable T was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable U was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable V was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable SH2O was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable SMC was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable STC was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable AERO_ACPREC was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable AERO_CUPREC was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable AERO_DEPDRY was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable AERO_OPT_R was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable DRE_SW_TOA was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable DRE_SW_SFC was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable DRE_LW_TOA was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable DRE_LW_SFC was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable ENG_SW_SFC was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable ADRYDEP was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable WETDEP was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable PH_NO2 was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable HSUM was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable POLR was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable aerosol_optical_depth_dim was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable aerosol_optical_depth was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable satellite_AOD_dim was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable satellite_AOD was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable aerosol_loading_dim was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable aerosol_loading was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable clear_sky_AOD_dim was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable clear_sky_AOD was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable layer_thickness was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable mid_layer_pressure was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable interface_pressure was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable relative_humidity was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable mid_layer_height was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable mid_layer_height_agl was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable air_density was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable dry_pm10_mass was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable dry_pm2p5_mass was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable QC was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable QR was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable QS was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable QG was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable aero_dust_001 was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable aero_dust_002 was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable aero_dust_003 was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable aero_dust_004 was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable aero_dust_005 was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable aero_dust_006 was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable aero_dust_007 was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable aero_dust_008 was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable aero_ssa_001 was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable aero_ssa_002 was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable aero_ssa_003 was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable aero_ssa_004 was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable aero_ssa_005 was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable aero_ssa_006 was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable aero_ssa_007 was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable aero_ssa_008 was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable aero_om_001 was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable aero_om_002 was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable aero_om_003 was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable aero_om_004 was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable aero_om_005 was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable aero_om_006 was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable aero_bc_001 was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable aero_bc_002 was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable aero_so4_001 was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable aero_no3_001 was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable aero_no3_002 was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable aero_no3_003 was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable aero_nh4_001 was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable aero_unsp_001 was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable aero_unsp_002 was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable aero_unsp_003 was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable aero_unsp_004 was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable aero_unsp_005 was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable NO2 was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable NO was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable O3 was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable NO3 was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable N2O5 was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable HNO3 was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable HONO was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable PNA was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable H2O2 was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable NTR was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable ROOH was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable FORM was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable ALD2 was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable ALDX was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable PAR was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable CO was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable MEPX was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable MEOH was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable FACD was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable PAN was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable PACD was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable AACD was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable PANX was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable OLE was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable ETH was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable IOLE was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable TOL was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable CRES was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable OPEN was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable MGLY was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable XYL was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable ISOP was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable ISPD was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable TERP was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable SO2 was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable SULF was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable ETOH was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable ETHA was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable CL2 was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable HOCL was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable FMCL was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable HCL was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable BENZENE was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable SESQ was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable NH3 was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable DMS was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable SOAP_I was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable SOAP_T was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable SOAP_F was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable SOAP_A was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable O was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable O1D was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable OH was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable HO2 was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable XO2 was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable XO2N was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable MEO2 was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable HCO3 was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable C2O3 was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable CXO3 was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable ROR was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable TO2 was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable TOLRO2 was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable CRO was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable XYLRO2 was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable ISOPRXN was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable TRPRXN was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable SULRXN was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable CL was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable CLO was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable TOLNRXN was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable TOLHRXN was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable XYLNRXN was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable XYLHRXN was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable BENZRO2 was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable BNZNRXN was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable BNZHRXN was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable SESQRXN was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable aerosol_extinction_dim was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable aerosol_extinction_DUST_1 was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable aerosol_extinction_DUST_2 was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable aerosol_extinction_DUST_3 was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable aerosol_extinction_DUST_4 was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable aerosol_extinction_DUST_5 was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable aerosol_extinction_DUST_6 was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable aerosol_extinction_DUST_7 was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable aerosol_extinction_DUST_8 was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable aerosol_extinction_SALT_total was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable aerosol_extinction_OM_total was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable aerosol_extinction_BC_total was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable aerosol_extinction_SO4_total was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable aerosol_extinction_NO3_total was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable aerosol_extinction_NH4_total was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable aerosol_extinction_UNSPC_1 was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable aerosol_extinction_UNSPC_2 was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable aerosol_extinction_UNSPC_3 was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable aerosol_extinction_UNSPC_4 was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n",
+      "/esarchive/scratch/avilanova/software/NES/nes/nc_projections/default_nes.py:2353: UserWarning: WARNING!!! Variable aerosol_extinction_UNSPC_5 was not loaded. It will not be written.\n",
+      "  warnings.warn(msg)\n"
+     ]
+    }
+   ],
+   "source": [
+    "grid.to_netcdf('grid_with_tz.nc')"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 12,
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "grid_with_tz = open_netcdf('grid_with_tz.nc')"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 13,
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "grid_with_tz.keep_vars('tz')"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 14,
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "grid_with_tz.load()"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 15,
+   "metadata": {},
+   "outputs": [
+    {
+     "data": {
+      "text/plain": [
+       "{'data': array([[[['nan', 'nan', 'nan', ..., 'Asia/Riyadh', 'Asia/Riyadh',\n",
+       "           'Asia/Aden'],\n",
+       "          ['nan', 'nan', 'nan', ..., 'Asia/Riyadh', 'Asia/Riyadh',\n",
+       "           'Asia/Riyadh'],\n",
+       "          ['nan', 'nan', 'nan', ..., 'Asia/Riyadh', 'Asia/Riyadh',\n",
+       "           'Asia/Riyadh'],\n",
+       "          ...,\n",
+       "          ['America/Pangnirtung', 'America/Pangnirtung',\n",
+       "           'America/Pangnirtung', ..., 'Asia/Tomsk', 'Asia/Tomsk',\n",
+       "           'Asia/Tomsk'],\n",
+       "          ['America/Pangnirtung', 'America/Pangnirtung',\n",
+       "           'America/Pangnirtung', ..., 'Asia/Tomsk', 'Asia/Tomsk',\n",
+       "           'Asia/Tomsk'],\n",
+       "          ['America/Toronto', 'America/Iqaluit', 'America/Pangnirtung',\n",
+       "           ..., 'Asia/Tomsk', 'Asia/Tomsk', 'Asia/Krasnoyarsk']]]],\n",
+       "       dtype=object),\n",
+       " 'dimensions': ('rlat', 'rlon', 'strlen'),\n",
+       " 'grid_mapping': 'rotated_pole',\n",
+       " 'coordinates': 'lat lon'}"
+      ]
+     },
+     "execution_count": 15,
+     "metadata": {},
+     "output_type": "execute_result"
+    }
+   ],
+   "source": [
+    "grid_with_tz.variables['tz']"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "### Add ozone data of last timestep"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 16,
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "grid.keep_vars('O3')"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 17,
+   "metadata": {},
+   "outputs": [],
    "source": [
     "grid.load()"
    ]
   },
   {
    "cell_type": "code",
-   "execution_count": 10,
+   "execution_count": 18,
    "metadata": {},
    "outputs": [],
    "source": [
@@ -251,7 +804,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 11,
+   "execution_count": 19,
    "metadata": {},
    "outputs": [
     {
@@ -390,7 +943,7 @@
        "[95121 rows x 3 columns]"
       ]
      },
-     "execution_count": 11,
+     "execution_count": 19,
      "metadata": {},
      "output_type": "execute_result"
     }
@@ -408,7 +961,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 12,
+   "execution_count": 20,
    "metadata": {},
    "outputs": [
     {
@@ -444,7 +997,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 13,
+   "execution_count": 21,
    "metadata": {},
    "outputs": [
     {
@@ -478,7 +1031,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 14,
+   "execution_count": 22,
    "metadata": {},
    "outputs": [
     {
@@ -514,7 +1067,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 15,
+   "execution_count": 23,
    "metadata": {},
    "outputs": [],
    "source": [
@@ -536,7 +1089,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 16,
+   "execution_count": 24,
    "metadata": {},
    "outputs": [],
    "source": [
@@ -549,7 +1102,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 17,
+   "execution_count": 25,
    "metadata": {},
    "outputs": [
     {
@@ -614,54 +1167,54 @@
        "      <td>...</td>\n",
        "    </tr>\n",
        "    <tr>\n",
-       "      <th>22495</th>\n",
-       "      <td>POLYGON ((18.00000 64.80000, 18.20000 64.80000...</td>\n",
-       "      <td>Europe/Stockholm</td>\n",
+       "      <th>14995</th>\n",
+       "      <td>POLYGON ((18.00000 54.80000, 18.20000 54.80000...</td>\n",
+       "      <td>Europe/Warsaw</td>\n",
        "    </tr>\n",
        "    <tr>\n",
-       "      <th>22496</th>\n",
-       "      <td>POLYGON ((18.20000 64.80000, 18.40000 64.80000...</td>\n",
-       "      <td>Europe/Stockholm</td>\n",
+       "      <th>14996</th>\n",
+       "      <td>POLYGON ((18.20000 54.80000, 18.40000 54.80000...</td>\n",
+       "      <td>Europe/Warsaw</td>\n",
        "    </tr>\n",
        "    <tr>\n",
-       "      <th>22497</th>\n",
-       "      <td>POLYGON ((18.40000 64.80000, 18.60000 64.80000...</td>\n",
-       "      <td>Europe/Stockholm</td>\n",
+       "      <th>14997</th>\n",
+       "      <td>POLYGON ((18.40000 54.80000, 18.60000 54.80000...</td>\n",
+       "      <td>Europe/Warsaw</td>\n",
        "    </tr>\n",
        "    <tr>\n",
-       "      <th>22498</th>\n",
-       "      <td>POLYGON ((18.60000 64.80000, 18.80000 64.80000...</td>\n",
-       "      <td>Europe/Stockholm</td>\n",
+       "      <th>14998</th>\n",
+       "      <td>POLYGON ((18.60000 54.80000, 18.80000 54.80000...</td>\n",
+       "      <td>Europe/Warsaw</td>\n",
        "    </tr>\n",
        "    <tr>\n",
-       "      <th>22499</th>\n",
-       "      <td>POLYGON ((18.80000 64.80000, 19.00000 64.80000...</td>\n",
-       "      <td>Europe/Stockholm</td>\n",
+       "      <th>14999</th>\n",
+       "      <td>POLYGON ((18.80000 54.80000, 19.00000 54.80000...</td>\n",
+       "      <td>NaN</td>\n",
        "    </tr>\n",
        "  </tbody>\n",
        "</table>\n",
-       "<p>22500 rows × 2 columns</p>\n",
+       "<p>15000 rows × 2 columns</p>\n",
        "</div>"
       ],
       "text/plain": [
-       "                                                geometry              tzid\n",
-       "FID                                                                       \n",
-       "0      POLYGON ((-11.00000 35.00000, -10.80000 35.000...               NaN\n",
-       "1      POLYGON ((-10.80000 35.00000, -10.60000 35.000...               NaN\n",
-       "2      POLYGON ((-10.60000 35.00000, -10.40000 35.000...               NaN\n",
-       "3      POLYGON ((-10.40000 35.00000, -10.20000 35.000...               NaN\n",
-       "4      POLYGON ((-10.20000 35.00000, -10.00000 35.000...               NaN\n",
-       "...                                                  ...               ...\n",
-       "22495  POLYGON ((18.00000 64.80000, 18.20000 64.80000...  Europe/Stockholm\n",
-       "22496  POLYGON ((18.20000 64.80000, 18.40000 64.80000...  Europe/Stockholm\n",
-       "22497  POLYGON ((18.40000 64.80000, 18.60000 64.80000...  Europe/Stockholm\n",
-       "22498  POLYGON ((18.60000 64.80000, 18.80000 64.80000...  Europe/Stockholm\n",
-       "22499  POLYGON ((18.80000 64.80000, 19.00000 64.80000...  Europe/Stockholm\n",
+       "                                                geometry           tzid\n",
+       "FID                                                                    \n",
+       "0      POLYGON ((-11.00000 35.00000, -10.80000 35.000...            NaN\n",
+       "1      POLYGON ((-10.80000 35.00000, -10.60000 35.000...            NaN\n",
+       "2      POLYGON ((-10.60000 35.00000, -10.40000 35.000...            NaN\n",
+       "3      POLYGON ((-10.40000 35.00000, -10.20000 35.000...            NaN\n",
+       "4      POLYGON ((-10.20000 35.00000, -10.00000 35.000...            NaN\n",
+       "...                                                  ...            ...\n",
+       "14995  POLYGON ((18.00000 54.80000, 18.20000 54.80000...  Europe/Warsaw\n",
+       "14996  POLYGON ((18.20000 54.80000, 18.40000 54.80000...  Europe/Warsaw\n",
+       "14997  POLYGON ((18.40000 54.80000, 18.60000 54.80000...  Europe/Warsaw\n",
+       "14998  POLYGON ((18.60000 54.80000, 18.80000 54.80000...  Europe/Warsaw\n",
+       "14999  POLYGON ((18.80000 54.80000, 19.00000 54.80000...            NaN\n",
        "\n",
-       "[22500 rows x 2 columns]"
+       "[15000 rows x 2 columns]"
       ]
      },
-     "execution_count": 17,
+     "execution_count": 25,
      "metadata": {},
      "output_type": "execute_result"
     }
@@ -679,7 +1232,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 18,
+   "execution_count": 26,
    "metadata": {},
    "outputs": [
     {