[Dart-dev] DART/releases Revision: 10753

[dart at ucar.edu]

nancy at ucar.edu
2016-11-15 11:05:19 -0700 (Tue, 15 Nov 2016)
289
minor bug fix in the AIRS converter to allow it to ignore
data which are marked with a good QC but have a missing value
(-9999) as the lowest valid starting level.  

this commit updates the trunk, the rma_trunk, and the
lanai release.  also updated the change log for the lanai release.




Modified: DART/branches/rma_trunk/observations/AIRS/airs_obs_mod.f90
===================================================================
--- DART/branches/rma_trunk/observations/AIRS/airs_obs_mod.f90	2016-11-15 04:44:47 UTC (rev 10752)
+++ DART/branches/rma_trunk/observations/AIRS/airs_obs_mod.f90	2016-11-15 18:05:19 UTC (rev 10753)
@@ -238,8 +238,11 @@
       obs_time = base_time + set_time(int(granule%Time(icol, irow)))
       call get_time(obs_time, seconds, days)
      
+      ! avoid -9999 in the starting vertical level
+      istart = granule%nBestStd(icol, irow)
+      if (istart < 1) goto 100   ! skip vert_T_loop
 
-      vert_T_loop: do ivert=granule%nBestStd(icol, irow), temperature_top_index
+      vert_T_loop: do ivert=istart, temperature_top_index
 
          tqc = 0   ! if we get here, the quality control is 'best' == 0
 
@@ -281,7 +284,12 @@
  
       enddo vert_T_loop
 
+100   continue
+
+      ! avoid -9999 in the starting vertical layer
       istart = granule%nSurfStd(icol, irow)
+      if (istart < 1) goto 200  !skip vert_Q_loop
+
       vert_Q_loop:  do ivert=istart,humidity_top_index
 
          if (granule%Qual_H2O(icol, irow) > 0) exit vert_Q_loop
@@ -330,6 +338,8 @@
  
       enddo vert_Q_loop
 
+200   continue
+
    enddo colloop
 enddo rowloop
 

Modified: DART/releases/Lanai/index.html
===================================================================
--- DART/releases/Lanai/index.html	2016-11-15 04:44:47 UTC (rev 10752)
+++ DART/releases/Lanai/index.html	2016-11-15 18:05:19 UTC (rev 10753)
@@ -201,6 +201,23 @@
 </TR>
 -->
 
+<TR><TD> 15 Nov 2016 </TD>
+    <TD> 10753 </TD>
+    <TD> Minor code change to the AIRS observation converter to allow it
+to convert input HDF files which have a missing value in the "lowest valid
+level" field.  Does not affect results for files which were able to be
+converted by the original code.
+    </TD>
+</TR>
+
+<TR><TD> 1 Nov 2016 </TD>
+    <TD> 10700 </TD>
+    <TD> Code tree moved to a new subversion server as an intermediate
+step to converting to GIT.  Revision numbers continue sequentially
+without a gap.
+    </TD>
+</TR>
+
 <TR><TD> 1 Aug 2016 </TD>
     <TD> 10604 </TD>
     <TD> Documentation change only.  Minor updates to the 

Modified: DART/releases/Lanai/observations/AIRS/airs_obs_mod.f90
===================================================================
--- DART/releases/Lanai/observations/AIRS/airs_obs_mod.f90	2016-11-15 04:44:47 UTC (rev 10752)
+++ DART/releases/Lanai/observations/AIRS/airs_obs_mod.f90	2016-11-15 18:05:19 UTC (rev 10753)
@@ -238,8 +238,11 @@
       obs_time = base_time + set_time(int(granule%Time(icol, irow)))
       call get_time(obs_time, seconds, days)
      
+      ! avoid -9999 in the starting vertical level
+      istart = granule%nBestStd(icol, irow)
+      if (istart < 1) goto 100   ! skip vert_T_loop
 
-      vert_T_loop: do ivert=granule%nBestStd(icol, irow), temperature_top_index
+      vert_T_loop: do ivert=istart, temperature_top_index
 
          tqc = 0   ! if we get here, the quality control is 'best' == 0
 
@@ -281,7 +284,12 @@
  
       enddo vert_T_loop
 
+100   continue
+
+      ! avoid -9999 in the starting vertical layer
       istart = granule%nSurfStd(icol, irow)
+      if (istart < 1) goto 200  !skip vert_Q_loop
+
       vert_Q_loop:  do ivert=istart,humidity_top_index
 
          if (granule%Qual_H2O(icol, irow) > 0) exit vert_Q_loop
@@ -330,6 +338,8 @@