Merge pull request #6 from JiliDong-NOAA/lambert_update

Lambert update; update for jet and write component coordinate change; sync with the latest develop branch
hafs-community · Jan 8, 2020 · 6e6ea73 · 6e6ea73
2 parents baef932 + 4c247c8
commit 6e6ea73
Show file tree

Hide file tree

Showing 21 changed files with 122 additions and 86 deletions.
diff --git a/modulefiles/post/lib-v8.0.0-cray-intel → modulefiles/post/lib-cray-intel b/modulefiles/post/lib-v8.0.0-cray-intel → modulefiles/post/lib-cray-intel
@@ -33,15 +33,18 @@ module use -a /gpfs/hps/nco/ops/nwprod/lib/modulefiles
 # Loading Intel Compiler Suite
 module load PrgEnv-intel
 module load craype-sandybridge
-module switch intel intel/15.0.3.187
+#module switch intel intel/15.0.3.187
+module switch intel intel/16.3.210
 module load craype/2.3.0
 
 # Loading nceplibs modules
 module load jasper-gnu-sandybridge/1.900.1
 module load png-gnu-sandybridge/1.2.49
 module load zlib-gnu-sandybridge/1.2.7
 module load g2-intel/3.1.0
-module load g2tmpl-intel/1.5.0
+module use -a /usrx/local/nceplibs/NCEPLIBS/modulefiles
+module load g2tmpl/1.6.0
+#module load g2tmpl-intel/1.5.0
 #module load xmlparse/v2.0.0
 
 module load w3emc-intel/2.2.0

diff --git a/modulefiles/post/lib-v8.0.0-hera → modulefiles/post/lib-hera b/modulefiles/post/lib-v8.0.0-hera → modulefiles/post/lib-hera
@@ -39,7 +39,7 @@ module load jasper/1.900.1
 module load png/1.2.44
 module load z/1.2.11
 module load g2/3.1.0
-module load g2tmpl/1.5.0
+module load g2tmpl/1.6.0
 #module load xmlparse/v2.0.0
 
 module load w3emc/2.3.1

diff --git a/modulefiles/post/lib-v8.0.0-wcoss → modulefiles/post/lib-wcoss b/modulefiles/post/lib-v8.0.0-wcoss → modulefiles/post/lib-wcoss
diff --git a/modulefiles/post/lib-v8.0.0-wcoss_dell_p3 → modulefiles/post/lib-wcoss_dell_p3 b/modulefiles/post/lib-v8.0.0-wcoss_dell_p3 → modulefiles/post/lib-wcoss_dell_p3
@@ -37,7 +37,8 @@ module load jasper/1.900.1
 module load libpng/1.2.59
 module load zlib/1.2.11
 module load g2/3.1.0
-module load g2tmpl/1.5.0
+module use -a /usrx/local/nceplibs/dev/NCEPLIBS/modulefiles
+module load g2tmpl/1.6.0
 #module load xmlparse/v2.0.0
 
 module load w3emc/2.3.0

diff --git a/modulefiles/post/v8.0.0-cray-intel b/modulefiles/post/v8.0.0-cray-intel
@@ -19,7 +19,8 @@ module use -a /usrx/local/prod/modulefiles
 module use -a /gpfs/hps/nco/ops/nwprod/lib/modulefiles
 module load PrgEnv-intel
 module load craype-sandybridge
-module switch intel intel/15.0.3.187
+#module switch intel intel/15.0.3.187
+module switch intel intel/16.3.210
 module load craype/2.3.0
 module load prod_util/1.0.33
 module load cray-libsci/13.0.3
@@ -47,7 +48,9 @@ module load cray-netcdf/4.3.2
 # module load NetCDF-cray-sandybridge/3.6.3
 #### module load cray-netcdf/4.3.2
 #module load g2tmpl-intel/1.4.0
-module load g2tmpl-intel/1.5.0
+module use -a /usrx/local/nceplibs/NCEPLIBS/modulefiles
+module load g2tmpl/1.6.0
+#module load g2tmpl-intel/1.5.0
 module load wrfio-intel/1.1.1
 
 setenv myFC ftn

