Ufo geovals2

modulefiles/gsi_orion.lua

@@ -19,6 +19,16 @@ load(pathJoin("hpc-impi", hpc_impi_ver))
 load(pathJoin("cmake", cmake_ver))
 
 load("gsi_common")
+setenv("crtm_ROOT","/work/noaa/da/eliu/JEDI-GDAS/crtm_v2.4.1-jedi.1-intel2022/build")
+setenv("crtm_VERSION","2.4.1-jedi.1")
+setenv("CRTM_INC","/work/noaa/da/eliu/JEDI-GDAS/crtm_v2.4.1-jedi.1-intel2022/build/module")
+setenv("CRTM_LIB","/work/noaa/da/eliu/JEDI-GDAS/crtm_v2.4.1-jedi.1-intel2022/build/lib/libcrtm_static.a")
+setenv("CRTM_FIX","/work/noaa/da/eliu/JEDI-GDAS/crtm_v2.4.1-jedi.1-fix_gdasapp/fix")
+whatis("Name: crtm")
+whatis("Version: 2.4.1-jedi.1")
+whatis("Category: library")
+whatis("Description: crtm library")
+
 
 load(pathJoin("prod_util", prod_util_ver))
 

src/gsi/CMakeLists.txt

@@ -146,6 +146,7 @@ target_link_libraries(gsi_fortran_obj PUBLIC nemsio::nemsio)
 target_link_libraries(gsi_fortran_obj PUBLIC ncio::ncio)
 target_link_libraries(gsi_fortran_obj PUBLIC w3emc::w3emc_d)
 target_link_libraries(gsi_fortran_obj PUBLIC sp::sp_d)
+add_library(crtm::crtm ALIAS crtm)
 target_link_libraries(gsi_fortran_obj PUBLIC bufr::bufr_d)
 target_link_libraries(gsi_fortran_obj PUBLIC crtm::crtm)
 if(GSI_MODE MATCHES "Regional")

src/gsi/crtm_interface.f90

@@ -85,6 +85,7 @@ module crtm_interface
 public destroy_crtm         ! Subroutine destroys initialization for crtm
 public sensorindex
 public surface
+public atmosphere
 public isatid               ! = 1  index of satellite id
 public itime                ! = 2  index of analysis relative obs time
 public ilon                 ! = 3  index of grid relative obs location (x)
@@ -125,6 +126,8 @@ module crtm_interface
 public idtw                 ! = 35/37 index of d(Tw)
 public idtc                 ! = 36/38 index of d(Tc)
 public itz_tr               ! = 37/39 index of d(Tz)/d(Tr)
+public n_clouds_fwd_wk
+public n_absorbers
 
 ! For TMI and GMI
 public iedge_log            ! = 32  ! index, if obs is to be obleted beause of locating near scan edges.
@@ -202,6 +205,7 @@ module crtm_interface
   logical        ,save :: mixed_use
   logical        ,save :: use_gfdl_qsat 
   integer(i_kind), parameter :: min_n_absorbers = 2
+  integer(i_kind) :: n_absorbers
 
   integer(i_kind),save :: iedge_log
   integer(i_kind),save :: ilzen_ang2,ilazi_ang2,iscan_ang2,iszen_ang2,isazi_ang2
