Feature/obs error assignment tool (#147)

* added source codes along with input and output example files

* Create README.md

* Update README.md

* Update README.md

* Updated some comments

* Update run_error_model_estimate.sh

* Update README.md

* Applied Linting

* Update README.md

* Update README.md

* Update README.md

* Fixed  pycodestyle violation already in develop

* Modified the ploting part of the python code

* fixed coding norms

* ddressed reviewer comments. added docstring in each function and modified paths

fix/cloudy_radiance_info.txt

@@ -0,0 +1,62 @@
+radiance_mod_instr_input::
+!obsname    obsloc    ex_obserr  ex_biascor  cld_effect
+ gmi        sea       ex_obserr3   .true.      .false.
+ amsua      sea       ex_obserr1   .true.      .true.
+ atms       sea       ex_obserr1   .true.      .true.
+::
+
+obs_amsua::
+!  Parameters for the observation error model
+!  cclr [kg/m2] & ccld [kg/m2]: range of cloud amounts 
+!  over which the main increase in error take place
+!  ch   cclr   ccld
+    1   0.050  0.600
+    2   0.030  0.450
+    3   0.030  0.400
+    4   0.020  0.450
+    5   0.000  1.000
+    6   0.100  1.500
+   15   0.030  0.200
+::        
+
+obs_atms::
+!  Parameters for the observation error model
+!  cclr [kg/m2] & ccld [kg/m2]: range of cloud amounts
+!  over which the main increase in error take place
+!  ch   cclr   ccld
+    1   0.030  0.350
+    2   0.030  0.380
+    3   0.030  0.400
+    4   0.020  0.450
+    5   0.030  0.500
+    6   0.080  1.000
+    7   0.150  1.000
+   16   0.020  0.350
+   17   0.030  0.500
+   18   0.030  0.500
+   19   0.030  0.500
+   20   0.030  0.500
+   21   0.050  0.500
+   22   0.100  0.500
+::
+
+obs_gmi::
+!  Parameters for the observation error model
+!  cclr [kg/m2] & ccld [kg/m2]: range of cloud amounts
+!  over which the main increase in error take place
+!  ch   cclr   ccld  
+    1   0.050  0.200 
+    2   0.050  0.200 
+    3   0.050  0.200 
+    4   0.050  0.200 
+    5   0.050  0.200 
+    6   0.050  0.200 
+    7   0.050  0.200 
+    8   0.050  0.200 
+    9   0.050  0.200 
+   10   0.050  0.300 
+   11   0.050  0.200 
+   12   0.050  0.300 
+   13   0.050  0.300 
+::
+

scripts/plot_proftrace.py

@@ -297,13 +297,13 @@ def collect_statistics(setdict):
             # entries in a list
             # if they are not a list, assert them as a list
             it = sd['it']
-            it = [it] if type(it) != list else it
+            it = [it] if type(it) is not list else it
             use = sd['use']
-            use = [use] if type(use) != list else use
+            use = [use] if type(use) is not list else use
             typ = sd['typ']
-            typ = [typ] if type(typ) != list else typ
+            typ = [typ] if type(typ) is not list else typ
             styp = sd['styp']
-            styp = [styp] if type(styp) != list else styp
+            styp = [styp] if type(styp) is not list else styp
             # extract variable from gdas
             stat = gdas.extract(var)  # t, uv, q, etc.
             #             date         it           obs          use