Add lat_lon_native diagnostic set for native grid visualization

Author: chengzhuzhang

analysis_data_preprocess/create_AOD550_climo.py

@@ -209,15 +209,15 @@
 ncatted -O -h -a cell_methods,AOD_550,o,c,"time: mean within years time: mean over years" $OBSFILE $OBSFILE
 ncatted -O -h -a cell_methods,AOD_550_ann,o,c,"time: mean over years" $OBSFILE $OBSFILE
 # Climatology files should contain a "season" attribute, which should
-# be either a three-letter season (e.g., ANN, DJF, MAR) or a two-digit
+# be either a three-letter season (e.g., ANN, DJF, MAR) or a two-digi
 # month (e.g. 01, 02, ... 12)/mo
 ncatted -O -h -a season,global,o,c,$SEASON $OBSFILE
 
 # Variable cleanup
 ncatted -O -a standard_name,lat,o,c,"latitude" $OBSFILE $OBSFILE
 ncatted -O -a standard_name,lon,o,c,"longitude" $OBSFILE $OBSFILE
 ncatted -O -a long_name,time,o,c,"time" -a standard_name,time,o,c,"time" $OBSFILE $OBSFILE
-# Flip FLNS/FLNSC signs to match model output
+# Flip FLNS/FLNSC signs to match model outpu
 #
 # Note: Simply reversing valid_min/valid_max was
 # throwing an error in ncview. By doing the below

analysis_data_preprocess/create_GPCP_climo.py

@@ -1,5 +1,5 @@
 # flake8: noqa
-#!/usr/bin/env cdat
+#!/usr/bin/env cda
 """
 This script takes an observational dataset that contains monthly data
 and creates the appropriate climo files. Aspects of this code have been taken
@@ -24,13 +24,13 @@
 Modified by Jill Zhang ([email protected]) 08-15-2018
 """
 
-# Upload the following libraries to call in this script
+# Upload the following libraries to call in this scrip
 import argparse
 import code
 import datetime
 import gc
 import glob
-import inspect
+import inspec
 import os
 import re
 import string
@@ -67,16 +67,16 @@
 output_hostpath="/p/user_pub/e3sm/zhang40/analysis_data_e3sm_diags/"+data_name+"/climatology/" #location to save climatology files
 if data_name =='GPCP':
     rootname='GPCP_v2.3' #name of the climo files
-    startyear='1979'  #start year of the dataset
+    startyear='1979'  #start year of the datase
     startmonth='01'   #start month of the dataset, typically 01
-    endyear='2017'    #end year of the dataset
+    endyear='2017'    #end year of the datase
     endmonth='12'     #end month of the dataset, typically 12
 
 if data_name =='GPCP_v2.2':
     rootname='GPCP_v2.2' #name of the climo files
-    startyear='1979'  #start year of the dataset
+    startyear='1979'  #start year of the datase
     startmonth='01'   #start month of the dataset, typically 01
-    endyear='2014'    #end year of the dataset
+    endyear='2014'    #end year of the datase
     endmonth='12'     #end month of the dataset, typically 12
 
 #^^^^^^^^^^^^^^^  ^^^^^^^^^^^^^^^ ^^^^^^^^^^^^^^^ ^^^^^^^^^^^^^^^  ^^^^^^^^^^^^^^^^^
@@ -150,7 +150,7 @@
 f_in=cdm.open(input_hostANDfilename)
 varlist=f_in.listvariables()
 
-print varlist
+print varlis
 print f_in.listdimension()
 
 ignore_variables=['bounds_lat','bounds_time','bounds_lon','bnds_lat','bnds_time','bnds_lon','date','lon_bnds','lat_bnds','time_bnds']  #variables to ignore when making climatology means
@@ -176,7 +176,7 @@
             outputdattable=eval(".".join(['cdu',fi,'climatology(dattable, criteriaarg = [0.99,None])'])) #only full DJF ranges are included in the calculation
             outputdattable.cell_methods="time: mean over years"
 
-            #specify the time axis long_name and the units - uses the startyear-1-1 as the starting point
+            #specify the time axis long_name and the units - uses the startyear-1-1 as the starting poin
             time_axis=outputdattable.getAxis(0)
             time_axis.climatology="climatology_bnds"
             time_axis.long_name="time"

analysis_data_preprocess/create_OAFLUX_climo.py

@@ -1,5 +1,5 @@
 # flake8: noqa
-#!/usr/bin/env cdat
+#!/usr/bin/env cda
 """
 This script takes an observational dataset that contains monthly data
 and creates the appropriate climo files. Aspects of this code have been taken
@@ -25,13 +25,13 @@
 Modified by Jill Zhang ([email protected]) 08-27-2018
 """
 
-# Upload the following libraries to call in this script
+# Upload the following libraries to call in this scrip
 import argparse
 import code
 import datetime
 import gc
 import glob
-import inspect
+import inspec
 import os
 import re
 import string
@@ -80,9 +80,9 @@
 )  # string of location+filename of input file
 
 rootname = "OAFlux"  # name of the climo files
-startyear = "1980"  # start year of the dataset
+startyear = "1980"  # start year of the datase
 startmonth = "01"  # start month of the dataset, typically 01