@@ -356,7 +360,6 @@ subroutine init_crtm(init_pass,mype_diaghdr,mype,nchanl,nreal,isis,obstype,radmo
 ! ...all "additional absorber" variables
   integer(i_kind) :: j,icount
   integer(i_kind) :: ig
-  integer(i_kind) :: n_absorbers
   logical quiet
   logical print_verbose
 

src/gsi/genstats_gps.f90

@@ -64,6 +64,15 @@ module m_gpsStats
       integer(i_kind) :: idv,iob	      ! device id and obs index for sorting
       real   (r_kind) :: elat, elon      ! earth lat-lon for redistribution
       !real   (r_kind) :: dlat, dlon      ! earth lat-lon for redistribution
+      !
+      real(r_kind),dimension(:),pointer :: tsenges => NULL()
+      real(r_kind),dimension(:),pointer :: tvirges => NULL()
+      real(r_kind),dimension(:),pointer :: sphmges => NULL()
+      real(r_kind),dimension(:),pointer :: hgtlges => NULL()
+      real(r_kind),dimension(:),pointer :: hgtiges => NULL()
+      real(r_kind),dimension(:),pointer :: prslges => NULL()
+      real(r_kind),dimension(:),pointer :: prsiges => NULL()
+      !
     end type gps_all_ob_type
 
     type gps_all_ob_head
@@ -427,7 +436,10 @@ subroutine genstats_gps(bwork,awork,toss_gps_sub,conv_diagsave,mype)
         end do
      END DO
      if(icnt > 0)then
-        nreal =22
+!
+!       nreal =22
+        nreal =34
+!
         ioff  =nreal
         if (lobsdiagsave) nreal=nreal+4*miter+1
         if (save_jacobian) then
@@ -439,8 +451,6 @@ subroutine genstats_gps(bwork,awork,toss_gps_sub,conv_diagsave,mype)
      end if
   endif
 
-
-
 ! Loop over data to apply final qc, superobs factors, accumulate
 ! statistics and (optionally) load diagnostic output arrays
   icnt=0
@@ -534,10 +544,12 @@ subroutine genstats_gps(bwork,awork,toss_gps_sub,conv_diagsave,mype)
           elev         = gps_allptr%rdiag(7)
           dobs         = gps_allptr%rdiag(17)
           if  (toss_gps(kprof) > zero .and. (dobs == toss_gps(kprof) .or. elev < dobs_height(kprof))) then ! SR from obs
+
               if(ratio_errors*data_ier > tiny_r_kind) then ! obs was good
                  if (luse) then
                     if(conv_diagsave) then
-                      gps_allptr%rdiag(10) = four
+                      !gps_allptr%rdiag(10) = four
+                      gps_allptr%rdiag(10) = 8
                       gps_allptr%rdiag(12) = -one
                       gps_allptr%rdiag(16) = zero
                       if(lobsdiagsave) gps_allptr%rdiag(mreal+jiter) = -one
@@ -704,7 +716,6 @@ subroutine genstats_gps(bwork,awork,toss_gps_sub,conv_diagsave,mype)
     endif
   endif
 
-
 ! Destroy arrays holding gps data
   call destroy_genstats_gps
 contains
@@ -760,19 +771,25 @@ subroutine contents_netcdf_diag_
 ! Observation class
   character(7),parameter     :: obsclass = '    gps'
 
-           call nc_diag_metadata("Station_ID",                            gps_allptr%cdiag             )
-           call nc_diag_metadata("Observation_Class",                     obsclass                     )
+           call nc_diag_metadata("Station_ID",                                 gps_allptr%cdiag             )
+           call nc_diag_metadata("Observation_Class",                          obsclass                     )
            obstype    = gps_allptr%rdiag(1) 
            obssubtype = gps_allptr%rdiag(2)
            call nc_diag_metadata("Observation_Type",                      obstype                      )
-           call nc_diag_metadata("Observation_Subtype",                   obssubtype                   )
+! rename it to be consistent with GMAO gsi_ncdiag
+!          call nc_diag_metadata("Observation_Subtype",                   obssubtype                   )
+           call nc_diag_metadata("record_number",                         obssubtype                   )
            call nc_diag_metadata_to_single("Latitude",                    gps_allptr%rdiag(3)          )
            call nc_diag_metadata_to_single("Longitude",                   gps_allptr%rdiag(4)          )
            call nc_diag_metadata_to_single("Incremental_Bending_Angle",   gps_allptr%rdiag(5)          )
-           call nc_diag_metadata_to_single("Pressure",                    gps_allptr%rdiag(6)          )
-           call nc_diag_metadata_to_single("Height",                      gps_allptr%rdiag(7)          )
+!          call nc_diag_metadata_to_single("Pressure",                    gps_allptr%rdiag(6)          )  !orig
+           call nc_diag_metadata_to_single("Pressure",                    gps_allptr%rdiag(6)*100.0    )
+! rename the variable as impact_height
+!          call nc_diag_metadata_to_single("Height",                      gps_allptr%rdiag(7)          )
+           call nc_diag_metadata_to_single("Impact_Height",               gps_allptr%rdiag(7)          )
            call nc_diag_metadata_to_single("Time",                        gps_allptr%rdiag(8)          )
-           call nc_diag_metadata_to_single("Model_Elevation",             gps_allptr%rdiag(9)          )
+!          call nc_diag_metadata_to_single("Model_Elevation",             gps_allptr%rdiag(9)          )  !orig
+           call nc_diag_metadata_to_single("surface_geopotential_height", gps_allptr%rdiag(9)          )
            call nc_diag_metadata_to_single("Setup_QC_Mark",               gps_allptr%rdiag(10)         )
            call nc_diag_metadata_to_single("Prep_Use_Flag",               gps_allptr%rdiag(11)         )
            call nc_diag_metadata_to_single("Analysis_Use_Flag",           gps_allptr%rdiag(12)         )
@@ -788,14 +805,33 @@ subroutine contents_netcdf_diag_
            call nc_diag_metadata_to_single("Temperature_at_Obs_Location", gps_allptr%rdiag(18)         )
            call nc_diag_metadata_to_single("Specific_Humidity_at_Obs_Location",gps_allptr%rdiag(21)    )
 
+           call nc_diag_metadata("impact_parameter",                           sngl(gps_allptr%rdiag(23))   )
+           call nc_diag_metadata("pccf",                                       sngl(gps_allptr%rdiag(24))   )
+           call nc_diag_metadata("reference_sat_id",                           int(gps_allptr%rdiag(25))   )
+           call nc_diag_metadata("earth_radius_of_curvature",                  sngl(gps_allptr%rdiag(26))   )
+           call nc_diag_metadata("geoid_height_above_reference_ellipsoid",     sngl(gps_allptr%rdiag(27))   )
+           call nc_diag_metadata("qfro",                                       int(gps_allptr%rdiag(28))   )
+           call nc_diag_metadata("ascending_flag",                              int(gps_allptr%rdiag(29))   )
+           call nc_diag_metadata("sensor_azimuth_angle",                       sngl(gps_allptr%rdiag(30))   )
+           call nc_diag_metadata("sat_constellation",                           int(gps_allptr%rdiag(31))   )
+           call nc_diag_metadata("occulting_sat",                               int(gps_allptr%rdiag(32))   )
+           call nc_diag_metadata("process_center",                              int(gps_allptr%rdiag(33))   )
+           call nc_diag_metadata("atmospheric_refractivity",                   sngl(gps_allptr%rdiag(34))   )
+           call nc_diag_metadata("surface_altitude",                           sngl(gps_allptr%rdiag(9))    )
+
            if (save_jacobian) then
               call readarray(dhx_dx, gps_allptr%rdiag(ioff+1:nreal))
               call nc_diag_data2d("Observation_Operator_Jacobian_stind", dhx_dx%st_ind(1:dhx_dx%nind))
               call nc_diag_data2d("Observation_Operator_Jacobian_endind", dhx_dx%end_ind(1:dhx_dx%nind))
               call nc_diag_data2d("Observation_Operator_Jacobian_val", real(dhx_dx%val(1:dhx_dx%nnz),r_single))
            endif
-
-
+              call nc_diag_data2d("atmosphere_pressure_coordinate", sngl(gps_allptr%prslges))
+              call nc_diag_data2d("atmosphere_pressure_coordinate_interface", sngl(gps_allptr%prsiges))
+              call nc_diag_data2d("air_temperature", sngl(gps_allptr%tsenges))
+              call nc_diag_data2d("virtual_temperature", sngl(gps_allptr%tvirges))
+              call nc_diag_data2d("specific_humidity", sngl(gps_allptr%sphmges))
+              call nc_diag_data2d("geopotential_height", sngl(gps_allptr%hgtlges))
+              call nc_diag_data2d("geopotential_height_levels", sngl(gps_allptr%hgtiges))
 
 !           call nc_diag_data2d("T_Jacobian",                              gps_allptr%mmpoint%jac_t             )
            if (lobsdiagsave) then

src/gsi/gesinfo.F90

@@ -588,6 +588,15 @@ subroutine gesinfo
   time_offset=time_offset+fha(3)/r60
 #endif
 
+  write(6,*)'emily time: nhr_assimilation = ', nhr_assimilation 
+  write(6,*)'emily time: idate nmin_an    = ', idate , nmin_an
+  write(6,*)'emily time: ibdate           = ', ibdate
+  write(6,*)'emily time: iadatebgn        = ', iadateend 
+  write(6,*)'emily time: iedate           = ', iedate
+  write(6,*)'emily time: iadateend        = ', iedate
+  write(6,*)'emily time: ianldate         = ', ianldate 
+  write(6,*)'emily time: iwinbgn          = ', iwinbgn 
+
 ! Get information about date/time and number of guess files
   if (regional) then
      if(wrf_nmm_regional) then