diff --git a/src/pypromice/process/join_l3.py b/src/pypromice/process/join_l3.py
index 10a2c769..25d28e75 100644
--- a/src/pypromice/process/join_l3.py
+++ b/src/pypromice/process/join_l3.py
@@ -211,7 +211,7 @@ def join_l3(config_folder, site, folder_l3, folder_gcnet, outpath, variables, me
             filepath = os.path.join(folder_gcnet, stid+'.csv')
             isNead = True
         if not os.path.isfile(filepath):
-            logger.info(stid+' is from an project '+folder_l3+' or '+folder_gcnet)
+            logger.info(stid+' was listed as station but could not be found in '+folder_l3+' nor '+folder_gcnet)
             continue
 
         l3, _ = loadArr(filepath, isNead)            
diff --git a/src/pypromice/process/write.py b/src/pypromice/process/write.py
index e8d9e6a7..89ea2441 100644
--- a/src/pypromice/process/write.py
+++ b/src/pypromice/process/write.py
@@ -63,7 +63,11 @@ def prepare_and_write(dataset, outpath, vars_df=None, meta_dict=None, time='60mi
     d2 = roundValues(d2, vars_df)
 
     # Get variable names to write out
-    col_names = getColNames(vars_df, d2, remove_nan_fields=True)
+    if d2.attrs['level'] == 'L3':
+        remove_nan_fields = True
+    else:
+        remove_nan_fields = False
+    col_names = getColNames(vars_df, d2, remove_nan_fields=remove_nan_fields)
 
     # Define filename based on resample rate
     t = int(pd.Timedelta((d2['time'][1] - d2['time'][0]).values).total_seconds())