00001 00002 !===================SUBROUTINE DELLWF===================================== 00003 00004 00005 subroutine dellwf(sib,sib_loc) 00006 00007 use kinds 00008 use sibtype 00009 00010 00011 !======================================================================== 00012 ! 00013 ! Calculation of partial derivatives of canopy and ground radiative 00014 ! heat fluxes with respect to Tc, Tg 00015 ! Here we are doing only the long wave radiative loss, which is the 00016 ! only radiative quantity we are trying to bring to the next time step. 00017 ! 00018 !======================================================================== 00019 00020 !------------------------------INPUT is coming from Netrad------------- 00021 ! 00022 ! dtc4, dtg4, dts4, which are the derivatives of the LW loss 00023 ! 00024 !++++++++++++++++++++++++++++++OUTPUT+++++++++++++++++++++++++++++++++++ 00025 ! 00026 ! LCDTC dLC/dTC 00027 ! LCDTG dLC/dTG 00028 ! LCDTS dLC/dTS 00029 ! LGDTG dLG/dTG 00030 ! LGDTC dLG/dTC 00031 ! LSDTS dLS/dTS 00032 ! LSDTC dLS/dTC 00033 ! 00034 !+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 00035 00036 IMPLICIT none 00037 00038 !---------------------------------------------------------------------- 00039 00040 type(sib_t), intent(inout) :: sib 00041 00042 type(sib_local_vars) ,intent(inout) :: sib_loc 00043 ! variables local to SiB 00044 00045 !---------------------------------------------------------------------- 00046 ! 00047 ! canopy leaves: 00048 sib_loc%lcdtc = 2 * sib_loc%dtc4 * sib_loc%fac1 00049 sib_loc%lcdtg = - sib_loc%dtg4 * sib_loc%fac1 00050 sib_loc%lcdts = - sib_loc%dts4 * sib_loc%fac1 00051 ! 00052 ! ground: 00053 ! 00054 sib_loc%lgdtg = sib_loc%dtg4 00055 sib_loc%lgdtc = - sib_loc%dtc4 * sib_loc%fac1 00056 ! 00057 ! snow: 00058 ! 00059 sib_loc%lsdts = sib_loc%dts4 00060 sib_loc%lsdtc = - sib_loc%dtc4 * sib_loc%fac1 00061 00062 00063 end subroutine dellwf
1.7.1