From e8f0f2ea1edde6aea3112ee3088fb4aca3f0dce5 Mon Sep 17 00:00:00 2001 From: WenMeng-NOAA <48260754+WenMeng-NOAA@users.noreply.github.com> Date: Fri, 12 Jun 2020 15:22:20 -0400 Subject: [PATCH] Post ceilling (#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. --- modulefiles/post/v8.0.0-hera | 6 - parm/global_1x1_paramlist_g2 | 10 + parm/global_1x1_paramlist_g2.f000 | 10 + parm/post_avblflds.xml | 6 +- parm/postcntrl_gfs.xml | 60 ++++ parm/postcntrl_gfs_f00.xml | 60 ++++ parm/postcntrl_gfs_f00_two.xml | 60 ++++ parm/postcntrl_gfs_goes.xml | 20 -- parm/postcntrl_gfs_two.xml | 60 ++++ parm/postxconfig-NT-GEFS.txt | 6 +- parm/postxconfig-NT-GFS-F00-TWO.txt | 335 +++++++++++++++++- parm/postxconfig-NT-GFS-F00.txt | 335 +++++++++++++++++- parm/postxconfig-NT-GFS-FLUX.txt | 6 +- parm/postxconfig-NT-GFS-GOES.txt | 150 +------- parm/postxconfig-NT-GFS-TWO.txt | 347 ++++++++++++++++++- parm/postxconfig-NT-GFS.txt | 341 +++++++++++++++++- sorc/ncep_post.fd/CLDRAD.f | 18 +- sorc/ncep_post.fd/INITPOST_GFS_NETCDF.f | 10 +- sorc/ncep_post.fd/INITPOST_GFS_NETCDF_PARA.f | 10 +- sorc/ncep_post.fd/MDLFLD.f | 4 +- 20 files changed, 1638 insertions(+), 216 deletions(-) diff --git a/modulefiles/post/v8.0.0-hera b/modulefiles/post/v8.0.0-hera index a46e9e9976..265c55ea2e 100644 --- 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 index 25f4e136b0..52eafa7f0e 100644 --- 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 index 57dfdfa24e..cb5beccdce 100644 --- 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 index de6f72922c..0d931adaf8 100755 --- a/parm/post_avblflds.xml +++ b/parm/post_avblflds.xml @@ -2556,7 +2556,7 @@ 300 AVE_TCDC_ON_LOW_CLOUD_LYR tmpl4_8 - TCDC + LCDC AVE low_cloud_lyr 3.0 @@ -2566,7 +2566,7 @@ 301 AVE_TCDC_ON_MID_CLOUD_LYR tmpl4_8 - TCDC + MCDC AVE mid_cloud_lyr 3.0 @@ -2576,7 +2576,7 @@ 302 AVE_TCDC_ON_HIGH_CLOUD_LYR tmpl4_8 - TCDC + HCDC AVE high_cloud_lyr 3.0 diff --git a/parm/postcntrl_gfs.xml b/parm/postcntrl_gfs.xml index 1ee82d8aa7..3b6c3b675e 100755 --- a/parm/postcntrl_gfs.xml +++ b/parm/postcntrl_gfs.xml @@ -1090,6 +1090,66 @@ 5.0 + + HGT_ON_CLOUD_CEILING + 6.0 + + + + INST_TCDC_ON_ENTIRE_ATMOS + 3.0 + + + + LCDC_ON_LOW_CLOUD_LYR + LCDC + 3.0 + + + + MCDC_ON_MID_CLOUD_LYR + MCDC + 3.0 + + + + HCDC_ON_HIGH_CLOUD_LYR + HCDC + 3.0 + + + + REFD_ON_SPEC_HGT_LVL_ABOVE_GRND + REFD + NCEP + 4000. 1000. + 4.0 + + + + REFD_ON_HYBRID_LVL + REFD + NCEP + 1. 2. + 4.0 + + + + MIXED_LAYER_CAPE_ON_SPEC_PRES_ABOVE_GRND + CAPE + 9000. + 0. + 4.0 + + + + MIXED_LAYER_CIN_ON_SPEC_PRES_ABOVE_GRND + CIN + 9000. + 0. + 4.0 + + diff --git a/parm/postcntrl_gfs_f00.xml b/parm/postcntrl_gfs_f00.xml index 7093006d83..6d73b0c75c 100755 --- a/parm/postcntrl_gfs_f00.xml +++ b/parm/postcntrl_gfs_f00.xml @@ -838,6 +838,66 @@ 5.0 + + HGT_ON_CLOUD_CEILING + 6.0 + + + + INST_TCDC_ON_ENTIRE_ATMOS + 3.0 + + + + LCDC_ON_LOW_CLOUD_LYR + LCDC + 3.0 + + + + MCDC_ON_MID_CLOUD_LYR + MCDC + 3.0 + + + + HCDC_ON_HIGH_CLOUD_LYR + HCDC + 3.0 + + + + REFD_ON_SPEC_HGT_LVL_ABOVE_GRND + REFD + NCEP + 4000. 1000. + 4.0 + + + + REFD_ON_HYBRID_LVL + REFD + NCEP + 1. 2. + 4.0 + + + + MIXED_LAYER_CAPE_ON_SPEC_PRES_ABOVE_GRND + CAPE + 9000. + 0. + 4.0 + + + + MIXED_LAYER_CIN_ON_SPEC_PRES_ABOVE_GRND + CIN + 9000. + 0. + 4.0 + + diff --git a/parm/postcntrl_gfs_f00_two.xml b/parm/postcntrl_gfs_f00_two.xml index fc9c06e009..497d1582ba 100755 --- a/parm/postcntrl_gfs_f00_two.xml +++ b/parm/postcntrl_gfs_f00_two.xml @@ -838,6 +838,66 @@ 5.0 + + HGT_ON_CLOUD_CEILING + 6.0 + + + + INST_TCDC_ON_ENTIRE_ATMOS + 3.0 + + + + LCDC_ON_LOW_CLOUD_LYR + LCDC + 3.0 + + + + MCDC_ON_MID_CLOUD_LYR + MCDC + 3.0 + + + + HCDC_ON_HIGH_CLOUD_LYR + HCDC + 3.0 + + + + REFD_ON_SPEC_HGT_LVL_ABOVE_GRND + REFD + NCEP + 4000. 1000. + 4.0 + + + + REFD_ON_HYBRID_LVL + REFD + NCEP + 1. 2. + 4.0 + + + + MIXED_LAYER_CAPE_ON_SPEC_PRES_ABOVE_GRND + CAPE + 9000. + 0. + 4.0 + + + + MIXED_LAYER_CIN_ON_SPEC_PRES_ABOVE_GRND + CIN + 9000. + 0. + 4.0 + + diff --git a/parm/postcntrl_gfs_goes.xml b/parm/postcntrl_gfs_goes.xml index 512a1e751d..16d5bbf9ec 100755 --- a/parm/postcntrl_gfs_goes.xml +++ b/parm/postcntrl_gfs_goes.xml @@ -19,26 +19,6 @@ fltng_pnt lossless - - SBT122_ON_TOP_OF_ATMOS - 4.0 - - - - SBT123_ON_TOP_OF_ATMOS - 4.0 - - - - SBT124_ON_TOP_OF_ATMOS - 4.0 - - - - SBT126_ON_TOP_OF_ATMOS - 4.0 - - SBTAGR8_ON_TOP_OF_ATMOS 4.0 diff --git a/parm/postcntrl_gfs_two.xml b/parm/postcntrl_gfs_two.xml index 4be9239d28..f42bb5b9f9 100755 --- a/parm/postcntrl_gfs_two.xml +++ b/parm/postcntrl_gfs_two.xml @@ -1090,6 +1090,66 @@ 5.0 + + HGT_ON_CLOUD_CEILING + 6.0 + + + + INST_TCDC_ON_ENTIRE_ATMOS + 3.0 + + + + LCDC_ON_LOW_CLOUD_LYR + LCDC + 3.0 + + + + MCDC_ON_MID_CLOUD_LYR + MCDC + 3.0 + + + + HCDC_ON_HIGH_CLOUD_LYR + HCDC + 3.0 + + + + REFD_ON_SPEC_HGT_LVL_ABOVE_GRND + REFD + NCEP + 4000. 1000. + 4.0 + + + + REFD_ON_HYBRID_LVL + REFD + NCEP + 1. 2. + 4.0 + + + + MIXED_LAYER_CAPE_ON_SPEC_PRES_ABOVE_GRND + CAPE + 9000. + 0. + 4.0 + + + + MIXED_LAYER_CIN_ON_SPEC_PRES_ABOVE_GRND + CIN + 9000. + 0. + 4.0 + + diff --git a/parm/postxconfig-NT-GEFS.txt b/parm/postxconfig-NT-GEFS.txt index 21e1f6f11b..553c254884 100644 --- a/parm/postxconfig-NT-GEFS.txt +++ b/parm/postxconfig-NT-GEFS.txt @@ -1983,7 +1983,7 @@ AVE_TCDC_ON_LOW_CLOUD_LYR ? 1 tmpl4_8 -TCDC +LCDC ? AVE low_cloud_lyr @@ -2020,7 +2020,7 @@ AVE_TCDC_ON_MID_CLOUD_LYR ? 1 tmpl4_8 -TCDC +MCDC ? AVE mid_cloud_lyr @@ -2057,7 +2057,7 @@ AVE_TCDC_ON_HIGH_CLOUD_LYR ? 1 tmpl4_8 -TCDC +HCDC ? AVE high_cloud_lyr diff --git a/parm/postxconfig-NT-GFS-F00-TWO.txt b/parm/postxconfig-NT-GFS-F00-TWO.txt index ade3b21ece..f0b1ea1fe9 100644 --- a/parm/postxconfig-NT-GFS-F00-TWO.txt +++ b/parm/postxconfig-NT-GFS-F00-TWO.txt @@ -1,6 +1,6 @@ 2 45 -140 +149 GFSPRS 0 ncep_nco @@ -5197,6 +5197,339 @@ surface ? ? ? +260 +HGT_ON_CLOUD_CEILING +? +1 +tmpl4_0 +HGT +? +? +cloud_ceilng +0 +? +0 +? +? +0 +? +0 +? +? +? +0 +0.0 +0 +0.0 +? +0 +0.0 +0 +0.0 +1 +6.0 +0 +0 +0 +? +? +? +161 +INST_TCDC_ON_ENTIRE_ATMOS +? +1 +tmpl4_0 +TCDC +? +? +entire_atmos +0 +? +0 +? +? +0 +? +0 +? +? +? +0 +0.0 +0 +0.0 +? +0 +0.0 +0 +0.0 +1 +3.0 +0 +0 +0 +? +? +? +37 +LCDC_ON_LOW_CLOUD_LYR +? +1 +tmpl4_0 +LCDC +? +? +low_cloud_lyr +0 +? +0 +? +? +0 +? +0 +? +? +? +0 +0.0 +0 +0.0 +? +0 +0.0 +0 +0.0 +1 +3.0 +0 +0 +0 +? +? +? +38 +MCDC_ON_MID_CLOUD_LYR +? +1 +tmpl4_0 +MCDC +? +? +mid_cloud_lyr +0 +? +0 +? +? +0 +? +0 +? +? +? +0 +0.0 +0 +0.0 +? +0 +0.0 +0 +0.0 +1 +3.0 +0 +0 +0 +? +? +? +39 +HCDC_ON_HIGH_CLOUD_LYR +? +1 +tmpl4_0 +HCDC +? +? +high_cloud_lyr +0 +? +0 +? +? +0 +? +0 +? +? +? +0 +0.0 +0 +0.0 +? +0 +0.0 +0 +0.0 +1 +3.0 +0 +0 +0 +? +? +? +253 +REFD_ON_SPEC_HGT_LVL_ABOVE_GRND +? +1 +tmpl4_0 +REFD +NCEP +? +spec_hgt_lvl_above_grnd +0 +? +2 +4000. 1000. +? +0 +? +0 +? +? +? +0 +0.0 +0 +0.0 +? +0 +0.0 +0 +0.0 +1 +4.0 +0 +0 +0 +? +? +? +250 +REFD_ON_HYBRID_LVL +? +1 +tmpl4_0 +REFD +NCEP +? +hybrid_lvl +0 +? +2 +1. 2. +? +0 +? +0 +? +? +? +0 +0.0 +0 +0.0 +? +0 +0.0 +0 +0.0 +1 +4.0 +0 +0 +0 +? +? +? +582 +MIXED_LAYER_CAPE_ON_SPEC_PRES_ABOVE_GRND +? +1 +tmpl4_0 +CAPE +? +? +spec_pres_above_grnd +0 +? +1 +9000. +spec_pres_above_grnd +0 +? +1 +0. +? +? +0 +0.0 +0 +0.0 +? +0 +0.0 +0 +0.0 +1 +4.0 +0 +0 +0 +? +? +? +583 +MIXED_LAYER_CIN_ON_SPEC_PRES_ABOVE_GRND +? +1 +tmpl4_0 +CIN +? +? +spec_pres_above_grnd +0 +? +1 +9000. +spec_pres_above_grnd +0 +? +1 +0. +? +? +0 +0.0 +0 +0.0 +? +0 +0.0 +0 +0.0 +1 +4.0 +0 +0 +0 +? +? +? GFSFLX 0 ncep_nco diff --git a/parm/postxconfig-NT-GFS-F00.txt b/parm/postxconfig-NT-GFS-F00.txt index eb9f711312..f20f9435d5 100644 --- a/parm/postxconfig-NT-GFS-F00.txt +++ b/parm/postxconfig-NT-GFS-F00.txt @@ -1,5 +1,5 @@ 1 -140 +149 GFSPRS 0 ncep_nco @@ -5196,3 +5196,336 @@ surface ? ? ? +260 +HGT_ON_CLOUD_CEILING +? +1 +tmpl4_0 +HGT +? +? +cloud_ceilng +0 +? +0 +? +? +0 +? +0 +? +? +? +0 +0.0 +0 +0.0 +? +0 +0.0 +0 +0.0 +1 +6.0 +0 +0 +0 +? +? +? +161 +INST_TCDC_ON_ENTIRE_ATMOS +? +1 +tmpl4_0 +TCDC +? +? +entire_atmos +0 +? +0 +? +? +0 +? +0 +? +? +? +0 +0.0 +0 +0.0 +? +0 +0.0 +0 +0.0 +1 +3.0 +0 +0 +0 +? +? +? +37 +LCDC_ON_LOW_CLOUD_LYR +? +1 +tmpl4_0 +LCDC +? +? +low_cloud_lyr +0 +? +0 +? +? +0 +? +0 +? +? +? +0 +0.0 +0 +0.0 +? +0 +0.0 +0 +0.0 +1 +3.0 +0 +0 +0 +? +? +? +38 +MCDC_ON_MID_CLOUD_LYR +? +1 +tmpl4_0 +MCDC +? +? +mid_cloud_lyr +0 +? +0 +? +? +0 +? +0 +? +? +? +0 +0.0 +0 +0.0 +? +0 +0.0 +0 +0.0 +1 +3.0 +0 +0 +0 +? +? +? +39 +HCDC_ON_HIGH_CLOUD_LYR +? +1 +tmpl4_0 +HCDC +? +? +high_cloud_lyr +0 +? +0 +? +? +0 +? +0 +? +? +? +0 +0.0 +0 +0.0 +? +0 +0.0 +0 +0.0 +1 +3.0 +0 +0 +0 +? +? +? +253 +REFD_ON_SPEC_HGT_LVL_ABOVE_GRND +? +1 +tmpl4_0 +REFD +NCEP +? +spec_hgt_lvl_above_grnd +0 +? +2 +4000. 1000. +? +0 +? +0 +? +? +? +0 +0.0 +0 +0.0 +? +0 +0.0 +0 +0.0 +1 +4.0 +0 +0 +0 +? +? +? +250 +REFD_ON_HYBRID_LVL +? +1 +tmpl4_0 +REFD +NCEP +? +hybrid_lvl +0 +? +2 +1. 2. +? +0 +? +0 +? +? +? +0 +0.0 +0 +0.0 +? +0 +0.0 +0 +0.0 +1 +4.0 +0 +0 +0 +? +? +? +582 +MIXED_LAYER_CAPE_ON_SPEC_PRES_ABOVE_GRND +? +1 +tmpl4_0 +CAPE +? +? +spec_pres_above_grnd +0 +? +1 +9000. +spec_pres_above_grnd +0 +? +1 +0. +? +? +0 +0.0 +0 +0.0 +? +0 +0.0 +0 +0.0 +1 +4.0 +0 +0 +0 +? +? +? +583 +MIXED_LAYER_CIN_ON_SPEC_PRES_ABOVE_GRND +? +1 +tmpl4_0 +CIN +? +? +spec_pres_above_grnd +0 +? +1 +9000. +spec_pres_above_grnd +0 +? +1 +0. +? +? +0 +0.0 +0 +0.0 +? +0 +0.0 +0 +0.0 +1 +4.0 +0 +0 +0 +? +? +? diff --git a/parm/postxconfig-NT-GFS-FLUX.txt b/parm/postxconfig-NT-GFS-FLUX.txt index 07ebd6025a..8bef03763b 100644 --- a/parm/postxconfig-NT-GFS-FLUX.txt +++ b/parm/postxconfig-NT-GFS-FLUX.txt @@ -798,7 +798,7 @@ AVE_TCDC_ON_LOW_CLOUD_LYR ? 1 tmpl4_8 -TCDC +LCDC ? AVE low_cloud_lyr @@ -835,7 +835,7 @@ AVE_TCDC_ON_MID_CLOUD_LYR ? 1 tmpl4_8 -TCDC +MCDC ? AVE mid_cloud_lyr @@ -872,7 +872,7 @@ AVE_TCDC_ON_HIGH_CLOUD_LYR ? 1 tmpl4_8 -TCDC +HCDC ? AVE high_cloud_lyr diff --git a/parm/postxconfig-NT-GFS-GOES.txt b/parm/postxconfig-NT-GFS-GOES.txt index e23e1f5d9a..510347582f 100644 --- a/parm/postxconfig-NT-GFS-GOES.txt +++ b/parm/postxconfig-NT-GFS-GOES.txt @@ -1,5 +1,5 @@ 1 -8 +4 GFSPRS 0 ncep_nco @@ -16,154 +16,6 @@ complex_packing_spatial_diff 2nd_ord_sptdiff fltng_pnt lossless -327 -SBT122_ON_TOP_OF_ATMOS -? -1 -tmpl4_0 -SBT122 -? -? -top_of_atmos -0 -? -0 -? -? -0 -? -0 -? -? -? -0 -0.0 -0 -0.0 -? -0 -0.0 -0 -0.0 -1 -4.0 -0 -0 -0 -? -? -? -328 -SBT123_ON_TOP_OF_ATMOS -? -1 -tmpl4_0 -SBT123 -? -? -top_of_atmos -0 -? -0 -? -? -0 -? -0 -? -? -? -0 -0.0 -0 -0.0 -? -0 -0.0 -0 -0.0 -1 -4.0 -0 -0 -0 -? -? -? -329 -SBT124_ON_TOP_OF_ATMOS -? -1 -tmpl4_0 -SBT124 -? -? -top_of_atmos -0 -? -0 -? -? -0 -? -0 -? -? -? -0 -0.0 -0 -0.0 -? -0 -0.0 -0 -0.0 -1 -4.0 -0 -0 -0 -? -? -? -330 -SBT126_ON_TOP_OF_ATMOS -? -1 -tmpl4_0 -SBT126 -? -? -top_of_atmos -0 -? -0 -? -? -0 -? -0 -? -? -? -0 -0.0 -0 -0.0 -? -0 -0.0 -0 -0.0 -1 -4.0 -0 -0 -0 -? -? -? 959 SBTAGR8_ON_TOP_OF_ATMOS ? diff --git a/parm/postxconfig-NT-GFS-TWO.txt b/parm/postxconfig-NT-GFS-TWO.txt index c895c9338e..cad2999588 100644 --- a/parm/postxconfig-NT-GFS-TWO.txt +++ b/parm/postxconfig-NT-GFS-TWO.txt @@ -1,6 +1,6 @@ 2 104 -187 +196 GFSPRS 0 ncep_nco @@ -2464,7 +2464,7 @@ AVE_TCDC_ON_LOW_CLOUD_LYR ? 1 tmpl4_8 -TCDC +LCDC ? AVE low_cloud_lyr @@ -2501,7 +2501,7 @@ AVE_TCDC_ON_MID_CLOUD_LYR ? 1 tmpl4_8 -TCDC +MCDC ? AVE mid_cloud_lyr @@ -2538,7 +2538,7 @@ AVE_TCDC_ON_HIGH_CLOUD_LYR ? 1 tmpl4_8 -TCDC +HCDC ? AVE high_cloud_lyr @@ -6936,6 +6936,339 @@ surface ? ? ? +260 +HGT_ON_CLOUD_CEILING +? +1 +tmpl4_0 +HGT +? +? +cloud_ceilng +0 +? +0 +? +? +0 +? +0 +? +? +? +0 +0.0 +0 +0.0 +? +0 +0.0 +0 +0.0 +1 +6.0 +0 +0 +0 +? +? +? +161 +INST_TCDC_ON_ENTIRE_ATMOS +? +1 +tmpl4_0 +TCDC +? +? +entire_atmos +0 +? +0 +? +? +0 +? +0 +? +? +? +0 +0.0 +0 +0.0 +? +0 +0.0 +0 +0.0 +1 +3.0 +0 +0 +0 +? +? +? +37 +LCDC_ON_LOW_CLOUD_LYR +? +1 +tmpl4_0 +LCDC +? +? +low_cloud_lyr +0 +? +0 +? +? +0 +? +0 +? +? +? +0 +0.0 +0 +0.0 +? +0 +0.0 +0 +0.0 +1 +3.0 +0 +0 +0 +? +? +? +38 +MCDC_ON_MID_CLOUD_LYR +? +1 +tmpl4_0 +MCDC +? +? +mid_cloud_lyr +0 +? +0 +? +? +0 +? +0 +? +? +? +0 +0.0 +0 +0.0 +? +0 +0.0 +0 +0.0 +1 +3.0 +0 +0 +0 +? +? +? +39 +HCDC_ON_HIGH_CLOUD_LYR +? +1 +tmpl4_0 +HCDC +? +? +high_cloud_lyr +0 +? +0 +? +? +0 +? +0 +? +? +? +0 +0.0 +0 +0.0 +? +0 +0.0 +0 +0.0 +1 +3.0 +0 +0 +0 +? +? +? +253 +REFD_ON_SPEC_HGT_LVL_ABOVE_GRND +? +1 +tmpl4_0 +REFD +NCEP +? +spec_hgt_lvl_above_grnd +0 +? +2 +4000. 1000. +? +0 +? +0 +? +? +? +0 +0.0 +0 +0.0 +? +0 +0.0 +0 +0.0 +1 +4.0 +0 +0 +0 +? +? +? +250 +REFD_ON_HYBRID_LVL +? +1 +tmpl4_0 +REFD +NCEP +? +hybrid_lvl +0 +? +2 +1. 2. +? +0 +? +0 +? +? +? +0 +0.0 +0 +0.0 +? +0 +0.0 +0 +0.0 +1 +4.0 +0 +0 +0 +? +? +? +582 +MIXED_LAYER_CAPE_ON_SPEC_PRES_ABOVE_GRND +? +1 +tmpl4_0 +CAPE +? +? +spec_pres_above_grnd +0 +? +1 +9000. +spec_pres_above_grnd +0 +? +1 +0. +? +? +0 +0.0 +0 +0.0 +? +0 +0.0 +0 +0.0 +1 +4.0 +0 +0 +0 +? +? +? +583 +MIXED_LAYER_CIN_ON_SPEC_PRES_ABOVE_GRND +? +1 +tmpl4_0 +CIN +? +? +spec_pres_above_grnd +0 +? +1 +9000. +spec_pres_above_grnd +0 +? +1 +0. +? +? +0 +0.0 +0 +0.0 +? +0 +0.0 +0 +0.0 +1 +4.0 +0 +0 +0 +? +? +? GFSFLX 0 ncep_nco @@ -7734,7 +8067,7 @@ AVE_TCDC_ON_LOW_CLOUD_LYR ? 1 tmpl4_8 -TCDC +LCDC ? AVE low_cloud_lyr @@ -7771,7 +8104,7 @@ AVE_TCDC_ON_MID_CLOUD_LYR ? 1 tmpl4_8 -TCDC +MCDC ? AVE mid_cloud_lyr @@ -7808,7 +8141,7 @@ AVE_TCDC_ON_HIGH_CLOUD_LYR ? 1 tmpl4_8 -TCDC +HCDC ? AVE high_cloud_lyr diff --git a/parm/postxconfig-NT-GFS.txt b/parm/postxconfig-NT-GFS.txt index e3cfa6cb9e..afa1652bea 100644 --- a/parm/postxconfig-NT-GFS.txt +++ b/parm/postxconfig-NT-GFS.txt @@ -1,5 +1,5 @@ 1 -187 +196 GFSPRS 0 ncep_nco @@ -2463,7 +2463,7 @@ AVE_TCDC_ON_LOW_CLOUD_LYR ? 1 tmpl4_8 -TCDC +LCDC ? AVE low_cloud_lyr @@ -2500,7 +2500,7 @@ AVE_TCDC_ON_MID_CLOUD_LYR ? 1 tmpl4_8 -TCDC +MCDC ? AVE mid_cloud_lyr @@ -2537,7 +2537,7 @@ AVE_TCDC_ON_HIGH_CLOUD_LYR ? 1 tmpl4_8 -TCDC +HCDC ? AVE high_cloud_lyr @@ -6935,3 +6935,336 @@ surface ? ? ? +260 +HGT_ON_CLOUD_CEILING +? +1 +tmpl4_0 +HGT +? +? +cloud_ceilng +0 +? +0 +? +? +0 +? +0 +? +? +? +0 +0.0 +0 +0.0 +? +0 +0.0 +0 +0.0 +1 +6.0 +0 +0 +0 +? +? +? +161 +INST_TCDC_ON_ENTIRE_ATMOS +? +1 +tmpl4_0 +TCDC +? +? +entire_atmos +0 +? +0 +? +? +0 +? +0 +? +? +? +0 +0.0 +0 +0.0 +? +0 +0.0 +0 +0.0 +1 +3.0 +0 +0 +0 +? +? +? +37 +LCDC_ON_LOW_CLOUD_LYR +? +1 +tmpl4_0 +LCDC +? +? +low_cloud_lyr +0 +? +0 +? +? +0 +? +0 +? +? +? +0 +0.0 +0 +0.0 +? +0 +0.0 +0 +0.0 +1 +3.0 +0 +0 +0 +? +? +? +38 +MCDC_ON_MID_CLOUD_LYR +? +1 +tmpl4_0 +MCDC +? +? +mid_cloud_lyr +0 +? +0 +? +? +0 +? +0 +? +? +? +0 +0.0 +0 +0.0 +? +0 +0.0 +0 +0.0 +1 +3.0 +0 +0 +0 +? +? +? +39 +HCDC_ON_HIGH_CLOUD_LYR +? +1 +tmpl4_0 +HCDC +? +? +high_cloud_lyr +0 +? +0 +? +? +0 +? +0 +? +? +? +0 +0.0 +0 +0.0 +? +0 +0.0 +0 +0.0 +1 +3.0 +0 +0 +0 +? +? +? +253 +REFD_ON_SPEC_HGT_LVL_ABOVE_GRND +? +1 +tmpl4_0 +REFD +NCEP +? +spec_hgt_lvl_above_grnd +0 +? +2 +4000. 1000. +? +0 +? +0 +? +? +? +0 +0.0 +0 +0.0 +? +0 +0.0 +0 +0.0 +1 +4.0 +0 +0 +0 +? +? +? +250 +REFD_ON_HYBRID_LVL +? +1 +tmpl4_0 +REFD +NCEP +? +hybrid_lvl +0 +? +2 +1. 2. +? +0 +? +0 +? +? +? +0 +0.0 +0 +0.0 +? +0 +0.0 +0 +0.0 +1 +4.0 +0 +0 +0 +? +? +? +582 +MIXED_LAYER_CAPE_ON_SPEC_PRES_ABOVE_GRND +? +1 +tmpl4_0 +CAPE +? +? +spec_pres_above_grnd +0 +? +1 +9000. +spec_pres_above_grnd +0 +? +1 +0. +? +? +0 +0.0 +0 +0.0 +? +0 +0.0 +0 +0.0 +1 +4.0 +0 +0 +0 +? +? +? +583 +MIXED_LAYER_CIN_ON_SPEC_PRES_ABOVE_GRND +? +1 +tmpl4_0 +CIN +? +? +spec_pres_above_grnd +0 +? +1 +9000. +spec_pres_above_grnd +0 +? +1 +0. +? +? +0 +0.0 +0 +0.0 +? +0 +0.0 +0 +0.0 +1 +4.0 +0 +0 +0 +? +? +? diff --git a/sorc/ncep_post.fd/CLDRAD.f b/sorc/ncep_post.fd/CLDRAD.f index 51ab9f1110..d2a2e5b453 100644 --- a/sorc/ncep_post.fd/CLDRAD.f +++ b/sorc/ncep_post.fd/CLDRAD.f @@ -968,7 +968,7 @@ SUBROUTINE CLDRAD ENDIF ! nmmb_clds1: IF ((MODELNAME=='NMM' .AND. GRIDTYPE=='B') .OR. & - MODELNAME=='FV3R') THEN + MODELNAME=='FV3R'.OR. MODELNAME=='GFS') THEN ! !-- Initialize low, middle, high, and total cloud cover; ! also a method for cloud ceiling height @@ -1274,18 +1274,8 @@ SUBROUTINE CLDRAD ! TOTAL CLOUD FRACTION (INSTANTANEOUS). IF ((IGET(161) > 0) .OR. (IGET(260) > 0)) THEN ! GRID1=SPVAL - IF(MODELNAME == 'GFS')THEN - IF (IGET(161) > 0) THEN + IF(MODELNAME=='NCAR' .OR. MODELNAME=='RAPR')THEN !$omp parallel do private(i,j) - DO J=JSTA,JEND - DO I=1,IM - GRID1(i,j) = SPVAL - EGRID1(i,j) = SPVAL - TCLD(i,j) = SPVAL - ENDDO - ENDDO - ENDIF - ELSE IF(MODELNAME .EQ. 'NCAR' .OR. MODELNAME == 'RAPR')THEN DO J=JSTA,JEND DO I=1,IM GRID1(i,j) = SPVAL @@ -1296,13 +1286,15 @@ SUBROUTINE CLDRAD ENDDO ENDDO - ELSE IF (MODELNAME.EQ.'NMM'.OR.MODELNAME.EQ.'FV3R')THEN + ELSE IF (MODELNAME=='NMM'.OR.MODELNAME=='FV3R' & + .OR. MODELNAME=='GFS')THEN DO J=JSTA,JEND DO I=1,IM ! EGRID1(I,J)=AMAX1(CFRACL(I,J), ! 1 AMAX1(CFRACM(I,J),CFRACH(I,J))) ! EGRID1(I,J)=1.-(1.-CFRACL(I,J))*(1.-CFRACM(I,J))* & ! & (1.-CFRACH(I,J)) + GRID1(i,j)=SPVAL EGRID1(I,J)=TCLD(I,J) ENDDO ENDDO diff --git a/sorc/ncep_post.fd/INITPOST_GFS_NETCDF.f b/sorc/ncep_post.fd/INITPOST_GFS_NETCDF.f index 114b175258..eaf6c83c62 100644 --- a/sorc/ncep_post.fd/INITPOST_GFS_NETCDF.f +++ b/sorc/ncep_post.fd/INITPOST_GFS_NETCDF.f @@ -458,6 +458,7 @@ SUBROUTINE INITPOST_GFS_NETCDF(ncid3d) end do lonstart = nint(glon1d(1)*gdsdegr) lonlast = nint(glon1d(im)*gdsdegr) + dxval = nint(abs(glon1d(1)-glon1d(2))*gdsdegr) else if(numDims==2)then Status=nf90_get_var(ncid3d,varid,dummy) if(maxval(abs(dummy))<2.0*pi)convert_rad_to_deg=.true. @@ -477,9 +478,11 @@ SUBROUTINE INITPOST_GFS_NETCDF(ncid3d) if(convert_rad_to_deg)then lonstart = nint(dummy(1,1)*gdsdegr)*180./pi lonlast = nint(dummy(im,jm)*gdsdegr)*180./pi + dxval = nint(abs(dummy(1,1)-dummy(2,1))*gdsdegr)*180./pi else lonstart = nint(dummy(1,1)*gdsdegr) lonlast = nint(dummy(im,jm)*gdsdegr) + dxval = nint(abs(dummy(1,1)-dummy(2,1))*gdsdegr) end if ! Jili Dong add support for regular lat lon (2019/03/22) start @@ -494,7 +497,7 @@ SUBROUTINE INITPOST_GFS_NETCDF(ncid3d) ! Jili Dong add support for regular lat lon (2019/03/22) end end if - print*,'lonstart,lonlast ',lonstart,lonlast + print*,'lonstart,lonlast,dxval ',lonstart,lonlast,dxval ! get latitude Status=nf90_inq_varid(ncid3d,'grid_yt',varid) Status=nf90_inquire_variable(ncid3d,varid,ndims = numDims) @@ -508,6 +511,7 @@ SUBROUTINE INITPOST_GFS_NETCDF(ncid3d) end do latstart = nint(glat1d(1)*gdsdegr) latlast = nint(glat1d(jm)*gdsdegr) + dyval = nint(abs(glat1d(1)-glat1d(2))*gdsdegr) else if(numDims==2)then Status=nf90_get_var(ncid3d,varid,dummy) if(maxval(abs(dummy))