Post ceilling (NOAA-EMC#136)

* Add the changes of outputing ceiling height for global FV3. * Output low,middle,high cloud fraction; radar reflectivity at model layer 1, 2; 1 and 4 km radar reflectivity. * Add the bug fix of initializing DBZI from Ruiyu. * 1)correct time averaged cloud fraction at low/mid/high cloud layer grib2 name from "TCDC" to "LCDC/MCDC/HCDC"; 2)output mixed layer CAPE/CIN; 3)remove simulated GOES-12 brightness temperature. * Clean up v8.0.0-hera.
lgannoaa · Jul 20, 2020 · e8f0f2e · e8f0f2e
1 parent 8802461
commit e8f0f2e
Show file tree

Hide file tree

Showing 20 changed files with 1,638 additions and 216 deletions.
diff --git a/modulefiles/post/v8.0.0-hera b/modulefiles/post/v8.0.0-hera
@@ -15,12 +15,9 @@ module-whatis "post"
 set ver v8.0.0
 
 # Loading Intel Compiler Suite
-#module load intel/19.0.4.243
-#module load impi/2019.0.4
 module load intel/18.0.5.274
 module load impi/2018.0.4
 
-#module use /contrib/modulefiles
 module use -a /scratch2/NCEPDEV/nwprod/NCEPLIBS/modulefiles
 
 # Loding nceplibs modules
@@ -38,15 +35,12 @@ module load ip/3.0.2
 module load sp/2.0.3
 module load w3nco/2.0.7
 module load crtm/2.3.0
-#module load netcdf/3.6.3
-#module load netcdf/4.7.0
 module load g2tmpl/1.6.0
 module load wrfio/1.1.1
 module load hdf5_parallel/1.10.6
 module load netcdf_parallel/4.7.4
 module load w3emc_para/2.4.0
 
-#setenv WRFPATH /scratch3/NCEPDEV/nwprod/sorc/wrf_shared.v1.1.0
 setenv myFC mpiifort
 setenv OPENMP "-qopenmp"
 setenv myFCFLAGS "-O3 -convert big_endian -traceback -g -fp-model source -qopenmp -fpp"

diff --git a/parm/global_1x1_paramlist_g2 b/parm/global_1x1_paramlist_g2
@@ -48,11 +48,13 @@ APTMP:2 m above ground
 var discipline=0 master_table=2 parmcat=0 parm=21:2 m above ground
 CAPE:180-0 mb above ground
 CAPE:255-0 mb above ground
+CAPE:90-0 mb above ground
 CAPE:surface
 CFRZR:surface
 CICEP:surface
 CIN:180-0 mb above ground
 CIN:255-0 mb above ground
+CIN:90-0 mb above ground
 CIN:surface
 CLWMR:50 mb
 CLWMR:1000 mb
@@ -114,6 +116,7 @@ GRLE:950 mb
 GRLE:975 mb
 GRLE:1 hybrid level
 GUST:surface
+HCDC:high cloud layer
 HGT:0C isotherm
 HGT:1000 mb
 HGT:100 mb
@@ -156,6 +159,7 @@ HGT:0.01 mb
 HGT:0.02 mb
 HGT:0.04 mb
 HGT:0.07 mb
+HGT:cloud ceiling
 HGT:highest tropospheric freezing level
 HGT:max wind
 HGT:PV=-2e-06 (Km^2/kg/s) surface
@@ -213,6 +217,7 @@ ICMR:950 mb
 ICMR:975 mb
 ICMR:1 hybrid level
 LAND:surface
+LCDC:low cloud layer
 LFTX:surface
 LHTFL:surface
 MSLET:mean sea level
@@ -257,6 +262,7 @@ O3MR:0.01 mb
 O3MR:0.02 mb
 O3MR:0.04 mb
 O3MR:0.07 mb
+MCDC:middle cloud layer
 PEVPR:surface
 PLPL:255-0 mb above ground
 POT:0.995 sigma level
@@ -278,6 +284,10 @@ PRES:tropopause
 PRMSL:mean sea level
 PWAT:entire atmosphere
 REFC:entire atmosphere
