[Dart-dev] DART/branches Revision: 12600

[dart at ucar.edu]

thoar at ucar.edu
2018-05-24 13:16:16 -0600 (Thu, 24 May 2018)
85
Adding support for a 'biomass' forward operator.

The algorithm came from Andy Fox.




Modified: DART/branches/cesm_clm/models/clm/model_mod.f90
===================================================================
--- DART/branches/cesm_clm/models/clm/model_mod.f90	2018-05-24 17:24:27 UTC (rev 12599)
+++ DART/branches/cesm_clm/models/clm/model_mod.f90	2018-05-24 19:16:16 UTC (rev 12600)
@@ -68,6 +68,8 @@
                              QTY_SNOWCOVER_FRAC,         &
                              QTY_SNOW_THICKNESS,         &
                              QTY_LEAF_CARBON,            &
+                             QTY_LIVE_STEM_CARBON,       &
+                             QTY_DEAD_STEM_CARBON,       &
                              QTY_LEAF_AREA_INDEX,        &
                              QTY_WATER_TABLE_DEPTH,      &
                              QTY_GEOPOTENTIAL_HEIGHT,    &
@@ -704,6 +706,8 @@
 ! of the state vector index in question ... as well as the actual surface area.
 ! This MUST stride through the state vector the same way the state vector is filled.
 
+!@todo BAD BAD BAD
+
 allocate(lonixy(model_size), latjxy(model_size), levels(model_size), landarea(model_size))
 
 ! Initialize all levels to surface. If there is a level, we will explicitly specify it.
@@ -2059,7 +2063,8 @@
           QTY_WATER_TABLE_DEPTH, QTY_VEGETATION_TEMPERATURE, &
           QTY_FRAC_PHOTO_AVAIL_RADIATION, &
           QTY_FPAR_SUNLIT_DIRECT, QTY_FPAR_SUNLIT_DIFFUSE, &
-          QTY_FPAR_SHADED_DIRECT, QTY_FPAR_SHADED_DIFFUSE)
+          QTY_FPAR_SHADED_DIRECT, QTY_FPAR_SHADED_DIFFUSE, &
+          QTY_LIVE_STEM_CARBON,   QTY_DEAD_STEM_CARBON)
 
       call compute_gridcell_value(state_handle, ens_size, location, obs_kind, expected_obs, istatus)
 

Modified: DART/branches/cesm_clm/observations/forward_operators/obs_def_tower_mod.f90
===================================================================
--- DART/branches/cesm_clm/observations/forward_operators/obs_def_tower_mod.f90	2018-05-24 17:24:27 UTC (rev 12599)
+++ DART/branches/cesm_clm/observations/forward_operators/obs_def_tower_mod.f90	2018-05-24 19:16:16 UTC (rev 12600)
@@ -15,8 +15,10 @@
 !MODIS_SNOWCOVER_FRAC,           QTY_SNOWCOVER_FRAC,             COMMON_CODE
 !MODIS_LEAF_AREA_INDEX,          QTY_LEAF_AREA_INDEX,            COMMON_CODE
 !MODIS_FPAR,                     QTY_FRACTION_ABSORBED_PAR,      COMMON_CODE
-!BIOMASS,                        QTY_BIOMASS,                    COMMON_CODE
+!BIOMASS,                        QTY_BIOMASS
 !LEAF_CARBON,                    QTY_LEAF_CARBON,                COMMON_CODE
+!LIVE_STEM_CARBON,               QTY_LIVE_STEM_CARBON,           COMMON_CODE
+!DEAD_STEM_CARBON,               QTY_DEAD_STEM_CARBON,           COMMON_CODE
 !LEAF_AREA_INDEX,                QTY_LEAF_AREA_INDEX,            COMMON_CODE
 !LEAF_NITROGEN,                  QTY_LEAF_NITROGEN,              COMMON_CODE
 !TOWER_AIR_TEMPERATURE,          QTY_TEMPERATURE,                COMMON_CODE
@@ -38,7 +40,8 @@
 !-----------------------------------------------------------------------------
 ! BEGIN DART PREPROCESS USE OF SPECIAL OBS_DEF MODULE
 !  use obs_def_tower_mod, only : get_scalar_from_history, &
-!                                calculate_albedo
+!                                calculate_albedo, &
+!                                calculate_biomass
 ! END DART PREPROCESS USE OF SPECIAL OBS_DEF MODULE
 !-----------------------------------------------------------------------------
 
@@ -64,6 +67,8 @@
 !                    copy_indices, location, obs_time, expected_obs, istatus)
 !  case(SURFACE_ALBEDO)
 !     call calculate_albedo(state_handle, ens_size, location, expected_obs, istatus)
+!  case(BIOMASS)
+!     call calculate_biomass(state_handle, ens_size, location, expected_obs, istatus)
 ! END DART PREPROCESS GET_EXPECTED_OBS_FROM_DEF
 !-----------------------------------------------------------------------------
 
@@ -75,7 +80,8 @@
 !         TOWER_GPP_FLUX, &
 !         TOWER_ER_FLUX, &
 !         SOIL_RESPIRATION_FLUX, &
-!         SURFACE_ALBEDO)
+!         SURFACE_ALBEDO, &
+!         BIOMASS)
 !       continue
 ! END DART PREPROCESS READ_OBS_DEF
 !-----------------------------------------------------------------------------
@@ -88,7 +94,8 @@
 !         TOWER_GPP_FLUX, &
 !         TOWER_ER_FLUX, &
 !         SOIL_RESPIRATION_FLUX, &
-!         SURFACE_ALBEDO)
+!         SURFACE_ALBEDO, &
+!         BIOMASS)
 !       continue
 ! END DART PREPROCESS WRITE_OBS_DEF
 !-----------------------------------------------------------------------------
@@ -101,7 +108,8 @@
 !         TOWER_GPP_FLUX, &
 !         TOWER_ER_FLUX, &
 !         SOIL_RESPIRATION_FLUX, &
-!         SURFACE_ALBEDO)
+!         SURFACE_ALBEDO, &
+!         BIOMASS)
 !       continue
 ! END DART PREPROCESS INTERACTIVE_OBS_DEF
 !-----------------------------------------------------------------------------
@@ -135,7 +143,10 @@
 use     obs_kind_mod, only : QTY_RADIATION_VISIBLE_DOWN, &
                              QTY_RADIATION_NEAR_IR_DOWN, &