[Dart-dev] DART/branches Revision: 12024

[dart at ucar.edu]

nancy at ucar.edu
2017-10-27 11:38:27 -0600 (Fri, 27 Oct 2017)
196
add a comment about how we should be handing the vortex
forward operators, and fix a syntax error in the TES nadir
forward operator.  i'm not sure it's correct yet, but at
least it compiles now.




Modified: DART/branches/rma_trunk/observations/forward_operators/obs_def_TES_nadir_mod.f90
===================================================================
--- DART/branches/rma_trunk/observations/forward_operators/obs_def_TES_nadir_mod.f90	2017-10-27 15:51:28 UTC (rev 12023)
+++ DART/branches/rma_trunk/observations/forward_operators/obs_def_TES_nadir_mod.f90	2017-10-27 17:38:27 UTC (rev 12024)
@@ -782,14 +782,16 @@
 !   coordinate fashion where we have N_layers evenly spaced in log pressure
 !   over the pressure difference between p_surf and p_top.
 ! p([0:N]) = psfc * exp( -( log(psfc/ptop)/N ) * [0:N] )
-where (istatus == 0) then
-   ret_p(1:N_layers, :) = 0.0_r8
-   ret_p(0, :) = psfc(:)
-   delz(:) = log( psfc(:) / ptop ) / real(N_layers,r8)
-elsewhere
-   ret_p(1:N_layers, :) = missing_r8
-   delz(:) = missing_r8
-endwhere
+do imem=1, ens_size
+   if (istatus(imem) == 0) then
+      ret_p(1:N_layers, imem) = 0.0_r8
+      ret_p(0, imem) = psfc(imem)
+      delz(imem) = log( psfc(imem) / ptop ) / real(N_layers,r8)
+   else
+      ret_p(1:N_layers, imem) = missing_r8
+      delz(imem) = missing_r8
+   endif
+enddo
 
 where (istatus == 0)
       ret_p(k, :) = psfc * exp( -delz * real(k,r8) )

Modified: DART/branches/rma_trunk/observations/forward_operators/obs_def_vortex_mod.f90
===================================================================
--- DART/branches/rma_trunk/observations/forward_operators/obs_def_vortex_mod.f90	2017-10-27 15:51:28 UTC (rev 12023)
+++ DART/branches/rma_trunk/observations/forward_operators/obs_def_vortex_mod.f90	2017-10-27 17:38:27 UTC (rev 12024)
@@ -4,6 +4,12 @@
 !
 ! $Id$
 
+!>@todo  make a namelist-defined grid (lower left corner, nx, ny, dx, dy)
+!> allocate, call interp for t, u, v, p and fill it.  store the location
+!> of the obs being interpolated.  compute the values.  wait for subsequent
+!> calls for this same location.
+!>
+
 ! BEGIN DART PREPROCESS KIND LIST
 ! VORTEX_LAT, QTY_VORTEX_LAT
 ! VORTEX_LON, QTY_VORTEX_LON