From 8389f813d2dfa96c4db00a23601e5ccb77a0bfbf Mon Sep 17 00:00:00 2001
From: Kyle <kyle.c.klenk@gmail.com>
Date: Fri, 9 Sep 2022 18:59:12 +0000
Subject: [PATCH] added new derivatives
---
build/source/dshare/get_ixname.f90 | 3 +++
build/source/dshare/popMetadat.f90 | 5 ++++-
build/source/dshare/var_lookup.f90 | 15 +++++++++------
3 files changed, 16 insertions(+), 7 deletions(-)
diff --git a/build/source/dshare/get_ixname.f90 b/build/source/dshare/get_ixname.f90
index c5817b9..966b413 100755
--- a/build/source/dshare/get_ixname.f90
+++ b/build/source/dshare/get_ixname.f90
@@ -784,6 +784,9 @@ contains
case('dVolHtCapBulk_dCanWat' ); get_ixderiv = iLookDERIV%dVolHtCapBulk_dCanWat ! derivative in bulk heat capacity w.r.t. volumetric water content
case('dVolHtCapBulk_dTk' ); get_ixderiv = iLookDERIV%dVolHtCapBulk_dTk ! derivative in bulk heat capacity w.r.t. temperature
case('dVolHtCapBulk_dTkCanopy' ); get_ixderiv = iLookDERIV%dVolHtCapBulk_dTkCanopy ! derivative in bulk heat capacity w.r.t. temperature
+ ! derivatives in time
+ case( 'mLayerdTemp_dt' ); get_ixderiv = iLookDERIV%mLayerdTemp_dt ! timestep change in layer temperature
+ case( 'scalarCanopydTemp_dt' ); get_ixderiv = iLookDERIV%scalarCanopydTemp_dt ! timestep change in canopy temperature
case default
get_ixderiv = integerMissing
diff --git a/build/source/dshare/popMetadat.f90 b/build/source/dshare/popMetadat.f90
index a099df9..14dc59b 100755
--- a/build/source/dshare/popMetadat.f90
+++ b/build/source/dshare/popMetadat.f90
@@ -616,7 +616,10 @@ subroutine popMetadat(err,message)
deriv_meta(iLookDERIV%dVolHtCapBulk_dCanWat) = var_info('dVolHtCapBulk_dCanWat' , 'derivative in bulk heat capacity w.r.t. volumetric water content' , 'J m-3 K-1' , get_ixVarType('scalarv'), iMissVec, iMissVec, .false.)
deriv_meta(iLookDERIV%dVolHtCapBulk_dTk) = var_info('dVolHtCapBulk_dTk' , 'derivative in bulk heat capacity w.r.t. temperature' , 'J m-3 K-2' , get_ixVarType('midToto'), iMissVec, iMissVec, .false.)
deriv_meta(iLookDERIV%dVolHtCapBulk_dTkCanopy) = var_info('dVolHtCapBulk_dTkCanopy' , 'derivative in bulk heat capacity w.r.t. temperature' , 'J m-3 K-2' , get_ixVarType('scalarv'), iMissVec, iMissVec, .false.)
-
+ ! derivatives in time
+ deriv_meta(iLookDERIV%mLayerdTemp_dt) = var_info('mLayerdTemp_dt' , 'timestep change in layer temperature' , 'K' ,get_ixVarType('midToto'), iMissVec, iMissVec, .false.)
+ deriv_meta(iLookDERIV%scalarCanopydTemp_dt) = var_info('scalarCanopydTemp_dt' , 'timestep change in canopy temperature' , 'K' ,get_ixVarType('scalarv'), iMissVec, iMissVec, .false.)
+
! -----
! * basin-wide runoff and aquifer fluxes...
! -----------------------------------------
diff --git a/build/source/dshare/var_lookup.f90 b/build/source/dshare/var_lookup.f90
index df8c155..f47a3f1 100755
--- a/build/source/dshare/var_lookup.f90
+++ b/build/source/dshare/var_lookup.f90
@@ -553,10 +553,10 @@ MODULE var_lookup
integer(i4b) :: scalarNetRadiation = integerMissing ! net radiation (W m-2)
endtype iLook_flux
- ! ***********************************************************************************************************
- ! (9) define derivatives
- ! ***********************************************************************************************************
- type, public :: iLook_deriv
+! ***********************************************************************************************************
+! (9) define derivatives
+! ***********************************************************************************************************
+type, public :: iLook_deriv
! derivatives in net vegetation energy fluxes w.r.t. relevant state variables
integer(i4b) :: dCanairNetFlux_dCanairTemp = integerMissing ! derivative in net canopy air space flux w.r.t. canopy air temperature (W m-2 K-1)
integer(i4b) :: dCanairNetFlux_dCanopyTemp = integerMissing ! derivative in net canopy air space flux w.r.t. canopy temperature (W m-2 K-1)
@@ -637,7 +637,10 @@ MODULE var_lookup
integer(i4b) :: dVolHtCapBulk_dCanWat = integerMissing ! derivative in bulk heat capacity w.r.t. volumetric water content
integer(i4b) :: dVolHtCapBulk_dTk = integerMissing ! derivative in bulk heat capacity w.r.t. temperature
integer(i4b) :: dVolHtCapBulk_dTkCanopy = integerMissing ! derivative in bulk heat capacity w.r.t. temperature
- endtype iLook_deriv
+ ! derivatives in time
+ integer(i4b) :: mLayerdTemp_dt = integerMissing ! timestep change in layer temperature
+ integer(i4b) :: scalarCanopydTemp_dt = integerMissing ! timestep change in canopy temperature
+endtype iLook_deriv
! ***********************************************************************************************************
! (10) define model indices
@@ -882,7 +885,7 @@ MODULE var_lookup
31, 32, 33, 34, 35, 36, 37, 38, 39, 40,&
41, 42, 43, 44, 45, 46, 47, 48, 49, 50,&
51, 52, 53, 54, 55, 56, 57, 58, 59, 60,&
- 61, 62, 63, 64, 65)
+ 61, 62, 63, 64, 65, 66, 67)
! named variables: model indices
type(iLook_index), public,parameter :: iLookINDEX =ilook_index ( 1, 2, 3, 4, 5, 6, 7, 8, 9, 10,&
--
GitLab