Merge pull request #109 from ScottSieron-NOAA/feature/gmi

call_crtm bug fix for gmi
NOAA-EMC · Feb 19, 2021 · a944725 · a944725
2 parents 9eb537e + 69f2251
commit a944725
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/gsi/crtm_interface.f90 b/src/gsi/crtm_interface.f90
@@ -2074,7 +2074,7 @@ subroutine call_crtm(obstype,obstime,data_s,nchanl,nreal,ich, &
      hwp_total = sum(hwp_guess(:))
      theta_700 = atmosphere(1)%temperature(idx700)*(r1000/atmosphere(1)%pressure(idx700))**rd_over_cp
      theta_sfc = data_s(itsavg)*(r100/ps)**rd_over_cp
-     stability = theta_700 - theta_sfc
+     if (present(stability)) stability = theta_700 - theta_sfc
   endif
 
 ! Set clouds for CRTM