diff --git a/modulefiles/post/v8.0.0-hera b/modulefiles/post/v8.0.0-hera
@@ -15,8 +15,10 @@ 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/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
@@ -39,7 +41,7 @@ module load w3nco/2.0.7
 module load crtm/2.2.5
 #module load netcdf/3.6.3
 module load netcdf/4.7.0
-module load g2tmpl/1.5.1
+module load g2tmpl/1.6.0
 module load wrfio/1.1.1
 
 #setenv WRFPATH /scratch3/NCEPDEV/nwprod/sorc/wrf_shared.v1.1.0

diff --git a/modulefiles/post/v8.0.0-jet b/modulefiles/post/v8.0.0-jet
@@ -22,46 +22,50 @@ module load impi
 module use /contrib/modulefiles
 
 # Loding nceplibs modules
-module use -a /mnt/lfs3/projects/hfv3gfs/nwprod/lib/modulefiles
-module load sigio/v2.0.1
+#module use -a /mnt/lfs3/projects/hfv3gfs/nwprod/lib/modulefiles
+module use -a /mnt/lfs3/projects/hfv3gfs/nwprod/NCEPLIBS/modulefiles
+module load sigio/v2.1.0
 module load jasper/v1.900.1
 module load png/v1.2.44
 module load z/v1.2.6
 module load sfcio/v1.0.0
-module load nemsio/v2.2.2
-module load bacio/v2.0.1
+module load nemsio/v2.2.3
+module load bacio/v2.0.2
 #module load g2/v2.5.2
 module load xmlparse/v2.0.0
 module load gfsio/v1.1.0
-module load ip/v3.0.0
+module load ip/v3.0.1
 module load sp/v2.0.2
-module load w3emc/v2.2.0
+module load w3emc/v2.3.0
 module load w3nco/v2.0.6
-module load crtm/v2.2.3
+module load crtm/v2.2.6
+module load g2/v3.1.0
+module load g2tmpl/v1.6.0
+module load wrfio/v1.1.1
 #module load netcdf/3.6.3
 module load netcdf/4.6.1
 
-setenv NCEPLIBS /mnt/lfs3/projects/hfv3gfs/gwv/ljtjet/lib
+#setenv NCEPLIBS /mnt/lfs3/projects/hfv3gfs/gwv/ljtjet/lib
 
-module use /mnt/lfs3/projects/hfv3gfs/gwv/ljtjet/lib/modulefiles
-module load g2tmpl-intel/1.5.0
+#module use /mnt/lfs3/projects/hfv3gfs/gwv/ljtjet/lib/modulefiles
+#module load g2tmpl-intel/1.5.0
 
 #module use /mnt/lfs3/projects/hfv3gfs/gwv/ljtjet/lib/wrf.post.lib/modulefiles
 #module load wrf-io-v1.1.1
 #
-setenv WRFIO_LIB /mnt/lfs3/projects/hfv3gfs/gwv/ljtjet/lib/wrf.post.lib/v1.1.1/lib/wrf_io/libwrfio_nf.a
+#setenv WRFIO_LIB /mnt/lfs3/projects/hfv3gfs/gwv/ljtjet/lib/wrf.post.lib/v1.1.1/lib/wrf_io/libwrfio_nf.a
 
 
-set dlib   /mnt/lfs3/projects/hfv3gfs/gwv/ltmp2/lib/g2/v3.1.0/
-set bname  "G2"
+#set dlib   /mnt/lfs3/projects/hfv3gfs/gwv/ltmp2/lib/g2/v3.1.0/
+#set bname  "G2"
 #
 ### Export environment variables
-setenv ${bname}_SRC  $dlib/src
-setenv ${bname}_INC4 $dlib/intel/include/g2_v3.1.0_4
-setenv ${bname}_INCd $dlib/intel/include/g2_v3.1.0_d
-setenv ${bname}_LIB4 $dlib/intel/libg2_v3.1.0_4.a
-setenv ${bname}_LIBd $dlib/intel/libg2_v3.1.0_d.a
-setenv ${bname}_VER v3.1.0
+#setenv ${bname}_SRC  $dlib/src
+#setenv ${bname}_INC4 $dlib/intel/include/g2_v3.1.0_4
+#setenv ${bname}_INCd $dlib/intel/include/g2_v3.1.0_d
+#setenv ${bname}_LIB4 $dlib/intel/libg2_v3.1.0_4.a
+#setenv ${bname}_LIBd $dlib/intel/libg2_v3.1.0_d.a
+#setenv ${bname}_VER v3.1.0
 
 #setenv WRFPATH /mnt/lfs3/projects/hfv3gfs/nwprod/wrf_shared.v1.1.0/ 
 setenv myFC mpiifort

