Merge pull request #143 from arulalant/master

Fixed TIGGE : OLR negative sign bug

Author: arulalant

bsubScripts/ncumeps_global_post/ncumeps_global_vsdb_anl.bash

@@ -20,7 +20,7 @@ localTable_relative_dir="$DIR/../../tables/local/ncmr/v1/"
 localTable_absolute_dir="$( cd "$localTable_relative_dir" && pwd )"
 localTable=$localTable_absolute_dir/ncmr_grib2_local_table
 
-export UMRIDER_STARTDATE=20160920
+
 export GRIB2TABLE=$localTable
 echo "export GRIB2TABLE="$GRIB2TABLE
 

bsubScripts/ncumeps_global_tigge/ncumeps_global_tigge_2ds1_fcst_00Z.bash

@@ -8,7 +8,6 @@
 #BSUB -x                      # exclusive mode
 #BSUB -R span[ptile=16]       # task per node
 #BSUB -q ensemble             	  # queue
-#BSUB -L /bin/bash
 #BSUB -e /gpfs3/home/umeps/UMRiderLogs/tigge/bsub/um2grb2.tig2d1.err.%J.%I.hybrid     # error file name in which %J is replaced by the job ID
 #BSUB -o /gpfs3/home/umeps/UMRiderLogs/tigge/bsub/um2grb2.tig2d1.out.%J.%I.hybrid     # output file name in which %J is replaced by the job ID
 

bsubScripts/ncumeps_global_tigge/tigge_create_tarball_g2files.py

@@ -14,7 +14,7 @@
 ## Arulalan.T
 ## 04-Mar-2016.
 
-import os, subprocess, datetime, getopt, sys, glob
+import os, subprocess, datetime, getopt, sys, glob, time
 
 pbzip2 = '/gpfs1/home/Libs/GNU/ZIPUTIL/pbzip2'
 pigz = '/gpfs1/home/Libs/GNU/ZIPUTIL/pigz'
@@ -84,15 +84,16 @@ def createTarBalls(path, today, member):
     print "gzip_cmd = ", gzip_cmd
     subprocess.call(gzip_cmd, shell=True)
     time.sleep(5)
-            
-    # remove yesterday directory!!!
-    yDayPath = os.path.join(path, '../../%s' % yDay)
-    if os.path.exists(yDayPath):    
-        cmd = "rm -rf %s" % yDayPath
-        print cmd
-        subprocess.call(cmd, shell=True)
-    # end of if os.path.exists(yDayPath):     
-     
+    print os.getcwd(), member    
+    if member == '000':        
+        # remove today directory!!!
+        print "path",  path
+        if os.path.exists(path):    
+            cmd = "rm -rf %s" % path
+            print cmd
+            subprocess.call(cmd, shell=True)
+         # end of if os.path.exists(yDayPath):     
+    # end of if member == '000': 
     os.chdir(cdir)  
 # end of def createTarBalls(path, today, ...):
 

g2utils/um2grb2tigge.py

@@ -2059,6 +2059,16 @@ def packEnsembles(arg, **kwarg):
             print regdCube.data.min(), regdCube.data.max()
         # end of if varName.endswith('_flux'):
 
+        if (varName, varSTASH) == ('toa_outgoing_longwave_flux', 'm01s02i205'):
+            # https://software.ecmwf.int/wiki/display/TIGGE/Rules+for+data+encoding+and+exchange
+            # as per the tigge statndard, The flux sign convention will be positive downwards.
+            # So, here toa-olr shoule be negative, but ncum model just gives the magnitude. So 
+            # lets fix this by ourself.
+            if regdCube.data.max() > 0 and regdCube.data.min() > 0:
+                # convert to negative if only this data sign is positive.
+                regdCube.data *= -1 # multiply with -1 to indicate this is upward flux.
+        # end of if (varName, varSTASH) == ('toa_outgoing_longwave_flux', 'm01s02i205'):
+        
         if (varName, varSTASH) in _precipVars_:
             # Since we are not using 'mask' option for extrapolate while 
             # doing linear regrid, which bring -ve values after regrid in