+REFD:1 hybrid level
+REFD:2 hybrid level
+REFD:1000 m above ground
+REFD:4000 m above ground
 RH:0.33-1 sigma layer
 RH:0.44-0.72 sigma layer
 RH:0.44-1 sigma layer

diff --git a/parm/global_1x1_paramlist_g2.f000 b/parm/global_1x1_paramlist_g2.f000
@@ -45,11 +45,13 @@ APTMP:2 m above ground
 var discipline=0 master_table=2 parmcat=0 parm=21:2 m above ground
 CAPE:180-0 mb above ground
 CAPE:255-0 mb above ground
+CAPE:90-0 mb above ground
 CAPE:surface
 CFRZR:surface
 CICEP:surface
 CIN:180-0 mb above ground
 CIN:255-0 mb above ground
+CIN:90-0 mb above ground
 CIN:surface
 CLWMR:50 mb
 CLWMR:1000 mb
@@ -106,6 +108,7 @@ GRLE:950 mb
 GRLE:975 mb
 GRLE:1 hybrid level
 GUST:surface
+HCDC:high cloud layer
 HGT:0C isotherm
 HGT:1000 mb
 HGT:100 mb
@@ -148,6 +151,7 @@ HGT:0.01 mb
 HGT:0.02 mb
 HGT:0.04 mb
 HGT:0.07 mb
+HGT:cloud ceiling
 HGT:highest tropospheric freezing level
 HGT:max wind
 HGT:PV=-2e-06 (Km^2/kg/s) surface
@@ -186,6 +190,7 @@ ICMR:950 mb
 ICMR:975 mb
 ICMR:1 hybrid level
 LAND:surface
+LCDC:low cloud layer
 LFTX:surface
 MSLET:mean sea level
 O3MR:1000 mb
@@ -229,6 +234,7 @@ O3MR:0.01 mb
 O3MR:0.02 mb
 O3MR:0.04 mb
 O3MR:0.07 mb
+MCDC:middle cloud layer
 PLPL:255-0 mb above ground
 POT:0.995 sigma level
 PRATE:surface
@@ -241,6 +247,10 @@ PRES:tropopause
 PRMSL:mean sea level
 PWAT:entire atmosphere
 REFC:entire atmosphere
+REFD:1 hybrid level
+REFD:2 hybrid level
+REFD:1000 m above ground
+REFD:4000 m above ground
 RH:0.33-1 sigma layer
 RH:0.44-0.72 sigma layer
 RH:0.44-1 sigma layer

diff --git a/parm/post_avblflds.xml b/parm/post_avblflds.xml
@@ -2556,7 +2556,7 @@
        <post_avblfldidx>300</post_avblfldidx>
        <shortname>AVE_TCDC_ON_LOW_CLOUD_LYR</shortname>
        <pdstmpl>tmpl4_8</pdstmpl>	   
-       <pname>TCDC</pname>
+       <pname>LCDC</pname>
        <stats_proc>AVE</stats_proc>
        <fixed_sfc1_type>low_cloud_lyr</fixed_sfc1_type>
        <scale>3.0</scale>
@@ -2566,7 +2566,7 @@
        <post_avblfldidx>301</post_avblfldidx>
        <shortname>AVE_TCDC_ON_MID_CLOUD_LYR</shortname>
        <pdstmpl>tmpl4_8</pdstmpl>	   
-       <pname>TCDC</pname>
+       <pname>MCDC</pname>
        <stats_proc>AVE</stats_proc>
        <fixed_sfc1_type>mid_cloud_lyr</fixed_sfc1_type>
        <scale>3.0</scale>
@@ -2576,7 +2576,7 @@
        <post_avblfldidx>302</post_avblfldidx>
        <shortname>AVE_TCDC_ON_HIGH_CLOUD_LYR</shortname>
        <pdstmpl>tmpl4_8</pdstmpl>	   