diff --git a/modulefiles/post/v8.0.0-wcoss_dell_p3 b/modulefiles/post/v8.0.0-wcoss_dell_p3
@@ -32,9 +32,10 @@ module load sp/2.0.2
 module load w3emc/2.3.0          
 module load w3nco/2.0.6
 module load NetCDF/4.5.0
-module load g2tmpl/1.5.0          
 module load crtm/2.2.6
 module load wrfio/1.1.1
+module use -a /usrx/local/nceplibs/dev/NCEPLIBS/modulefiles
+module load g2tmpl/1.6.0          
 
 setenv myFC mpiifort
 setenv OPENMP "-qopenmp"

diff --git a/parm/hafs.xml b/parm/hafs.xml
@@ -13,7 +13,7 @@
     <gen_proc_type>fcst</gen_proc_type>
     <time_range_unit>hour</time_range_unit>
     <orig_center>nws_ncep</orig_center>
-    <gen_proc>hur_hafs</gen_proc>
+    <gen_proc>hur_weather_res_and_fcst</gen_proc>
     <packing_method>complex_packing_spatial_diff</packing_method>
     <order_of_sptdiff>2nd_ord_sptdiff</order_of_sptdiff>
     <field_datatype>fltng_pnt</field_datatype>

diff --git a/parm/post_avblflds.xml b/parm/post_avblflds.xml
@@ -6821,7 +6821,6 @@
        <post_avblfldidx>958</post_avblfldidx>
        <shortname>SBTAGR7_ON_TOP_OF_ATMOS</shortname>
        <pname>SBTAGR7</pname>
-       <table_info>NCEP</table_info>
        <fixed_sfc1_type>top_of_atmos</fixed_sfc1_type>
        <scale>4.0</scale>
     </param>
@@ -6830,7 +6829,6 @@
        <post_avblfldidx>959</post_avblfldidx>
        <shortname>SBTAGR8_ON_TOP_OF_ATMOS</shortname>
        <pname>SBTAGR8</pname>
-       <table_info>NCEP</table_info>
        <fixed_sfc1_type>top_of_atmos</fixed_sfc1_type>
        <scale>4.0</scale>
     </param>
@@ -6839,7 +6837,6 @@
        <post_avblfldidx>960</post_avblfldidx>
        <shortname>SBTAGR9_ON_TOP_OF_ATMOS</shortname>
        <pname>SBTAGR9</pname>
-       <table_info>NCEP</table_info>
        <fixed_sfc1_type>top_of_atmos</fixed_sfc1_type>
        <scale>4.0</scale>
     </param>
@@ -6848,7 +6845,6 @@
        <post_avblfldidx>961</post_avblfldidx>
        <shortname>SBTAGR10_ON_TOP_OF_ATMOS</shortname>
        <pname>SBTAGR10</pname>
-       <table_info>NCEP</table_info>
        <fixed_sfc1_type>top_of_atmos</fixed_sfc1_type>
        <scale>4.0</scale>
     </param>
@@ -6857,7 +6853,6 @@
        <post_avblfldidx>962</post_avblfldidx>
        <shortname>SBTAGR11_ON_TOP_OF_ATMOS</shortname>
        <pname>SBTAGR11</pname>
-       <table_info>NCEP</table_info>
        <fixed_sfc1_type>top_of_atmos</fixed_sfc1_type>
        <scale>4.0</scale>
     </param>
@@ -6866,7 +6861,6 @@
        <post_avblfldidx>963</post_avblfldidx>
        <shortname>SBTAGR12_ON_TOP_OF_ATMOS</shortname>
        <pname>SBTAGR12</pname>
-       <table_info>NCEP</table_info>
        <fixed_sfc1_type>top_of_atmos</fixed_sfc1_type>
        <scale>4.0</scale>
     </param>
@@ -6875,7 +6869,6 @@
        <post_avblfldidx>964</post_avblfldidx>
        <shortname>SBTAGR13_ON_TOP_OF_ATMOS</shortname>
        <pname>SBTAGR13</pname>
