[Dart-dev] [9607] DART/trunk/obs_def/obs_def_COSMOS_mod.f90: this obs_def was bypassing the assim_model layer by

nancy at ucar.edu nancy at ucar.edu
Tue Jan 26 13:19:33 MST 2016


Revision: 9607
Author:   nancy
Date:     2016-01-26 13:19:32 -0700 (Tue, 26 Jan 2016)
Log Message:
-----------
this obs_def was bypassing the assim_model layer by
calling the model_mod directly instead of going through
the assim_model code like all the other obs_defs do.

Modified Paths:
--------------
    DART/trunk/obs_def/obs_def_COSMOS_mod.f90

-------------- next part --------------
Modified: DART/trunk/obs_def/obs_def_COSMOS_mod.f90
===================================================================
--- DART/trunk/obs_def/obs_def_COSMOS_mod.f90	2016-01-26 20:17:49 UTC (rev 9606)
+++ DART/trunk/obs_def/obs_def_COSMOS_mod.f90	2016-01-26 20:19:32 UTC (rev 9607)
@@ -76,7 +76,7 @@
                              vert_is_height,   VERTISHEIGHT,            &
                              vert_is_level,    VERTISLEVEL,             &
                              set_location_missing
-use        model_mod, only : model_interpolate
+use  assim_model_mod, only : interpolate
 use     obs_kind_mod, only : KIND_GEOPOTENTIAL_HEIGHT, KIND_SOIL_MOISTURE
 
 use typesizes
@@ -507,7 +507,7 @@
 nlevels = 0
 COUNTLEVELS : do i = 1,maxlayers
    loc = set_location(loc_lon, loc_lat, real(i,r8), VERTISLEVEL)
-   call model_interpolate(state,loc,KIND_GEOPOTENTIAL_HEIGHT,loc_value,istatus)
+   call interpolate(state,loc,KIND_GEOPOTENTIAL_HEIGHT,loc_value,istatus)
    if (istatus /= 0) exit COUNTLEVELS
    nlevels = nlevels + 1
 enddo COUNTLEVELS
@@ -529,11 +529,11 @@
 
 FINDLEVELS : do i = 1,nlevels
    loc = set_location(loc_lon, loc_lat, real(i,r8), VERTISLEVEL)
-   call model_interpolate(state,loc,KIND_GEOPOTENTIAL_HEIGHT,layerz(i),istatus)
+   call interpolate(state,loc,KIND_GEOPOTENTIAL_HEIGHT,layerz(i),istatus)
    ! not checking this error code because it worked just a few lines earlier
 
    loc = set_location(loc_lon, loc_lat, layerz(i), VERTISHEIGHT)
-   call model_interpolate(state,loc,KIND_SOIL_MOISTURE,loc_value,istatus)
+   call interpolate(state,loc,KIND_SOIL_MOISTURE,loc_value,istatus)
 
    if (istatus /= 0) then
       write(string1,*) 'FAILED to determine soil moisture for layer',i


More information about the Dart-dev mailing list