-       <pname>TCDC</pname>
+       <pname>HCDC</pname>
        <stats_proc>AVE</stats_proc>
        <fixed_sfc1_type>high_cloud_lyr</fixed_sfc1_type>
        <scale>3.0</scale>

diff --git a/parm/postcntrl_gfs.xml b/parm/postcntrl_gfs.xml
@@ -1090,6 +1090,66 @@
        <scale>5.0</scale>
     </param>
 
+    <param>
+       <shortname>HGT_ON_CLOUD_CEILING</shortname>
+       <scale>6.0</scale>
+    </param>
+
+    <param>
+       <shortname>INST_TCDC_ON_ENTIRE_ATMOS</shortname>
+       <scale>3.0</scale>
+    </param>
+
+    <param>
+       <shortname>LCDC_ON_LOW_CLOUD_LYR</shortname>
+       <pname>LCDC</pname>
+       <scale>3.0</scale>
+    </param>
+
+    <param>
+       <shortname>MCDC_ON_MID_CLOUD_LYR</shortname>
+       <pname>MCDC</pname>
+       <scale>3.0</scale>
+    </param>
+
+    <param>
+       <shortname>HCDC_ON_HIGH_CLOUD_LYR</shortname>
+       <pname>HCDC</pname>
+       <scale>3.0</scale>
+    </param>
+
+    <param>
+       <shortname>REFD_ON_SPEC_HGT_LVL_ABOVE_GRND</shortname>
+       <pname>REFD</pname>
+       <table_info>NCEP</table_info>
+       <level>4000. 1000.</level>
+       <scale>4.0</scale>
+    </param>
+
+    <param>
+       <shortname>REFD_ON_HYBRID_LVL</shortname>
+       <pname>REFD</pname>
+       <table_info>NCEP</table_info>
+       <level>1. 2.</level>
+       <scale>4.0</scale>
+    </param>
+
+    <param>
+        <shortname>MIXED_LAYER_CAPE_ON_SPEC_PRES_ABOVE_GRND</shortname>
+        <pname>CAPE</pname>
+        <level>9000.</level>
+        <level2>0.</level2>
+        <scale>4.0</scale>
+    </param>
+
+    <param>
+        <shortname>MIXED_LAYER_CIN_ON_SPEC_PRES_ABOVE_GRND</shortname>
+        <pname>CIN</pname>
+        <level>9000.</level>
+        <level2>0.</level2>
+        <scale>4.0</scale>
+    </param>
+
     <!--  comments wherever I like -->
   </paramset>
 </postxml>
diff --git a/parm/postcntrl_gfs_f00.xml b/parm/postcntrl_gfs_f00.xml
@@ -838,6 +838,66 @@
        <scale>5.0</scale>
     </param>
 
+    <param>
+       <shortname>HGT_ON_CLOUD_CEILING</shortname>
+       <scale>6.0</scale>
+    </param>
+
+    <param>
+       <shortname>INST_TCDC_ON_ENTIRE_ATMOS</shortname>
+       <scale>3.0</scale>
+    </param>
+
+    <param>
+       <shortname>LCDC_ON_LOW_CLOUD_LYR</shortname>
+       <pname>LCDC</pname>
+       <scale>3.0</scale>
+    </param>
+
+    <param>
+       <shortname>MCDC_ON_MID_CLOUD_LYR</shortname>
+       <pname>MCDC</pname>
+       <scale>3.0</scale>
+    </param>
+
+    <param>
+       <shortname>HCDC_ON_HIGH_CLOUD_LYR</shortname>
+       <pname>HCDC</pname>
+       <scale>3.0</scale>
+    </param>
+
+    <param>
+       <shortname>REFD_ON_SPEC_HGT_LVL_ABOVE_GRND</shortname>
+       <pname>REFD</pname>
+       <table_info>NCEP</table_info>
+       <level>4000. 1000.</level>
+       <scale>4.0</scale>
+    </param>
+
+    <param>
+       <shortname>REFD_ON_HYBRID_LVL</shortname>
+       <pname>REFD</pname>
+       <table_info>NCEP</table_info>
+       <level>1. 2.</level>
+       <scale>4.0</scale>
+    </param>
+
+    <param>
+        <shortname>MIXED_LAYER_CAPE_ON_SPEC_PRES_ABOVE_GRND</shortname>
+        <pname>CAPE</pname>
+        <level>9000.</level>
+        <level2>0.</level2>
+        <scale>4.0</scale>
+    </param>
+
+    <param>
+        <shortname>MIXED_LAYER_CIN_ON_SPEC_PRES_ABOVE_GRND</shortname>
+        <pname>CIN</pname>
+        <level>9000.</level>
+        <level2>0.</level2>
+        <scale>4.0</scale>
+    </param>
+
     <!--  comments wherever I like -->
   </paramset>
 </postxml>