-       <table_info>NCEP</table_info>
        <fixed_sfc1_type>top_of_atmos</fixed_sfc1_type>
        <scale>4.0</scale>
     </param>
@@ -6884,7 +6877,6 @@
        <post_avblfldidx>965</post_avblfldidx>
        <shortname>SBTAGR14_ON_TOP_OF_ATMOS</shortname>
        <pname>SBTAGR14</pname>
-       <table_info>NCEP</table_info>
        <fixed_sfc1_type>top_of_atmos</fixed_sfc1_type>
        <scale>4.0</scale>
     </param>
@@ -6893,7 +6885,6 @@
        <post_avblfldidx>966</post_avblfldidx>
        <shortname>SBTAGR15_ON_TOP_OF_ATMOS</shortname>
        <pname>SBTAGR15</pname>
-       <table_info>NCEP</table_info>
        <fixed_sfc1_type>top_of_atmos</fixed_sfc1_type>
        <scale>4.0</scale>
     </param>
@@ -6902,7 +6893,6 @@
        <post_avblfldidx>967</post_avblfldidx>
        <shortname>SBTAGR16_ON_TOP_OF_ATMOS</shortname>
        <pname>SBTAGR16</pname>
-       <table_info>NCEP</table_info>
        <fixed_sfc1_type>top_of_atmos</fixed_sfc1_type>
        <scale>4.0</scale>
     </param>

diff --git a/parm/postxconfig-NT-GFS-GOES.txt b/parm/postxconfig-NT-GFS-GOES.txt
@@ -170,7 +170,7 @@ SBTAGR8_ON_TOP_OF_ATMOS
 1
 tmpl4_0
 SBTAGR8
-NCEP
+?
 ?
 top_of_atmos
 0
@@ -207,7 +207,7 @@ SBTAGR9_ON_TOP_OF_ATMOS
 1
 tmpl4_0
 SBTAGR9
-NCEP
+?
 ?
 top_of_atmos
 0
@@ -244,7 +244,7 @@ SBTAGR10_ON_TOP_OF_ATMOS
 1
 tmpl4_0
 SBTAGR10
-NCEP
+?
 ?
 top_of_atmos
 0
@@ -281,7 +281,7 @@ SBTAGR13_ON_TOP_OF_ATMOS
 1
 tmpl4_0
 SBTAGR13
-NCEP
+?
 ?
 top_of_atmos
 0

diff --git a/parm/postxconfig-NT-hafs.txt b/parm/postxconfig-NT-hafs.txt
@@ -11,7 +11,7 @@ fcst
 fcst
 hour
 nws_ncep
-hur_hafs
+hur_weather_res_and_fcst
 complex_packing_spatial_diff
 2nd_ord_sptdiff
 fltng_pnt

diff --git a/sorc/ncep_post.fd/CALRAD_WCLOUD_newcrtm.f b/sorc/ncep_post.fd/CALRAD_WCLOUD_newcrtm.f
@@ -358,19 +358,19 @@ SUBROUTINE CALRAD_WCLOUD
               czen(i,j)=cos(sun_zenith_rad)
 	   end do
 	end do
-        if(jj>=jsta .and. jj<=jend)                                  &
+        if(jj>=jsta .and. jj<=jend.and.debugprint)                                  &
             print*,'sample GFS zenith angle=',acos(czen(ii,jj))*rtd   
 !     end if	       
      ! Initialize CRTM.  Load satellite sensor array.
      ! The optional arguments Process_ID and Output_Process_ID limit
      ! generation of runtime informative output to mpi task
      ! Output_Process_ID (which here is set to be task 0)
-     print*,'success in CALRAD= ',success
+     if(me==0)print*,'success in CALRAD= ',success
      allocate( channelinfo(n_sensors))
 
      error_status = crtm_init(sensorlist_local,channelinfo,   &
           Process_ID=0,Output_Process_ID=0 )
-     print*, 'channelinfo after init= ',channelinfo(1)%sensor_id, &
+     if(me==0)print*, 'channelinfo after init= ',channelinfo(1)%sensor_id, &
               channelinfo(2)%sensor_id
      if (error_status /= 0_i_kind)                                  &
          write(6,*)'ERROR*** crtm_init error_status=',error_status