-endyear = "2005"  # end year of the dataset
+endyear = "2005"  # end year of the datase
 endmonth = "12"  # end month of the dataset, typically 12
 # ^^^^^^^^^^^^^^^  ^^^^^^^^^^^^^^^ ^^^^^^^^^^^^^^^ ^^^^^^^^^^^^^^^  ^^^^^^^^^^^^^^^^^
 
@@ -250,7 +250,7 @@
             )  # only full DJF ranges are included in the calculation
             outputdattable.cell_methods = "time: mean over years"
 
-            # specify the time axis long_name and the units - uses the startyear-1-1 as the starting point
+            # specify the time axis long_name and the units - uses the startyear-1-1 as the starting poin
             time_axis = outputdattable.getAxis(0)
             time_axis.climatology = "climatology_bnds"
             time_axis.long_name = "time"

analysis_data_preprocess/create_PminusE_GPCP-OAFlux_timeseries.py

@@ -21,8 +21,8 @@
 fin1 = cdms2.open(input_path1 + filename1)
 fin2 = cdms2.open(input_path2 + filename2)
 
-startyear = "1979"  # start year of the dataset
-endyear = "2013"  # end year of the dataset
+startyear = "1979"  # start year of the datase
+endyear = "2013"  # end year of the datase
 
 prect = fin1("PRECT", time=("1979-01-15", "2013-12-15", "ccb"))
 qflux = fin2("QFLX", time=("1979-01-15", "2013-12-15", "ccb"))

analysis_data_preprocess/fix_units_1.5hr_time_shift_forTRMM.py

@@ -25,7 +25,7 @@
     # The shift may be due to the assumption that ncclimo makes that times are encoded using the CESM/E3SM convention where times refer to the end of the timestep
     print('var.time',var.getTime()[:])
     for ilat in range(nlat):
-        print,'ilat= ',ilat
+        print,'ilat= ',ila
         var[0,ilat,:] = var0[7,ilat,:]
         var[1:8,ilat,:] = var0[0:7,ilat,:]
     newtime=cdms2.createAxis(MV2.zeros(8))

analysis_data_preprocess/process_GPCP_dataset.py

@@ -1,5 +1,5 @@
 # flake8: noqa
-#!/usr/bin/env cdat
+#!/usr/bin/env cda
 """
 This script takes the GPCP dataset with monthly data on it from
 "http://www.esrl.noaa.gov/psd/data/gridded/data.gpcp.html
@@ -15,7 +15,7 @@
 import datetime
 import gc
 import glob
-import inspect
+import inspec
 import os
 import re
 import string
@@ -70,7 +70,7 @@
 
 
 filecount = 0
-fisc=varlist
+fisc=varlis
 f_index=0
 for fi in fisc:
     print "##########################"
@@ -129,7 +129,7 @@
 
     f_out.write(dattable) ;
     print "".join(["** Finished processing: ",fi," **"])
-    filecount = filecount + 1; filecount_s = '%06d' % filecount
+    filecount = filecount + 1; filecount_s = '%06d' % filecoun
     f_index=f_index + 1
 
 f_out.close()

analysis_data_preprocess/process_OAFLUX_dataset.py

@@ -1,5 +1,5 @@
 # flake8: noqa
-#!/usr/bin/env cdat
+#!/usr/bin/env cda
 """
 Based on Chris Terai's original script for GPCP
 Modified for OAFlux by Jill Zhang ([email protected]) 2018-09-15
@@ -9,7 +9,7 @@
 import datetime
 import gc
 import glob
-import inspect
+import inspec
 import os
 import re
 import string

analysis_data_preprocess/split_vars.py

@@ -24,7 +24,7 @@
     years = "200101_201512"
 datapath = (
     "/p/user_pub/e3sm/zhang40/analysis_data_e3sm_diags/CERES-EBAF/"
-    + dataset
+    + datase
     + "/time_series/"
 )
 filename = datapath + data + "_" + years + ".nc"

auxiliary_tools/aerosol_budget.py

@@ -1,4 +1,4 @@
-# NOTE: This module uses the deprecated e3sm_diags.driver.utils.dataset.Dataset
+# NOTE: This module uses the deprecated e3sm_diags.driver.utils.dataset.Datase
 # class, which was replaced by e3sm_diags.driver.utils.dataset_xr.Dataset.
 import csv
 import json
@@ -78,7 +78,7 @@ def generate_metrics_dic(data, aerosol, season):
         "Burden (Tg)": f"{burden_total:.3f}",
         "Lifetime (Days)": f"{burden_total/sink*365:.3f}",
     }
-    return metrics_dict
+    return metrics_dic
 
 
 param = CoreParameter()

auxiliary_tools/cdat_regression_testing/562-index-html/run_script.py

@@ -1,5 +1,5 @@
-# python -m auxiliary_tools.cdat_regression_testing.792-lat-lon-run-script.792_lat_lon_run_script
-from auxiliary_tools.cdat_regression_testing.base_run_script import run_set
+# python -m auxiliary_tools.cdat_regression_testing.792-lat-lon-run-script.792_lat_lon_run_scrip
+from auxiliary_tools.cdat_regression_testing.base_run_script import run_se
 
 SET_NAME = "polar"
 SET_DIR = "562-index-html"