diff --git a/parm/postcntrl_gfs_f00_two.xml b/parm/postcntrl_gfs_f00_two.xml
@@ -838,6 +838,66 @@
        <scale>5.0</scale>
     </param>
 
+    <param>
+       <shortname>HGT_ON_CLOUD_CEILING</shortname>
+       <scale>6.0</scale>
+    </param>
+
+    <param>
+       <shortname>INST_TCDC_ON_ENTIRE_ATMOS</shortname>
+       <scale>3.0</scale>
+    </param>
+
+    <param>
+       <shortname>LCDC_ON_LOW_CLOUD_LYR</shortname>
+       <pname>LCDC</pname>
+       <scale>3.0</scale>
+    </param>
+
+    <param>
+       <shortname>MCDC_ON_MID_CLOUD_LYR</shortname>
+       <pname>MCDC</pname>
+       <scale>3.0</scale>
+    </param>
+
+    <param>
+       <shortname>HCDC_ON_HIGH_CLOUD_LYR</shortname>
+       <pname>HCDC</pname>
+       <scale>3.0</scale>
+    </param>
+
+    <param>
+       <shortname>REFD_ON_SPEC_HGT_LVL_ABOVE_GRND</shortname>
+       <pname>REFD</pname>
+       <table_info>NCEP</table_info>
+       <level>4000. 1000.</level>
+       <scale>4.0</scale>
+    </param>
+
+    <param>
+       <shortname>REFD_ON_HYBRID_LVL</shortname>
+       <pname>REFD</pname>
+       <table_info>NCEP</table_info>
+       <level>1. 2.</level>
+       <scale>4.0</scale>
+    </param>
+
+    <param>
+        <shortname>MIXED_LAYER_CAPE_ON_SPEC_PRES_ABOVE_GRND</shortname>
+        <pname>CAPE</pname>
+        <level>9000.</level>
+        <level2>0.</level2>
+        <scale>4.0</scale>
+    </param>
+
+    <param>
+        <shortname>MIXED_LAYER_CIN_ON_SPEC_PRES_ABOVE_GRND</shortname>
+        <pname>CIN</pname>
+        <level>9000.</level>
+        <level2>0.</level2>
+        <scale>4.0</scale>
+    </param>
+
     <!--  comments wherever I like -->
   </paramset>
 

diff --git a/parm/postcntrl_gfs_goes.xml b/parm/postcntrl_gfs_goes.xml
@@ -19,26 +19,6 @@
     <field_datatype>fltng_pnt</field_datatype>
     <comprs_type>lossless</comprs_type>
 
-    <param>
-       <shortname>SBT122_ON_TOP_OF_ATMOS</shortname>
-       <scale>4.0</scale>
-    </param>
-
-    <param>
-       <shortname>SBT123_ON_TOP_OF_ATMOS</shortname>
-       <scale>4.0</scale>
-    </param>
-
-    <param>
-       <shortname>SBT124_ON_TOP_OF_ATMOS</shortname>
-       <scale>4.0</scale>
-    </param>
-
-    <param>
-       <shortname>SBT126_ON_TOP_OF_ATMOS</shortname>
-       <scale>4.0</scale>
-    </param>
-
     <param>
        <shortname>SBTAGR8_ON_TOP_OF_ATMOS</shortname>
        <scale>4.0</scale>