@@ -477,7 +477,7 @@ SUBROUTINE CALRAD_WCLOUD
      ! Loop over data types to process    
      sensordo: do isat=1,n_sensors
 
-        print*,'n_sensor,obstype,isis',isat,obslist(isat),sensorlist(isat)
+        if(me==0)print*,'n_sensor,obstype,isis',isat,obslist(isat),sensorlist(isat)
 
         obstype=obslist(isat) 
         isis=trim(sensorlist(isat))
@@ -510,7 +510,7 @@ SUBROUTINE CALRAD_WCLOUD
              (isis=='abi_g17'  .and. post_abig17) .OR. &
              (isis=='abi_gr'   .and. post_abigr) .OR. &
              (isis=='seviri_m10' .and. iget(876)>0) )then
-           print*,'obstype, isis= ',obstype,isis
+           if(me==0)print*,'obstype, isis= ',obstype,isis
            !       isis='amsua_n15'
 
            ! Initialize logical flags for satellite platform
@@ -680,7 +680,7 @@ SUBROUTINE CALRAD_WCLOUD
                          .and. geometryinfo(1)%sensor_zenith_angle >= 0.0_r_kind)THEN
                        geometryinfo(1)%source_zenith_angle = acos(czen(i,j))*rtd ! solar zenith angle
                        geometryinfo(1)%sensor_scan_angle   = 0. ! scan angle, assuming nadir
-                       if(i==ii.and.j==jj)print*,'sample geometry ',                   &
+                       if(i==ii.and.j==jj.and.debugprint)print*,'sample geometry ',                   &
                                geometryinfo(1)%sensor_zenith_angle                     &
                               ,geometryinfo(1)%source_zenith_angle                     &
                               ,czen(i,j)*rtd 
@@ -705,7 +705,7 @@ SUBROUTINE CALRAD_WCLOUD
                           else
                              snoeqv=0.
                           end if
-                          if(i==ii.and.j==jj)print*,'sno,itype,ivgtyp B cing snfrc = ',  &
+                          if(i==ii.and.j==jj.and.debugprint)print*,'sno,itype,ivgtyp B cing snfrc = ',  &
                                                      snoeqv,itype,IVGTYP(I,J)
                           if(sm(i,j) > 0.1)then
                              sfcpct(4)=0.
@@ -869,7 +869,7 @@ SUBROUTINE CALRAD_WCLOUD
                           if(surface(1)%snow_depth<0. .or.  surface(1)%snow_depth>10000.) &
                              print*,'bad snow_depth'
                        end if
-                       if(i==ii.and.j==jj)print*,'sample surface in CALRAD=', &
+                       if(i==ii.and.j==jj.and.debugprint)print*,'sample surface in CALRAD=', &
                              i,j,surface(1)%wind_speed,surface(1)%water_coverage,       &
                              surface(1)%land_coverage,surface(1)%ice_coverage,          &
                              surface(1)%snow_coverage,surface(1)%land_temperature,      &
@@ -882,7 +882,7 @@ SUBROUTINE CALRAD_WCLOUD
 
                        !       Load atmosphere profiles into RTM model layers
                        !       CRTM counts from top down just as post does
-                       if(i==ii.and.j==jj)print*,'TOA= ',atmosphere(1)%level_pressure(0)
+                       if(i==ii.and.j==jj.and.debugprint)print*,'TOA= ',atmosphere(1)%level_pressure(0)
                        do k = 1,lm
                           atmosphere(1)%level_pressure(k) = pint(i,j,k+1)/r100
                           atmosphere(1)%pressure(k)       = pmid(i,j,k)/r100
@@ -910,7 +910,7 @@ SUBROUTINE CALRAD_WCLOUD
                              !     &      atmosphere(1)%absorber(k,1)>1.)  &
                              !     &      print*,'bad atmosphere o3'
                           end if
-                          if(i==ii.and.j==jj)print*,'sample atmosphere in CALRAD=',  &
+                          if(i==ii.and.j==jj.and.debugprint)print*,'sample atmosphere in CALRAD=',  &
    	                        i,j,k,atmosphere(1)%level_pressure(k),atmosphere(1)%pressure(k),  &
                                 atmosphere(1)%temperature(k),atmosphere(1)%absorber(k,1),  &
                                 atmosphere(1)%absorber(k,2)