653 AH_MIN, AHProf_24hr, AH_SLOPE_Cooling, & ! input&inout in alphabetical order
655 alb, AlbMax_DecTr, AlbMax_EveTr, AlbMax_Grass, &
656 AlbMin_DecTr, AlbMin_EveTr, AlbMin_Grass, &
657 alpha_bioCO2, alpha_enh_bioCO2, alt, kdown, avRh, avU1, BaseT, BaseTe, &
658 beta_bioCO2, beta_enh_bioCO2, bldgH, CapMax_dec, CapMin_dec, &
659 chAnOHM, CO2PointSource, cpAnOHM, CRWmax, CRWmin, DayWat, DayWatPer, &
660 DecTreeH, DiagMethod, Diagnose, DRAINRT, &
661 dt_since_start, dqndt, qn_av, dqnsdt, qn_s_av, &
662 EF_umolCO2perJ, emis, EmissionsMethod, EnEF_v_Jkm, endDLS, EveTreeH, FAIBldg, &
663 FAIDecTree, FAIEveTree, FAIMethod, Faut, FcEF_v_kgkm, fcld_obs, FlowChange, &
664 FrFossilFuel_Heat, FrFossilFuel_NonHeat, g_max, g_k, g_q_base, g_q_shape, g_t, g_sm, GDD_id, &
665 GDDFull, Gridiv, gsModel, H_maintain, HDD_id, HumActivity_24hr, &
666 IceFrac, id, Ie_a, Ie_end, Ie_m, Ie_start, imin, &
667 InternalWaterUse_h, &
668 IrrFracPaved, IrrFracBldgs, &
669 IrrFracEveTr, IrrFracDecTr, IrrFracGrass, &
670 IrrFracBSoil, IrrFracWater, &
672 iy, kkAnOHM, Kmax, LAI_id, LAIMax, LAIMin, LAI_obs, &
673 LAIPower, LAIType, lat, lenDay_id, ldown_obs, lng, MaxConductance, MaxFCMetab, MaxQFMetab, &
674 SnowWater, MinFCMetab, MinQFMetab, min_res_bioCO2, &
675 NARP_EMIS_SNOW, NARP_TRANS_SITE, NetRadiationMethod, &
677 n_vegetation_region_urban, &
678 n_stream_sw_urban, n_stream_lw_urban, &
679 sw_dn_direct_frac, air_ext_sw, air_ssa_sw, &
680 veg_ssa_sw, air_ext_lw, air_ssa_lw, veg_ssa_lw, &
681 veg_fsd_const, veg_contact_fraction_const, &
682 ground_albedo_dir_mult_fact, use_sw_direct_albedo, & !input
683 height, building_frac, veg_frac, building_scale, veg_scale, & !input: SPARTACUS
684 alb_roof, emis_roof, alb_wall, emis_wall, &
685 roof_albedo_dir_mult_fact, wall_specular_frac, &
686 OHM_coef, OHMIncQF, OHM_threshSW, &
687 OHM_threshWD, PipeCapacity, PopDensDaytime, &
688 PopDensNighttime, PopProf_24hr, PorMax_dec, PorMin_dec, &
689 Precip, PrecipLimit, PrecipLimitAlb, Press_hPa, &
690 QF0_BEU, Qf_A, Qf_B, Qf_C, &
691 qn1_obs, qs_obs, qf_obs, &
692 RadMeltFact, RAINCOVER, RainMaxRes, resp_a, resp_b, &
693 RoughLenHeatMethod, RoughLenMomMethod, RunoffToWater, S1, S2, &
694 SatHydraulicConduct, SDDFull, SDD_id, SMDMethod, SnowAlb, SnowAlbMax, &
695 SnowAlbMin, SnowPackLimit, SnowDens, SnowDensMax, SnowDensMin, SnowfallCum, SnowFrac, &
696 SnowLimBldg, SnowLimPaved, snowFrac_obs, SnowPack, SnowProf_24hr, SnowUse, SoilDepth, &
697 StabilityMethod, startDLS, &
698 soilstore_surf, SoilStoreCap_surf, state_surf, StateLimit_surf, WetThresh_surf, &
699 soilstore_roof, SoilStoreCap_roof, state_roof, StateLimit_roof, WetThresh_roof, &
700 soilstore_wall, SoilStoreCap_wall, state_wall, StateLimit_wall, WetThresh_wall, &
701 StorageHeatMethod, StoreDrainPrm, SurfaceArea, Tair_av, tau_a, tau_f, tau_r, &
703 BaseT_Cooling, BaseT_Heating, Temp_C, TempMeltFact, TH, &
704 theta_bioCO2, timezone, TL, TrafficRate, TrafficUnits, &
706 tsfc_roof, tsfc_wall, tsfc_surf, &
707 temp_roof, temp_wall, temp_surf, &
708 tin_roof, tin_wall, tin_surf, &
709 k_roof, k_wall, k_surf, &
710 cp_roof, cp_wall, cp_surf, &
711 dz_roof, dz_wall, dz_surf, &
712 TraffProf_24hr, Ts5mindata_ir, tstep, tstep_prev, veg_type, &
713 WaterDist, WaterUseMethod, wu_m3, &
714 WUDay_id, DecidCap_id, albDecTr_id, albEveTr_id, albGrass_id, porosity_id, &
715 WUProfA_24hr, WUProfM_24hr, xsmd, Z, z0m_in, zdm_in, &
727 INTEGER,
PARAMETER :: AerodynamicResistanceMethod = 2
728 INTEGER,
PARAMETER :: BaseTMethod = 2
729 INTEGER,
INTENT(IN) :: Diagnose
730 INTEGER,
PARAMETER :: DiagQN = 0
731 INTEGER,
PARAMETER :: DiagQS = 0
732 INTEGER,
INTENT(IN) :: startDLS
733 INTEGER,
INTENT(IN) :: endDLS
734 INTEGER,
INTENT(IN) :: EmissionsMethod
735 INTEGER,
INTENT(IN) :: Gridiv
736 INTEGER,
INTENT(IN) :: nlayer
737 INTEGER,
INTENT(IN) :: gsModel
738 INTEGER,
INTENT(IN) :: id
739 INTEGER,
INTENT(IN) :: Ie_end
740 INTEGER,
INTENT(IN) :: Ie_start
741 INTEGER,
INTENT(IN) :: isec
742 INTEGER,
INTENT(IN) :: imin
743 INTEGER,
INTENT(IN) :: it
744 INTEGER,
PARAMETER :: EvapMethod = 2
745 INTEGER,
INTENT(IN) :: iy
746 INTEGER,
PARAMETER :: LAImethod = 1
747 INTEGER,
INTENT(IN) :: NetRadiationMethod
748 INTEGER,
INTENT(IN) :: OHMIncQF
749 INTEGER,
INTENT(IN) :: RoughLenHeatMethod
750 INTEGER,
INTENT(IN) :: RoughLenMomMethod
751 INTEGER,
INTENT(IN) :: FAIMethod
752 INTEGER,
INTENT(IN) :: SMDMethod
753 INTEGER,
INTENT(IN) :: SnowUse
754 INTEGER,
INTENT(IN) :: StabilityMethod
755 INTEGER,
INTENT(IN) :: StorageHeatMethod
756 INTEGER,
INTENT(in) :: DiagMethod
757 INTEGER,
INTENT(IN) :: tstep
758 INTEGER,
INTENT(IN) :: tstep_prev
759 INTEGER,
INTENT(in) :: dt_since_start
760 INTEGER,
INTENT(IN) :: veg_type
761 INTEGER,
INTENT(IN) :: WaterUseMethod
763 REAL(KIND(1D0)),
INTENT(IN) :: AlbMax_DecTr
764 REAL(KIND(1D0)),
INTENT(IN) :: AlbMax_EveTr
765 REAL(KIND(1D0)),
INTENT(IN) :: AlbMax_Grass
766 REAL(KIND(1D0)),
INTENT(IN) :: AlbMin_DecTr
767 REAL(KIND(1D0)),
INTENT(IN) :: AlbMin_EveTr
768 REAL(KIND(1D0)),
INTENT(IN) :: AlbMin_Grass
769 REAL(KIND(1D0)),
INTENT(IN) :: alt
770 REAL(KIND(1D0)),
INTENT(IN) :: kdown
771 REAL(KIND(1D0)),
INTENT(IN) :: avRh
772 REAL(KIND(1D0)),
INTENT(IN) :: avU1
773 REAL(KIND(1D0)),
PARAMETER :: BaseT_HC = 18.2
774 REAL(KIND(1D0)),
INTENT(IN) :: bldgH
775 REAL(KIND(1D0)),
INTENT(IN) :: CapMax_dec
776 REAL(KIND(1D0)),
INTENT(IN) :: CapMin_dec
777 REAL(KIND(1D0)),
INTENT(IN) :: CO2PointSource
778 REAL(KIND(1D0)),
INTENT(IN) :: CRWmax
779 REAL(KIND(1D0)),
INTENT(IN) :: CRWmin
780 REAL(KIND(1D0)),
INTENT(IN) :: DecTreeH
781 REAL(KIND(1D0)),
INTENT(IN) :: DRAINRT
782 REAL(KIND(1D0)),
INTENT(IN) :: EF_umolCO2perJ
783 REAL(KIND(1D0)),
INTENT(IN) :: EnEF_v_Jkm
784 REAL(KIND(1D0)),
INTENT(IN) :: EveTreeH
785 REAL(KIND(1D0)),
INTENT(IN) :: FAIBldg
786 REAL(KIND(1D0)),
INTENT(IN) :: FAIDecTree
787 REAL(KIND(1D0)),
INTENT(IN) :: FAIEveTree
788 REAL(KIND(1D0)),
INTENT(IN) :: Faut
789 REAL(KIND(1D0)),
INTENT(IN) :: fcld_obs
790 REAL(KIND(1D0)),
INTENT(IN) :: FlowChange
791 REAL(KIND(1D0)),
INTENT(IN) :: FrFossilFuel_Heat
792 REAL(KIND(1D0)),
INTENT(IN) :: FrFossilFuel_NonHeat
793 REAL(KIND(1D0)),
INTENT(IN) :: g_max
794 REAL(KIND(1D0)),
INTENT(IN) :: g_k
795 REAL(KIND(1D0)),
INTENT(IN) :: g_q_base
796 REAL(KIND(1D0)),
INTENT(IN) :: g_q_shape
797 REAL(KIND(1D0)),
INTENT(IN) :: g_t
798 REAL(KIND(1D0)),
INTENT(IN) :: g_sm
799 REAL(KIND(1D0)),
INTENT(IN) :: H_maintain
800 REAL(KIND(1D0)),
INTENT(IN) :: InternalWaterUse_h
801 REAL(KIND(1D0)),
INTENT(IN) :: IrrFracPaved
802 REAL(KIND(1D0)),
INTENT(IN) :: IrrFracBldgs
803 REAL(KIND(1D0)),
INTENT(IN) :: IrrFracDecTr
804 REAL(KIND(1D0)),
INTENT(IN) :: IrrFracEveTr
805 REAL(KIND(1D0)),
INTENT(IN) :: IrrFracGrass
806 REAL(KIND(1D0)),
INTENT(IN) :: IrrFracBSoil
807 REAL(KIND(1D0)),
INTENT(IN) :: IrrFracWater
808 REAL(KIND(1D0)),
INTENT(IN) :: Kmax
809 REAL(KIND(1D0)),
INTENT(IN) :: LAI_obs
810 REAL(KIND(1D0)),
INTENT(IN) :: lat
811 REAL(KIND(1D0)),
INTENT(IN) :: ldown_obs
812 REAL(KIND(1D0)),
INTENT(IN) :: lng
813 REAL(KIND(1D0)),
INTENT(IN) :: MaxFCMetab
814 REAL(KIND(1D0)),
INTENT(IN) :: MaxQFMetab
815 REAL(KIND(1D0)),
INTENT(IN) :: MinFCMetab
816 REAL(KIND(1D0)),
INTENT(IN) :: MinQFMetab
817 REAL(KIND(1D0)),
INTENT(IN) :: NARP_EMIS_SNOW
818 REAL(KIND(1D0)),
INTENT(IN) :: NARP_TRANS_SITE
819 REAL(KIND(1D0)),
INTENT(IN) :: PipeCapacity
820 REAL(KIND(1D0)),
INTENT(IN) :: PopDensNighttime
821 REAL(KIND(1D0)),
INTENT(IN) :: PorMax_dec
822 REAL(KIND(1D0)),
INTENT(IN) :: PorMin_dec
823 REAL(KIND(1D0)),
INTENT(IN) :: Precip
824 REAL(KIND(1D0)),
INTENT(IN) :: PrecipLimit
825 REAL(KIND(1D0)),
INTENT(IN) :: PrecipLimitAlb
826 REAL(KIND(1D0)),
INTENT(IN) :: Press_hPa
827 REAL(KIND(1D0)),
INTENT(IN) :: qn1_obs
828 REAL(KIND(1D0)),
INTENT(IN) :: qs_obs
829 REAL(KIND(1D0)),
INTENT(IN) :: qf_obs
830 REAL(KIND(1D0)),
INTENT(IN) :: RadMeltFact
831 REAL(KIND(1D0)),
INTENT(IN) :: RAINCOVER
832 REAL(KIND(1D0)),
INTENT(IN) :: RainMaxRes
833 REAL(KIND(1D0)),
INTENT(IN) :: RunoffToWater
834 REAL(KIND(1D0)),
INTENT(IN) :: S1
835 REAL(KIND(1D0)),
INTENT(IN) :: S2
836 REAL(KIND(1D0)),
INTENT(IN) :: SnowAlbMax
837 REAL(KIND(1D0)),
INTENT(IN) :: SnowAlbMin
838 REAL(KIND(1D0)),
INTENT(IN) :: SnowDensMax
839 REAL(KIND(1D0)),
INTENT(IN) :: SnowDensMin
840 REAL(KIND(1D0)),
INTENT(IN) :: SnowLimBldg
841 REAL(KIND(1D0)),
INTENT(IN) :: SnowLimPaved
842 REAL(KIND(1D0)),
INTENT(IN) :: snowFrac_obs
843 REAL(KIND(1D0)),
INTENT(IN) :: SurfaceArea
844 REAL(KIND(1D0)),
INTENT(IN) :: tau_a
845 REAL(KIND(1D0)),
INTENT(IN) :: tau_f
846 REAL(KIND(1D0)),
INTENT(IN) :: tau_r
847 REAL(KIND(1D0)),
INTENT(IN) :: Temp_C
848 REAL(KIND(1D0)),
INTENT(IN) :: TempMeltFact
849 REAL(KIND(1D0)),
INTENT(IN) :: TH
850 REAL(KIND(1D0)),
INTENT(IN) :: timezone
851 REAL(KIND(1D0)),
INTENT(IN) :: TL
852 REAL(KIND(1D0)),
INTENT(IN) :: TrafficUnits
853 REAL(KIND(1D0)),
INTENT(IN) :: wu_m3
854 REAL(KIND(1D0)),
INTENT(IN) :: xsmd
855 REAL(KIND(1D0)),
INTENT(IN) :: Z
856 REAL(KIND(1D0)),
INTENT(IN) :: z0m_in
857 REAL(KIND(1D0)),
INTENT(IN) :: zdm_in
859 INTEGER,
DIMENSION(NVEGSURF),
INTENT(IN) :: LAIType
861 REAL(KIND(1D0)),
DIMENSION(2),
INTENT(IN) :: AH_MIN
862 REAL(KIND(1D0)),
DIMENSION(2),
INTENT(IN) :: AH_SLOPE_Cooling
863 REAL(KIND(1D0)),
DIMENSION(2),
INTENT(IN) :: AH_SLOPE_Heating
864 REAL(KIND(1D0)),
DIMENSION(2),
INTENT(IN) :: FcEF_v_kgkm
865 REAL(KIND(1D0)),
DIMENSION(2),
INTENT(IN) :: QF0_BEU
866 REAL(KIND(1D0)),
DIMENSION(2),
INTENT(IN) :: Qf_A
867 REAL(KIND(1D0)),
DIMENSION(2),
INTENT(IN) :: Qf_B
868 REAL(KIND(1D0)),
DIMENSION(2),
INTENT(IN) :: Qf_C
869 REAL(KIND(1D0)),
DIMENSION(2),
INTENT(IN) :: PopDensDaytime
870 REAL(KIND(1D0)),
DIMENSION(2),
INTENT(IN) :: BaseT_Cooling
871 REAL(KIND(1D0)),
DIMENSION(2),
INTENT(IN) :: BaseT_Heating
872 REAL(KIND(1D0)),
DIMENSION(2),
INTENT(IN) :: TrafficRate
873 REAL(KIND(1D0)),
DIMENSION(3),
INTENT(IN) :: Ie_a
874 REAL(KIND(1D0)),
DIMENSION(3),
INTENT(IN) :: Ie_m
875 REAL(KIND(1D0)),
DIMENSION(3),
INTENT(IN) :: MaxConductance
876 REAL(KIND(1D0)),
DIMENSION(7),
INTENT(IN) :: DayWat
877 REAL(KIND(1D0)),
DIMENSION(7),
INTENT(IN) :: DayWatPer
878 REAL(KIND(1D0)),
DIMENSION(nsurf + 1),
INTENT(IN) :: OHM_threshSW
879 REAL(KIND(1D0)),
DIMENSION(nsurf + 1),
INTENT(IN) :: OHM_threshWD
880 REAL(KIND(1D0)),
DIMENSION(NSURF),
INTENT(IN) :: chAnOHM
881 REAL(KIND(1D0)),
DIMENSION(NSURF),
INTENT(IN) :: cpAnOHM
882 REAL(KIND(1D0)),
DIMENSION(NSURF),
INTENT(IN) :: emis
883 REAL(KIND(1D0)),
DIMENSION(NSURF),
INTENT(IN) :: kkAnOHM
884 REAL(KIND(1D0)),
DIMENSION(NSURF),
INTENT(IN) :: SatHydraulicConduct
885 REAL(KIND(1D0)),
DIMENSION(NSURF),
INTENT(IN) :: sfr_surf
886 REAL(KIND(1D0)),
DIMENSION(NSURF),
INTENT(IN) :: SnowPackLimit
887 REAL(KIND(1D0)),
DIMENSION(NSURF),
INTENT(IN) :: SoilDepth
888 REAL(KIND(1D0)),
DIMENSION(NSURF),
INTENT(IN) :: SoilStoreCap_surf
889 REAL(KIND(1D0)),
DIMENSION(NSURF),
INTENT(IN) :: StateLimit_surf
890 REAL(KIND(1D0)),
DIMENSION(NSURF),
INTENT(IN) :: WetThresh_surf
891 REAL(KIND(1D0)),
DIMENSION(NVEGSURF),
INTENT(IN) :: alpha_bioCO2
892 REAL(KIND(1D0)),
DIMENSION(NVEGSURF),
INTENT(IN) :: alpha_enh_bioCO2
893 REAL(KIND(1D0)),
DIMENSION(NVEGSURF),
INTENT(IN) :: BaseT
894 REAL(KIND(1D0)),
DIMENSION(NVEGSURF),
INTENT(IN) :: BaseTe
895 REAL(KIND(1D0)),
DIMENSION(NVEGSURF),
INTENT(IN) :: beta_bioCO2
896 REAL(KIND(1D0)),
DIMENSION(NVEGSURF),
INTENT(IN) :: beta_enh_bioCO2
897 REAL(KIND(1D0)),
DIMENSION(NVEGSURF),
INTENT(IN) :: GDDFull
898 REAL(KIND(1D0)),
DIMENSION(NVEGSURF),
INTENT(IN) :: LAIMax
899 REAL(KIND(1D0)),
DIMENSION(NVEGSURF),
INTENT(IN) :: LAIMin
900 REAL(KIND(1D0)),
DIMENSION(NVEGSURF),
INTENT(IN) :: min_res_bioCO2
901 REAL(KIND(1D0)),
DIMENSION(NVEGSURF),
INTENT(IN) :: resp_a
902 REAL(KIND(1D0)),
DIMENSION(NVEGSURF),
INTENT(IN) :: resp_b
903 REAL(KIND(1D0)),
DIMENSION(NVEGSURF),
INTENT(IN) :: SDDFull
904 REAL(KIND(1D0)),
DIMENSION(0:23, 2),
INTENT(IN) :: SnowProf_24hr
905 REAL(KIND(1D0)),
DIMENSION(NVEGSURF),
INTENT(IN) :: theta_bioCO2
906 REAL(KIND(1D0)),
DIMENSION(4, NVEGSURF),
INTENT(IN) :: LAIPower
907 REAL(KIND(1D0)),
DIMENSION(nsurf + 1, 4, 3),
INTENT(IN) :: OHM_coef
908 REAL(KIND(1D0)),
DIMENSION(NSURF + 1, NSURF - 1),
INTENT(IN) :: WaterDist
909 REAL(KIND(1D0)),
DIMENSION(:),
INTENT(IN) :: Ts5mindata_ir
910 REAL(KIND(1D0)),
DIMENSION(10) :: MetForcingData_grid
913 REAL(KIND(1D0)),
DIMENSION(0:23, 2),
INTENT(IN) :: AHProf_24hr
914 REAL(KIND(1D0)),
DIMENSION(0:23, 2),
INTENT(IN) :: HumActivity_24hr
915 REAL(KIND(1D0)),
DIMENSION(0:23, 2),
INTENT(IN) :: PopProf_24hr
916 REAL(KIND(1D0)),
DIMENSION(0:23, 2),
INTENT(IN) :: TraffProf_24hr
917 REAL(KIND(1D0)),
DIMENSION(0:23, 2),
INTENT(IN) :: WUProfA_24hr
918 REAL(KIND(1D0)),
DIMENSION(0:23, 2),
INTENT(IN) :: WUProfM_24hr
922 REAL(KIND(1D0)),
DIMENSION(nlayer),
INTENT(IN) :: SoilStoreCap_roof
923 REAL(KIND(1D0)),
DIMENSION(nlayer),
INTENT(IN) :: StateLimit_roof
924 REAL(KIND(1D0)),
DIMENSION(nlayer),
INTENT(IN) :: wetthresh_roof
925 REAL(KIND(1D0)),
DIMENSION(nlayer),
INTENT(INOUT) :: soilstore_roof
926 REAL(KIND(1D0)),
DIMENSION(nlayer),
INTENT(INOUT) :: state_roof
928 REAL(KIND(1D0)),
DIMENSION(nlayer),
INTENT(IN) :: SoilStoreCap_wall
929 REAL(KIND(1D0)),
DIMENSION(nlayer),
INTENT(IN) :: StateLimit_wall
930 REAL(KIND(1D0)),
DIMENSION(nlayer),
INTENT(IN) :: wetthresh_wall
931 REAL(KIND(1D0)),
DIMENSION(nlayer),
INTENT(INOUT) :: soilstore_wall
932 REAL(KIND(1D0)),
DIMENSION(nlayer),
INTENT(INOUT) :: state_wall
939 REAL(KIND(1D0)),
INTENT(INOUT) :: qn_av
940 REAL(KIND(1D0)),
INTENT(INOUT) :: dqndt
941 REAL(KIND(1D0)),
INTENT(INOUT) :: qn_s_av
942 REAL(KIND(1D0)),
INTENT(INOUT) :: dqnsdt
945 REAL(KIND(1D0)),
INTENT(INOUT) :: SnowfallCum
946 REAL(KIND(1D0)),
INTENT(INOUT) :: SnowAlb
947 REAL(KIND(1D0)),
DIMENSION(NSURF),
INTENT(INOUT) :: IceFrac
948 REAL(KIND(1D0)),
DIMENSION(NSURF),
INTENT(INOUT) :: SnowWater
949 REAL(KIND(1D0)),
DIMENSION(NSURF),
INTENT(INOUT) :: SnowDens
950 REAL(KIND(1D0)),
DIMENSION(NSURF),
INTENT(INOUT) :: SnowFrac
951 REAL(KIND(1D0)),
DIMENSION(NSURF),
INTENT(INOUT) :: SnowPack
954 REAL(KIND(1D0)),
DIMENSION(NSURF),
INTENT(INOUT) :: soilstore_surf
955 REAL(KIND(1D0)),
DIMENSION(NSURF),
INTENT(INOUT) :: state_surf
956 REAL(KIND(1D0)),
DIMENSION(6, NSURF),
INTENT(INOUT) :: StoreDrainPrm
959 REAL(KIND(1D0)),
DIMENSION(NSURF),
INTENT(INOUT) :: alb
960 REAL(KIND(1D0)),
DIMENSION(nvegsurf),
INTENT(INOUT) :: GDD_id
961 REAL(KIND(1D0)),
DIMENSION(nvegsurf),
INTENT(INout) :: SDD_id
962 REAL(KIND(1D0)),
DIMENSION(nvegsurf),
INTENT(INOUT) :: LAI_id
963 REAL(KIND(1D0)),
INTENT(INout) :: Tmin_id
964 REAL(KIND(1D0)),
INTENT(INout) :: Tmax_id
965 REAL(KIND(1D0)),
INTENT(INout) :: lenDay_id
966 REAL(KIND(1D0)),
INTENT(INOUT) :: DecidCap_id
967 REAL(KIND(1D0)),
INTENT(INOUT) :: albDecTr_id
968 REAL(KIND(1D0)),
INTENT(INOUT) :: albEveTr_id
969 REAL(KIND(1D0)),
INTENT(INOUT) :: albGrass_id
970 REAL(KIND(1D0)),
INTENT(INOUT) :: porosity_id
973 REAL(KIND(1D0)),
DIMENSION(12),
INTENT(INOUT) :: HDD_id
976 REAL(KIND(1D0)),
DIMENSION(9),
INTENT(INOUT) :: WUDay_id
979 REAL(KIND(1D0)),
INTENT(INOUT) :: Tair_av
982 REAL(KIND(1D0)),
DIMENSION(nlayer, ndepth),
INTENT(INOUT) :: temp_roof
983 REAL(KIND(1D0)),
DIMENSION(nlayer, ndepth),
INTENT(INOUT) :: temp_wall
984 REAL(KIND(1D0)),
DIMENSION(nsurf, ndepth),
INTENT(INOUT) :: temp_surf
986 REAL(KIND(1D0)),
DIMENSION(nlayer),
INTENT(INOUT) :: tsfc_roof
987 REAL(KIND(1D0)),
DIMENSION(nlayer),
INTENT(INOUT) :: tsfc_wall
988 REAL(KIND(1D0)),
DIMENSION(nsurf),
INTENT(INOUT) :: tsfc_surf
991 INTEGER,
INTENT(IN) :: n_vegetation_region_urban
992 INTEGER,
INTENT(IN) :: n_stream_sw_urban
993 INTEGER,
INTENT(IN) :: n_stream_lw_urban
994 REAL(KIND(1D0)),
INTENT(IN) :: sw_dn_direct_frac
995 REAL(KIND(1D0)),
INTENT(IN) :: air_ext_sw
996 REAL(KIND(1D0)),
INTENT(IN) :: air_ssa_sw
997 REAL(KIND(1D0)),
INTENT(IN) :: veg_ssa_sw
998 REAL(KIND(1D0)),
INTENT(IN) :: air_ext_lw
999 REAL(KIND(1D0)),
INTENT(IN) :: air_ssa_lw
1000 REAL(KIND(1D0)),
INTENT(IN) :: veg_ssa_lw
1001 REAL(KIND(1D0)),
INTENT(IN) :: veg_fsd_const
1002 REAL(KIND(1D0)),
INTENT(IN) :: veg_contact_fraction_const
1003 REAL(KIND(1D0)),
INTENT(IN) :: ground_albedo_dir_mult_fact
1009 REAL(KIND(1D0)),
DIMENSION(5) :: datetimeLine
1010 REAL(KIND(1D0)),
DIMENSION(ncolumnsDataOutSUEWS - 5) :: dataOutLineSUEWS
1011 REAL(KIND(1D0)),
DIMENSION(ncolumnsDataOutSnow - 5) :: dataOutLineSnow
1012 REAL(KIND(1D0)),
DIMENSION(ncolumnsDataOutESTM - 5) :: dataOutLineESTM
1013 REAL(KIND(1D0)),
DIMENSION(ncolumnsDataOutEHC - 5) :: dataOutLineEHC
1014 REAL(KIND(1D0)),
DIMENSION(ncolumnsDataOutRSL - 5) :: dataoutLineRSL
1015 REAL(KIND(1D0)),
DIMENSION(ncolumnsDataOutBEERS - 5) :: dataOutLineBEERS
1016 REAL(KIND(1D0)),
DIMENSION(ncolumnsDataOutDebug - 5) :: dataOutLineDebug
1017 REAL(KIND(1D0)),
DIMENSION(ncolumnsDataOutSPARTACUS - 5) :: dataOutLineSPARTACUS
1018 REAL(KIND(1D0)),
DIMENSION(ncolumnsDataOutDailyState - 5) :: dataOutLineDailyState
1020 TYPE(
output_line),
INTENT(OUT) :: output_line_suews
1025 REAL(KIND(1D0)) :: a1
1026 REAL(KIND(1D0)) :: a2
1027 REAL(KIND(1D0)) :: a3
1028 REAL(KIND(1D0)) :: AdditionalWater
1029 REAL(KIND(1D0)) :: U10_ms
1030 REAL(KIND(1D0)) :: azimuth
1031 REAL(KIND(1D0)) :: chSnow_per_interval
1033 REAL(KIND(1D0)) :: dens_dry
1034 REAL(KIND(1D0)) :: deltaLAI
1035 REAL(KIND(1D0)) :: drain_per_tstep
1036 REAL(KIND(1D0)) :: Ea_hPa
1037 REAL(KIND(1D0)) :: QE_LUMPS
1038 REAL(KIND(1D0)) :: es_hPa
1039 REAL(KIND(1D0)) :: ev_per_tstep
1040 REAL(KIND(1D0)) :: wu_ext
1041 REAL(KIND(1D0)) :: Fc
1042 REAL(KIND(1D0)) :: Fc_anthro
1043 REAL(KIND(1D0)) :: Fc_biogen
1044 REAL(KIND(1D0)) :: Fc_build
1045 REAL(KIND(1D0)) :: fcld
1046 REAL(KIND(1D0)) :: Fc_metab
1047 REAL(KIND(1D0)) :: Fc_photo
1048 REAL(KIND(1D0)) :: Fc_point
1049 REAL(KIND(1D0)) :: Fc_respi
1050 REAL(KIND(1D0)) :: Fc_traff
1051 REAL(KIND(1D0)) :: gfunc
1052 REAL(KIND(1D0)) :: gsc
1053 REAL(KIND(1D0)) :: QH_LUMPS
1054 REAL(KIND(1D0)) :: wu_int
1055 REAL(KIND(1D0)) :: kclear
1056 REAL(KIND(1D0)) :: kup
1057 REAL(KIND(1D0)) :: ldown
1058 REAL(KIND(1D0)) :: lup
1059 REAL(KIND(1D0)) :: L_mod
1060 REAL(KIND(1D0)) :: mwh
1061 REAL(KIND(1D0)) :: mwstore
1062 REAL(KIND(1D0)) :: NWstate_per_tstep
1063 REAL(KIND(1D0)) :: FAI
1064 REAL(KIND(1D0)) :: PAI
1065 REAL(KIND(1D0)) :: zL
1066 REAL(KIND(1D0)) :: q2_gkg
1067 REAL(KIND(1D0)) :: qe
1068 REAL(KIND(1D0)) :: qf
1069 REAL(KIND(1D0)) :: QF_SAHP
1070 REAL(KIND(1D0)) :: qh
1071 REAL(KIND(1D0)) :: qh_residual
1072 REAL(KIND(1D0)) :: qh_resist
1073 REAL(KIND(1D0)) :: Qm
1074 REAL(KIND(1D0)) :: QmFreez
1075 REAL(KIND(1D0)) :: QmRain
1076 REAL(KIND(1D0)) :: qn
1077 REAL(KIND(1D0)) :: qn_snow
1078 REAL(KIND(1D0)) :: qn_snowfree
1079 REAL(KIND(1D0)) :: qs
1080 REAL(KIND(1D0)) :: RA_h
1081 REAL(KIND(1D0)) :: RS
1082 REAL(KIND(1D0)),
DIMENSION(NSURF) :: rss_surf
1083 REAL(KIND(1D0)) :: RH2
1084 REAL(KIND(1D0)) :: runoffAGveg
1085 REAL(KIND(1D0)) :: runoffAGimpervious
1086 REAL(KIND(1D0)) :: runoff_per_tstep
1087 REAL(KIND(1D0)) :: runoffPipes
1088 REAL(KIND(1D0)) :: runoffSoil_per_tstep
1089 REAL(KIND(1D0)) :: runoffwaterbody
1090 REAL(KIND(1D0)) :: smd
1091 REAL(KIND(1D0)) :: SoilState
1092 REAL(KIND(1D0)) :: state_per_tstep
1093 REAL(KIND(1D0)) :: surf_chang_per_tstep
1094 REAL(KIND(1D0)) :: swe
1095 REAL(KIND(1D0)) :: t2_C
1096 REAL(KIND(1D0)) :: TSfc_C
1097 REAL(KIND(1D0)) :: TempVeg
1098 REAL(KIND(1D0)) :: tot_chang_per_tstep
1099 REAL(KIND(1D0)) :: TStar
1100 REAL(KIND(1D0)) :: tsurf
1101 REAL(KIND(1D0)) :: UStar
1102 REAL(KIND(1D0)) :: VPD_Pa
1103 REAL(KIND(1D0)) :: z0m
1104 REAL(KIND(1D0)) :: zdm
1105 REAL(KIND(1D0)) :: ZENITH_deg
1106 REAL(KIND(1D0)) :: zH
1108 REAL(KIND(1D0)),
DIMENSION(2) :: SnowRemoval
1109 REAL(KIND(1D0)),
DIMENSION(NSURF) :: wu_surf
1111 REAL(KIND(1D0)),
DIMENSION(nsurf) :: kup_ind_snow
1116 REAL(KIND(1D0)),
DIMENSION(NSURF) :: qn_ind_snow
1118 REAL(KIND(1D0)),
DIMENSION(NSURF) :: runoffSoil
1119 REAL(KIND(1D0)),
DIMENSION(NSURF) :: smd_nsurf
1122 REAL(KIND(1D0)),
DIMENSION(nsurf) :: Tsurf_ind_snow
1125 INTEGER,
DIMENSION(3) :: dayofWeek_id
1128 REAL(KIND(1D0)) :: dq
1129 REAL(KIND(1D0)) :: lvS_J_kg
1130 REAL(KIND(1D0)) :: psyc_hPa
1131 REAL(KIND(1D0)) :: z0v
1132 REAL(KIND(1D0)) :: z0vSnow
1133 REAL(KIND(1D0)) :: RAsnow
1134 REAL(KIND(1D0)) :: RB
1135 REAL(KIND(1D0)) :: runoff_per_interval
1136 REAL(KIND(1D0)) :: s_hPa
1137 REAL(KIND(1D0)) :: sIce_hpa
1138 REAL(KIND(1D0)) :: SoilMoistCap
1140 REAL(KIND(1D0)) :: VegPhenLumps
1141 REAL(KIND(1D0)) :: VPd_hpa
1142 REAL(KIND(1D0)) :: vsmd
1143 REAL(KIND(1D0)) :: ZZD
1145 REAL(KIND(1D0)),
DIMENSION(NSURF) :: deltaQi
1146 REAL(KIND(1D0)),
DIMENSION(NSURF) :: drain_surf
1149 REAL(KIND(1D0)),
DIMENSION(NSURF) :: tsurf_ind
1154 REAL(KIND(1D0)) :: addImpervious = 0
1155 REAL(KIND(1D0)) :: addPipes = 0
1156 REAL(KIND(1D0)) :: addVeg = 0
1157 REAL(KIND(1D0)) :: addWaterBody = 0
1158 REAL(KIND(1D0)),
DIMENSION(NSURF) :: AddWater = 0
1159 REAL(KIND(1D0)),
DIMENSION(NSURF) :: frac_water2runoff = 0
1163 REAL(KIND(1D0)) :: nsh_real
1164 REAL(KIND(1D0)) :: tstep_real
1165 REAL(KIND(1D0)) :: dectime
1168 REAL(KIND(1D0)) :: VegFraction
1169 REAL(KIND(1D0)) :: ImpervFraction
1170 REAL(KIND(1D0)) :: PervFraction
1171 REAL(KIND(1D0)) :: NonWaterFraction
1174 REAL(KIND(1D0)) :: albedo_snow
1184 REAL(KIND(1D0)) :: qn_av_prev, qn_av_next
1185 REAL(KIND(1D0)) :: dqndt_prev, dqndt_next
1186 REAL(KIND(1D0)) :: qn_s_av_prev, qn_s_av_next
1187 REAL(KIND(1D0)) :: dqnsdt_prev, dqnsdt_next
1190 REAL(KIND(1D0)) :: SnowfallCum_prev, SnowfallCum_next
1191 REAL(KIND(1D0)) :: SnowAlb_prev, SnowAlb_next
1193 REAL(KIND(1D0)),
DIMENSION(NSURF) :: IceFrac_prev, IceFrac_next
1194 REAL(KIND(1D0)),
DIMENSION(NSURF) :: SnowWater_prev, SnowWater_next
1195 REAL(KIND(1D0)),
DIMENSION(NSURF) :: SnowDens_prev, SnowDens_next
1196 REAL(KIND(1D0)),
DIMENSION(NSURF) :: SnowFrac_prev, SnowFrac_next
1197 REAL(KIND(1D0)),
DIMENSION(NSURF) :: SnowPack_prev, SnowPack_next
1200 REAL(KIND(1D0)),
DIMENSION(NSURF) :: soilstore_surf_prev, soilstore_surf_next
1201 REAL(KIND(1D0)),
DIMENSION(nlayer) :: soilstore_roof_prev, soilstore_roof_next
1202 REAL(KIND(1D0)),
DIMENSION(nlayer) :: soilstore_wall_prev, soilstore_wall_next
1203 REAL(KIND(1D0)),
DIMENSION(NSURF) :: state_surf_prev, state_surf_next
1204 REAL(KIND(1D0)),
DIMENSION(nlayer) :: state_roof_prev, state_roof_next
1205 REAL(KIND(1D0)),
DIMENSION(nlayer) :: state_wall_prev, state_wall_next
1206 REAL(KIND(1D0)),
DIMENSION(NSURF) :: ev0_surf
1207 REAL(KIND(1D0)),
DIMENSION(NSURF) :: ev_surf
1208 REAL(KIND(1D0)),
DIMENSION(nlayer) :: ev_roof
1209 REAL(KIND(1D0)),
DIMENSION(nlayer) :: ev_wall
1210 REAL(KIND(1D0)),
DIMENSION(6, NSURF) :: StoreDrainPrm_prev, StoreDrainPrm_next
1213 REAL(KIND(1D0)),
DIMENSION(NSURF) :: alb_prev, alb_next
1214 REAL(KIND(1D0)),
DIMENSION(nvegsurf) :: GDD_id_prev, GDD_id_next
1215 REAL(KIND(1D0)),
DIMENSION(nvegsurf) :: LAI_id_prev, LAI_id_next
1216 REAL(KIND(1D0)),
DIMENSION(nvegsurf) :: SDD_id_prev, SDD_id_next
1218 REAL(KIND(1D0)) :: DecidCap_id_prev, DecidCap_id_next
1219 REAL(KIND(1D0)) :: albDecTr_id_prev, albDecTr_id_next
1220 REAL(KIND(1D0)) :: albEveTr_id_prev, albEveTr_id_next
1221 REAL(KIND(1D0)) :: albGrass_id_prev, albGrass_id_next
1222 REAL(KIND(1D0)) :: porosity_id_prev, porosity_id_next
1224 REAL(KIND(1D0)) :: Tmin_id_prev, Tmin_id_next
1225 REAL(KIND(1D0)) :: Tmax_id_prev, Tmax_id_next
1226 REAL(KIND(1D0)) :: lenDay_id_prev, lenDay_id_next
1229 REAL(KIND(1D0)),
DIMENSION(12) :: HDD_id_prev, HDD_id_next
1232 REAL(KIND(1D0)),
DIMENSION(9) :: WUDay_id_prev, WUDay_id_next
1234 REAL(KIND(1D0)) :: Tair_av_prev, Tair_av_next
1238 INTEGER,
PARAMETER :: nz = 90
1241 LOGICAL :: flag_converge
1242 REAL(KIND(1D0)) :: Ts_iter
1243 REAL(KIND(1D0)) :: dif_tsfc_iter
1244 REAL(KIND(1D0)) :: QH_Init
1251 REAL(KIND(1D0)),
DIMENSION(nlayer) :: tsfc_out_roof, tsfc0_out_roof
1252 REAL(KIND(1D0)),
DIMENSION(nlayer),
INTENT(in) :: tin_roof
1253 REAL(KIND(1D0)),
DIMENSION(nlayer) :: sfr_roof
1254 REAL(KIND(1D0)),
DIMENSION(nlayer, ndepth) :: temp_in_roof
1255 REAL(KIND(1D0)),
DIMENSION(nlayer, ndepth),
INTENT(in) :: k_roof
1256 REAL(KIND(1D0)),
DIMENSION(nlayer, ndepth),
INTENT(in) :: cp_roof
1257 REAL(KIND(1D0)),
DIMENSION(nlayer, ndepth),
INTENT(in) :: dz_roof
1259 REAL(KIND(1D0)),
DIMENSION(nlayer) :: tsfc_out_wall, tsfc0_out_wall
1260 REAL(KIND(1D0)),
DIMENSION(nlayer),
INTENT(in) :: tin_wall
1261 REAL(KIND(1D0)),
DIMENSION(nlayer) :: sfr_wall
1262 REAL(KIND(1D0)),
DIMENSION(nlayer, ndepth) :: temp_in_wall
1263 REAL(KIND(1D0)),
DIMENSION(nlayer, ndepth),
INTENT(in) :: k_wall
1264 REAL(KIND(1D0)),
DIMENSION(nlayer, ndepth),
INTENT(in) :: cp_wall
1265 REAL(KIND(1D0)),
DIMENSION(nlayer, ndepth),
INTENT(in) :: dz_wall
1267 REAL(KIND(1D0)),
DIMENSION(nsurf) :: tsfc_out_surf, tsfc0_out_surf
1268 REAL(KIND(1D0)),
DIMENSION(nsurf),
INTENT(in) :: tin_surf
1269 REAL(KIND(1D0)),
DIMENSION(nsurf, ndepth) :: temp_in_surf
1270 REAL(KIND(1D0)),
DIMENSION(nsurf, ndepth),
INTENT(in) :: k_surf
1271 REAL(KIND(1D0)),
DIMENSION(nsurf, ndepth),
INTENT(in) :: cp_surf
1272 REAL(KIND(1D0)),
DIMENSION(nsurf, ndepth),
INTENT(in) :: dz_surf
1278 REAL(KIND(1D0)),
DIMENSION(nlayer) :: QS_roof
1280 REAL(KIND(1D0)),
DIMENSION(nlayer, ndepth) :: temp_out_roof
1283 REAL(KIND(1D0)),
DIMENSION(nlayer) :: QG_roof
1284 REAL(KIND(1D0)),
DIMENSION(nlayer) :: QN_roof
1285 REAL(KIND(1D0)),
DIMENSION(nlayer) :: qe_roof
1286 REAL(KIND(1D0)),
DIMENSION(nlayer) :: qh_roof
1287 REAL(KIND(1D0)),
DIMENSION(nlayer) :: qh_resist_roof
1291 REAL(KIND(1D0)),
DIMENSION(nlayer) :: QS_wall
1293 REAL(KIND(1D0)),
DIMENSION(nlayer, ndepth) :: temp_out_wall
1296 REAL(KIND(1D0)),
DIMENSION(nlayer) :: QG_wall
1297 REAL(KIND(1D0)),
DIMENSION(nlayer) :: QN_wall
1298 REAL(KIND(1D0)),
DIMENSION(nlayer) :: qe_wall
1299 REAL(KIND(1D0)),
DIMENSION(nlayer) :: qh_wall
1300 REAL(KIND(1D0)),
DIMENSION(nlayer) :: qh_resist_wall
1304 REAL(KIND(1D0)),
DIMENSION(nsurf, ndepth) :: temp_out_surf
1307 REAL(KIND(1D0)),
DIMENSION(nsurf) :: QG_surf
1308 REAL(KIND(1D0)),
DIMENSION(nsurf) :: QN_surf
1309 REAL(KIND(1D0)),
DIMENSION(nsurf) :: qs_surf
1310 REAL(KIND(1D0)),
DIMENSION(nsurf) :: qe0_surf
1311 REAL(KIND(1D0)),
DIMENSION(nsurf) :: qe_surf
1312 REAL(KIND(1D0)),
DIMENSION(nsurf) :: qh_surf
1313 REAL(KIND(1D0)),
DIMENSION(nsurf) :: qh_resist_surf
1322 REAL(KIND(1D0)) :: ratio_iter
1324 LOGICAL,
INTENT(IN) :: use_sw_direct_albedo
1326 REAL(KIND(1D0)),
DIMENSION(nlayer + 1),
INTENT(IN) :: height
1327 REAL(KIND(1D0)),
DIMENSION(nlayer),
INTENT(IN) :: building_frac
1328 REAL(KIND(1D0)),
DIMENSION(nlayer),
INTENT(IN) :: veg_frac
1329 REAL(KIND(1D0)),
DIMENSION(nlayer),
INTENT(IN) :: building_scale
1330 REAL(KIND(1D0)),
DIMENSION(nlayer),
INTENT(IN) :: veg_scale
1331 REAL(KIND(1D0)),
DIMENSION(nlayer),
INTENT(IN) :: alb_roof
1332 REAL(KIND(1D0)),
DIMENSION(nlayer),
INTENT(IN) :: emis_roof
1333 REAL(KIND(1D0)),
DIMENSION(nlayer),
INTENT(IN) :: alb_wall
1334 REAL(KIND(1D0)),
DIMENSION(nlayer),
INTENT(IN) :: emis_wall
1335 REAL(KIND(1D0)),
DIMENSION(nspec, nlayer),
INTENT(IN) :: roof_albedo_dir_mult_fact
1336 REAL(KIND(1D0)),
DIMENSION(nspec, nlayer),
INTENT(IN) :: wall_specular_frac
1338 REAL(KIND(1D0)) :: g_kdown
1339 REAL(KIND(1D0)) :: g_dq
1340 REAL(KIND(1D0)) :: g_ta
1341 REAL(KIND(1D0)) :: g_smd
1342 REAL(KIND(1D0)) :: g_lai
1345 REAL(KIND(1D0)) :: FAIBldg_use
1346 REAL(KIND(1D0)) :: FAIEveTree_use
1347 REAL(KIND(1D0)) :: FAIDecTree_use
1354 chsnow_per_interval = 0.
1367 IF (snowuse == 0)
THEN
1380 qn_s_av_prev = qn_s_av
1381 dqnsdt_prev = dqnsdt
1382 snowfallcum_prev = snowfallcum
1383 snowalb_prev = snowalb
1384 icefrac_prev = icefrac
1385 snowwater_prev = snowwater
1386 snowdens_prev = snowdens
1387 snowfrac_prev = merge(snowfrac_obs, snowfrac, netradiationmethod == 0)
1388 snowpack_prev = snowpack
1389 state_surf_prev = state_surf
1390 soilstore_surf_prev = soilstore_surf
1391 IF (storageheatmethod == 5)
THEN
1392 state_roof_prev = state_roof
1393 state_wall_prev = state_wall
1394 soilstore_roof_prev = soilstore_roof
1395 soilstore_wall_prev = soilstore_wall
1397 tair_av_prev = tair_av
1398 lai_id_prev = lai_id
1399 gdd_id_prev = gdd_id
1400 sdd_id_prev = sdd_id
1401 tmin_id_prev = tmin_id
1402 tmax_id_prev = tmax_id
1403 lenday_id_prev = lenday_id
1404 storedrainprm_prev = storedrainprm
1405 decidcap_id_prev = decidcap_id
1406 porosity_id_prev = porosity_id
1408 albdectr_id_prev = albdectr_id
1409 albevetr_id_prev = albevetr_id
1410 albgrass_id_prev = albgrass_id
1411 hdd_id_prev = hdd_id
1412 wuday_id_prev = wuday_id
1416 IF (storageheatmethod == 5)
THEN
1417 temp_in_roof = temp_roof
1418 temp_in_wall = temp_wall
1419 temp_in_surf = temp_surf
1429 qn_s_av_next = qn_s_av
1430 dqnsdt_next = dqnsdt
1431 snowfallcum_next = snowfallcum
1432 snowalb_next = snowalb
1433 icefrac_next = icefrac
1434 snowwater_next = snowwater
1435 snowdens_next = snowdens
1436 snowfrac_next = snowfrac_prev
1437 snowpack_next = snowpack
1438 state_surf_next = state_surf
1439 soilstore_surf_next = soilstore_surf
1441 IF (storageheatmethod == 5)
THEN
1443 soilstore_roof_next = soilstore_roof
1444 soilstore_wall_next = soilstore_wall
1445 state_roof_next = state_roof
1446 state_wall_next = state_wall
1450 tair_av_next = tair_av
1451 lai_id_next = lai_id
1452 gdd_id_next = gdd_id
1453 sdd_id_next = sdd_id
1454 tmin_id_next = tmin_id
1455 tmax_id_next = tmax_id
1456 lenday_id_next = lenday_id
1457 storedrainprm_next = storedrainprm
1458 decidcap_id_next = decidcap_id
1459 porosity_id_next = porosity_id
1461 albdectr_id_next = albdectr_id
1462 albevetr_id_next = albevetr_id
1463 albgrass_id_next = albgrass_id
1464 hdd_id_next = hdd_id
1465 wuday_id_next = wuday_id
1468 dataoutlinesnow = -999.
1469 dataoutlineestm = -999.
1470 dataoutlineehc = -999.
1471 dataoutlinersl = -999.
1472 dataoutlinebeers = -999.
1473 dataoutlinedebug = -999.
1474 dataoutlinespartacus = -999.
1475 dataoutlinedailystate = -999.
1482 flag_converge = .false.
1485 tsfc_out_surf = tsfc_surf
1486 tsfc0_out_surf = tsfc_surf
1488 IF (storageheatmethod == 5 .OR. netradiationmethod > 1000)
THEN
1489 tsfc_out_roof = tsfc_roof
1490 tsfc0_out_roof = tsfc_roof
1491 tsfc_out_wall = tsfc_wall
1492 tsfc0_out_wall = tsfc_wall
1499 DO WHILE ((.NOT. flag_converge) .AND. i_iter < max_iter)
1505 id, it, imin, isec, &
1511 nsh, nsh_real, tstep_real)
1515 storageheatmethod, netradiationmethod, &
1517 building_frac, building_scale, height, &
1518 vegfraction, impervfraction, pervfraction, nonwaterfraction, &
1528 id, startdls, enddls, &
1532 tair_av_next =
cal_tair_av(tair_av_prev, dt_since_start, tstep, temp_c)
1537 IF (diagnose == 1)
WRITE (*, *)
'Calling SUEWS_cal_RoughnessParameters...'
1538 IF (diagnose == 1) print *,
'z0m_in =', z0m_in
1540 roughlenmommethod, faimethod, &
1543 bldgh, evetreeh, dectreeh, &
1544 porosity_id_prev, faibldg, faievetree, faidectree, &
1545 z0m_in, zdm_in, z, &
1551 IF (diagnose == 1)
WRITE (*, *)
'Calling NARP_cal_SunPosition...'
1553 REAL(iy, KIND(1D0)), &
1554 dectime - tstep/2/86400, &
1555 timezone, lat, lng, alt, &
1556 azimuth, zenith_deg)
1559 IF (diagnose == 1)
WRITE (*, *)
'Calling SUEWS_cal_DailyState...'
1561 iy, id, it, imin, isec, tstep, tstep_prev, dt_since_start, dayofweek_id, &
1562 tmin_id_prev, tmax_id_prev, lenday_id_prev, &
1564 waterusemethod, ie_start, ie_end, &
1565 laimethod, laitype, &
1566 nsh_real, kdown, temp_c, precip, baset_hc, &
1567 baset_heating, baset_cooling, &
1568 lat, faut, lai_obs, &
1569 albmax_dectr, albmax_evetr, albmax_grass, &
1570 albmin_dectr, albmin_evetr, albmin_grass, &
1571 capmax_dec, capmin_dec, pormax_dec, pormin_dec, &
1572 ie_a, ie_m, daywatper, daywat, &
1573 baset, basete, gddfull, sddfull, laimin, laimax, laipower, &
1574 decidcap_id_prev, storedrainprm_prev, lai_id_prev, gdd_id_prev, sdd_id_prev, &
1575 albdectr_id_prev, albevetr_id_prev, albgrass_id_prev, porosity_id_prev, &
1577 state_surf_prev, soilstore_surf_prev, soilstorecap_surf, h_maintain, &
1579 tmin_id_next, tmax_id_next, lenday_id_next, &
1580 albdectr_id_next, albevetr_id_next, albgrass_id_next, porosity_id_next, &
1581 decidcap_id_next, storedrainprm_next, lai_id_next, gdd_id_next, sdd_id_next, wuday_id_next)
1584 IF (diagnose == 1)
WRITE (*, *)
'Calling LUMPS_cal_AtmMoist...'
1586 temp_c, press_hpa, avrh, dectime, &
1588 es_hpa, ea_hpa, vpd_hpa, vpd_pa, dq, dens_dry,
avcp,
avdens)
1591 IF (diagnose == 1)
WRITE (*, *)
'Calling SUEWS_update_SoilMoist...'
1594 soilstorecap_surf, sfr_surf, soilstore_surf_prev, &
1595 soilmoistcap, soilstate, &
1598 IF (diagnose == 1)
WRITE (*, *)
'Calling SUEWS_cal_WaterUse...'
1602 wu_m3, surfacearea, sfr_surf, &
1603 irrfracpaved, irrfracbldgs, &
1604 irrfracevetr, irrfracdectr, irrfracgrass, &
1605 irrfracbsoil, irrfracwater, &
1606 dayofweek_id, wuprofa_24hr, wuprofm_24hr, &
1607 internalwateruse_h, hdd_id_next, wuday_id_next, &
1608 waterusemethod, nsh, it, imin, dls, &
1609 wu_surf, wu_int, wu_ext)
1613 ah_min, ahprof_24hr, ah_slope_cooling, ah_slope_heating, co2pointsource, &
1614 dayofweek_id, dls, ef_umolco2perj, emissionsmethod, enef_v_jkm, &
1615 fcef_v_kgkm, frfossilfuel_heat, frfossilfuel_nonheat, hdd_id_next, humactivity_24hr, &
1616 imin, it, maxfcmetab, maxqfmetab, minfcmetab, minqfmetab, &
1617 popdensdaytime, popdensnighttime, popprof_24hr, qf, qf0_beu, qf_a, qf_b, qf_c, &
1618 qf_obs, qf_sahp, surfacearea, baset_cooling, baset_heating, &
1619 temp_c, trafficrate, trafficunits, traffprof_24hr, &
1620 fc_anthro, fc_build, fc_metab, fc_point, fc_traff)
1631 storageheatmethod, netradiationmethod, snowuse, &
1632 tstep, nlayer, snowpack_prev, tau_a, tau_f, snowalbmax, snowalbmin, &
1633 diagnose, ldown_obs, fcld_obs, &
1634 dectime, zenith_deg, ts_iter, kdown, temp_c, avrh, ea_hpa, qn1_obs, &
1635 snowalb_prev, snowfrac_prev, diagqn, &
1636 narp_trans_site, narp_emis_snow, icefrac_prev, &
1637 sfr_surf, sfr_roof, sfr_wall, &
1638 tsfc_out_surf, tsfc_out_roof, tsfc_out_wall, &
1639 emis, alb_prev, albdectr_id_next, albevetr_id_next, albgrass_id_next, &
1641 n_vegetation_region_urban, &
1642 n_stream_sw_urban, n_stream_lw_urban, &
1643 sw_dn_direct_frac, air_ext_sw, air_ssa_sw, &
1644 veg_ssa_sw, air_ext_lw, air_ssa_lw, veg_ssa_lw, &
1645 veg_fsd_const, veg_contact_fraction_const, &
1646 ground_albedo_dir_mult_fact, use_sw_direct_albedo, &
1647 height, building_frac, veg_frac, building_scale, veg_scale, &
1648 alb_roof, emis_roof, alb_wall, emis_wall, &
1649 roof_albedo_dir_mult_fact, wall_specular_frac, &
1650 alb_next, ldown, fcld, &
1651 qn_surf, qn_roof, qn_wall, &
1652 qn, qn_snowfree, qn_snow, kclear, kup, lup, tsurf, &
1653 qn_ind_snow, kup_ind_snow, tsurf_ind_snow, tsurf_ind, &
1654 albedo_snow, snowalb_next, &
1655 dataoutlinespartacus)
1668 IF (i_iter == 1)
THEN
1669 qg_surf = 0.1*qn_surf
1670 qg_roof = 0.1*qn_roof
1671 qg_wall = 0.1*qn_wall
1673 qg_surf = qn_surf + qf - (qh_surf + qe_surf)
1674 qg_roof = qn_roof + qf - (qh_roof + qe_roof)
1675 qg_wall = qn_wall + qf - (qh_wall + qe_wall)
1689 storageheatmethod, qs_obs, ohmincqf, gridiv, &
1690 id, tstep, dt_since_start, diagnose, &
1692 qg_surf, qg_roof, qg_wall, &
1693 tsfc_out_roof, tin_roof, temp_in_roof, k_roof, cp_roof, dz_roof, sfr_roof, &
1694 tsfc_out_wall, tin_wall, temp_in_wall, k_wall, cp_wall, dz_wall, sfr_wall, &
1695 tsfc_out_surf, tin_surf, temp_in_surf, k_surf, cp_surf, dz_surf, sfr_surf, &
1696 ohm_coef, ohm_threshsw, ohm_threshwd, &
1697 soilstore_surf_prev, soilstorecap_surf, state_surf_prev, snowuse, snowfrac_prev, diagqs, &
1698 hdd_id, metforcingdata_grid, ts5mindata_ir, qf, qn, &
1699 kdown, avu1, temp_c, zenith_deg, avrh, press_hpa, ldown, &
1700 bldgh, alb, emis, cpanohm, kkanohm, chanohm, emissionsmethod, &
1701 tair_av, qn_av_prev, dqndt_prev, qn_s_av_prev, dqnsdt_prev, &
1703 qn_snow, dataoutlineestm, qs, &
1704 qn_av_next, dqndt_next, qn_s_av_next, dqnsdt_next, &
1705 deltaqi, a1, a2, a3, &
1706 temp_out_roof, qs_roof, &
1707 temp_out_wall, qs_wall, &
1708 temp_out_surf, qs_surf)
1737 IF (diagnose == 1)
WRITE (*, *)
'Calling MeltHeat'
1740 IF (diagnose == 1)
WRITE (*, *)
'Calling LUMPS_cal_QHQE...'
1741 IF (i_iter == 1)
THEN
1745 snowuse, qn, qf, qs, temp_c, vegfraction,
avcp, press_hpa,
lv_j_kg, &
1746 tstep_real, drainrt, nsh_real, &
1747 precip, rainmaxres, raincover, sfr_surf, lai_id_next, laimax, laimin, &
1749 qe_lumps, psyc_hpa, s_hpa, sice_hpa, tempveg, vegphenlumps)
1761 snowuse, nonwaterfraction, addpipes, addimpervious, addveg, addwaterbody, &
1762 state_surf_prev, sfr_surf, storedrainprm_next, waterdist, nsh_real, &
1764 drain_surf, frac_water2runoff, &
1765 additionalwater, runoffpipes, runoff_per_interval, &
1772 diagnose, aerodynamicresistancemethod, roughlenheatmethod, snowuse, &
1773 id, it, gsmodel, smdmethod, &
1775 avu1, temp_c, vegfraction, kdown, &
1777 g_max, g_k, g_q_base, g_q_shape, &
1778 g_t, g_sm, s1, s2, &
1780 dq, xsmd, vsmd, maxconductance, laimax, lai_id_next, snowfrac_prev, sfr_surf, &
1781 g_kdown, g_dq, g_ta, g_smd, g_lai, &
1782 ustar, tstar, l_mod, &
1783 zl, gsc, rs, ra_h, rasnow, rb, z0v, z0vsnow)
1788 IF (snowuse == 1)
THEN
1793 tstep, imin, it, evapmethod, dayofweek_id, crwmin, crwmax, &
1795 rasnow, psyc_hpa, sice_hpa, tau_r, &
1796 radmeltfact, tempmeltfact, snowalbmax, preciplimit, preciplimitalb, &
1797 qn_ind_snow, kup_ind_snow, deltaqi, tsurf_ind_snow, &
1799 pervfraction, vegfraction, addimpervious, qn_snowfree, qf, qs, vpd_hpa, s_hpa, &
1800 rs, ra_h, rb, snowdensmax, snowdensmin, precip, pipecapacity, runofftowater, &
1801 addveg, snowlimpaved, snowlimbldg, &
1802 flowchange, drain_surf, wetthresh_surf, soilstorecap_surf, &
1803 tsurf_ind, sfr_surf, &
1804 addwater, frac_water2runoff, storedrainprm_next, snowpacklimit, snowprof_24hr, &
1805 snowpack_prev, snowfrac_prev, snowwater_prev, icefrac_prev, snowdens_prev, &
1806 snowfallcum_prev, state_surf_prev, soilstore_surf_prev, &
1809 snowpack_next, snowfrac_next, snowwater_next, icefrac_next, snowdens_next, &
1810 snowfallcum_next, state_surf_next, soilstore_surf_next, &
1811 state_per_tstep, nwstate_per_tstep, &
1812 qe, qe_surf, qe_roof, qe_wall, &
1814 swe, chsnow_per_interval, ev_per_tstep, runoff_per_tstep, &
1815 surf_chang_per_tstep, runoffpipes, mwstore, runoffwaterbody, &
1816 runoffagveg, runoffagimpervious, rss_surf, &
1823 diagnose, storageheatmethod, nlayer, &
1830 qf, vpd_hpa, s_hpa, rs, ra_h, rb, &
1831 precip, pipecapacity, runofftowater, &
1832 nonwaterfraction, wu_surf, addveg, addwaterbody, addwater, &
1833 flowchange, drain_surf, &
1834 frac_water2runoff, storedrainprm_next, &
1835 sfr_surf, statelimit_surf, soilstorecap_surf, wetthresh_surf, &
1836 state_surf_prev, soilstore_surf_prev, qn_surf, qs_surf, &
1837 sfr_roof, statelimit_roof, soilstorecap_roof, wetthresh_roof, &
1838 state_roof_prev, soilstore_roof_prev, qn_roof, qs_roof, &
1839 sfr_wall, statelimit_wall, soilstorecap_wall, wetthresh_wall, &
1840 state_wall_prev, soilstore_wall_prev, qn_wall, qs_wall, &
1841 state_surf_next, soilstore_surf_next, ev_surf, &
1842 state_roof_next, soilstore_roof_next, ev_roof, &
1843 state_wall_next, soilstore_wall_next, ev_wall, &
1844 state_per_tstep, nwstate_per_tstep, &
1845 ev0_surf, qe0_surf, &
1846 qe, qe_surf, qe_roof, qe_wall, &
1847 ev_per_tstep, runoff_per_tstep, &
1848 surf_chang_per_tstep, runoffpipes, &
1850 runoffagveg, runoffagimpervious, rss_surf)
1853 IF (diagnose == 1) print *,
'before SUEWS_cal_SoilState soilstore_id = ', soilstore_surf_next
1857 IF (diagnose == 1)
WRITE (*, *)
'Calling SUEWS_cal_HorizontalSoilWater...'
1860 soilstorecap_surf, &
1862 sathydraulicconduct, &
1866 soilstore_surf_next, &
1868 runoffsoil_per_tstep &
1872 IF (diagnose == 1)
WRITE (*, *)
'Calling SUEWS_cal_SoilState...'
1874 smdmethod, xsmd, nonwaterfraction, soilmoistcap, &
1875 soilstorecap_surf, surf_chang_per_tstep, &
1876 soilstore_surf_next, soilstore_surf_prev, sfr_surf, &
1877 smd, smd_nsurf, tot_chang_per_tstep, soilstate)
1880 IF (diagnose == 1)
WRITE (*, *)
'Calling SUEWS_cal_QH...'
1882 1, nlayer, storageheatmethod, &
1883 qn, qf, qmrain, qe, qs, qmfreez, qm,
avdens,
avcp, &
1884 sfr_surf, sfr_roof, sfr_wall, &
1885 tsfc_out_surf, tsfc_out_roof, tsfc_out_wall, &
1888 qh, qh_residual, qh_resist, &
1889 qh_resist_surf, qh_resist_roof, qh_resist_wall)
1918 tsfc0_out_surf = min(tsfc_out_surf, temp_c + 50)
1919 tsfc0_out_roof = min(tsfc_out_roof, temp_c + 50)
1920 tsfc0_out_wall = min(tsfc_out_wall, temp_c + 50)
1922 qh_surf = qn_surf + qf - qs_surf - qe_surf
1923 qh_roof = qn_roof + qf - qs_roof - qe_roof
1924 qh_wall = qn_wall + qf - qs_wall - qe_wall
1926 IF (diagnose == 1)
THEN
1927 print *,
'qn_surf before QH back env.:', qn_surf
1928 print *,
'qf before QH back env.:', qf
1929 print *,
'qs_surf before QH back env.:', qs_surf
1930 print *,
'qe_surf before QH back env.:', qe_surf
1931 print *,
'qh_surf before QH back env.:', qh_surf
1933 print *,
'qn_roof before QH back env.:', qn_roof
1934 print *,
'qs_roof before QH back env.:', qs_roof
1935 print *,
'qe_roof before QH back env.:', qe_roof
1936 print *,
'qh_roof before QH back env.:', qh_roof
1939 DO i_surf = 1,
nsurf
1951 DO i_surf = 1, nlayer
1956 IF (diagnose == 1) print *,
'tsfc_surf after QH back env.:', tsfc_out_surf
1958 IF (diagnose == 1) print *, &
1959 'tsfc_surf abs. diff.:', maxval(abs(tsfc_out_surf - tsfc0_out_surf)), maxloc(abs(tsfc_out_surf - tsfc0_out_surf))
1960 dif_tsfc_iter = maxval(abs(tsfc_out_surf - tsfc0_out_surf))
1961 IF (storageheatmethod == 5)
THEN
1962 IF (diagnose == 1) print *, &
1963 'tsfc_roof abs. diff.:', maxval(abs(tsfc_out_roof - tsfc0_out_roof)), maxloc(abs(tsfc_out_roof - tsfc0_out_roof))
1964 dif_tsfc_iter = max(maxval(abs(tsfc_out_roof - tsfc0_out_roof)), dif_tsfc_iter)
1965 IF (diagnose == 1) print *, &
1966 'tsfc_wall abs. diff.:', maxval(abs(tsfc_out_wall - tsfc0_out_wall)), maxloc(abs(tsfc_out_wall - tsfc0_out_wall))
1967 dif_tsfc_iter = max(maxval(abs(tsfc0_out_wall - tsfc_out_wall)), dif_tsfc_iter)
1974 tsfc_out_surf = (tsfc0_out_surf*(1 - ratio_iter) + tsfc_out_surf*ratio_iter)
1975 tsfc_out_roof = (tsfc0_out_roof*(1 - ratio_iter) + tsfc_out_roof*ratio_iter)
1976 tsfc_out_wall = (tsfc0_out_wall*(1 - ratio_iter) + tsfc_out_wall*ratio_iter)
1997 IF (dif_tsfc_iter > .1)
THEN
1998 flag_converge = .false.
2000 flag_converge = .true.
2010 IF (diagnose == 1 .AND. i_iter == max_iter)
THEN
2030 IF (diagnose == 1)
WRITE (*, *)
'Calling RSLProfile...'
2033 zh, z0m, zdm, z0v, &
2034 l_mod, sfr_surf, fai, pai, &
2035 stabilitymethod, ra_h, &
2037 avu1, temp_c, avrh, press_hpa, z, qh, qe, &
2038 t2_c, q2_gkg, u10_ms, rh2, &
2043 alpha_bioco2, alpha_enh_bioco2, kdown, avrh, beta_bioco2, beta_enh_bioco2, &
2044 dectime, diagnose, emissionsmethod, fc_anthro, g_max, g_k, g_q_base, g_q_shape, &
2045 g_t, g_sm, gfunc, gsmodel, id, it, kmax, lai_id_next, laimin, &
2046 laimax, maxconductance, min_res_bioco2, press_hpa, resp_a, &
2047 resp_b, s1, s2, sfr_surf, smdmethod, snowfrac, t2_c, temp_c, theta_bioco2, th, tl, vsmd, xsmd, &
2048 fc, fc_biogen, fc_photo, fc_respi)
2057 qn_s_av = qn_s_av_next
2058 dqnsdt = dqnsdt_next
2059 snowfallcum = snowfallcum_next
2060 snowalb = snowalb_next
2061 icefrac = icefrac_next
2062 snowwater = snowwater_next
2063 snowdens = snowdens_next
2064 snowfrac = snowfrac_next
2065 snowpack = snowpack_next
2067 soilstore_surf = soilstore_surf_next
2068 state_surf = state_surf_next
2070 gdd_id = gdd_id_next
2071 sdd_id = sdd_id_next
2072 lai_id = lai_id_next
2073 decidcap_id = decidcap_id_next
2074 albdectr_id = albdectr_id_next
2075 albevetr_id = albevetr_id_next
2076 albgrass_id = albgrass_id_next
2077 porosity_id = porosity_id_next
2078 storedrainprm = storedrainprm_next
2079 tair_av = tair_av_next
2080 tmin_id = tmin_id_next
2081 tmax_id = tmax_id_next
2082 lenday_id = lenday_id_next
2083 hdd_id = hdd_id_next
2084 wuday_id = wuday_id_next
2086 IF (storageheatmethod == 5)
THEN
2088 temp_roof = temp_out_roof
2089 temp_wall = temp_out_wall
2090 temp_surf = temp_out_surf
2091 tsfc_roof = tsfc_out_roof
2092 tsfc_wall = tsfc_out_wall
2093 tsfc_surf = tsfc_out_surf
2095 soilstore_roof = soilstore_roof_next
2096 state_roof = state_roof_next
2097 soilstore_wall = soilstore_wall_next
2098 state_wall = state_wall_next
2112 pai = sfr_surf(2)/sum(sfr_surf(1:2))
2113 CALL beers_cal_main(iy, id, dectime, pai, fai, kdown, ldown, temp_c, avrh, &
2114 press_hpa, tsfc_c, lat, lng, alt, timezone, zenith_deg, azimuth, &
2115 alb(1), alb(2), emis(1), emis(2), &
2120 dataoutlinebeers =
set_nan(dataoutlinebeers)
2125 additionalwater, alb, kdown, u10_ms, azimuth, &
2126 chsnow_per_interval, dectime, &
2127 drain_per_tstep, qe_lumps, ev_per_tstep, wu_ext, fc, fc_build, fcld, &
2128 fc_metab, fc_photo, fc_respi, fc_point, fc_traff, flowchange, &
2129 qh_lumps, id, imin, wu_int, it, iy, &
2130 kup, lai_id, ldown, l_mod, lup, mwh, &
2132 nsh_real, nwstate_per_tstep, precip, q2_gkg, &
2133 qe, qf, qh, qh_resist, qm, qmfreez, &
2134 qmrain, qn, qn_snow, qn_snowfree, qs, ra_h, &
2135 rs, rh2, runoffagimpervious, runoffagveg, &
2136 runoff_per_tstep, runoffpipes, runoffsoil_per_tstep, &
2137 runoffwaterbody, sfr_surf, smd, smd_nsurf, snowalb, snowremoval, &
2138 state_surf_next, state_per_tstep, surf_chang_per_tstep, swe, t2_c, tsfc_c, &
2139 tot_chang_per_tstep, tsurf, ustar, &
2141 z0m, zdm, zenith_deg, &
2142 datetimeline, dataoutlinesuews)
2145 iy, id, it, imin, dectime, nlayer, &
2146 tsfc_out_surf, qs_surf, &
2161 datetimeline, dataoutlineehc)
2165 it, imin, nsh_real, &
2166 gdd_id, hdd_id, lai_id, &
2168 tmin_id, tmax_id, lenday_id, &
2176 snowalb, snowdens, &
2178 dataoutlinedailystate)
2182 dataoutlinedebug = &
2184 qn_surf, qs_surf, qe0_surf, qe_surf, qh_surf, &
2185 wu_surf, ev0_surf, ev_surf, drain_surf, state_surf_prev, state_surf_next, soilstore_surf_prev, soilstore_surf_next, &
2186 rs, ra_h, rb, rasnow, rss_surf, &
2187 vsmd, s1/g_sm + s2, g_sm, g_sm*(vsmd - s1/g_sm + s2), &
2188 g_kdown, g_dq, g_ta, g_smd, g_lai, &
2191 faibldg_use, faievetree_use, faidectree_use, fai, &
2196 output_line_suews%datetimeLine = datetimeline
2197 output_line_suews%dataOutLineSUEWS = [datetimeline, dataoutlinesuews]
2198 output_line_suews%dataOutLineEHC = [datetimeline, dataoutlineehc]
2199 output_line_suews%dataOutLineDailyState = [datetimeline, dataoutlinedailystate]
2200 output_line_suews%dataOutLineBEERS = [datetimeline, dataoutlinebeers]
2201 output_line_suews%dataOutLineDebug = [datetimeline, dataoutlinedebug]
2202 output_line_suews%dataOutLineSPARTACUS = [datetimeline, dataoutlinespartacus]
2203 output_line_suews%dataOutLineSnow = [datetimeline, dataoutlinesnow]
2204 output_line_suews%dataoutLineRSL = [datetimeline, dataoutlinersl]
2205 output_line_suews%dataOutLineESTM = [datetimeline, dataoutlineestm]
2210 AH_MIN, AHProf_24hr, AH_SLOPE_Cooling, AH_SLOPE_Heating, & ! input&inout in alphabetical order
2211 alb, AlbMax_DecTr, AlbMax_EveTr, AlbMax_Grass, &
2212 AlbMin_DecTr, AlbMin_EveTr, AlbMin_Grass, &
2213 alpha_bioCO2, alpha_enh_bioCO2, alt, kdown, avRh, avU1, BaseT, BaseTe, &
2214 beta_bioCO2, beta_enh_bioCO2, bldgH, CapMax_dec, CapMin_dec, &
2215 chAnOHM, CO2PointSource, cpAnOHM, CRWmax, CRWmin, DayWat, DayWatPer, &
2216 DecTreeH, DiagMethod, Diagnose, DRAINRT, &
2217 dt_since_start, dqndt, qn_av, dqnsdt, qn_s_av, &
2218 EF_umolCO2perJ, emis, EmissionsMethod, EnEF_v_Jkm, endDLS, EveTreeH, FAIBldg, &
2219 FAIDecTree, FAIEveTree, FAIMethod, Faut, FcEF_v_kgkm, fcld_obs, FlowChange, &
2220 FrFossilFuel_Heat, FrFossilFuel_NonHeat, g_max, g_k, g_q_base, g_q_shape, g_t, g_sm, GDD_id, &
2221 GDDFull, Gridiv, gsModel, H_maintain, HDD_id, HumActivity_24hr, &
2222 IceFrac, id, Ie_a, Ie_end, Ie_m, Ie_start, imin, &
2223 InternalWaterUse_h, &
2224 IrrFracPaved, IrrFracBldgs, &
2225 IrrFracEveTr, IrrFracDecTr, IrrFracGrass, &
2226 IrrFracBSoil, IrrFracWater, &
2228 iy, kkAnOHM, Kmax, LAI_id, LAIMax, LAIMin, LAI_obs, &
2229 LAIPower, LAIType, lat, lenDay_id, ldown_obs, lng, MaxConductance, MaxFCMetab, MaxQFMetab, &
2230 SnowWater, MinFCMetab, MinQFMetab, min_res_bioCO2, &
2231 NARP_EMIS_SNOW, NARP_TRANS_SITE, NetRadiationMethod, &
2233 n_vegetation_region_urban, &
2234 n_stream_sw_urban, n_stream_lw_urban, &
2235 sw_dn_direct_frac, air_ext_sw, air_ssa_sw, &
2236 veg_ssa_sw, air_ext_lw, air_ssa_lw, veg_ssa_lw, &
2237 veg_fsd_const, veg_contact_fraction_const, &
2238 ground_albedo_dir_mult_fact, use_sw_direct_albedo, & !input
2239 height, building_frac, veg_frac, building_scale, veg_scale, & !input: SPARTACUS
2240 alb_roof, emis_roof, alb_wall, emis_wall, &
2241 roof_albedo_dir_mult_fact, wall_specular_frac, &
2242 OHM_coef, OHMIncQF, OHM_threshSW, &
2243 OHM_threshWD, PipeCapacity, PopDensDaytime, &
2244 PopDensNighttime, PopProf_24hr, PorMax_dec, PorMin_dec, &
2245 Precip, PrecipLimit, PrecipLimitAlb, Press_hPa, &
2246 QF0_BEU, Qf_A, Qf_B, Qf_C, &
2247 qn1_obs, qs_obs, qf_obs, &
2248 RadMeltFact, RAINCOVER, RainMaxRes, resp_a, resp_b, &
2249 RoughLenHeatMethod, RoughLenMomMethod, RunoffToWater, S1, S2, &
2250 SatHydraulicConduct, SDDFull, SDD_id, SMDMethod, SnowAlb, SnowAlbMax, &
2251 SnowAlbMin, SnowPackLimit, SnowDens, SnowDensMax, SnowDensMin, SnowfallCum, SnowFrac, &
2252 SnowLimBldg, SnowLimPaved, snowFrac_obs, SnowPack, SnowProf_24hr, SnowUse, SoilDepth, &
2253 StabilityMethod, startDLS, &
2254 soilstore_surf, SoilStoreCap_surf, state_surf, StateLimit_surf, WetThresh_surf, &
2255 soilstore_roof, SoilStoreCap_roof, state_roof, StateLimit_roof, WetThresh_roof, &
2256 soilstore_wall, SoilStoreCap_wall, state_wall, StateLimit_wall, WetThresh_wall, &
2257 StorageHeatMethod, StoreDrainPrm, SurfaceArea, Tair_av, tau_a, tau_f, tau_r, &
2259 BaseT_Cooling, BaseT_Heating, Temp_C, TempMeltFact, TH, &
2260 theta_bioCO2, timezone, TL, TrafficRate, TrafficUnits, &
2262 tsfc_roof, tsfc_wall, tsfc_surf, &
2263 temp_roof, temp_wall, temp_surf, &
2264 tin_roof, tin_wall, tin_surf, &
2265 k_roof, k_wall, k_surf, &
2266 cp_roof, cp_wall, cp_surf, &
2267 dz_roof, dz_wall, dz_surf, &
2268 TraffProf_24hr, Ts5mindata_ir, tstep, tstep_prev, veg_type, &
2269 WaterDist, WaterUseMethod, wu_m3, &
2270 WUDay_id, DecidCap_id, albDecTr_id, albEveTr_id, albGrass_id, porosity_id, &
2271 WUProfA_24hr, WUProfM_24hr, xsmd, Z, z0m_in, zdm_in, &
2288 REAL(KIND(1D0)),
DIMENSION(:),
INTENT(IN) :: Ts5mindata_ir
2289 REAL(KIND(1D0)),
DIMENSION(10) :: MetForcingData_grid
2291 INTEGER,
INTENT(IN) :: nlayer
2295 REAL(KIND(1D0)),
INTENT(IN) :: lat
2296 REAL(KIND(1D0)),
INTENT(IN) :: lng
2297 REAL(KIND(1D0)),
INTENT(IN) :: alt
2298 INTEGER,
INTENT(IN) :: Gridiv
2299 REAL(KIND(1D0)),
INTENT(IN) :: timezone
2300 REAL(KIND(1D0)),
INTENT(IN) :: SurfaceArea
2301 REAL(KIND(1D0)),
INTENT(IN) :: Z
2302 REAL(KIND(1D0)),
INTENT(IN) :: z0m_in
2303 REAL(KIND(1D0)),
INTENT(IN) :: zdm_in
2304 REAL(KIND(1D0)),
INTENT(IN) :: PipeCapacity
2305 REAL(KIND(1D0)),
INTENT(IN) :: RunoffToWater
2306 REAL(KIND(1D0)),
INTENT(IN) :: NARP_TRANS_SITE
2307 REAL(KIND(1D0)),
INTENT(IN) :: CO2PointSource
2308 REAL(KIND(1D0)),
INTENT(IN) :: FlowChange
2312 REAL(KIND(1D0)),
INTENT(IN) :: kdown
2313 REAL(KIND(1D0)),
INTENT(IN) :: ldown_obs
2314 REAL(KIND(1D0)),
INTENT(IN) :: avRh
2315 REAL(KIND(1D0)),
INTENT(IN) :: Press_hPa
2316 REAL(KIND(1D0)),
INTENT(IN) :: avU1
2317 REAL(KIND(1D0)),
INTENT(IN) :: Precip
2318 REAL(KIND(1D0)),
INTENT(IN) :: fcld_obs
2319 REAL(KIND(1D0)),
INTENT(IN) :: LAI_obs
2320 REAL(KIND(1D0)),
INTENT(IN) :: snowFrac_obs
2321 REAL(KIND(1D0)),
INTENT(IN) :: xsmd
2322 REAL(KIND(1D0)),
INTENT(IN) :: qn1_obs
2323 REAL(KIND(1D0)),
INTENT(IN) :: qs_obs
2324 REAL(KIND(1D0)),
INTENT(IN) :: qf_obs
2325 REAL(KIND(1D0)),
INTENT(IN) :: Temp_C
2326 REAL(KIND(1D0)),
INTENT(IN) :: wu_m3
2328 REAL(KIND(1D0)),
INTENT(INOUT) :: Tair_av
2332 INTEGER,
INTENT(IN) :: id
2333 INTEGER,
INTENT(IN) :: imin
2334 INTEGER,
INTENT(IN) :: isec
2335 INTEGER,
INTENT(IN) :: it
2336 INTEGER,
INTENT(IN) :: iy
2337 INTEGER,
INTENT(IN) :: tstep
2338 INTEGER,
INTENT(IN) :: tstep_prev
2339 INTEGER,
INTENT(in) :: dt_since_start
2343 INTEGER,
INTENT(IN) :: Diagnose
2344 INTEGER,
INTENT(in) :: DiagMethod
2345 INTEGER,
INTENT(IN) :: EmissionsMethod
2346 INTEGER,
INTENT(IN) :: RoughLenHeatMethod
2347 INTEGER,
INTENT(IN) :: RoughLenMomMethod
2348 INTEGER,
INTENT(IN) :: FAIMethod
2349 INTEGER,
INTENT(IN) :: SMDMethod
2350 INTEGER,
INTENT(IN) :: WaterUseMethod
2351 INTEGER,
INTENT(IN) :: NetRadiationMethod
2352 INTEGER,
INTENT(IN) :: StabilityMethod
2353 INTEGER,
INTENT(IN) :: StorageHeatMethod
2354 INTEGER,
INTENT(IN) :: SnowUse
2355 LOGICAL,
INTENT(IN) :: use_sw_direct_albedo
2356 INTEGER,
INTENT(IN) :: OHMIncQF
2360 REAL(KIND(1D0)),
INTENT(IN) :: RAINCOVER
2361 REAL(KIND(1D0)),
INTENT(IN) :: RainMaxRes
2362 REAL(KIND(1D0)),
INTENT(IN) :: DRAINRT
2363 INTEGER,
INTENT(IN) :: veg_type
2367 REAL(KIND(1D0)),
DIMENSION(nlayer),
INTENT(IN) :: SoilStoreCap_roof
2368 REAL(KIND(1D0)),
DIMENSION(nlayer),
INTENT(IN) :: StateLimit_roof
2369 REAL(KIND(1D0)),
DIMENSION(nlayer),
INTENT(IN) :: wetthresh_roof
2370 REAL(KIND(1D0)),
DIMENSION(nlayer),
INTENT(IN) :: SoilStoreCap_wall
2371 REAL(KIND(1D0)),
DIMENSION(nlayer),
INTENT(IN) :: StateLimit_wall
2372 REAL(KIND(1D0)),
DIMENSION(nlayer),
INTENT(IN) :: wetthresh_wall
2373 REAL(KIND(1D0)),
DIMENSION(nlayer),
INTENT(in) :: tin_roof
2374 REAL(KIND(1D0)),
DIMENSION(nlayer, ndepth),
INTENT(in) :: k_roof
2375 REAL(KIND(1D0)),
DIMENSION(nlayer, ndepth),
INTENT(in) :: cp_roof
2376 REAL(KIND(1D0)),
DIMENSION(nlayer, ndepth),
INTENT(in) :: dz_roof
2377 REAL(KIND(1D0)),
DIMENSION(nlayer),
INTENT(in) :: tin_wall
2378 REAL(KIND(1D0)),
DIMENSION(nlayer, ndepth),
INTENT(in) :: k_wall
2379 REAL(KIND(1D0)),
DIMENSION(nlayer, ndepth),
INTENT(in) :: cp_wall
2380 REAL(KIND(1D0)),
DIMENSION(nlayer, ndepth),
INTENT(in) :: dz_wall
2381 REAL(KIND(1D0)),
DIMENSION(nsurf),
INTENT(in) :: tin_surf
2382 REAL(KIND(1D0)),
DIMENSION(nsurf, ndepth),
INTENT(in) :: k_surf
2383 REAL(KIND(1D0)),
DIMENSION(nsurf, ndepth),
INTENT(in) :: cp_surf
2384 REAL(KIND(1D0)),
DIMENSION(nsurf, ndepth),
INTENT(in) :: dz_surf
2388 REAL(KIND(1D0)),
INTENT(IN) :: air_ext_lw
2389 REAL(KIND(1D0)),
INTENT(IN) :: air_ext_sw
2390 REAL(KIND(1D0)),
INTENT(IN) :: air_ssa_lw
2391 REAL(KIND(1D0)),
INTENT(IN) :: air_ssa_sw
2392 REAL(KIND(1D0)),
INTENT(IN) :: veg_ssa_lw
2393 REAL(KIND(1D0)),
INTENT(IN) :: veg_ssa_sw
2394 REAL(KIND(1D0)),
DIMENSION(nlayer + 1),
INTENT(IN) :: height
2395 REAL(KIND(1D0)),
INTENT(IN) :: ground_albedo_dir_mult_fact
2396 INTEGER,
INTENT(IN) :: n_stream_lw_urban
2397 INTEGER,
INTENT(IN) :: n_stream_sw_urban
2398 INTEGER,
INTENT(IN) :: n_vegetation_region_urban
2399 REAL(KIND(1D0)),
INTENT(IN) :: sw_dn_direct_frac
2400 REAL(KIND(1D0)),
INTENT(IN) :: veg_contact_fraction_const
2401 REAL(KIND(1D0)),
INTENT(IN) :: veg_fsd_const
2405 REAL(KIND(1D0)),
DIMENSION(nlayer),
INTENT(IN) :: building_frac
2406 REAL(KIND(1D0)),
DIMENSION(nlayer),
INTENT(IN) :: building_scale
2407 REAL(KIND(1D0)),
DIMENSION(nlayer),
INTENT(IN) :: veg_frac
2408 REAL(KIND(1D0)),
DIMENSION(nlayer),
INTENT(IN) :: veg_scale
2409 REAL(KIND(1D0)),
DIMENSION(nlayer),
INTENT(IN) :: alb_roof
2410 REAL(KIND(1D0)),
DIMENSION(nlayer),
INTENT(IN) :: emis_roof
2411 REAL(KIND(1D0)),
DIMENSION(nlayer),
INTENT(IN) :: alb_wall
2412 REAL(KIND(1D0)),
DIMENSION(nlayer),
INTENT(IN) :: emis_wall
2413 REAL(KIND(1D0)),
DIMENSION(nspec, nlayer),
INTENT(IN) :: roof_albedo_dir_mult_fact
2414 REAL(KIND(1D0)),
DIMENSION(nspec, nlayer),
INTENT(IN) :: wall_specular_frac
2418 INTEGER,
INTENT(IN) :: startDLS
2419 INTEGER,
INTENT(IN) :: endDLS
2420 REAL(KIND(1D0)),
DIMENSION(2),
INTENT(IN) :: QF0_BEU
2421 REAL(KIND(1D0)),
DIMENSION(2),
INTENT(IN) :: Qf_A
2422 REAL(KIND(1D0)),
DIMENSION(2),
INTENT(IN) :: Qf_B
2423 REAL(KIND(1D0)),
DIMENSION(2),
INTENT(IN) :: Qf_C
2424 REAL(KIND(1D0)),
DIMENSION(2),
INTENT(IN) :: BaseT_Cooling
2425 REAL(KIND(1D0)),
DIMENSION(2),
INTENT(IN) :: BaseT_Heating
2426 REAL(KIND(1D0)),
DIMENSION(2),
INTENT(IN) :: PopDensDaytime
2427 REAL(KIND(1D0)),
INTENT(IN) :: PopDensNighttime
2428 REAL(KIND(1D0)),
DIMENSION(0:23, 2),
INTENT(IN) :: PopProf_24hr
2429 REAL(KIND(1D0)),
DIMENSION(2),
INTENT(IN) :: AH_MIN
2430 REAL(KIND(1D0)),
DIMENSION(0:23, 2),
INTENT(IN) :: AHProf_24hr
2431 REAL(KIND(1D0)),
DIMENSION(2),
INTENT(IN) :: AH_SLOPE_Cooling
2432 REAL(KIND(1D0)),
DIMENSION(2),
INTENT(IN) :: AH_SLOPE_Heating
2433 REAL(KIND(1D0)),
INTENT(IN) :: EF_umolCO2perJ
2434 REAL(KIND(1D0)),
INTENT(IN) :: EnEF_v_Jkm
2435 REAL(KIND(1D0)),
INTENT(IN) :: FrFossilFuel_Heat
2436 REAL(KIND(1D0)),
INTENT(IN) :: FrFossilFuel_NonHeat
2437 REAL(KIND(1D0)),
DIMENSION(2),
INTENT(IN) :: FcEF_v_kgkm
2438 REAL(KIND(1D0)),
DIMENSION(0:23, 2),
INTENT(IN) :: HumActivity_24hr
2439 REAL(KIND(1D0)),
INTENT(IN) :: MaxFCMetab
2440 REAL(KIND(1D0)),
INTENT(IN) :: MaxQFMetab
2441 REAL(KIND(1D0)),
INTENT(IN) :: MinFCMetab
2442 REAL(KIND(1D0)),
INTENT(IN) :: MinQFMetab
2443 REAL(KIND(1D0)),
DIMENSION(2),
INTENT(IN) :: TrafficRate
2444 REAL(KIND(1D0)),
INTENT(IN) :: TrafficUnits
2445 REAL(KIND(1D0)),
DIMENSION(0:23, 2),
INTENT(IN) :: TraffProf_24hr
2449 REAL(KIND(1D0)),
INTENT(IN) :: H_maintain
2450 REAL(KIND(1D0)),
INTENT(IN) :: Faut
2451 REAL(KIND(1D0)),
DIMENSION(3),
INTENT(IN) :: Ie_a
2452 REAL(KIND(1D0)),
DIMENSION(3),
INTENT(IN) :: Ie_m
2453 INTEGER,
INTENT(IN) :: Ie_end
2454 INTEGER,
INTENT(IN) :: Ie_start
2455 REAL(KIND(1D0)),
INTENT(IN) :: InternalWaterUse_h
2456 REAL(KIND(1D0)),
DIMENSION(7),
INTENT(IN) :: DayWat
2457 REAL(KIND(1D0)),
DIMENSION(7),
INTENT(IN) :: DayWatPer
2458 REAL(KIND(1D0)),
DIMENSION(0:23, 2),
INTENT(IN) :: WUProfA_24hr
2459 REAL(KIND(1D0)),
DIMENSION(0:23, 2),
INTENT(IN) :: WUProfM_24hr
2463 REAL(KIND(1D0)),
INTENT(IN) :: CRWmax
2464 REAL(KIND(1D0)),
INTENT(IN) :: CRWmin
2465 REAL(KIND(1D0)),
INTENT(IN) :: NARP_EMIS_SNOW
2466 REAL(KIND(1D0)),
INTENT(IN) :: PrecipLimit
2467 REAL(KIND(1D0)),
INTENT(IN) :: PrecipLimitAlb
2468 REAL(KIND(1D0)),
INTENT(IN) :: SnowAlbMax
2469 REAL(KIND(1D0)),
INTENT(IN) :: SnowAlbMin
2470 REAL(KIND(1D0)),
INTENT(IN) :: SnowDensMax
2471 REAL(KIND(1D0)),
INTENT(IN) :: SnowDensMin
2472 REAL(KIND(1D0)),
INTENT(IN) :: SnowLimBldg
2473 REAL(KIND(1D0)),
INTENT(IN) :: SnowLimPaved
2474 REAL(KIND(1D0)),
DIMENSION(NSURF),
INTENT(IN) :: SnowPackLimit
2475 REAL(KIND(1D0)),
DIMENSION(0:23, 2),
INTENT(IN) :: SnowProf_24hr
2476 REAL(KIND(1D0)),
INTENT(IN) :: tau_a
2477 REAL(KIND(1D0)),
INTENT(IN) :: tau_f
2478 REAL(KIND(1D0)),
INTENT(IN) :: tau_r
2479 REAL(KIND(1D0)),
INTENT(IN) :: TempMeltFact
2480 REAL(KIND(1D0)),
INTENT(IN) :: RadMeltFact
2484 REAL(KIND(1D0)),
INTENT(IN) :: g_max
2485 REAL(KIND(1D0)),
INTENT(IN) :: g_k
2486 REAL(KIND(1D0)),
INTENT(IN) :: g_q_base
2487 REAL(KIND(1D0)),
INTENT(IN) :: g_q_shape
2488 REAL(KIND(1D0)),
INTENT(IN) :: g_t
2489 REAL(KIND(1D0)),
INTENT(IN) :: g_sm
2490 REAL(KIND(1D0)),
INTENT(IN) :: Kmax
2491 INTEGER,
INTENT(IN) :: gsModel
2492 REAL(KIND(1D0)),
INTENT(IN) :: S1
2493 REAL(KIND(1D0)),
INTENT(IN) :: S2
2494 REAL(KIND(1D0)),
INTENT(IN) :: TH
2495 REAL(KIND(1D0)),
INTENT(IN) :: TL
2506 REAL(KIND(1D0)),
DIMENSION(NSURF),
INTENT(IN) :: sfr_surf
2507 REAL(KIND(1D0)),
DIMENSION(NSURF),
INTENT(IN) :: emis
2508 REAL(KIND(1D0)),
DIMENSION(NSURF),
INTENT(IN) :: chAnOHM
2509 REAL(KIND(1D0)),
DIMENSION(NSURF),
INTENT(IN) :: cpAnOHM
2510 REAL(KIND(1D0)),
DIMENSION(NSURF),
INTENT(IN) :: kkAnOHM
2511 REAL(KIND(1D0)),
DIMENSION(nsurf + 1),
INTENT(IN) :: OHM_threshSW
2512 REAL(KIND(1D0)),
DIMENSION(nsurf + 1),
INTENT(IN) :: OHM_threshWD
2513 REAL(KIND(1D0)),
DIMENSION(nsurf + 1, 4, 3),
INTENT(IN) :: OHM_coef
2514 REAL(KIND(1D0)),
DIMENSION(NSURF),
INTENT(IN) :: SoilDepth
2515 REAL(KIND(1D0)),
DIMENSION(NSURF),
INTENT(IN) :: SoilStoreCap_surf
2516 REAL(KIND(1D0)),
DIMENSION(NSURF),
INTENT(IN) :: SatHydraulicConduct
2517 REAL(KIND(1D0)),
DIMENSION(NSURF),
INTENT(IN) :: StateLimit_surf
2518 REAL(KIND(1D0)),
DIMENSION(NSURF),
INTENT(IN) :: WetThresh_surf
2519 REAL(KIND(1D0)),
DIMENSION(NSURF + 1, NSURF - 1),
INTENT(IN) :: WaterDist
2521 REAL(KIND(1D0)),
INTENT(IN) :: IrrFracPaved
2522 REAL(KIND(1D0)),
INTENT(IN) :: IrrFracBldgs
2523 REAL(KIND(1D0)),
INTENT(IN) :: IrrFracDecTr
2524 REAL(KIND(1D0)),
INTENT(IN) :: IrrFracEveTr
2525 REAL(KIND(1D0)),
INTENT(IN) :: IrrFracGrass
2526 REAL(KIND(1D0)),
INTENT(IN) :: IrrFracBSoil
2527 REAL(KIND(1D0)),
INTENT(IN) :: IrrFracWater
2529 REAL(KIND(1D0)),
INTENT(IN) :: bldgH
2530 REAL(KIND(1D0)),
INTENT(IN) :: FAIBldg
2532 REAL(KIND(1D0)),
INTENT(IN) :: DecTreeH
2533 REAL(KIND(1D0)),
INTENT(IN) :: FAIDecTree
2534 REAL(KIND(1D0)),
INTENT(IN) :: CapMax_dec
2535 REAL(KIND(1D0)),
INTENT(IN) :: CapMin_dec
2536 REAL(KIND(1D0)),
INTENT(IN) :: PorMax_dec
2537 REAL(KIND(1D0)),
INTENT(IN) :: PorMin_dec
2538 REAL(KIND(1D0)),
INTENT(IN) :: AlbMax_DecTr
2539 REAL(KIND(1D0)),
INTENT(IN) :: AlbMin_DecTr
2541 REAL(KIND(1D0)),
INTENT(IN) :: EveTreeH
2542 REAL(KIND(1D0)),
INTENT(IN) :: FAIEveTree
2543 REAL(KIND(1D0)),
INTENT(IN) :: AlbMax_EveTr
2544 REAL(KIND(1D0)),
INTENT(IN) :: AlbMin_EveTr
2546 REAL(KIND(1D0)),
INTENT(IN) :: AlbMax_Grass
2547 REAL(KIND(1D0)),
INTENT(IN) :: AlbMin_Grass
2549 REAL(KIND(1D0)),
DIMENSION(NVEGSURF),
INTENT(IN) :: beta_bioCO2
2550 REAL(KIND(1D0)),
DIMENSION(NVEGSURF),
INTENT(IN) :: beta_enh_bioCO2
2551 REAL(KIND(1D0)),
DIMENSION(NVEGSURF),
INTENT(IN) :: alpha_bioCO2
2552 REAL(KIND(1D0)),
DIMENSION(NVEGSURF),
INTENT(IN) :: alpha_enh_bioCO2
2553 REAL(KIND(1D0)),
DIMENSION(NVEGSURF),
INTENT(IN) :: resp_a
2554 REAL(KIND(1D0)),
DIMENSION(NVEGSURF),
INTENT(IN) :: resp_b
2555 REAL(KIND(1D0)),
DIMENSION(NVEGSURF),
INTENT(IN) :: min_res_bioCO2
2556 REAL(KIND(1D0)),
DIMENSION(NVEGSURF),
INTENT(IN) :: theta_bioCO2
2558 REAL(KIND(1D0)),
DIMENSION(NVEGSURF),
INTENT(IN) :: BaseT
2559 REAL(KIND(1D0)),
DIMENSION(NVEGSURF),
INTENT(IN) :: BaseTe
2560 REAL(KIND(1D0)),
DIMENSION(NVEGSURF),
INTENT(IN) :: SDDFull
2561 REAL(KIND(1D0)),
DIMENSION(NVEGSURF),
INTENT(IN) :: GDDFull
2562 REAL(KIND(1D0)),
DIMENSION(NVEGSURF),
INTENT(IN) :: LAIMax
2563 REAL(KIND(1D0)),
DIMENSION(NVEGSURF),
INTENT(IN) :: LAIMin
2564 REAL(KIND(1D0)),
DIMENSION(4, NVEGSURF),
INTENT(IN) :: LAIPower
2565 INTEGER,
DIMENSION(NVEGSURF),
INTENT(IN) :: LAIType
2567 REAL(KIND(1D0)),
DIMENSION(3),
INTENT(IN) :: MaxConductance
2572 REAL(KIND(1D0)),
DIMENSION(12),
INTENT(INOUT) :: HDD_id
2576 REAL(KIND(1D0)),
DIMENSION(nlayer),
INTENT(INOUT) :: soilstore_roof
2577 REAL(KIND(1D0)),
DIMENSION(nlayer),
INTENT(INOUT) :: state_roof
2578 REAL(KIND(1D0)),
DIMENSION(nlayer),
INTENT(INOUT) :: soilstore_wall
2579 REAL(KIND(1D0)),
DIMENSION(nlayer),
INTENT(INOUT) :: state_wall
2580 REAL(KIND(1D0)),
DIMENSION(NSURF),
INTENT(INOUT) :: soilstore_surf
2581 REAL(KIND(1D0)),
DIMENSION(NSURF),
INTENT(INOUT) :: state_surf
2582 REAL(KIND(1D0)),
DIMENSION(9),
INTENT(INOUT) :: WUDay_id
2586 REAL(KIND(1D0)),
DIMENSION(nlayer, ndepth),
INTENT(INOUT) :: temp_roof
2587 REAL(KIND(1D0)),
DIMENSION(nlayer),
INTENT(INOUT) :: tsfc_roof
2588 REAL(KIND(1D0)),
DIMENSION(nlayer, ndepth),
INTENT(INOUT) :: temp_wall
2589 REAL(KIND(1D0)),
DIMENSION(nlayer),
INTENT(INOUT) :: tsfc_wall
2590 REAL(KIND(1D0)),
DIMENSION(nsurf),
INTENT(INOUT) :: tsfc_surf
2591 REAL(KIND(1D0)),
DIMENSION(nsurf, ndepth),
INTENT(INOUT) :: temp_surf
2595 REAL(KIND(1D0)),
INTENT(INOUT) :: qn_av
2596 REAL(KIND(1D0)),
INTENT(INOUT) :: dqndt
2597 REAL(KIND(1D0)),
INTENT(INOUT) :: qn_s_av
2598 REAL(KIND(1D0)),
INTENT(INOUT) :: dqnsdt
2602 REAL(KIND(1D0)),
INTENT(INOUT) :: SnowfallCum
2603 REAL(KIND(1D0)),
INTENT(INOUT) :: SnowAlb
2604 REAL(KIND(1D0)),
DIMENSION(NSURF),
INTENT(INOUT) :: IceFrac
2605 REAL(KIND(1D0)),
DIMENSION(NSURF),
INTENT(INOUT) :: SnowWater
2606 REAL(KIND(1D0)),
DIMENSION(NSURF),
INTENT(INOUT) :: SnowDens
2607 REAL(KIND(1D0)),
DIMENSION(NSURF),
INTENT(INOUT) :: SnowFrac
2608 REAL(KIND(1D0)),
DIMENSION(NSURF),
INTENT(INOUT) :: SnowPack
2612 REAL(KIND(1D0)),
DIMENSION(NSURF),
INTENT(INOUT) :: alb
2613 REAL(KIND(1D0)),
DIMENSION(nvegsurf),
INTENT(INOUT) :: GDD_id
2614 REAL(KIND(1D0)),
DIMENSION(nvegsurf),
INTENT(INout) :: SDD_id
2615 REAL(KIND(1D0)),
DIMENSION(nvegsurf),
INTENT(INOUT) :: LAI_id
2616 REAL(KIND(1D0)),
INTENT(INout) :: Tmin_id
2617 REAL(KIND(1D0)),
INTENT(INout) :: Tmax_id
2618 REAL(KIND(1D0)),
INTENT(INout) :: lenDay_id
2619 REAL(KIND(1D0)),
INTENT(INOUT) :: DecidCap_id
2620 REAL(KIND(1D0)),
INTENT(INOUT) :: albDecTr_id
2621 REAL(KIND(1D0)),
INTENT(INOUT) :: albEveTr_id
2622 REAL(KIND(1D0)),
INTENT(INOUT) :: albGrass_id
2623 REAL(KIND(1D0)),
INTENT(INOUT) :: porosity_id
2624 REAL(KIND(1D0)),
DIMENSION(6, NSURF),
INTENT(INOUT) :: StoreDrainPrm
2628 REAL(KIND(1D0)),
DIMENSION(5) :: datetimeLine
2629 REAL(KIND(1D0)),
DIMENSION(ncolumnsDataOutSUEWS - 5) :: dataOutLineSUEWS
2630 REAL(KIND(1D0)),
DIMENSION(ncolumnsDataOutSnow - 5) :: dataOutLineSnow
2631 REAL(KIND(1D0)),
DIMENSION(ncolumnsDataOutESTM - 5) :: dataOutLineESTM
2632 REAL(KIND(1D0)),
DIMENSION(ncolumnsDataOutEHC - 5) :: dataOutLineEHC
2633 REAL(KIND(1D0)),
DIMENSION(ncolumnsDataOutRSL - 5) :: dataoutLineRSL
2634 REAL(KIND(1D0)),
DIMENSION(ncolumnsDataOutBEERS - 5) :: dataOutLineBEERS
2635 REAL(KIND(1D0)),
DIMENSION(ncolumnsDataOutDebug - 5) :: dataOutLineDebug
2636 REAL(KIND(1D0)),
DIMENSION(ncolumnsDataOutSPARTACUS - 5) :: dataOutLineSPARTACUS
2637 REAL(KIND(1D0)),
DIMENSION(ncolumnsDataOutDailyState - 5) :: dataOutLineDailyState
2639 TYPE(
output_line),
INTENT(OUT) :: output_line_suews
2644 REAL(KIND(1D0)) :: a1
2645 REAL(KIND(1D0)) :: a2
2646 REAL(KIND(1D0)) :: a3
2647 REAL(KIND(1D0)) :: AdditionalWater = 0
2648 REAL(KIND(1D0)) :: U10_ms
2649 REAL(KIND(1D0)) :: azimuth
2650 REAL(KIND(1D0)) :: chSnow_per_interval
2652 REAL(KIND(1D0)) :: dens_dry
2653 REAL(KIND(1D0)) :: deltaLAI
2654 REAL(KIND(1D0)) :: drain_per_tstep
2655 REAL(KIND(1D0)) :: Ea_hPa
2656 REAL(KIND(1D0)) :: QE_LUMPS
2657 REAL(KIND(1D0)) :: es_hPa
2658 REAL(KIND(1D0)) :: ev_per_tstep
2659 REAL(KIND(1D0)) :: wu_ext
2660 REAL(KIND(1D0)) :: Fc
2661 REAL(KIND(1D0)) :: Fc_anthro
2662 REAL(KIND(1D0)) :: Fc_biogen
2663 REAL(KIND(1D0)) :: Fc_build
2664 REAL(KIND(1D0)) :: fcld
2665 REAL(KIND(1D0)) :: Fc_metab
2666 REAL(KIND(1D0)) :: Fc_photo
2667 REAL(KIND(1D0)) :: Fc_point
2668 REAL(KIND(1D0)) :: Fc_respi
2669 REAL(KIND(1D0)) :: Fc_traff
2670 REAL(KIND(1D0)) :: gfunc
2671 REAL(KIND(1D0)) :: gsc
2672 REAL(KIND(1D0)) :: QH_LUMPS
2673 REAL(KIND(1D0)) :: wu_int
2674 REAL(KIND(1D0)) :: kclear
2675 REAL(KIND(1D0)) :: kup
2676 REAL(KIND(1D0)) :: ldown
2677 REAL(KIND(1D0)) :: lup
2678 REAL(KIND(1D0)) :: L_mod
2679 REAL(KIND(1D0)) :: mwh
2680 REAL(KIND(1D0)) :: mwstore
2681 REAL(KIND(1D0)) :: NWstate_per_tstep
2682 REAL(KIND(1D0)) :: FAI
2683 REAL(KIND(1D0)) :: PAI
2684 REAL(KIND(1D0)) :: zL
2685 REAL(KIND(1D0)) :: q2_gkg
2686 REAL(KIND(1D0)) :: qe
2687 REAL(KIND(1D0)) :: qf
2688 REAL(KIND(1D0)) :: QF_SAHP
2689 REAL(KIND(1D0)) :: qh
2690 REAL(KIND(1D0)) :: qh_residual
2691 REAL(KIND(1D0)) :: qh_resist
2692 REAL(KIND(1D0)) :: Qm
2693 REAL(KIND(1D0)) :: QmFreez
2694 REAL(KIND(1D0)) :: QmRain
2695 REAL(KIND(1D0)) :: qn
2696 REAL(KIND(1D0)) :: qn_snow
2697 REAL(KIND(1D0)) :: qn_snowfree
2698 REAL(KIND(1D0)) :: qs
2699 REAL(KIND(1D0)) :: RA_h
2700 REAL(KIND(1D0)) :: RS
2701 REAL(KIND(1D0)),
DIMENSION(NSURF) :: rss_surf
2702 REAL(KIND(1D0)) :: RH2
2703 REAL(KIND(1D0)) :: runoffAGveg
2704 REAL(KIND(1D0)) :: runoffAGimpervious
2705 REAL(KIND(1D0)) :: runoff_per_tstep
2706 REAL(KIND(1D0)) :: runoffPipes
2707 REAL(KIND(1D0)) :: runoffSoil_per_tstep
2708 REAL(KIND(1D0)) :: runoffwaterbody
2709 REAL(KIND(1D0)) :: smd
2710 REAL(KIND(1D0)) :: SoilState
2711 REAL(KIND(1D0)) :: state_per_tstep
2712 REAL(KIND(1D0)) :: surf_chang_per_tstep
2713 REAL(KIND(1D0)) :: swe
2714 REAL(KIND(1D0)) :: t2_C
2715 REAL(KIND(1D0)) :: TSfc_C
2716 REAL(KIND(1D0)) :: TempVeg
2717 REAL(KIND(1D0)) :: tot_chang_per_tstep
2718 REAL(KIND(1D0)) :: TStar
2719 REAL(KIND(1D0)) :: tsurf
2720 REAL(KIND(1D0)) :: UStar
2721 REAL(KIND(1D0)) :: VPD_Pa
2722 REAL(KIND(1D0)) :: z0m
2723 REAL(KIND(1D0)) :: zdm
2724 REAL(KIND(1D0)) :: ZENITH_deg
2725 REAL(KIND(1D0)) :: zH
2727 REAL(KIND(1D0)),
DIMENSION(2) :: SnowRemoval
2728 REAL(KIND(1D0)),
DIMENSION(NSURF) :: wu_surf
2730 REAL(KIND(1D0)),
DIMENSION(nsurf) :: kup_ind_snow
2735 REAL(KIND(1D0)),
DIMENSION(NSURF) :: qn_ind_snow
2737 REAL(KIND(1D0)),
DIMENSION(NSURF) :: runoffSoil
2738 REAL(KIND(1D0)),
DIMENSION(NSURF) :: smd_nsurf
2741 REAL(KIND(1D0)),
DIMENSION(nsurf) :: Tsurf_ind_snow
2744 INTEGER,
DIMENSION(3) :: dayofWeek_id
2747 REAL(KIND(1D0)) :: dq
2748 REAL(KIND(1D0)) :: lvS_J_kg
2749 REAL(KIND(1D0)) :: psyc_hPa
2750 REAL(KIND(1D0)) :: z0v
2751 REAL(KIND(1D0)) :: z0vSnow
2752 REAL(KIND(1D0)) :: RAsnow
2753 REAL(KIND(1D0)) :: RB
2754 REAL(KIND(1D0)) :: runoff_per_interval
2755 REAL(KIND(1D0)) :: s_hPa
2756 REAL(KIND(1D0)) :: sIce_hpa
2757 REAL(KIND(1D0)) :: SoilMoistCap
2759 REAL(KIND(1D0)) :: VegPhenLumps
2760 REAL(KIND(1D0)) :: VPd_hpa
2761 REAL(KIND(1D0)) :: vsmd
2762 REAL(KIND(1D0)) :: ZZD
2764 REAL(KIND(1D0)),
DIMENSION(NSURF) :: deltaQi
2765 REAL(KIND(1D0)),
DIMENSION(NSURF) :: drain_surf
2768 REAL(KIND(1D0)),
DIMENSION(NSURF) :: tsurf_ind
2773 REAL(KIND(1D0)) :: addImpervious = 0
2774 REAL(KIND(1D0)) :: addPipes = 0
2775 REAL(KIND(1D0)) :: addVeg = 0
2776 REAL(KIND(1D0)) :: addWaterBody = 0
2777 REAL(KIND(1D0)),
DIMENSION(NSURF) :: AddWater = 0
2778 REAL(KIND(1D0)),
DIMENSION(NSURF) :: frac_water2runoff = 0
2782 REAL(KIND(1D0)) :: nsh_real
2783 REAL(KIND(1D0)) :: tstep_real
2784 REAL(KIND(1D0)) :: dectime
2787 REAL(KIND(1D0)) :: VegFraction
2788 REAL(KIND(1D0)) :: ImpervFraction
2789 REAL(KIND(1D0)) :: PervFraction
2790 REAL(KIND(1D0)) :: NonWaterFraction
2793 REAL(KIND(1D0)) :: albedo_snow
2803 TYPE(
ohm_state) :: ohmState_prev, ohmState_next
2806 TYPE(
snow_state) :: snowState_prev, snowState_next
2809 TYPE(
hydro_state) :: hydroState_prev, hydroState_next
2810 REAL(KIND(1D0)),
DIMENSION(NSURF) :: ev0_surf
2811 REAL(KIND(1D0)),
DIMENSION(NSURF) :: ev_surf
2812 REAL(KIND(1D0)),
DIMENSION(nlayer) :: ev_roof
2813 REAL(KIND(1D0)),
DIMENSION(nlayer) :: ev_wall
2821 REAL(KIND(1D0)) :: Tair_av_prev, Tair_av_next
2824 LOGICAL :: flag_converge
2825 REAL(KIND(1D0)) :: Ts_iter
2826 REAL(KIND(1D0)) :: dif_tsfc_iter
2827 REAL(KIND(1D0)) :: QH_Init
2834 TYPE(
heat_state) :: heatState_in, heatState_out
2835 REAL(KIND(1D0)),
DIMENSION(nlayer) :: sfr_roof
2836 REAL(KIND(1D0)),
DIMENSION(nlayer) :: sfr_wall
2837 REAL(KIND(1D0)),
DIMENSION(nsurf) :: tsfc0_out_roof
2838 REAL(KIND(1D0)),
DIMENSION(nsurf) :: tsfc0_out_wall
2839 REAL(KIND(1D0)),
DIMENSION(nsurf) :: tsfc0_out_surf
2845 REAL(KIND(1D0)),
DIMENSION(nlayer) :: QS_roof
2847 REAL(KIND(1D0)),
DIMENSION(nlayer, ndepth) :: temp_out_roof
2850 REAL(KIND(1D0)),
DIMENSION(nlayer) :: QG_roof
2851 REAL(KIND(1D0)),
DIMENSION(nlayer) :: QN_roof
2852 REAL(KIND(1D0)),
DIMENSION(nlayer) :: qe_roof
2853 REAL(KIND(1D0)),
DIMENSION(nlayer) :: qh_roof
2854 REAL(KIND(1D0)),
DIMENSION(nlayer) :: qh_resist_roof
2858 REAL(KIND(1D0)),
DIMENSION(nlayer) :: QS_wall
2860 REAL(KIND(1D0)),
DIMENSION(nlayer, ndepth) :: temp_out_wall
2863 REAL(KIND(1D0)),
DIMENSION(nlayer) :: QG_wall
2864 REAL(KIND(1D0)),
DIMENSION(nlayer) :: QN_wall
2865 REAL(KIND(1D0)),
DIMENSION(nlayer) :: qe_wall
2866 REAL(KIND(1D0)),
DIMENSION(nlayer) :: qh_wall
2867 REAL(KIND(1D0)),
DIMENSION(nlayer) :: qh_resist_wall
2871 REAL(KIND(1D0)),
DIMENSION(nsurf, ndepth) :: temp_out_surf
2874 REAL(KIND(1D0)),
DIMENSION(nsurf) :: QG_surf
2875 REAL(KIND(1D0)),
DIMENSION(nsurf) :: QN_surf
2876 REAL(KIND(1D0)),
DIMENSION(nsurf) :: qs_surf
2877 REAL(KIND(1D0)),
DIMENSION(nsurf) :: qe0_surf
2878 REAL(KIND(1D0)),
DIMENSION(nsurf) :: qe_surf
2879 REAL(KIND(1D0)),
DIMENSION(nsurf) :: qh_surf
2880 REAL(KIND(1D0)),
DIMENSION(nsurf) :: qh_resist_surf
2889 REAL(KIND(1D0)) :: ratio_iter
2891 REAL(KIND(1D0)) :: g_kdown
2892 REAL(KIND(1D0)) :: g_dq
2893 REAL(KIND(1D0)) :: g_ta
2894 REAL(KIND(1D0)) :: g_smd
2895 REAL(KIND(1D0)) :: g_lai
2899 INTEGER,
PARAMETER :: nz = 90
2900 INTEGER,
PARAMETER :: AerodynamicResistanceMethod = 2
2901 INTEGER,
PARAMETER :: BaseTMethod = 2
2902 INTEGER,
PARAMETER :: DiagQN = 0
2903 INTEGER,
PARAMETER :: DiagQS = 0
2904 INTEGER,
PARAMETER :: EvapMethod = 2
2905 INTEGER,
PARAMETER :: LAImethod = 1
2906 REAL(KIND(1D0)),
PARAMETER :: BaseT_HC = 18.2
2909 REAL(KIND(1D0)) :: FAIBldg_use
2910 REAL(KIND(1D0)) :: FAIEveTree_use
2911 REAL(KIND(1D0)) :: FAIDecTree_use
2914 ALLOCATE (hydrostate_prev%soilstore_roof(nlayer))
2915 ALLOCATE (hydrostate_prev%state_roof(nlayer))
2916 ALLOCATE (hydrostate_prev%soilstore_wall(nlayer))
2917 ALLOCATE (hydrostate_prev%state_wall(nlayer))
2918 ALLOCATE (hydrostate_next%soilstore_roof(nlayer))
2919 ALLOCATE (hydrostate_next%state_roof(nlayer))
2920 ALLOCATE (hydrostate_next%soilstore_wall(nlayer))
2921 ALLOCATE (hydrostate_next%state_wall(nlayer))
2923 ALLOCATE (heatstate_in%temp_roof(nlayer,
ndepth))
2924 ALLOCATE (heatstate_in%temp_wall(nlayer,
ndepth))
2925 ALLOCATE (heatstate_in%tsfc_roof(nlayer))
2926 ALLOCATE (heatstate_in%tsfc_wall(nlayer))
2927 ALLOCATE (heatstate_in%tsfc_surf(
nsurf))
2930 ALLOCATE (heatstate_out%temp_roof(nlayer,
ndepth))
2931 ALLOCATE (heatstate_out%temp_wall(nlayer,
ndepth))
2932 ALLOCATE (heatstate_out%tsfc_roof(nlayer))
2933 ALLOCATE (heatstate_out%tsfc_wall(nlayer))
2934 ALLOCATE (heatstate_out%tsfc_surf(
nsurf))
2941 siteinfo%gridiv = gridiv
2942 siteinfo%timezone = timezone
2943 siteinfo%surfacearea = surfacearea
2945 siteinfo%z0m_in = z0m_in
2946 siteinfo%zdm_in = zdm_in
2947 siteinfo%pipecapacity = pipecapacity
2948 siteinfo%runofftowater = runofftowater
2949 siteinfo%narp_trans_site = narp_trans_site
2950 siteinfo%CO2PointSource = co2pointsource
2951 siteinfo%flowchange = flowchange
2953 forcing%kdown = kdown
2954 forcing%ldown = ldown_obs
2956 forcing%pres = press_hpa
2958 forcing%rain = precip
2960 forcing%fcld = fcld_obs
2961 forcing%LAI_obs = lai_obs
2962 forcing%snowfrac = snowfrac_obs
2964 forcing%qn1_obs = qn1_obs
2965 forcing%qs_obs = qs_obs
2966 forcing%qf_obs = qf_obs
2967 forcing%Tair = tair_av
2968 forcing%temp_c = temp_c
2976 timer%tstep_prev = tstep_prev
2977 timer%dt_since_start = dt_since_start
2979 methodprm%Diagnose = diagnose
2980 methodprm%DiagMethod = diagmethod
2981 methodprm%EmissionsMethod = emissionsmethod
2982 methodprm%RoughLenHeatMethod = roughlenheatmethod
2983 methodprm%RoughLenMomMethod = roughlenmommethod
2984 methodprm%FAIMethod = faimethod
2985 methodprm%SMDMethod = smdmethod
2986 methodprm%WaterUseMethod = waterusemethod
2987 methodprm%NetRadiationMethod = netradiationmethod
2988 methodprm%StabilityMethod = stabilitymethod
2989 methodprm%StorageHeatMethod = storageheatmethod
2990 methodprm%SnowUse = snowuse
2991 methodprm%use_sw_direct_albedo = use_sw_direct_albedo
2992 methodprm%ohmIncQF = ohmincqf
2994 lumpsprm%raincover = raincover
2995 lumpsprm%rainmaxres = rainmaxres
2996 lumpsprm%drainrt = drainrt
2997 lumpsprm%veg_type = veg_type
3000 ALLOCATE (ehcprm%soil_storecap_roof(nlayer))
3001 ALLOCATE (ehcprm%soil_storecap_wall(nlayer))
3002 ALLOCATE (ehcprm%state_limit_roof(nlayer))
3003 ALLOCATE (ehcprm%state_limit_wall(nlayer))
3004 ALLOCATE (ehcprm%wet_thresh_roof(nlayer))
3005 ALLOCATE (ehcprm%wet_thresh_wall(nlayer))
3006 ALLOCATE (ehcprm%tin_roof(nlayer))
3007 ALLOCATE (ehcprm%tin_wall(nlayer))
3008 ALLOCATE (ehcprm%tin_surf(nlayer))
3009 ALLOCATE (ehcprm%k_roof(nlayer,
ndepth))
3010 ALLOCATE (ehcprm%k_wall(nlayer,
ndepth))
3011 ALLOCATE (ehcprm%k_surf(nlayer,
ndepth))
3012 ALLOCATE (ehcprm%cp_roof(nlayer,
ndepth))
3013 ALLOCATE (ehcprm%cp_wall(nlayer,
ndepth))
3014 ALLOCATE (ehcprm%cp_surf(nlayer,
ndepth))
3015 ALLOCATE (ehcprm%dz_roof(nlayer,
ndepth))
3016 ALLOCATE (ehcprm%dz_wall(nlayer,
ndepth))
3017 ALLOCATE (ehcprm%dz_surf(nlayer,
ndepth))
3018 ehcprm%soil_storecap_roof = soilstorecap_roof
3019 ehcprm%soil_storecap_wall = soilstorecap_wall
3020 ehcprm%state_limit_roof = statelimit_roof
3021 ehcprm%state_limit_wall = statelimit_wall
3022 ehcprm%wet_thresh_roof = wetthresh_roof
3023 ehcprm%wet_thresh_wall = wetthresh_wall
3024 ehcprm%tin_roof = tin_roof
3025 ehcprm%tin_wall = tin_wall
3026 ehcprm%tin_surf = tin_surf
3027 ehcprm%k_roof = k_roof
3028 ehcprm%k_wall = k_wall
3029 ehcprm%k_surf = k_surf
3030 ehcprm%cp_roof = cp_roof
3031 ehcprm%cp_wall = cp_wall
3032 ehcprm%cp_surf = cp_surf
3033 ehcprm%dz_roof = dz_roof
3034 ehcprm%dz_wall = dz_wall
3035 ehcprm%dz_surf = dz_surf
3037 ALLOCATE (spartacusprm%height(nlayer + 1))
3038 spartacusprm%air_ext_lw = air_ext_lw
3039 spartacusprm%air_ext_sw = air_ext_sw
3040 spartacusprm%air_ssa_lw = air_ssa_lw
3041 spartacusprm%air_ssa_sw = air_ssa_sw
3042 spartacusprm%veg_ssa_lw = veg_ssa_lw
3043 spartacusprm%veg_ssa_sw = veg_ssa_sw
3044 spartacusprm%height = height
3045 spartacusprm%ground_albedo_dir_mult_fact = ground_albedo_dir_mult_fact
3046 spartacusprm%n_stream_lw_urban = n_stream_lw_urban
3047 spartacusprm%n_stream_sw_urban = n_stream_sw_urban
3048 spartacusprm%n_vegetation_region_urban = n_vegetation_region_urban
3049 spartacusprm%sw_dn_direct_frac = sw_dn_direct_frac
3050 spartacusprm%veg_contact_fraction_const = veg_contact_fraction_const
3051 spartacusprm%veg_fsd_const = veg_fsd_const
3053 ALLOCATE (spartacuslayerprm%building_frac(nlayer))
3054 ALLOCATE (spartacuslayerprm%building_scale(nlayer))
3055 ALLOCATE (spartacuslayerprm%veg_frac(nlayer))
3056 ALLOCATE (spartacuslayerprm%veg_scale(nlayer))
3057 ALLOCATE (spartacuslayerprm%alb_roof(nlayer))
3058 ALLOCATE (spartacuslayerprm%emis_roof(nlayer))
3059 ALLOCATE (spartacuslayerprm%alb_wall(nlayer))
3060 ALLOCATE (spartacuslayerprm%emis_wall(nlayer))
3061 ALLOCATE (spartacuslayerprm%roof_albedo_dir_mult_fact(
nspec, nlayer))
3062 ALLOCATE (spartacuslayerprm%wall_specular_frac(
nspec, nlayer))
3063 spartacuslayerprm%building_frac = building_frac
3064 spartacuslayerprm%building_scale = building_scale
3065 spartacuslayerprm%veg_frac = veg_frac
3066 spartacuslayerprm%veg_scale = veg_scale
3067 spartacuslayerprm%alb_roof = alb_roof
3068 spartacuslayerprm%emis_roof = emis_roof
3069 spartacuslayerprm%alb_wall = alb_wall
3070 spartacuslayerprm%emis_wall = emis_wall
3071 spartacuslayerprm%roof_albedo_dir_mult_fact = roof_albedo_dir_mult_fact
3072 spartacuslayerprm%wall_specular_frac = wall_specular_frac
3074 ahemisprm%startdls = startdls
3075 ahemisprm%enddls = enddls
3076 ahemisprm%anthroheat%qf0_beu_working = qf0_beu(1)
3077 ahemisprm%anthroheat%qf0_beu_holiday = qf0_beu(2)
3078 ahemisprm%anthroheat%qf_a_working = qf_a(1)
3079 ahemisprm%anthroheat%qf_a_holiday = qf_a(2)
3080 ahemisprm%anthroheat%qf_b_working = qf_b(1)
3081 ahemisprm%anthroheat%qf_b_holiday = qf_b(2)
3082 ahemisprm%anthroheat%qf_c_working = qf_c(1)
3083 ahemisprm%anthroheat%qf_c_holiday = qf_c(2)
3084 ahemisprm%anthroheat%baset_cooling_working = baset_cooling(1)
3085 ahemisprm%anthroheat%baset_cooling_holiday = baset_cooling(2)
3086 ahemisprm%anthroheat%baset_heating_working = baset_heating(1)
3087 ahemisprm%anthroheat%baset_heating_holiday = baset_heating(2)
3088 ahemisprm%anthroheat%popdensdaytime_working = popdensdaytime(1)
3089 ahemisprm%anthroheat%popdensdaytime_holiday = popdensdaytime(2)
3090 ahemisprm%anthroheat%popdensnighttime = popdensnighttime
3091 ahemisprm%anthroheat%popprof_24hr_working = popprof_24hr(:, 1)
3092 ahemisprm%anthroheat%popprof_24hr_holiday = popprof_24hr(:, 2)
3093 ahemisprm%anthroheat%ah_min_working = ah_min(1)
3094 ahemisprm%anthroheat%ah_min_holiday = ah_min(2)
3095 ahemisprm%anthroheat%ahprof_24hr_working = ahprof_24hr(:, 1)
3096 ahemisprm%anthroheat%ahprof_24hr_holiday = ahprof_24hr(:, 2)
3097 ahemisprm%anthroheat%ah_slope_cooling_working = ah_slope_cooling(1)
3098 ahemisprm%anthroheat%ah_slope_cooling_holiday = ah_slope_cooling(2)
3099 ahemisprm%anthroheat%ah_slope_heating_working = ah_slope_heating(1)
3100 ahemisprm%anthroheat%ah_slope_heating_holiday = ah_slope_heating(2)
3101 ahemisprm%EF_umolCO2perJ = ef_umolco2perj
3102 ahemisprm%EnEF_v_Jkm = enef_v_jkm
3103 ahemisprm%FrFossilFuel_Heat = frfossilfuel_heat
3104 ahemisprm%FrFossilFuel_NonHeat = frfossilfuel_nonheat
3105 ahemisprm%FcEF_v_kgkm = fcef_v_kgkm
3106 ahemisprm%HumActivity_24hr_working = humactivity_24hr(:, 1)
3107 ahemisprm%HumActivity_24hr_holiday = humactivity_24hr(:, 2)
3108 ahemisprm%MaxFCMetab = maxfcmetab
3109 ahemisprm%MaxQFMetab = maxqfmetab
3110 ahemisprm%MinFCMetab = minfcmetab
3111 ahemisprm%MinQFMetab = minqfmetab
3112 ahemisprm%TrafficRate_working = trafficrate(1)
3113 ahemisprm%TrafficRate_holiday = trafficrate(2)
3114 ahemisprm%TrafficUnits = trafficunits
3115 ahemisprm%TraffProf_24hr_working = traffprof_24hr(:, 1)
3116 ahemisprm%TraffProf_24hr_holiday = traffprof_24hr(:, 2)
3118 irrprm%h_maintain = h_maintain
3122 irrprm%ie_start = ie_start
3123 irrprm%ie_end = ie_end
3124 irrprm%internalwateruse_h = internalwateruse_h
3125 irrprm%irr_daywater%monday_flag = daywat(1)
3126 irrprm%irr_daywater%monday_percent = daywatper(1)
3127 irrprm%irr_daywater%tuesday_flag = daywat(2)
3128 irrprm%irr_daywater%tuesday_percent = daywatper(2)
3129 irrprm%irr_daywater%wednesday_flag = daywat(3)
3130 irrprm%irr_daywater%wednesday_percent = daywatper(3)
3131 irrprm%irr_daywater%thursday_flag = daywat(4)
3132 irrprm%irr_daywater%thursday_percent = daywatper(4)
3133 irrprm%irr_daywater%friday_flag = daywat(5)
3134 irrprm%irr_daywater%friday_percent = daywatper(5)
3135 irrprm%irr_daywater%saturday_flag = daywat(6)
3136 irrprm%irr_daywater%saturday_percent = daywatper(6)
3137 irrprm%irr_daywater%sunday_flag = daywat(7)
3138 irrprm%irr_daywater%sunday_percent = daywatper(7)
3139 irrprm%wuprofa_24hr_working = wuprofa_24hr(:, 1)
3140 irrprm%wuprofa_24hr_holiday = wuprofa_24hr(:, 2)
3141 irrprm%wuprofm_24hr_working = wuprofm_24hr(:, 1)
3142 irrprm%wuprofm_24hr_holiday = wuprofm_24hr(:, 2)
3144 snowprm%crwmax = crwmax
3145 snowprm%crwmin = crwmin
3146 snowprm%narp_emis_snow = narp_emis_snow
3147 snowprm%preciplimit = preciplimit
3148 snowprm%preciplimitalb = preciplimitalb
3149 snowprm%snowalbmax = snowalbmax
3150 snowprm%snowalbmin = snowalbmin
3151 snowprm%snowdensmax = snowdensmax
3152 snowprm%snowdensmin = snowdensmin
3153 snowprm%snowlimbldg = snowlimbldg
3154 snowprm%snowlimpaved = snowlimpaved
3155 snowprm%snowpacklimit = snowpacklimit
3156 snowprm%snowprof_24hr_working = snowprof_24hr(:, 1)
3157 snowprm%snowprof_24hr_holiday = snowprof_24hr(:, 2)
3158 snowprm%tau_a = tau_a
3159 snowprm%tau_f = tau_f
3160 snowprm%tau_r = tau_r
3161 snowprm%tempmeltfact = tempmeltfact
3162 snowprm%radmeltfact = radmeltfact
3164 conductanceprm%g_max = g_max
3165 conductanceprm%g_k = g_k
3166 conductanceprm%g_q_base = g_q_base
3167 conductanceprm%g_q_shape = g_q_shape
3168 conductanceprm%g_t = g_t
3169 conductanceprm%g_sm = g_sm
3170 conductanceprm%kmax = kmax
3171 conductanceprm%gsmodel = gsmodel
3172 conductanceprm%s1 = s1
3173 conductanceprm%s2 = s2
3174 conductanceprm%TH = th
3175 conductanceprm%TL = tl
3177 pavedprm%sfr = sfr_surf(
pavsurf)
3179 pavedprm%ohm%chanohm = chanohm(
pavsurf)
3180 pavedprm%ohm%cpanohm = cpanohm(
pavsurf)
3181 pavedprm%ohm%kkanohm = kkanohm(
pavsurf)
3182 pavedprm%ohm%ohm_threshsw = ohm_threshsw(
pavsurf)
3183 pavedprm%ohm%ohm_threshwd = ohm_threshwd(
pavsurf)
3187 pavedprm%ohm%ohm_coef_lc(1)%summer_wet = ohm_coef(
pavsurf, 1, 1)
3189 pavedprm%ohm%ohm_coef_lc(1)%summer_dry = ohm_coef(
pavsurf, 2, 1)
3191 pavedprm%ohm%ohm_coef_lc(1)%winter_wet = ohm_coef(
pavsurf, 3, 1)
3193 pavedprm%ohm%ohm_coef_lc(1)%winter_dry = ohm_coef(
pavsurf, 4, 1)
3196 pavedprm%ohm%ohm_coef_lc(2)%summer_wet = ohm_coef(
pavsurf, 1, 2)
3197 pavedprm%ohm%ohm_coef_lc(2)%summer_dry = ohm_coef(
pavsurf, 2, 2)
3198 pavedprm%ohm%ohm_coef_lc(2)%winter_wet = ohm_coef(
pavsurf, 3, 2)
3199 pavedprm%ohm%ohm_coef_lc(2)%winter_dry = ohm_coef(
pavsurf, 4, 2)
3201 pavedprm%ohm%ohm_coef_lc(3)%summer_wet = ohm_coef(
pavsurf, 1, 3)
3202 pavedprm%ohm%ohm_coef_lc(3)%summer_dry = ohm_coef(
pavsurf, 2, 3)
3203 pavedprm%ohm%ohm_coef_lc(3)%winter_wet = ohm_coef(
pavsurf, 3, 3)
3204 pavedprm%ohm%ohm_coef_lc(3)%winter_dry = ohm_coef(
pavsurf, 4, 3)
3207 pavedprm%soil%soildepth = soildepth(
pavsurf)
3208 pavedprm%soil%soilstorecap = soilstorecap_surf(
pavsurf)
3209 pavedprm%soil%sathydraulicconduct = sathydraulicconduct(
pavsurf)
3210 pavedprm%statelimit = statelimit_surf(
pavsurf)
3211 pavedprm%irrfracpaved = irrfracpaved
3212 pavedprm%wetthresh = wetthresh_surf(
pavsurf)
3219 pavedprm%waterdist%to_paved = waterdist(1,
pavsurf)
3220 pavedprm%waterdist%to_bldg = waterdist(2,
pavsurf)
3221 pavedprm%waterdist%to_evetr = waterdist(3,
pavsurf)
3222 pavedprm%waterdist%to_dectr = waterdist(4,
pavsurf)
3223 pavedprm%waterdist%to_grass = waterdist(5,
pavsurf)
3224 pavedprm%waterdist%to_bsoil = waterdist(6,
pavsurf)
3225 pavedprm%waterdist%to_water = waterdist(7,
pavsurf)
3226 pavedprm%waterdist%to_soilstore = waterdist(8,
pavsurf)
3229 bldgprm%faibldg = faibldg
3230 bldgprm%bldgh = bldgh
3232 bldgprm%ohm%chanohm = chanohm(
bldgsurf)
3233 bldgprm%ohm%cpanohm = cpanohm(
bldgsurf)
3234 bldgprm%ohm%kkanohm = kkanohm(
bldgsurf)
3235 bldgprm%ohm%ohm_threshsw = ohm_threshsw(
bldgsurf)
3236 bldgprm%ohm%ohm_threshwd = ohm_threshwd(
bldgsurf)
3237 bldgprm%ohm%ohm_coef_lc(1)%summer_wet = ohm_coef(
bldgsurf, 1, 1)
3238 bldgprm%ohm%ohm_coef_lc(1)%summer_dry = ohm_coef(
bldgsurf, 2, 1)
3239 bldgprm%ohm%ohm_coef_lc(1)%winter_wet = ohm_coef(
bldgsurf, 3, 1)
3240 bldgprm%ohm%ohm_coef_lc(1)%winter_dry = ohm_coef(
bldgsurf, 4, 1)
3243 bldgprm%ohm%ohm_coef_lc(2)%summer_wet = ohm_coef(
bldgsurf, 1, 2)
3244 bldgprm%ohm%ohm_coef_lc(2)%summer_dry = ohm_coef(
bldgsurf, 2, 2)
3245 bldgprm%ohm%ohm_coef_lc(2)%winter_wet = ohm_coef(
bldgsurf, 3, 2)
3246 bldgprm%ohm%ohm_coef_lc(2)%winter_dry = ohm_coef(
bldgsurf, 4, 2)
3249 bldgprm%ohm%ohm_coef_lc(3)%summer_wet = ohm_coef(
bldgsurf, 1, 3)
3250 bldgprm%ohm%ohm_coef_lc(3)%summer_dry = ohm_coef(
bldgsurf, 2, 3)
3251 bldgprm%ohm%ohm_coef_lc(3)%winter_wet = ohm_coef(
bldgsurf, 3, 3)
3252 bldgprm%ohm%ohm_coef_lc(3)%winter_dry = ohm_coef(
bldgsurf, 4, 3)
3255 bldgprm%soil%soildepth = soildepth(
bldgsurf)
3256 bldgprm%soil%soilstorecap = soilstorecap_surf(
bldgsurf)
3257 bldgprm%soil%sathydraulicconduct = sathydraulicconduct(
bldgsurf)
3258 bldgprm%statelimit = statelimit_surf(
bldgsurf)
3259 bldgprm%irrfracbldgs = irrfracbldgs
3260 bldgprm%wetthresh = wetthresh_surf(
bldgsurf)
3267 bldgprm%waterdist%to_paved = waterdist(1,
bldgsurf)
3268 bldgprm%waterdist%to_bldg = waterdist(2,
bldgsurf)
3269 bldgprm%waterdist%to_evetr = waterdist(3,
bldgsurf)
3270 bldgprm%waterdist%to_dectr = waterdist(4,
bldgsurf)
3271 bldgprm%waterdist%to_grass = waterdist(5,
bldgsurf)
3272 bldgprm%waterdist%to_bsoil = waterdist(6,
bldgsurf)
3273 bldgprm%waterdist%to_water = waterdist(7,
bldgsurf)
3274 bldgprm%waterdist%to_soilstore = waterdist(8,
bldgsurf)
3278 dectrprm%faidectree = faidectree
3279 dectrprm%dectreeh = dectreeh
3280 dectrprm%pormin_dec = pormin_dec
3281 dectrprm%pormax_dec = pormax_dec
3282 dectrprm%alb_min = albmin_dectr
3283 dectrprm%alb_max = albmax_dectr
3284 dectrprm%ohm%chanohm = chanohm(
decidsurf)
3285 dectrprm%ohm%cpanohm = cpanohm(
decidsurf)
3286 dectrprm%ohm%kkanohm = kkanohm(
decidsurf)
3287 dectrprm%ohm%ohm_threshsw = ohm_threshsw(
decidsurf)
3288 dectrprm%ohm%ohm_threshwd = ohm_threshwd(
decidsurf)
3290 dectrprm%ohm%ohm_coef_lc(1)%summer_wet = ohm_coef(
decidsurf, 1, 1)
3291 dectrprm%ohm%ohm_coef_lc(1)%summer_dry = ohm_coef(
decidsurf, 2, 1)
3292 dectrprm%ohm%ohm_coef_lc(1)%winter_wet = ohm_coef(
decidsurf, 3, 1)
3293 dectrprm%ohm%ohm_coef_lc(1)%winter_dry = ohm_coef(
decidsurf, 4, 1)
3296 dectrprm%ohm%ohm_coef_lc(2)%summer_wet = ohm_coef(
decidsurf, 1, 2)
3297 dectrprm%ohm%ohm_coef_lc(2)%summer_dry = ohm_coef(
decidsurf, 2, 2)
3298 dectrprm%ohm%ohm_coef_lc(2)%winter_wet = ohm_coef(
decidsurf, 3, 2)
3299 dectrprm%ohm%ohm_coef_lc(2)%winter_dry = ohm_coef(
decidsurf, 4, 2)
3302 dectrprm%ohm%ohm_coef_lc(3)%summer_wet = ohm_coef(
decidsurf, 1, 3)
3303 dectrprm%ohm%ohm_coef_lc(3)%summer_dry = ohm_coef(
decidsurf, 2, 3)
3304 dectrprm%ohm%ohm_coef_lc(3)%winter_wet = ohm_coef(
decidsurf, 3, 3)
3305 dectrprm%ohm%ohm_coef_lc(3)%winter_dry = ohm_coef(
decidsurf, 4, 3)
3308 dectrprm%soil%soildepth = soildepth(
decidsurf)
3309 dectrprm%soil%soilstorecap = soilstorecap_surf(
decidsurf)
3310 dectrprm%soil%sathydraulicconduct = sathydraulicconduct(
decidsurf)
3312 dectrprm%capmax_dec = capmax_dec
3313 dectrprm%capmin_dec = capmin_dec
3314 dectrprm%irrfracdectr = irrfracdectr
3315 dectrprm%wetthresh = wetthresh_surf(
decidsurf)
3316 dectrprm%bioco2%beta_bioco2 = beta_bioco2(
ivdecid)
3317 dectrprm%bioco2%beta_enh_bioco2 = beta_enh_bioco2(
ivdecid)
3318 dectrprm%bioco2%alpha_bioco2 = alpha_bioco2(
ivdecid)
3319 dectrprm%bioco2%alpha_enh_bioco2 = alpha_enh_bioco2(
ivdecid)
3320 dectrprm%bioco2%resp_a = resp_a(
ivdecid)
3321 dectrprm%bioco2%resp_b = resp_b(
ivdecid)
3322 dectrprm%bioco2%min_res_bioCO2 = min_res_bioco2(
ivdecid)
3323 dectrprm%bioco2%theta_bioco2 = theta_bioco2(
ivdecid)
3324 dectrprm%maxconductance = maxconductance(
ivdecid)
3325 dectrprm%lai%baset = baset(
ivdecid)
3326 dectrprm%lai%gddfull = gddfull(
ivdecid)
3327 dectrprm%lai%basete = basete(
ivdecid)
3328 dectrprm%lai%sddfull = sddfull(
ivdecid)
3329 dectrprm%lai%laimin = laimin(
ivdecid)
3330 dectrprm%lai%laimax = laimax(
ivdecid)
3331 dectrprm%lai%laipower = laipower(:,
ivdecid)
3332 dectrprm%lai%laitype = laitype(
ivdecid)
3339 dectrprm%waterdist%to_paved = waterdist(1,
decidsurf)
3340 dectrprm%waterdist%to_bldg = waterdist(2,
decidsurf)
3341 dectrprm%waterdist%to_evetr = waterdist(3,
decidsurf)
3342 dectrprm%waterdist%to_dectr = waterdist(4,
decidsurf)
3343 dectrprm%waterdist%to_grass = waterdist(5,
decidsurf)
3344 dectrprm%waterdist%to_bsoil = waterdist(6,
decidsurf)
3345 dectrprm%waterdist%to_water = waterdist(7,
decidsurf)
3346 dectrprm%waterdist%to_soilstore = waterdist(8,
decidsurf)
3350 evetrprm%faievetree = faievetree
3351 evetrprm%evetreeh = evetreeh
3352 evetrprm%alb_min = albmin_evetr
3353 evetrprm%alb_max = albmax_evetr
3354 evetrprm%ohm%chanohm = chanohm(
conifsurf)
3355 evetrprm%ohm%cpanohm = cpanohm(
conifsurf)
3356 evetrprm%ohm%kkanohm = kkanohm(
conifsurf)
3357 evetrprm%ohm%ohm_threshsw = ohm_threshsw(
conifsurf)
3358 evetrprm%ohm%ohm_threshwd = ohm_threshwd(
conifsurf)
3359 evetrprm%ohm%ohm_coef_lc(1)%summer_wet = ohm_coef(
conifsurf, 1, 1)
3360 evetrprm%ohm%ohm_coef_lc(1)%summer_dry = ohm_coef(
conifsurf, 2, 1)
3361 evetrprm%ohm%ohm_coef_lc(1)%winter_wet = ohm_coef(
conifsurf, 3, 1)
3362 evetrprm%ohm%ohm_coef_lc(1)%winter_dry = ohm_coef(
conifsurf, 4, 1)
3365 evetrprm%ohm%ohm_coef_lc(2)%summer_wet = ohm_coef(
conifsurf, 1, 2)
3366 evetrprm%ohm%ohm_coef_lc(2)%summer_dry = ohm_coef(
conifsurf, 2, 2)
3367 evetrprm%ohm%ohm_coef_lc(2)%winter_wet = ohm_coef(
conifsurf, 3, 2)
3368 evetrprm%ohm%ohm_coef_lc(2)%winter_dry = ohm_coef(
conifsurf, 4, 2)
3371 evetrprm%ohm%ohm_coef_lc(3)%summer_wet = ohm_coef(
conifsurf, 1, 3)
3372 evetrprm%ohm%ohm_coef_lc(3)%summer_dry = ohm_coef(
conifsurf, 2, 3)
3373 evetrprm%ohm%ohm_coef_lc(3)%winter_wet = ohm_coef(
conifsurf, 3, 3)
3374 evetrprm%ohm%ohm_coef_lc(3)%winter_dry = ohm_coef(
conifsurf, 4, 3)
3377 evetrprm%soil%soildepth = soildepth(
conifsurf)
3378 evetrprm%soil%soilstorecap = soilstorecap_surf(
conifsurf)
3379 evetrprm%soil%sathydraulicconduct = sathydraulicconduct(
conifsurf)
3380 evetrprm%statelimit = statelimit_surf(
conifsurf)
3381 evetrprm%irrfracevetr = irrfracevetr
3382 evetrprm%wetthresh = wetthresh_surf(
conifsurf)
3383 evetrprm%bioco2%beta_bioco2 = beta_bioco2(
ivconif)
3384 evetrprm%bioco2%beta_enh_bioco2 = beta_enh_bioco2(
ivconif)
3385 evetrprm%bioco2%alpha_bioco2 = alpha_bioco2(
ivconif)
3386 evetrprm%bioco2%alpha_enh_bioco2 = alpha_enh_bioco2(
ivconif)
3387 evetrprm%bioco2%resp_a = resp_a(
ivconif)
3388 evetrprm%bioco2%resp_b = resp_b(
ivconif)
3389 evetrprm%bioco2%min_res_bioCO2 = min_res_bioco2(
ivconif)
3390 evetrprm%bioco2%theta_bioco2 = theta_bioco2(
ivconif)
3391 evetrprm%maxconductance = maxconductance(
ivconif)
3392 evetrprm%lai%baset = baset(
ivconif)
3393 evetrprm%lai%gddfull = gddfull(
ivconif)
3394 evetrprm%lai%basete = basete(
ivconif)
3395 evetrprm%lai%sddfull = sddfull(
ivconif)
3396 evetrprm%lai%laimin = laimin(
ivconif)
3397 evetrprm%lai%laimax = laimax(
ivconif)
3398 evetrprm%lai%laipower = laipower(:,
ivconif)
3399 evetrprm%lai%laitype = laitype(
ivconif)
3406 evetrprm%waterdist%to_paved = waterdist(1,
conifsurf)
3407 evetrprm%waterdist%to_bldg = waterdist(2,
conifsurf)
3408 evetrprm%waterdist%to_evetr = waterdist(3,
conifsurf)
3409 evetrprm%waterdist%to_dectr = waterdist(4,
conifsurf)
3410 evetrprm%waterdist%to_grass = waterdist(5,
conifsurf)
3411 evetrprm%waterdist%to_bsoil = waterdist(6,
conifsurf)
3412 evetrprm%waterdist%to_water = waterdist(7,
conifsurf)
3413 evetrprm%waterdist%to_soilstore = waterdist(8,
conifsurf)
3417 grassprm%alb_min = albmin_grass
3418 grassprm%alb_max = albmax_grass
3419 grassprm%ohm%chanohm = chanohm(
grasssurf)
3420 grassprm%ohm%cpanohm = cpanohm(
grasssurf)
3421 grassprm%ohm%kkanohm = kkanohm(
grasssurf)
3422 grassprm%ohm%ohm_threshsw = ohm_threshsw(
grasssurf)
3423 grassprm%ohm%ohm_threshwd = ohm_threshwd(
grasssurf)
3424 grassprm%ohm%ohm_coef_lc(1)%summer_wet = ohm_coef(
grasssurf, 1, 1)
3425 grassprm%ohm%ohm_coef_lc(1)%summer_dry = ohm_coef(
grasssurf, 2, 1)
3426 grassprm%ohm%ohm_coef_lc(1)%winter_wet = ohm_coef(
grasssurf, 3, 1)
3427 grassprm%ohm%ohm_coef_lc(1)%winter_dry = ohm_coef(
grasssurf, 4, 1)
3430 grassprm%ohm%ohm_coef_lc(2)%summer_wet = ohm_coef(
grasssurf, 1, 2)
3431 grassprm%ohm%ohm_coef_lc(2)%summer_dry = ohm_coef(
grasssurf, 2, 2)
3432 grassprm%ohm%ohm_coef_lc(2)%winter_wet = ohm_coef(
grasssurf, 3, 2)
3433 grassprm%ohm%ohm_coef_lc(2)%winter_dry = ohm_coef(
grasssurf, 4, 2)
3436 grassprm%ohm%ohm_coef_lc(3)%summer_wet = ohm_coef(
grasssurf, 1, 3)
3437 grassprm%ohm%ohm_coef_lc(3)%summer_dry = ohm_coef(
grasssurf, 2, 3)
3438 grassprm%ohm%ohm_coef_lc(3)%winter_wet = ohm_coef(
grasssurf, 3, 3)
3439 grassprm%ohm%ohm_coef_lc(3)%winter_dry = ohm_coef(
grasssurf, 4, 3)
3442 grassprm%soil%soildepth = soildepth(
grasssurf)
3443 grassprm%soil%soilstorecap = soilstorecap_surf(
grasssurf)
3444 grassprm%soil%sathydraulicconduct = sathydraulicconduct(
grasssurf)
3445 grassprm%statelimit = statelimit_surf(
grasssurf)
3446 grassprm%irrfracgrass = irrfracgrass
3447 grassprm%wetthresh = wetthresh_surf(
grasssurf)
3448 grassprm%bioco2%beta_bioco2 = beta_bioco2(
ivgrass)
3449 grassprm%bioco2%beta_enh_bioco2 = beta_enh_bioco2(
ivgrass)
3450 grassprm%bioco2%alpha_bioco2 = alpha_bioco2(
ivgrass)
3451 grassprm%bioco2%alpha_enh_bioco2 = alpha_enh_bioco2(
ivgrass)
3452 grassprm%bioco2%resp_a = resp_a(
ivgrass)
3453 grassprm%bioco2%resp_b = resp_b(
ivgrass)
3454 grassprm%bioco2%min_res_bioCO2 = min_res_bioco2(
ivgrass)
3455 grassprm%bioco2%theta_bioco2 = theta_bioco2(
ivgrass)
3456 grassprm%maxconductance = maxconductance(
ivgrass)
3457 grassprm%lai%baset = baset(
ivgrass)
3458 grassprm%lai%gddfull = gddfull(
ivgrass)
3459 grassprm%lai%basete = basete(
ivgrass)
3460 grassprm%lai%sddfull = sddfull(
ivgrass)
3461 grassprm%lai%laimin = laimin(
ivgrass)
3462 grassprm%lai%laimax = laimax(
ivgrass)
3463 grassprm%lai%laipower = laipower(:,
ivgrass)
3464 grassprm%lai%laitype = laitype(
ivgrass)
3471 grassprm%waterdist%to_paved = waterdist(1,
grasssurf)
3472 grassprm%waterdist%to_bldg = waterdist(2,
grasssurf)
3473 grassprm%waterdist%to_evetr = waterdist(3,
grasssurf)
3474 grassprm%waterdist%to_dectr = waterdist(4,
grasssurf)
3475 grassprm%waterdist%to_grass = waterdist(5,
grasssurf)
3476 grassprm%waterdist%to_bsoil = waterdist(6,
grasssurf)
3477 grassprm%waterdist%to_water = waterdist(7,
grasssurf)
3478 grassprm%waterdist%to_soilstore = waterdist(8,
grasssurf)
3482 bsoilprm%ohm%chanohm = chanohm(
bsoilsurf)
3483 bsoilprm%ohm%cpanohm = cpanohm(
bsoilsurf)
3484 bsoilprm%ohm%kkanohm = kkanohm(
bsoilsurf)
3485 bsoilprm%ohm%ohm_threshsw = ohm_threshsw(
bsoilsurf)
3486 bsoilprm%ohm%ohm_threshwd = ohm_threshwd(
bsoilsurf)
3487 bsoilprm%ohm%ohm_coef_lc(1)%summer_wet = ohm_coef(
bsoilsurf, 1, 1)
3488 bsoilprm%ohm%ohm_coef_lc(1)%summer_dry = ohm_coef(
bsoilsurf, 2, 1)
3489 bsoilprm%ohm%ohm_coef_lc(1)%winter_wet = ohm_coef(
bsoilsurf, 3, 1)
3490 bsoilprm%ohm%ohm_coef_lc(1)%winter_dry = ohm_coef(
bsoilsurf, 4, 1)
3493 bsoilprm%ohm%ohm_coef_lc(2)%summer_wet = ohm_coef(
bsoilsurf, 1, 2)
3494 bsoilprm%ohm%ohm_coef_lc(2)%summer_dry = ohm_coef(
bsoilsurf, 2, 2)
3495 bsoilprm%ohm%ohm_coef_lc(2)%winter_wet = ohm_coef(
bsoilsurf, 3, 2)
3496 bsoilprm%ohm%ohm_coef_lc(2)%winter_dry = ohm_coef(
bsoilsurf, 4, 2)
3499 bsoilprm%ohm%ohm_coef_lc(3)%summer_wet = ohm_coef(
bsoilsurf, 1, 3)
3500 bsoilprm%ohm%ohm_coef_lc(3)%summer_dry = ohm_coef(
bsoilsurf, 2, 3)
3501 bsoilprm%ohm%ohm_coef_lc(3)%winter_wet = ohm_coef(
bsoilsurf, 3, 3)
3502 bsoilprm%ohm%ohm_coef_lc(3)%winter_dry = ohm_coef(
bsoilsurf, 4, 3)
3505 bsoilprm%soil%soildepth = soildepth(
bsoilsurf)
3506 bsoilprm%soil%soilstorecap = soilstorecap_surf(
bsoilsurf)
3507 bsoilprm%soil%sathydraulicconduct = sathydraulicconduct(
bsoilsurf)
3508 bsoilprm%statelimit = statelimit_surf(
bsoilsurf)
3509 bsoilprm%irrfracbsoil = irrfracbsoil
3510 bsoilprm%wetthresh = wetthresh_surf(
bsoilsurf)
3517 bsoilprm%waterdist%to_paved = waterdist(1,
bsoilsurf)
3518 bsoilprm%waterdist%to_bldg = waterdist(2,
bsoilsurf)
3519 bsoilprm%waterdist%to_evetr = waterdist(3,
bsoilsurf)
3520 bsoilprm%waterdist%to_dectr = waterdist(4,
bsoilsurf)
3521 bsoilprm%waterdist%to_grass = waterdist(5,
bsoilsurf)
3522 bsoilprm%waterdist%to_bsoil = waterdist(6,
bsoilsurf)
3523 bsoilprm%waterdist%to_water = waterdist(7,
bsoilsurf)
3524 bsoilprm%waterdist%to_soilstore = waterdist(8,
bsoilsurf)
3528 waterprm%ohm%chanohm = chanohm(
watersurf)
3529 waterprm%ohm%cpanohm = cpanohm(
watersurf)
3530 waterprm%ohm%kkanohm = kkanohm(
watersurf)
3531 waterprm%ohm%ohm_threshsw = ohm_threshsw(
watersurf)
3532 waterprm%ohm%ohm_threshwd = ohm_threshwd(
watersurf)
3533 waterprm%ohm%ohm_coef_lc(1)%summer_wet = ohm_coef(
watersurf, 1, 1)
3534 waterprm%ohm%ohm_coef_lc(1)%summer_dry = ohm_coef(
watersurf, 2, 1)
3535 waterprm%ohm%ohm_coef_lc(1)%winter_wet = ohm_coef(
watersurf, 3, 1)
3536 waterprm%ohm%ohm_coef_lc(1)%winter_dry = ohm_coef(
watersurf, 4, 1)
3539 waterprm%ohm%ohm_coef_lc(2)%summer_wet = ohm_coef(
watersurf, 1, 2)
3540 waterprm%ohm%ohm_coef_lc(2)%summer_dry = ohm_coef(
watersurf, 2, 2)
3541 waterprm%ohm%ohm_coef_lc(2)%winter_wet = ohm_coef(
watersurf, 3, 2)
3542 waterprm%ohm%ohm_coef_lc(2)%winter_dry = ohm_coef(
watersurf, 4, 2)
3545 waterprm%ohm%ohm_coef_lc(3)%summer_wet = ohm_coef(
watersurf, 1, 3)
3546 waterprm%ohm%ohm_coef_lc(3)%summer_dry = ohm_coef(
watersurf, 2, 3)
3547 waterprm%ohm%ohm_coef_lc(3)%winter_wet = ohm_coef(
watersurf, 3, 3)
3548 waterprm%ohm%ohm_coef_lc(3)%winter_dry = ohm_coef(
watersurf, 4, 3)
3551 waterprm%soil%soildepth = soildepth(
watersurf)
3552 waterprm%soil%soilstorecap = soilstorecap_surf(
watersurf)
3553 waterprm%soil%sathydraulicconduct = sathydraulicconduct(
watersurf)
3554 waterprm%statelimit = statelimit_surf(
watersurf)
3555 waterprm%irrfracwater = irrfracwater
3565 anthroheatstate%HDD_id = hdd_id
3569 ALLOCATE (hydrostate%soilstore_roof(nlayer))
3570 ALLOCATE (hydrostate%state_roof(nlayer))
3571 ALLOCATE (hydrostate%soilstore_wall(nlayer))
3572 ALLOCATE (hydrostate%state_wall(nlayer))
3573 hydrostate%soilstore_roof = soilstore_roof
3574 hydrostate%state_roof = state_roof
3575 hydrostate%soilstore_wall = soilstore_wall
3576 hydrostate%state_wall = state_wall
3577 hydrostate%soilstore_surf = soilstore_surf
3578 hydrostate%state_surf = state_surf
3579 hydrostate%WUDay_id = wuday_id
3581 ALLOCATE (heatstate%temp_roof(nlayer,
ndepth))
3582 ALLOCATE (heatstate%temp_wall(nlayer,
ndepth))
3583 ALLOCATE (heatstate%tsfc_roof(nlayer))
3584 ALLOCATE (heatstate%tsfc_wall(nlayer))
3585 ALLOCATE (heatstate%tsfc_surf(
nsurf))
3587 heatstate%temp_roof = temp_roof
3588 heatstate%temp_wall = temp_wall
3589 heatstate%temp_surf = temp_surf
3590 heatstate%tsfc_roof = tsfc_roof
3591 heatstate%tsfc_wall = tsfc_wall
3592 heatstate%tsfc_surf = tsfc_surf
3593 heatstate%temp_surf = temp_surf
3596 ohmstate%qn_av = qn_av
3597 ohmstate%dqndt = dqndt
3598 ohmstate%qn_s_av = qn_s_av
3599 ohmstate%dqnsdt = dqnsdt
3602 snowstate%snowfallCum = snowfallcum
3603 snowstate%snowalb = snowalb
3604 snowstate%icefrac = icefrac
3605 snowstate%snowdens = snowdens
3606 snowstate%snowfrac = snowfrac
3607 snowstate%snowpack = snowpack
3608 snowstate%snowwater = snowwater
3612 phenstate%lai_id = lai_id
3613 phenstate%SDD_id = sdd_id
3614 phenstate%GDD_id = gdd_id
3615 phenstate%porosity_id = porosity_id
3616 phenstate%decidcap_id = decidcap_id
3617 phenstate%albDecTr_id = albdectr_id
3618 phenstate%albEveTr_id = albevetr_id
3619 phenstate%albGrass_id = albgrass_id
3620 phenstate%Tmin_id = tmin_id
3621 phenstate%Tmax_id = tmax_id
3622 phenstate%lenDay_id = lenday_id
3623 phenstate%StoreDrainPrm = storedrainprm
3630 chsnow_per_interval = 0.
3643 IF (methodprm%snowuse == 0)
THEN
3644 snowstate%SnowDens = 0.
3645 snowstate%SnowFrac = 0.
3646 snowstate%SnowWater = 0.
3647 snowstate%SnowAlb = 0.
3648 snowstate%IceFrac = 0.
3649 snowstate%SnowPack = 0.
3654 ohmstate_prev = ohmstate
3656 snowstate_prev = snowstate
3657 snowstate_prev%snowfrac = merge(forcing%snowfrac, snowstate%SnowFrac, methodprm%NetRadiationMethod == 0)
3659 hydrostate_prev = hydrostate
3667 tair_av_prev = forcing%Tair
3668 phenstate_prev = phenstate
3669 anthroheatstate_prev = anthroheatstate
3673 heatstate_in = heatstate
3685 ohmstate_next = ohmstate
3686 snowstate_next = snowstate
3687 hydrostate_next = hydrostate
3691 hydrostate_next = hydrostate
3702 tair_av_next = forcing%Tair
3703 phenstate_next = phenstate
3704 anthroheatstate_next = anthroheatstate
3707 dataoutlinesnow = -999.
3708 dataoutlineestm = -999.
3709 dataoutlineehc = -999.
3710 dataoutlinersl = -999.
3711 dataoutlinebeers = -999.
3712 dataoutlinedebug = -999.
3713 dataoutlinespartacus = -999.
3714 dataoutlinedailystate = -999.
3721 flag_converge = .false.
3724 heatstate_out = heatstate
3725 tsfc0_out_surf = heatstate%tsfc_surf
3729 IF (storageheatmethod == 5 .OR. netradiationmethod > 1000)
THEN
3730 tsfc0_out_roof = heatstate%tsfc_roof
3731 tsfc0_out_wall = heatstate%tsfc_wall
3744 DO WHILE ((.NOT. flag_converge) .AND. i_iter < max_iter)
3753 timer%id, timer%it, timer%imin, timer%isec, &
3762 nsh, nsh_real, tstep_real)
3772 methodprm%StorageHeatMethod, methodprm%NetRadiationMethod, &
3774 pavedprm%sfr, bldgprm%sfr, evetrprm%sfr, dectrprm%sfr, grassprm%sfr, bsoilprm%sfr, waterprm%sfr, &
3775 spartacuslayerprm%building_frac, spartacuslayerprm%building_scale, spartacusprm%height, &
3776 vegfraction, impervfraction, pervfraction, nonwaterfraction, &
3784 timer%iy, timer%id, siteinfo%lat, &
3792 timer%id, ahemisprm%startDLS, ahemisprm%endDLS, &
3797 tair_av_next =
cal_tair_av(tair_av_prev, timer%dt_since_start, timer%tstep, forcing%temp_c)
3802 IF (methodprm%Diagnose == 1)
WRITE (*, *)
'Calling SUEWS_cal_RoughnessParameters...'
3803 IF (methodprm%Diagnose == 1) print *,
'z0m_in =', z0m_in
3812 methodprm%RoughLenMomMethod, &
3813 methodprm%FAIMethod, &
3814 pavedprm%sfr, bldgprm%sfr, evetrprm%sfr, dectrprm%sfr, grassprm%sfr, bsoilprm%sfr, waterprm%sfr, &
3815 siteinfo%SurfaceArea, &
3816 bldgprm%bldgH, evetrprm%EveTreeH, dectrprm%DecTreeH, &
3817 phenstate_prev%porosity_id, bldgprm%FAIBldg, evetrprm%FAIEveTree, dectrprm%FAIDecTree, &
3818 siteinfo%z0m_in, siteinfo%zdm_in, siteinfo%Z, &
3819 faibldg_use, faievetree_use, faidectree_use, &
3825 IF (methodprm%Diagnose == 1)
WRITE (*, *)
'Calling NARP_cal_SunPosition...'
3832 REAL(timer%iy, KIND(1D0)), &
3833 dectime - timer%tstep/2/86400, &
3834 siteInfo%timezone, siteInfo%lat, siteInfo%lon, siteInfo%alt, &
3835 azimuth, zenith_deg)
3838 IF (methodprm%Diagnose == 1)
WRITE (*, *)
'Calling SUEWS_cal_DailyState...'
3864 timer%iy, timer%id, timer%it, timer%imin, timer%isec, &
3865 timer%tstep, timer%tstep_prev, timer%dt_since_start, dayofweek_id, &
3866 phenstate_prev%Tmin_id, phenstate_prev%Tmax_id, phenstate_prev%lenDay_id, &
3868 methodprm%WaterUseMethod, irrprm%Ie_start, irrprm%Ie_end, &
3870 evetrprm%lai%laitype, dectrprm%lai%laitype, grassprm%lai%laitype, &
3871 nsh_real, forcing%kdown, forcing%Temp_C, forcing%rain, baset_hc, &
3872 ahemisprm%anthroheat%BaseT_Heating_working, ahemisprm%anthroheat%BaseT_Heating_holiday, &
3873 ahemisprm%anthroheat%BaseT_Cooling_working, ahemisprm%anthroheat%BaseT_Cooling_holiday, &
3874 siteinfo%lat, irrprm%Faut, forcing%LAI_obs, &
3875 evetrprm%Alb_Max, dectrprm%Alb_Max, grassprm%Alb_Max, &
3876 evetrprm%Alb_Min, dectrprm%Alb_Min, grassprm%Alb_Min, &
3877 dectrprm%CapMax_dec, dectrprm%CapMin_dec, dectrprm%PorMax_dec, dectrprm%PorMin_dec, &
3878 irrprm%Ie_a, irrprm%Ie_m, &
3879 irrprm%irr_daywater%monday_percent, irrprm%irr_daywater%tuesday_percent, &
3880 irrprm%irr_daywater%wednesday_percent, irrprm%irr_daywater%thursday_percent, &
3881 irrprm%irr_daywater%friday_percent, irrprm%irr_daywater%saturday_percent, &
3882 irrprm%irr_daywater%sunday_percent, &
3883 irrprm%irr_daywater%monday_flag, irrprm%irr_daywater%tuesday_flag, irrprm%irr_daywater%wednesday_flag, &
3884 irrprm%irr_daywater%thursday_flag, irrprm%irr_daywater%friday_flag, irrprm%irr_daywater%saturday_flag, &
3885 irrprm%irr_daywater%sunday_flag, &
3886 evetrprm%lai%BaseT, dectrprm%lai%BaseT, grassprm%lai%BaseT, &
3887 evetrprm%lai%BaseTe, dectrprm%lai%BaseTe, grassprm%lai%BaseTe, &
3888 evetrprm%lai%gddfull, dectrprm%lai%gddfull, grassprm%lai%gddfull, &
3889 evetrprm%lai%sddfull, dectrprm%lai%sddfull, grassprm%lai%sddfull, &
3890 evetrprm%lai%laimin, dectrprm%lai%laimin, grassprm%lai%laimin, &
3891 evetrprm%lai%laimax, dectrprm%lai%laimax, grassprm%lai%laimax, &
3892 evetrprm%lai%laipower, dectrprm%lai%laipower, grassprm%lai%laipower, &
3893 phenstate_prev%DecidCap_id, phenstate_prev%StoreDrainPrm, &
3894 phenstate_prev%LAI_id, phenstate_prev%GDD_id, phenstate_prev%SDD_id, &
3895 phenstate_prev%albDecTr_id, phenstate_prev%albEveTr_id, phenstate_prev%albGrass_id, phenstate_prev%porosity_id, &
3896 anthroheatstate_prev%HDD_id, &
3897 hydrostate_prev%state_surf, hydrostate_prev%soilstore_surf, soilstorecap_surf, irrprm%H_maintain, &
3898 anthroheatstate_next%HDD_id, &
3899 phenstate_next%Tmin_id, phenstate_next%Tmax_id, phenstate_next%lenDay_id, &
3900 phenstate_next%albDecTr_id, phenstate_next%albEveTr_id, phenstate_next%albGrass_id, phenstate_next%porosity_id, &
3901 phenstate_next%DecidCap_id, phenstate_next%StoreDrainPrm, phenstate_next%LAI_id, &
3902 phenstate_next%GDD_id, phenstate_next%SDD_id, &
3903 hydrostate_next%WUDay_id)
3906 IF (methodprm%Diagnose == 1)
WRITE (*, *)
'Calling LUMPS_cal_AtmMoist...'
3912 forcing%Temp_C, forcing%pres, forcing%RH, dectime, &
3914 es_hpa, ea_hpa, vpd_hpa, vpd_pa, dq, dens_dry,
avcp,
avdens)
3917 IF (methodprm%Diagnose == 1)
WRITE (*, *)
'Calling SUEWS_update_SoilMoist...'
3925 pavedprm%sfr, bldgprm%sfr, evetrprm%sfr, dectrprm%sfr, grassprm%sfr, bsoilprm%sfr, waterprm%sfr, &
3926 pavedprm%soil%soilstorecap, bldgprm%soil%soilstorecap, &
3927 evetrprm%soil%soilstorecap, dectrprm%soil%soilstorecap, &
3928 grassprm%soil%soilstorecap, bsoilprm%soil%soilstorecap, waterprm%soil%soilstorecap, &
3929 hydrostate_prev%soilstore_surf, &
3930 soilmoistcap, soilstate, &
3933 IF (methodprm%Diagnose == 1)
WRITE (*, *)
'Calling SUEWS_cal_WaterUse...'
3947 wu_m3, siteinfo%SurfaceArea, &
3948 pavedprm%sfr, bldgprm%sfr, evetrprm%sfr, dectrprm%sfr, grassprm%sfr, bsoilprm%sfr, waterprm%sfr, &
3949 pavedprm%IrrFracPaved, bldgprm%IrrFracBldgs, &
3950 evetrprm%IrrFracEveTr, dectrprm%IrrFracDecTr, grassprm%IrrFracGrass, &
3951 bsoilprm%IrrFracBSoil, waterprm%IrrFracWater, &
3953 irrprm%wuprofa_24hr_working, irrprm%wuprofa_24hr_holiday, &
3954 irrprm%wuprofm_24hr_working, irrprm%wuprofm_24hr_holiday, &
3955 irrprm%InternalWaterUse_h, &
3956 anthroheatstate_next%HDD_id, hydrostate_next%WUDay_id, &
3957 methodprm%WaterUseMethod, nsh, timer%it, timer%imin, dls, &
3958 wu_surf, wu_int, wu_ext)
3972 ahemisprm%anthroheat%ah_min_working, ahemisprm%anthroheat%ah_min_holiday, &
3973 ahemisprm%anthroheat%ahprof_24hr_working, ahemisprm%anthroheat%ahprof_24hr_holiday, &
3974 ahemisprm%anthroheat%ah_slope_cooling_working, ahemisprm%anthroheat%ah_slope_cooling_holiday, &
3975 ahemisprm%anthroheat%ah_slope_heating_working, ahemisprm%anthroheat%ah_slope_heating_holiday, &
3976 siteinfo%CO2PointSource, &
3977 dayofweek_id, dls, ahemisprm%EF_umolCO2perJ, methodprm%EmissionsMethod, ahemisprm%EnEF_v_Jkm, &
3978 ahemisprm%FcEF_v_kgkm, ahemisprm%FrFossilFuel_Heat, ahemisprm%FrFossilFuel_NonHeat, &
3979 anthroheatstate_next%HDD_id, &
3980 ahemisprm%HumActivity_24hr_working, ahemisprm%HumActivity_24hr_holiday, &
3981 timer%imin, timer%it, ahemisprm%MaxFCMetab, ahemisprm%MaxQFMetab, ahemisprm%MinFCMetab, ahemisprm%MinQFMetab, &
3982 ahemisprm%anthroheat%popdensdaytime_working, ahemisprm%anthroheat%popdensdaytime_holiday, &
3983 ahemisprm%anthroheat%popdensnighttime, &
3984 ahemisprm%anthroheat%popprof_24hr_working, ahemisprm%anthroheat%popprof_24hr_holiday, &
3986 ahemisprm%anthroheat%qf0_beu_working, ahemisprm%anthroheat%qf0_beu_holiday, &
3987 ahemisprm%anthroheat%qf_a_working, ahemisprm%anthroheat%qf_a_holiday, &
3988 ahemisprm%anthroheat%qf_b_working, ahemisprm%anthroheat%qf_b_holiday, &
3989 ahemisprm%anthroheat%qf_c_working, ahemisprm%anthroheat%qf_c_holiday, &
3990 forcing%QF_obs, qf_sahp, siteinfo%SurfaceArea, &
3991 ahemisprm%anthroheat%baset_cooling_working, ahemisprm%anthroheat%baset_cooling_holiday, &
3992 ahemisprm%anthroheat%baset_heating_working, ahemisprm%anthroheat%baset_heating_holiday, &
3993 forcing%Temp_C, ahemisprm%TrafficRate_working, ahemisprm%TrafficRate_holiday, &
3994 ahemisprm%TrafficUnits, &
3995 ahemisprm%TraffProf_24hr_working, ahemisprm%TraffProf_24hr_holiday, &
3996 fc_anthro, fc_build, fc_metab, fc_point, fc_traff)
4034 methodprm%StorageHeatMethod, methodprm%NetRadiationMethod, methodprm%SnowUse, &
4035 timer%tstep, nlayer, snowstate_prev%SnowPack, snowprm%tau_a, snowprm%tau_f, snowprm%SnowAlbMax, snowprm%SnowAlbMin, &
4036 methodprm%Diagnose, forcing%ldown, forcing%fcld, &
4037 dectime, zenith_deg, ts_iter, forcing%kdown, forcing%Temp_C, forcing%RH, ea_hpa, forcing%qn1_obs, &
4038 snowstate_prev%snowalb, snowstate_prev%snowFrac, diagqn, &
4039 siteinfo%NARP_TRANS_SITE, snowprm%NARP_EMIS_SNOW, snowstate_prev%IceFrac, &
4040 pavedprm%sfr, bldgprm%sfr, evetrprm%sfr, dectrprm%sfr, grassprm%sfr, bsoilprm%sfr, waterprm%sfr, &
4041 sfr_roof, sfr_wall, &
4042 heatstate_out%tsfc_surf, heatstate_out%tsfc_roof, heatstate_out%tsfc_wall, &
4043 pavedprm%emis, bldgprm%emis, evetrprm%emis, dectrprm%emis, grassprm%emis, bsoilprm%emis, waterprm%emis, &
4044 phenstate_prev%alb, phenstate_next%albDecTr_id, phenstate_next%albEveTr_id, phenstate_next%albGrass_id, &
4046 spartacusprm%n_vegetation_region_urban, &
4047 spartacusprm%n_stream_sw_urban, spartacusprm%n_stream_lw_urban, &
4048 spartacusprm%sw_dn_direct_frac, spartacusprm%air_ext_sw, spartacusprm%air_ssa_sw, &
4049 spartacusprm%veg_ssa_sw, spartacusprm%air_ext_lw, spartacusprm%air_ssa_lw, spartacusprm%veg_ssa_lw, &
4050 spartacusprm%veg_fsd_const, spartacusprm%veg_contact_fraction_const, &
4051 spartacusprm%ground_albedo_dir_mult_fact, methodprm%use_sw_direct_albedo, &
4052 spartacusprm%height, spartacuslayerprm%building_frac, &
4053 spartacuslayerprm%veg_frac, spartacuslayerprm%building_scale, spartacuslayerprm%veg_scale, &
4054 spartacuslayerprm%alb_roof, spartacuslayerprm%emis_roof, spartacuslayerprm%alb_wall, spartacuslayerprm%emis_wall, &
4055 spartacuslayerprm%roof_albedo_dir_mult_fact, spartacuslayerprm%wall_specular_frac, &
4056 phenstate_next%alb, ldown, fcld, &
4057 qn_surf, qn_roof, qn_wall, &
4058 qn, qn_snowfree, qn_snow, kclear, kup, lup, tsurf, &
4059 qn_ind_snow, kup_ind_snow, tsurf_ind_snow, tsurf_ind, &
4060 albedo_snow, snowstate_next%SnowAlb, &
4061 dataoutlinespartacus)
4074 IF (i_iter == 1)
THEN
4075 qg_surf = 0.1*qn_surf
4076 qg_roof = 0.1*qn_roof
4077 qg_wall = 0.1*qn_wall
4079 qg_surf = qn_surf + qf - (qh_surf + qe_surf)
4080 qg_roof = qn_roof + qf - (qh_roof + qe_roof)
4081 qg_wall = qn_wall + qf - (qh_wall + qe_wall)
4116 methodprm%StorageHeatMethod, forcing%qs_obs, methodprm%OHMIncQF, siteinfo%Gridiv, &
4117 timer%id, timer%tstep, timer%dt_since_start, methodprm%Diagnose, &
4119 qg_surf, qg_roof, qg_wall, &
4120 heatstate_out%tsfc_roof, ehcprm%tin_roof, &
4121 heatstate_in%temp_roof, ehcprm%k_roof, &
4122 ehcprm%cp_roof, ehcprm%dz_roof, sfr_roof, &
4123 heatstate_out%tsfc_wall, ehcprm%tin_wall, &
4124 heatstate_in%temp_wall, ehcprm%k_wall, &
4125 ehcprm%cp_wall, ehcprm%dz_wall, sfr_wall, &
4126 heatstate_out%tsfc_surf, ehcprm%tin_surf, &
4127 heatstate_in%temp_surf, ehcprm%k_surf, &
4128 ehcprm%cp_surf, ehcprm%dz_surf, &
4129 pavedprm%sfr, bldgprm%sfr, evetrprm%sfr, dectrprm%sfr, grassprm%sfr, bsoilprm%sfr, waterprm%sfr, &
4130 pavedprm%ohm%ohm_coef_lc, bldgprm%ohm%ohm_coef_lc, evetrprm%ohm%ohm_coef_lc, &
4131 dectrprm%ohm%ohm_coef_lc, grassprm%ohm%ohm_coef_lc, bsoilprm%ohm%ohm_coef_lc, &
4132 waterprm%ohm%ohm_coef_lc, &
4133 pavedprm%ohm%ohm_threshsw, bldgprm%ohm%ohm_threshsw, evetrprm%ohm%ohm_threshsw, &
4134 dectrprm%ohm%ohm_threshsw, grassprm%ohm%ohm_threshsw, bsoilprm%ohm%ohm_threshsw, &
4135 waterprm%ohm%ohm_threshsw, &
4136 pavedprm%ohm%ohm_threshwd, bldgprm%ohm%ohm_threshwd, evetrprm%ohm%ohm_threshwd, &
4137 dectrprm%ohm%ohm_threshwd, grassprm%ohm%ohm_threshwd, bsoilprm%ohm%ohm_threshwd, &
4138 waterprm%ohm%ohm_threshwd, &
4139 hydrostate_prev%soilstore_surf, &
4140 pavedprm%soil%soilstorecap, bldgprm%soil%soilstorecap, &
4141 evetrprm%soil%soilstorecap, dectrprm%soil%soilstorecap, &
4142 grassprm%soil%soilstorecap, bsoilprm%soil%soilstorecap, waterprm%soil%soilstorecap, &
4143 hydrostate_prev%state_surf, methodprm%SnowUse, snowstate_prev%SnowFrac, diagqs, &
4144 anthroheatstate%HDD_id, metforcingdata_grid, ts5mindata_ir, qf, qn, &
4145 forcing%kdown, forcing%U, forcing%temp_c, zenith_deg, forcing%RH, forcing%pres, ldown, &
4146 bldgprm%bldgh, phenstate%alb, &
4147 pavedprm%emis, bldgprm%emis, evetrprm%emis, dectrprm%emis, grassprm%emis, bsoilprm%emis, waterprm%emis, &
4148 pavedprm%ohm%cpanohm, bldgprm%ohm%cpanohm, evetrprm%ohm%cpanohm, &
4149 dectrprm%ohm%cpanohm, grassprm%ohm%cpanohm, bsoilprm%ohm%cpanohm, waterprm%ohm%cpanohm, &
4150 pavedprm%ohm%kkanohm, bldgprm%ohm%kkanohm, evetrprm%ohm%kkanohm, &
4151 dectrprm%ohm%kkanohm, grassprm%ohm%kkanohm, bsoilprm%ohm%kkanohm, waterprm%ohm%kkanohm, &
4152 pavedprm%ohm%chanohm, bldgprm%ohm%chanohm, evetrprm%ohm%chanohm, &
4153 dectrprm%ohm%chanohm, grassprm%ohm%chanohm, bsoilprm%ohm%chanohm, waterprm%ohm%chanohm, &
4154 methodprm%EmissionsMethod, &
4155 forcing%Tair, ohmstate_prev%qn_av, ohmstate_prev%dqndt, ohmstate_prev%qn_s_av, ohmstate_prev%dqnsdt, &
4156 phenstate%StoreDrainPrm, &
4157 qn_snow, dataoutlineestm, qs, &
4158 ohmstate_next%qn_av, ohmstate_next%dqndt, ohmstate_next%qn_s_av, ohmstate_next%dqnsdt, &
4159 deltaqi, a1, a2, a3, &
4160 heatstate_out%temp_roof, qs_roof, &
4161 heatstate_out%temp_wall, qs_wall, &
4162 heatstate_out%temp_surf, qs_surf)
4190 IF (diagnose == 1)
WRITE (*, *)
'Calling MeltHeat'
4193 IF (diagnose == 1)
WRITE (*, *)
'Calling LUMPS_cal_QHQE...'
4194 IF (i_iter == 1)
THEN
4204 lumpsprm%veg_type, &
4205 methodprm%SnowUse, qn, qf, qs, forcing%Temp_C, vegfraction,
avcp, forcing%pres,
lv_j_kg, &
4206 tstep_real, lumpsprm%drainrt, nsh_real, &
4207 forcing%rain, lumpsprm%rainmaxres, lumpsprm%raincover, &
4208 pavedprm%sfr, bldgprm%sfr, evetrprm%sfr, dectrprm%sfr, grassprm%sfr, bsoilprm%sfr, waterprm%sfr, &
4209 phenstate_next%LAI_id, &
4210 evetrprm%lai%laimax, dectrprm%lai%laimax, grassprm%lai%laimax, &
4211 evetrprm%lai%laimin, dectrprm%lai%laimin, grassprm%lai%laimin, &
4213 qe_lumps, psyc_hpa, s_hpa, sice_hpa, tempveg, vegphenlumps)
4223 IF (diagnose == 1)
WRITE (*, *)
'Calling SUEWS_cal_Water...'
4234 methodprm%Diagnose, &
4235 methodprm%SnowUse, nonwaterfraction, addpipes, addimpervious, addveg, addwaterbody, &
4236 hydrostate_prev%state_surf, &
4237 pavedprm%sfr, bldgprm%sfr, evetrprm%sfr, dectrprm%sfr, grassprm%sfr, bsoilprm%sfr, waterprm%sfr, &
4238 phenstate_next%StoreDrainPrm, &
4239 pavedprm%waterdist%to_paved, pavedprm%waterdist%to_bldg, pavedprm%waterdist%to_evetr, &
4240 pavedprm%waterdist%to_dectr, pavedprm%waterdist%to_grass, pavedprm%waterdist%to_bsoil, pavedprm%waterdist%to_water, &
4241 pavedprm%waterdist%to_soilstore, &
4242 bldgprm%waterdist%to_paved, bldgprm%waterdist%to_bldg, bldgprm%waterdist%to_evetr, &
4243 bldgprm%waterdist%to_dectr, bldgprm%waterdist%to_grass, bldgprm%waterdist%to_bsoil, bldgprm%waterdist%to_water, &
4244 bldgprm%waterdist%to_soilstore, &
4245 evetrprm%waterdist%to_paved, evetrprm%waterdist%to_bldg, evetrprm%waterdist%to_evetr, &
4246 evetrprm%waterdist%to_dectr, evetrprm%waterdist%to_grass, evetrprm%waterdist%to_bsoil, evetrprm%waterdist%to_water, &
4247 evetrprm%waterdist%to_soilstore, &
4248 dectrprm%waterdist%to_paved, dectrprm%waterdist%to_bldg, dectrprm%waterdist%to_evetr, &
4249 dectrprm%waterdist%to_dectr, dectrprm%waterdist%to_grass, dectrprm%waterdist%to_bsoil, dectrprm%waterdist%to_water, &
4250 dectrprm%waterdist%to_soilstore, &
4251 grassprm%waterdist%to_paved, grassprm%waterdist%to_bldg, grassprm%waterdist%to_evetr, &
4252 grassprm%waterdist%to_dectr, grassprm%waterdist%to_grass, grassprm%waterdist%to_bsoil, grassprm%waterdist%to_water, &
4253 grassprm%waterdist%to_soilstore, &
4254 bsoilprm%waterdist%to_paved, bsoilprm%waterdist%to_bldg, bsoilprm%waterdist%to_evetr, &
4255 bsoilprm%waterdist%to_dectr, bsoilprm%waterdist%to_grass, bsoilprm%waterdist%to_bsoil, &
4256 bsoilprm%waterdist%to_water, &
4257 bsoilprm%waterdist%to_soilstore, &
4260 drain_surf, frac_water2runoff, &
4261 additionalwater, runoffpipes, runoff_per_interval, &
4266 IF (diagnose == 1)
WRITE (*, *)
'Calling SUEWS_cal_Resistance...'
4282 methodprm%StabilityMethod, &
4283 methodprm%Diagnose, aerodynamicresistancemethod, methodprm%RoughLenHeatMethod, methodprm%SnowUse, &
4284 timer%id, timer%it, &
4285 conductanceprm%gsModel, &
4286 methodprm%SMDMethod, &
4288 forcing%U, forcing%Temp_C, vegfraction, forcing%kdown, &
4289 conductanceprm%Kmax, conductanceprm%g_max, conductanceprm%g_k, conductanceprm%g_q_base, conductanceprm%g_q_shape, &
4290 conductanceprm%g_t, conductanceprm%g_sm, conductanceprm%s1, conductanceprm%s2, &
4291 conductanceprm%th, conductanceprm%tl, &
4292 dq, forcing%xsmd, vsmd, &
4293 evetrprm%maxconductance, dectrprm%maxconductance, grassprm%maxconductance, &
4294 evetrprm%lai%laimax, dectrprm%lai%laimax, grassprm%lai%laimax, &
4295 phenstate_next%LAI_id, snowstate_prev%SnowFrac, &
4296 pavedprm%sfr, bldgprm%sfr, evetrprm%sfr, dectrprm%sfr, grassprm%sfr, bsoilprm%sfr, waterprm%sfr, &
4297 g_kdown, g_dq, g_ta, g_smd, g_lai, &
4298 ustar, tstar, l_mod, &
4299 zl, gsc, rs, ra_h, rasnow, rb, z0v, z0vsnow)
4304 IF (snowuse == 1)
THEN
4335 methodprm%Diagnose, nlayer, &
4336 tstep, timer%imin, timer%it, evapmethod, dayofweek_id, snowprm%CRWmin, snowprm%CRWmax, &
4337 dectime,
avdens,
avcp,
lv_j_kg, lvs_j_kg, forcing%RH, forcing%Pres, forcing%Temp_C, &
4338 rasnow, psyc_hpa, sice_hpa, snowprm%tau_r, &
4339 radmeltfact, tempmeltfact, snowprm%SnowAlbMax, snowprm%PrecipLimit, snowprm%PrecipLimitAlb, &
4340 qn_ind_snow, kup_ind_snow, deltaqi, tsurf_ind_snow, &
4341 snowstate_next%SnowAlb, &
4342 pervfraction, vegfraction, addimpervious, qn_snowfree, qf, qs, vpd_hpa, s_hpa, &
4343 rs, ra_h, rb, snowprm%SnowDensMax, snowprm%snowdensmin, forcing%rain, siteinfo%PipeCapacity, siteinfo%RunoffToWater, &
4344 addveg, snowprm%SnowLimPaved, snowprm%SnowLimBldg, &
4345 siteinfo%FlowChange, drain_surf, &
4346 pavedprm%wetthresh, bldgprm%wetthresh, evetrprm%wetthresh, dectrprm%wetthresh, &
4347 grassprm%wetthresh, bsoilprm%wetthresh, waterprm%wetthresh, &
4348 pavedprm%soil%soilstorecap, bldgprm%soil%soilstorecap, &
4349 evetrprm%soil%soilstorecap, dectrprm%soil%soilstorecap, &
4350 grassprm%soil%soilstorecap, bsoilprm%soil%soilstorecap, waterprm%soil%soilstorecap, &
4352 pavedprm%sfr, bldgprm%sfr, evetrprm%sfr, dectrprm%sfr, grassprm%sfr, bsoilprm%sfr, waterprm%sfr, &
4353 addwater, frac_water2runoff, phenstate_next%StoreDrainPrm, snowprm%SnowPackLimit, &
4354 snowprm%snowprof_24hr_working, snowprm%snowprof_24hr_holiday, &
4355 snowstate_prev%SnowPack, snowstate_prev%snowFrac, snowstate_prev%SnowWater, &
4356 snowstate_prev%IceFrac, snowstate_prev%SnowDens, &
4357 snowstate_prev%SnowfallCum, hydrostate_prev%state_surf, hydrostate_prev%soilstore_surf, &
4360 snowstate_next%SnowPack, snowstate_next%SnowFrac, snowstate_next%SnowWater, &
4361 snowstate_next%iceFrac, snowstate_next%SnowDens, &
4362 snowstate_next%SnowfallCum, hydrostate_next%state_surf, hydrostate_next%soilstore_surf, &
4363 state_per_tstep, nwstate_per_tstep, &
4364 qe, qe_surf, qe_roof, qe_wall, &
4365 snowstate_next%SnowAlb, &
4366 swe, chsnow_per_interval, ev_per_tstep, runoff_per_tstep, &
4367 surf_chang_per_tstep, runoffpipes, mwstore, runoffwaterbody, &
4368 runoffagveg, runoffagimpervious, rss_surf, &
4373 IF (diagnose == 1)
WRITE (*, *)
'Calling SUEWS_cal_QE...'
4405 methodprm%Diagnose, methodprm%storageheatmethod, nlayer, &
4412 qf, vpd_hpa, s_hpa, rs, ra_h, rb, &
4413 forcing%rain, siteinfo%pipecapacity, siteinfo%runofftowater, &
4414 nonwaterfraction, wu_surf, addveg, addwaterbody, addwater, &
4415 siteinfo%flowchange, drain_surf, &
4416 frac_water2runoff, phenstate_next%StoreDrainPrm, &
4417 pavedprm%sfr, bldgprm%sfr, evetrprm%sfr, dectrprm%sfr, grassprm%sfr, bsoilprm%sfr, waterprm%sfr, &
4418 pavedprm%statelimit, bldgprm%statelimit, evetrprm%statelimit, &
4419 dectrprm%statelimit, grassprm%statelimit, bsoilprm%statelimit, waterprm%statelimit, &
4420 pavedprm%soil%soilstorecap, bldgprm%soil%soilstorecap, &
4421 evetrprm%soil%soilstorecap, dectrprm%soil%soilstorecap, &
4422 grassprm%soil%soilstorecap, bsoilprm%soil%soilstorecap, waterprm%soil%soilstorecap, &
4423 pavedprm%wetthresh, bldgprm%wetthresh, evetrprm%wetthresh, &
4424 dectrprm%wetthresh, grassprm%wetthresh, bsoilprm%wetthresh, waterprm%wetthresh, &
4425 hydrostate_prev%state_surf, hydrostate_prev%soilstore_surf, qn_surf, qs_surf, &
4426 sfr_roof, ehcprm%state_limit_roof, ehcprm%soil_storecap_roof, ehcprm%wet_thresh_roof, &
4427 hydrostate_prev%state_roof, hydrostate_prev%soilstore_roof, qn_roof, qs_roof, &
4428 sfr_wall, ehcprm%state_limit_wall, ehcprm%soil_storecap_wall, ehcprm%wet_thresh_wall, &
4429 hydrostate_prev%state_wall, hydrostate_prev%soilstore_wall, qn_wall, qs_wall, &
4430 hydrostate_next%state_surf, hydrostate_next%soilstore_surf, ev_surf, &
4431 hydrostate_next%state_roof, hydrostate_next%soilstore_roof, ev_roof, &
4432 hydrostate_next%state_wall, hydrostate_next%soilstore_wall, ev_wall, &
4433 state_per_tstep, nwstate_per_tstep, &
4434 ev0_surf, qe0_surf, &
4435 qe, qe_surf, qe_roof, qe_wall, &
4436 ev_per_tstep, runoff_per_tstep, &
4437 surf_chang_per_tstep, runoffpipes, &
4439 runoffagveg, runoffagimpervious, rss_surf)
4443 IF (diagnose == 1) print *,
'before SUEWS_cal_SoilState soilstore_id = ', hydrostate_next%soilstore_surf
4447 IF (diagnose == 1)
WRITE (*, *)
'Calling SUEWS_cal_HorizontalSoilWater...'
4461 pavedprm%sfr, bldgprm%sfr, evetrprm%sfr, dectrprm%sfr, grassprm%sfr, bsoilprm%sfr, waterprm%sfr, &
4462 pavedprm%soil%soilstorecap, bldgprm%soil%soilstorecap, evetrprm%soil%soilstorecap, &
4463 dectrprm%soil%soilstorecap, grassprm%soil%soilstorecap, bsoilprm%soil%soilstorecap, &
4464 waterprm%soil%soilstorecap, &
4465 pavedprm%soil%soildepth, bldgprm%soil%soildepth, evetrprm%soil%soildepth, &
4466 dectrprm%soil%soildepth, grassprm%soil%soildepth, bsoilprm%soil%soildepth, &
4467 waterprm%soil%soildepth, &
4468 pavedprm%soil%sathydraulicconduct, bldgprm%soil%sathydraulicconduct, evetrprm%soil%sathydraulicconduct, &
4469 dectrprm%soil%sathydraulicconduct, grassprm%soil%sathydraulicconduct, bsoilprm%soil%sathydraulicconduct, &
4470 waterprm%soil%sathydraulicconduct, &
4471 siteinfo%SurfaceArea, &
4474 hydrostate_next%soilstore_surf, &
4476 runoffsoil_per_tstep &
4480 IF (diagnose == 1)
WRITE (*, *)
'Calling SUEWS_cal_SoilState...'
4487 methodprm%SMDMethod, forcing%xsmd, nonwaterfraction, soilmoistcap, &
4488 pavedprm%soil%soilstorecap, bldgprm%soil%soilstorecap, evetrprm%soil%soilstorecap, &
4489 dectrprm%soil%soilstorecap, grassprm%soil%soilstorecap, bsoilprm%soil%soilstorecap, &
4490 waterprm%soil%soilstorecap, &
4491 surf_chang_per_tstep, &
4492 hydrostate_next%soilstore_surf, hydrostate_prev%soilstore_surf, &
4493 pavedprm%sfr, bldgprm%sfr, evetrprm%sfr, dectrprm%sfr, grassprm%sfr, bsoilprm%sfr, waterprm%sfr, &
4494 smd, smd_nsurf, tot_chang_per_tstep, soilstate)
4497 IF (diagnose == 1)
WRITE (*, *)
'Calling SUEWS_cal_QH...'
4509 1, nlayer, methodprm%StorageHeatMethod, &
4510 qn, qf, qmrain, qe, qs, qmfreez, qm,
avdens,
avcp, &
4511 pavedprm%sfr, bldgprm%sfr, evetrprm%sfr, dectrprm%sfr, grassprm%sfr, bsoilprm%sfr, waterprm%sfr, &
4512 sfr_roof, sfr_wall, &
4513 heatstate_out%tsfc_surf, heatstate_out%tsfc_roof, heatstate_out%tsfc_wall, &
4516 qh, qh_residual, qh_resist, &
4517 qh_resist_surf, qh_resist_roof, qh_resist_wall)
4546 tsfc0_out_surf = min(heatstate_out%tsfc_surf, forcing%Temp_C + 50)
4547 tsfc0_out_roof = min(heatstate_out%tsfc_roof, forcing%Temp_C + 50)
4548 tsfc0_out_wall = min(heatstate_out%tsfc_wall, forcing%Temp_C + 50)
4550 qh_surf = qn_surf + qf - qs_surf - qe_surf
4551 qh_roof = qn_roof + qf - qs_roof - qe_roof
4552 qh_wall = qn_wall + qf - qs_wall - qe_wall
4554 IF (methodprm%diagnose == 1)
THEN
4555 print *,
'qn_surf before QH back env.:', qn_surf
4556 print *,
'qf before QH back env.:', qf
4557 print *,
'qs_surf before QH back env.:', qs_surf
4558 print *,
'qe_surf before QH back env.:', qe_surf
4559 print *,
'qh_surf before QH back env.:', qh_surf
4561 print *,
'qn_roof before QH back env.:', qn_roof
4562 print *,
'qs_roof before QH back env.:', qs_roof
4563 print *,
'qe_roof before QH back env.:', qe_roof
4564 print *,
'qh_roof before QH back env.:', qh_roof
4567 DO i_surf = 1,
nsurf
4569 heatstate_out%tsfc_surf(i_surf) =
cal_tsfc(qh_surf(i_surf),
avdens,
avcp, ra_h, forcing%temp_c)
4579 DO i_surf = 1, nlayer
4580 heatstate_out%tsfc_roof(i_surf) =
cal_tsfc(qh_roof(i_surf),
avdens,
avcp, ra_h, forcing%temp_c)
4581 heatstate_out%tsfc_wall(i_surf) =
cal_tsfc(qh_wall(i_surf),
avdens,
avcp, ra_h, forcing%temp_c)
4584 IF (methodprm%diagnose == 1) print *,
'tsfc_surf after QH back env.:', heatstate_out%tsfc_surf
4586 IF (methodprm%diagnose == 1) print *, &
4587 'tsfc_surf abs. diff.:', maxval(abs(heatstate_out%tsfc_surf - tsfc0_out_surf)), &
4588 maxloc(abs(heatstate_out%tsfc_surf - tsfc0_out_surf))
4589 dif_tsfc_iter = maxval(abs(heatstate_out%tsfc_surf - tsfc0_out_surf))
4590 IF (methodprm%StorageHeatMethod == 5)
THEN
4591 IF (methodprm%diagnose == 1) print *, &
4592 'tsfc_roof abs. diff.:', maxval(abs(heatstate_out%tsfc_roof - tsfc0_out_roof)), &
4593 maxloc(abs(heatstate_out%tsfc_roof - tsfc0_out_roof))
4594 dif_tsfc_iter = max(maxval(abs(heatstate_out%tsfc_roof - tsfc0_out_roof)), dif_tsfc_iter)
4595 IF (methodprm%diagnose == 1) print *, &
4596 'tsfc_wall abs. diff.:', maxval(abs(heatstate_out%tsfc_wall - tsfc0_out_wall)), &
4597 maxloc(abs(heatstate_out%tsfc_wall - tsfc0_out_wall))
4598 dif_tsfc_iter = max(maxval(abs(tsfc0_out_wall - heatstate_out%tsfc_wall)), dif_tsfc_iter)
4605 heatstate_out%tsfc_surf = (tsfc0_out_surf*(1 - ratio_iter) + heatstate_out%tsfc_surf*ratio_iter)
4606 heatstate_out%tsfc_roof = (tsfc0_out_roof*(1 - ratio_iter) + heatstate_out%tsfc_roof*ratio_iter)
4607 heatstate_out%tsfc_wall = (tsfc0_out_wall*(1 - ratio_iter) + heatstate_out%tsfc_wall*ratio_iter)
4628 IF (dif_tsfc_iter > .1)
THEN
4629 flag_converge = .false.
4631 flag_converge = .true.
4641 IF (methodprm%Diagnose == 1 .AND. i_iter == max_iter)
THEN
4661 IF (methodprm%Diagnose == 1)
WRITE (*, *)
'Calling RSLProfile...'
4673 methodprm%DiagMethod, &
4674 zh, z0m, zdm, z0v, &
4676 pavedprm%sfr, bldgprm%sfr, evetrprm%sfr, dectrprm%sfr, grassprm%sfr, bsoilprm%sfr, waterprm%sfr, &
4678 methodprm%StabilityMethod, ra_h, &
4680 forcing%U, forcing%Temp_C, forcing%RH, forcing%pres, siteinfo%z, qh, qe, &
4681 t2_c, q2_gkg, u10_ms, rh2, &
4694 evetrprm%bioco2%alpha_bioco2, dectrprm%bioco2%alpha_bioco2, grassprm%bioco2%alpha_bioco2, &
4695 evetrprm%bioco2%alpha_enh_bioco2, dectrprm%bioco2%alpha_enh_bioco2, grassprm%bioco2%alpha_enh_bioco2, &
4696 forcing%kdown, forcing%RH, &
4697 evetrprm%bioco2%beta_bioCO2, dectrprm%bioco2%beta_bioCO2, grassprm%bioco2%beta_bioCO2, &
4698 evetrprm%bioco2%beta_enh_bioco2, dectrprm%bioco2%beta_enh_bioco2, grassprm%bioco2%beta_enh_bioco2, &
4699 dectime, methodprm%Diagnose, methodprm%EmissionsMethod, fc_anthro, &
4700 conductanceprm%g_max, conductanceprm%g_k, conductanceprm%g_q_base, conductanceprm%g_q_shape, &
4701 conductanceprm%g_t, conductanceprm%g_sm, gfunc, conductanceprm%gsmodel, &
4702 timer%id, timer%it, conductanceprm%Kmax, &
4703 phenstate_next%LAI_id, &
4704 evetrprm%lai%laimin, dectrprm%lai%laimin, grassprm%lai%laimin, &
4705 evetrprm%lai%laimax, dectrprm%lai%laimax, grassprm%lai%laimax, &
4706 evetrprm%MaxConductance, dectrprm%MaxConductance, grassprm%MaxConductance, &
4707 evetrprm%bioco2%min_res_bioCO2, dectrprm%bioco2%min_res_bioCO2, grassprm%bioco2%min_res_bioCO2, &
4709 evetrprm%bioco2%resp_a, dectrprm%bioco2%resp_a, grassprm%bioco2%resp_a, &
4710 evetrprm%bioco2%resp_b, dectrprm%bioco2%resp_b, grassprm%bioco2%resp_b, &
4711 conductanceprm%S1, conductanceprm%S2, &
4712 pavedprm%sfr, bldgprm%sfr, evetrprm%sfr, dectrprm%sfr, grassprm%sfr, bsoilprm%sfr, waterprm%sfr, &
4713 methodprm%SMDMethod, snowstate%SnowFrac, &
4714 t2_c, forcing%Temp_C, &
4715 evetrprm%bioco2%theta_bioCO2, dectrprm%bioco2%theta_bioCO2, grassprm%bioco2%theta_bioco2, &
4716 conductanceprm%TH, conductanceprm%TL, vsmd, forcing%xsmd, &
4717 fc, fc_biogen, fc_photo, fc_respi)
4724 ohmstate%qn_av = ohmstate_next%qn_av
4725 ohmstate%dqndt = ohmstate_next%dqndt
4726 ohmstate%qn_s_av = ohmstate_next%qn_s_av
4727 ohmstate%dqnsdt = ohmstate_next%dqnsdt
4728 snowstate%SnowfallCum = snowstate_next%SnowfallCum
4729 snowstate%SnowAlb = snowstate_next%SnowAlb
4730 snowstate%IceFrac = snowstate_next%IceFrac
4731 snowstate%SnowWater = snowstate_next%SnowWater
4732 snowstate%SnowDens = snowstate_next%SnowDens
4733 snowstate%SnowFrac = snowstate_next%SnowFrac
4734 snowstate%SnowPack = snowstate_next%SnowPack
4736 hydrostate%soilstore_surf = hydrostate_next%soilstore_surf
4737 hydrostate%state_surf = hydrostate_next%state_surf
4738 phenstate%alb = phenstate_next%alb
4739 phenstate%GDD_id = phenstate_next%GDD_id
4740 phenstate%SDD_id = phenstate_next%SDD_id
4741 phenstate%LAI_id = phenstate_next%LAI_id
4742 phenstate%DecidCap_id = phenstate_next%DecidCap_id
4743 phenstate%albDecTr_id = phenstate_next%albDecTr_id
4744 phenstate%albEveTr_id = phenstate_next%albEveTr_id
4745 phenstate%albGrass_id = phenstate_next%albGrass_id
4746 phenstate%porosity_id = phenstate_next%porosity_id
4747 phenstate%StoreDrainPrm = phenstate_next%StoreDrainPrm
4749 forcing%Tair = tair_av_next
4750 phenstate%Tmin_id = phenstate_next%Tmin_id
4751 phenstate%Tmax_id = phenstate_next%Tmax_id
4752 phenstate%lenday_id = phenstate_next%lenday_id
4753 anthroheatstate%HDD_id = anthroheatstate_next%HDD_id
4754 hydrostate%WUDay_id = hydrostate_next%WUDay_id
4756 IF (methodprm%StorageHeatMethod == 5)
THEN
4758 heatstate%temp_roof = heatstate_out%temp_roof
4759 heatstate%temp_wall = heatstate_out%temp_wall
4760 heatstate%temp_surf = heatstate_out%temp_surf
4761 heatstate%tsfc_roof = heatstate_out%tsfc_roof
4762 heatstate%tsfc_wall = heatstate_out%tsfc_wall
4763 heatstate%tsfc_surf = heatstate_out%tsfc_surf
4765 hydrostate%soilstore_roof = hydrostate_next%soilstore_roof
4766 hydrostate%state_roof = hydrostate_next%state_roof
4767 hydrostate%soilstore_wall = hydrostate_next%soilstore_wall
4768 hydrostate%state_wall = hydrostate_next%state_wall
4773 tair_av = tair_av_next
4775 qn_av = ohmstate_next%qn_av
4776 dqndt = ohmstate_next%dqndt
4777 qn_s_av = ohmstate_next%qn_s_av
4778 dqnsdt = ohmstate_next%dqnsdt
4780 snowfallcum = snowstate_next%SnowfallCum
4781 snowalb = snowstate_next%SnowAlb
4782 icefrac = snowstate_next%IceFrac
4783 snowwater = snowstate_next%SnowWater
4784 snowdens = snowstate_next%SnowDens
4785 snowfrac = snowstate_next%SnowFrac
4786 snowpack = snowstate_next%SnowPack
4788 soilstore_surf = hydrostate_next%soilstore_surf
4789 state_surf = hydrostate_next%state_surf
4790 wuday_id = hydrostate_next%WUDay_id
4792 alb = phenstate_next%alb
4793 gdd_id = phenstate_next%GDD_id
4794 sdd_id = phenstate_next%SDD_id
4795 lai_id = phenstate_next%LAI_id
4796 tmin_id = phenstate_next%Tmin_id
4797 tmax_id = phenstate_next%Tmax_id
4798 lenday_id = phenstate_next%lenday_id
4799 decidcap_id = phenstate_next%DecidCap_id
4800 albdectr_id = phenstate_next%albDecTr_id
4801 albevetr_id = phenstate_next%albEveTr_id
4802 albgrass_id = phenstate_next%albGrass_id
4803 porosity_id = phenstate_next%porosity_id
4804 storedrainprm = phenstate_next%StoreDrainPrm
4807 hdd_id = anthroheatstate_next%HDD_id
4809 IF (methodprm%StorageHeatMethod == 5)
THEN
4811 temp_roof = heatstate_out%temp_roof
4812 temp_wall = heatstate_out%temp_wall
4813 temp_surf = heatstate_out%temp_surf
4814 tsfc_roof = heatstate_out%tsfc_roof
4815 tsfc_wall = heatstate_out%tsfc_wall
4816 tsfc_surf = heatstate_out%tsfc_surf
4818 soilstore_roof = hydrostate_next%soilstore_roof
4819 state_roof = hydrostate_next%state_roof
4820 soilstore_wall = hydrostate_next%soilstore_wall
4821 state_wall = hydrostate_next%state_wall
4835 pai = sfr_surf(2)/sum(sfr_surf(1:2))
4840 CALL beers_cal_main(timer%iy, timer%id, dectime, pai, fai, forcing%kdown, ldown, forcing%Temp_C, forcing%RH, &
4841 forcing%Pres, tsfc_c, siteinfo%lat, siteinfo%lon, siteinfo%alt, siteinfo%timezone, zenith_deg, azimuth, &
4842 phenstate%alb(1), phenstate%alb(2), pavedprm%emis, bldgprm%emis, &
4847 dataoutlinebeers =
set_nan(dataoutlinebeers)
4871 additionalwater, phenstate%alb, forcing%kdown, u10_ms, azimuth, &
4872 chsnow_per_interval, dectime, &
4873 drain_per_tstep, qe_lumps, ev_per_tstep, wu_ext, fc, fc_build, fcld, &
4874 fc_metab, fc_photo, fc_respi, fc_point, fc_traff, siteinfo%FlowChange, &
4875 qh_lumps, timer%id, timer%imin, wu_int, timer%it, timer%iy, &
4876 kup, phenstate%LAI_id, ldown, l_mod, lup, mwh, &
4878 nsh_real, nwstate_per_tstep, precip, q2_gkg, &
4879 qe, qf, qh, qh_resist, qm, qmfreez, &
4880 qmrain, qn, qn_snow, qn_snowfree, qs, ra_h, &
4881 rs, rh2, runoffagimpervious, runoffagveg, &
4882 runoff_per_tstep, runoffpipes, runoffsoil_per_tstep, &
4885 smd, smd_nsurf, snowstate%SnowAlb, snowremoval, &
4886 hydrostate%state_surf, state_per_tstep, surf_chang_per_tstep, swe, t2_c, tsfc_c, &
4887 tot_chang_per_tstep, tsurf, ustar, &
4889 z0m, zdm, zenith_deg, &
4890 datetimeline, dataoutlinesuews)
4911 timer%iy, timer%id, timer%it, timer%imin, dectime, nlayer, &
4912 heatstate_out%tsfc_surf, qs_surf, &
4913 heatstate_out%tsfc_roof, &
4918 hydrostate%state_roof, &
4919 hydrostate%soilstore_roof, &
4920 heatstate_out%tsfc_wall, &
4925 hydrostate%state_wall, &
4926 hydrostate%soilstore_wall, &
4927 datetimeline, dataoutlineehc)
4946 timer%it, timer%imin, nsh_real, &
4947 phenstate%GDD_id, anthroheatstate%HDD_id, phenstate%LAI_id, &
4949 phenstate%Tmin_id, phenstate%Tmax_id, phenstate%lenday_id, &
4950 phenstate%DecidCap_id, &
4951 phenstate%albDecTr_id, &
4952 phenstate%albEveTr_id, &
4953 phenstate%albGrass_id, &
4954 phenstate%porosity_id, &
4955 hydrostate%WUDay_id, &
4957 snowstate%SnowAlb, snowstate%SnowDens, &
4959 dataoutlinedailystate)
4963 dataoutlinedebug = &
4965 qn_surf, qs_surf, qe0_surf, qe_surf, qh_surf, &
4966 wu_surf, ev0_surf, ev_surf, drain_surf, &
4967 hydrostate_prev%state_surf, hydrostate_next%state_surf, &
4968 hydrostate_prev%soilstore_surf, hydrostate_next%soilstore_surf, &
4969 rs, ra_h, rb, rasnow, rss_surf, &
4970 vsmd, s1/g_sm + s2, g_sm, g_sm*(vsmd - s1/g_sm + s2), &
4971 g_kdown, g_dq, g_ta, g_smd, g_lai, &
4974 faibldg_use, faievetree_use, faidectree_use, fai, &
4979 output_line_suews%datetimeLine = datetimeline
4980 output_line_suews%dataOutLineSUEWS = [datetimeline, dataoutlinesuews]
4981 output_line_suews%dataOutLineEHC = [datetimeline, dataoutlineehc]
4982 output_line_suews%dataOutLineDailyState = [datetimeline, dataoutlinedailystate]
4983 output_line_suews%dataOutLineBEERS = [datetimeline, dataoutlinebeers]
4984 output_line_suews%dataOutLineDebug = [datetimeline, dataoutlinedebug]
4985 output_line_suews%dataOutLineSPARTACUS = [datetimeline, dataoutlinespartacus]
4986 output_line_suews%dataOutLineSnow = [datetimeline, dataoutlinesnow]
4987 output_line_suews%dataoutLineRSL = [datetimeline, dataoutlinersl]
4988 output_line_suews%dataOutLineESTM = [datetimeline, dataoutlineestm]
10091 MetForcingBlock, len_sim, &
10092 AH_MIN, AHProf_24hr, AH_SLOPE_Cooling, & ! input&inout in alphabetical order
10093 AH_SLOPE_Heating, &
10094 alb, AlbMax_DecTr, AlbMax_EveTr, AlbMax_Grass, &
10095 AlbMin_DecTr, AlbMin_EveTr, AlbMin_Grass, &
10096 alpha_bioCO2, alpha_enh_bioCO2, alt, BaseT, BaseTe, &
10097 beta_bioCO2, beta_enh_bioCO2, bldgH, CapMax_dec, CapMin_dec, &
10098 chAnOHM, CO2PointSource, cpAnOHM, CRWmax, CRWmin, DayWat, DayWatPer, &
10099 DecTreeH, DiagMethod, Diagnose, DRAINRT, &
10100 dt_since_start, dqndt, qn_av, dqnsdt, qn_s_av, &
10101 EF_umolCO2perJ, emis, EmissionsMethod, EnEF_v_Jkm, endDLS, EveTreeH, FAIBldg, &
10102 FAIDecTree, FAIEveTree, FAIMethod, Faut, FcEF_v_kgkm, FlowChange, &
10103 FrFossilFuel_Heat, FrFossilFuel_NonHeat, G_max, G_k, G_q_base, G_q_shape, G_t, G_sm, GDD_id, &
10104 GDDFull, Gridiv, gsModel, H_maintain, HDD_id, HumActivity_24hr, &
10105 IceFrac, Ie_a, Ie_end, Ie_m, Ie_start, &
10106 InternalWaterUse_h, &
10107 IrrFracPaved, IrrFracBldgs, &
10108 IrrFracEveTr, IrrFracDecTr, IrrFracGrass, &
10109 IrrFracBSoil, IrrFracWater, &
10110 kkAnOHM, Kmax, LAI_id, LAIMax, LAIMin, &
10111 LAIPower, LAIType, lat, lng, MaxConductance, MaxFCMetab, MaxQFMetab, &
10112 SnowWater, MinFCMetab, MinQFMetab, min_res_bioCO2, &
10113 NARP_EMIS_SNOW, NARP_TRANS_SITE, NetRadiationMethod, &
10114 OHM_coef, OHMIncQF, OHM_threshSW, &
10115 OHM_threshWD, PipeCapacity, PopDensDaytime, &
10116 PopDensNighttime, PopProf_24hr, PorMax_dec, PorMin_dec, &
10117 PrecipLimit, PrecipLimitAlb, &
10118 QF0_BEU, Qf_A, Qf_B, Qf_C, &
10120 n_vegetation_region_urban, &
10121 n_stream_sw_urban, n_stream_lw_urban, &
10122 sw_dn_direct_frac, air_ext_sw, air_ssa_sw, &
10123 veg_ssa_sw, air_ext_lw, air_ssa_lw, veg_ssa_lw, &
10124 veg_fsd_const, veg_contact_fraction_const, &
10125 ground_albedo_dir_mult_fact, use_sw_direct_albedo, & !input
10126 height, building_frac, veg_frac, building_scale, veg_scale, & !input: SPARTACUS
10127 alb_roof, emis_roof, alb_wall, emis_wall, &
10128 roof_albedo_dir_mult_fact, wall_specular_frac, &
10129 RadMeltFact, RAINCOVER, RainMaxRes, resp_a, resp_b, &
10130 RoughLenHeatMethod, RoughLenMomMethod, RunoffToWater, S1, S2, &
10131 SatHydraulicConduct, SDDFull, SDD_id, SMDMethod, SnowAlb, SnowAlbMax, &
10132 SnowAlbMin, SnowPackLimit, SnowDens, SnowDensMax, SnowDensMin, SnowfallCum, SnowFrac, &
10133 SnowLimBldg, SnowLimPaved, SnowPack, SnowProf_24hr, SnowUse, SoilDepth, &
10134 StabilityMethod, startDLS, &
10135 soilstore_surf, SoilStoreCap_surf, state_surf, StateLimit_surf, WetThresh_surf, &
10136 soilstore_roof, SoilStoreCap_roof, state_roof, StateLimit_roof, WetThresh_roof, &
10137 soilstore_wall, SoilStoreCap_wall, state_wall, StateLimit_wall, WetThresh_wall, &
10138 StorageHeatMethod, StoreDrainPrm, SurfaceArea, Tair_av, tau_a, tau_f, tau_r, &
10139 BaseT_Cooling, BaseT_Heating, TempMeltFact, TH, &
10140 theta_bioCO2, timezone, TL, TrafficRate, TrafficUnits, &
10142 tsfc_roof, tsfc_wall, tsfc_surf, &
10143 temp_roof, temp_wall, temp_surf, &
10144 tin_roof, tin_wall, tin_surf, &
10145 k_wall, k_roof, k_surf, &
10146 cp_wall, cp_roof, cp_surf, &
10147 dz_wall, dz_roof, dz_surf, &
10148 Tmin_id, Tmax_id, lenday_id, &
10149 TraffProf_24hr, Ts5mindata_ir, tstep, tstep_prev, veg_type, &
10150 WaterDist, WaterUseMethod, &
10151 WUDay_id, DecidCap_id, albDecTr_id, albEveTr_id, albGrass_id, porosity_id, &
10152 WUProfA_24hr, WUProfM_24hr, Z, z0m_in, zdm_in, &
10153 output_block_suews)
10161 REAL(KIND(1D0)),
DIMENSION(len_sim, 24),
INTENT(IN) :: MetForcingBlock
10162 INTEGER,
INTENT(IN) :: len_sim
10164 INTEGER,
INTENT(IN) :: nlayer
10167 INTEGER,
INTENT(INOUT) :: Diagnose
10170 INTEGER,
INTENT(IN) :: startDLS
10171 INTEGER,
INTENT(IN) :: endDLS
10172 INTEGER,
INTENT(IN) :: EmissionsMethod
10173 INTEGER,
INTENT(IN) :: Gridiv
10174 INTEGER,
INTENT(IN) :: gsModel
10175 INTEGER,
INTENT(IN) :: Ie_end
10176 INTEGER,
INTENT(IN) :: Ie_start
10178 INTEGER,
INTENT(in) :: DiagMethod
10179 INTEGER,
INTENT(IN) :: NetRadiationMethod
10180 INTEGER,
INTENT(IN) :: OHMIncQF
10181 INTEGER,
INTENT(IN) :: RoughLenHeatMethod
10182 INTEGER,
INTENT(IN) :: RoughLenMomMethod
10183 INTEGER,
INTENT(IN) :: FAIMethod
10184 INTEGER,
INTENT(IN) :: SMDMethod
10185 INTEGER,
INTENT(IN) :: SnowUse
10186 INTEGER,
INTENT(IN) :: StabilityMethod
10187 INTEGER,
INTENT(IN) :: StorageHeatMethod
10188 INTEGER,
INTENT(IN) :: tstep
10189 INTEGER,
INTENT(IN) :: tstep_prev
10191 INTEGER,
INTENT(inout) :: dt_since_start
10192 INTEGER,
INTENT(IN) :: veg_type
10193 INTEGER,
INTENT(IN) :: WaterUseMethod
10195 INTEGER,
DIMENSION(NVEGSURF),
INTENT(IN) :: LAIType
10197 REAL(KIND(1D0)),
INTENT(IN) :: AlbMax_DecTr
10198 REAL(KIND(1D0)),
INTENT(IN) :: AlbMax_EveTr
10199 REAL(KIND(1D0)),
INTENT(IN) :: AlbMax_Grass
10200 REAL(KIND(1D0)),
INTENT(IN) :: AlbMin_DecTr
10201 REAL(KIND(1D0)),
INTENT(IN) :: AlbMin_EveTr
10202 REAL(KIND(1D0)),
INTENT(IN) :: AlbMin_Grass
10203 REAL(KIND(1D0)),
INTENT(IN) :: alt
10207 REAL(KIND(1D0)),
INTENT(IN) :: bldgH
10208 REAL(KIND(1D0)),
INTENT(IN) :: CapMax_dec
10209 REAL(KIND(1D0)),
INTENT(IN) :: CapMin_dec
10210 REAL(KIND(1D0)),
INTENT(IN) :: CO2PointSource
10211 REAL(KIND(1D0)),
INTENT(IN) :: CRWmax
10212 REAL(KIND(1D0)),
INTENT(IN) :: CRWmin
10213 REAL(KIND(1D0)),
INTENT(IN) :: DecTreeH
10214 REAL(KIND(1D0)),
INTENT(IN) :: DRAINRT
10215 REAL(KIND(1D0)),
INTENT(IN) :: EF_umolCO2perJ
10216 REAL(KIND(1D0)),
INTENT(IN) :: EnEF_v_Jkm
10217 REAL(KIND(1D0)),
INTENT(IN) :: EveTreeH
10218 REAL(KIND(1D0)),
INTENT(IN) :: FAIBldg
10219 REAL(KIND(1D0)),
INTENT(IN) :: FAIDecTree
10220 REAL(KIND(1D0)),
INTENT(IN) :: FAIEveTree
10221 REAL(KIND(1D0)),
INTENT(IN) :: Faut
10223 REAL(KIND(1D0)),
INTENT(IN) :: FlowChange
10224 REAL(KIND(1D0)),
INTENT(IN) :: FrFossilFuel_Heat
10225 REAL(KIND(1D0)),
INTENT(IN) :: FrFossilFuel_NonHeat
10226 REAL(KIND(1D0)),
INTENT(IN) :: G_max
10227 REAL(KIND(1D0)),
INTENT(IN) :: G_k
10228 REAL(KIND(1D0)),
INTENT(IN) :: G_q_base
10229 REAL(KIND(1D0)),
INTENT(IN) :: G_q_shape
10230 REAL(KIND(1D0)),
INTENT(IN) :: G_t
10231 REAL(KIND(1D0)),
INTENT(IN) :: G_sm
10232 REAL(KIND(1D0)),
INTENT(IN) :: H_maintain
10233 REAL(KIND(1D0)),
INTENT(IN) :: InternalWaterUse_h
10234 REAL(KIND(1D0)),
INTENT(IN) :: IrrFracPaved
10235 REAL(KIND(1D0)),
INTENT(IN) :: IrrFracBldgs
10236 REAL(KIND(1D0)),
INTENT(IN) :: IrrFracEveTr
10237 REAL(KIND(1D0)),
INTENT(IN) :: IrrFracDecTr
10238 REAL(KIND(1D0)),
INTENT(IN) :: IrrFracGrass
10239 REAL(KIND(1D0)),
INTENT(IN) :: IrrFracBSoil
10240 REAL(KIND(1D0)),
INTENT(IN) :: IrrFracWater
10241 REAL(KIND(1D0)),
INTENT(IN) :: Kmax
10243 REAL(KIND(1D0)),
INTENT(IN) :: lat
10245 REAL(KIND(1D0)),
INTENT(IN) :: lng
10246 REAL(KIND(1D0)),
INTENT(IN) :: MaxFCMetab
10247 REAL(KIND(1D0)),
INTENT(IN) :: MaxQFMetab
10248 REAL(KIND(1D0)),
INTENT(IN) :: MinFCMetab
10249 REAL(KIND(1D0)),
INTENT(IN) :: MinQFMetab
10250 REAL(KIND(1D0)),
INTENT(IN) :: NARP_EMIS_SNOW
10251 REAL(KIND(1D0)),
INTENT(IN) :: NARP_TRANS_SITE
10252 REAL(KIND(1D0)),
INTENT(IN) :: PipeCapacity
10253 REAL(KIND(1D0)),
INTENT(IN) :: PopDensNighttime
10254 REAL(KIND(1D0)),
INTENT(IN) :: PorMax_dec
10255 REAL(KIND(1D0)),
INTENT(IN) :: PorMin_dec
10257 REAL(KIND(1D0)),
INTENT(IN) :: PrecipLimit
10258 REAL(KIND(1D0)),
INTENT(IN) :: PrecipLimitAlb
10264 REAL(KIND(1D0)),
INTENT(IN) :: RadMeltFact
10265 REAL(KIND(1D0)),
INTENT(IN) :: RAINCOVER
10266 REAL(KIND(1D0)),
INTENT(IN) :: RainMaxRes
10267 REAL(KIND(1D0)),
INTENT(IN) :: RunoffToWater
10268 REAL(KIND(1D0)),
INTENT(IN) :: S1
10269 REAL(KIND(1D0)),
INTENT(IN) :: S2
10270 REAL(KIND(1D0)),
INTENT(IN) :: SnowAlbMax
10271 REAL(KIND(1D0)),
INTENT(IN) :: SnowAlbMin
10272 REAL(KIND(1D0)),
INTENT(IN) :: SnowDensMax
10273 REAL(KIND(1D0)),
INTENT(IN) :: SnowDensMin
10274 REAL(KIND(1D0)),
INTENT(IN) :: SnowLimBldg
10275 REAL(KIND(1D0)),
INTENT(IN) :: SnowLimPaved
10277 REAL(KIND(1D0)),
INTENT(IN) :: SurfaceArea
10278 REAL(KIND(1D0)),
INTENT(IN) :: tau_a
10279 REAL(KIND(1D0)),
INTENT(IN) :: tau_f
10280 REAL(KIND(1D0)),
INTENT(IN) :: tau_r
10282 REAL(KIND(1D0)),
INTENT(IN) :: TempMeltFact
10283 REAL(KIND(1D0)),
INTENT(IN) :: TH
10284 REAL(KIND(1D0)),
INTENT(IN) :: timezone
10285 REAL(KIND(1D0)),
INTENT(IN) :: TL
10286 REAL(KIND(1D0)),
INTENT(IN) :: TrafficUnits
10288 REAL(KIND(1D0)),
INTENT(IN) :: Z
10289 REAL(KIND(1D0)),
INTENT(IN) :: z0m_in
10290 REAL(KIND(1D0)),
INTENT(IN) :: zdm_in
10292 REAL(KIND(1D0)),
DIMENSION(2),
INTENT(IN) :: AH_MIN
10293 REAL(KIND(1D0)),
DIMENSION(2),
INTENT(IN) :: AH_SLOPE_Cooling
10294 REAL(KIND(1D0)),
DIMENSION(2),
INTENT(IN) :: AH_SLOPE_Heating
10295 REAL(KIND(1D0)),
DIMENSION(2),
INTENT(IN) :: FcEF_v_kgkm
10296 REAL(KIND(1D0)),
DIMENSION(2),
INTENT(IN) :: QF0_BEU
10297 REAL(KIND(1D0)),
DIMENSION(2),
INTENT(IN) :: Qf_A
10298 REAL(KIND(1D0)),
DIMENSION(2),
INTENT(IN) :: Qf_B
10299 REAL(KIND(1D0)),
DIMENSION(2),
INTENT(IN) :: Qf_C
10301 REAL(KIND(1D0)),
DIMENSION(2),
INTENT(IN) :: PopDensDaytime
10302 REAL(KIND(1D0)),
DIMENSION(2),
INTENT(IN) :: BaseT_Cooling
10303 REAL(KIND(1D0)),
DIMENSION(2),
INTENT(IN) :: BaseT_Heating
10304 REAL(KIND(1D0)),
DIMENSION(2),
INTENT(IN) :: TrafficRate
10305 REAL(KIND(1D0)),
DIMENSION(3),
INTENT(IN) :: Ie_a
10306 REAL(KIND(1D0)),
DIMENSION(3),
INTENT(IN) :: Ie_m
10307 REAL(KIND(1D0)),
DIMENSION(3),
INTENT(IN) :: MaxConductance
10308 REAL(KIND(1D0)),
DIMENSION(7),
INTENT(IN) :: DayWat
10309 REAL(KIND(1D0)),
DIMENSION(7),
INTENT(IN) :: DayWatPer
10310 REAL(KIND(1D0)),
DIMENSION(nsurf + 1),
INTENT(IN) :: OHM_threshSW
10311 REAL(KIND(1D0)),
DIMENSION(nsurf + 1),
INTENT(IN) :: OHM_threshWD
10312 REAL(KIND(1D0)),
DIMENSION(NSURF),
INTENT(IN) :: chAnOHM
10313 REAL(KIND(1D0)),
DIMENSION(NSURF),
INTENT(IN) :: cpAnOHM
10314 REAL(KIND(1D0)),
DIMENSION(NSURF),
INTENT(IN) :: emis
10315 REAL(KIND(1D0)),
DIMENSION(NSURF),
INTENT(IN) :: kkAnOHM
10316 REAL(KIND(1D0)),
DIMENSION(NSURF),
INTENT(IN) :: SatHydraulicConduct
10317 REAL(KIND(1D0)),
DIMENSION(NSURF),
INTENT(IN) :: sfr_surf
10318 REAL(KIND(1D0)),
DIMENSION(NSURF),
INTENT(IN) :: SnowPackLimit
10319 REAL(KIND(1D0)),
DIMENSION(NSURF),
INTENT(IN) :: SoilDepth
10320 REAL(KIND(1D0)),
DIMENSION(NSURF),
INTENT(IN) :: SoilStoreCap_surf
10321 REAL(KIND(1D0)),
DIMENSION(NSURF),
INTENT(IN) :: StateLimit_surf
10322 REAL(KIND(1D0)),
DIMENSION(NSURF),
INTENT(IN) :: WetThresh_surf
10323 REAL(KIND(1D0)),
DIMENSION(NVEGSURF),
INTENT(IN) :: alpha_bioCO2
10324 REAL(KIND(1D0)),
DIMENSION(NVEGSURF),
INTENT(IN) :: alpha_enh_bioCO2
10325 REAL(KIND(1D0)),
DIMENSION(NVEGSURF),
INTENT(IN) :: BaseT
10326 REAL(KIND(1D0)),
DIMENSION(NVEGSURF),
INTENT(IN) :: BaseTe
10327 REAL(KIND(1D0)),
DIMENSION(NVEGSURF),
INTENT(IN) :: beta_bioCO2
10328 REAL(KIND(1D0)),
DIMENSION(NVEGSURF),
INTENT(IN) :: beta_enh_bioCO2
10329 REAL(KIND(1D0)),
DIMENSION(NVEGSURF),
INTENT(IN) :: GDDFull
10330 REAL(KIND(1D0)),
DIMENSION(NVEGSURF),
INTENT(IN) :: LAIMax
10331 REAL(KIND(1D0)),
DIMENSION(NVEGSURF),
INTENT(IN) :: LAIMin
10332 REAL(KIND(1D0)),
DIMENSION(NVEGSURF),
INTENT(IN) :: min_res_bioCO2
10333 REAL(KIND(1D0)),
DIMENSION(NVEGSURF),
INTENT(IN) :: resp_a
10334 REAL(KIND(1D0)),
DIMENSION(NVEGSURF),
INTENT(IN) :: resp_b
10335 REAL(KIND(1D0)),
DIMENSION(NVEGSURF),
INTENT(IN) :: SDDFull
10336 REAL(KIND(1D0)),
DIMENSION(0:23, 2),
INTENT(IN) :: SnowProf_24hr
10337 REAL(KIND(1D0)),
DIMENSION(NVEGSURF),
INTENT(IN) :: theta_bioCO2
10338 REAL(KIND(1D0)),
DIMENSION(4, NVEGSURF),
INTENT(IN) :: LAIPower
10339 REAL(KIND(1D0)),
DIMENSION(nsurf + 1, 4, 3),
INTENT(IN) :: OHM_coef
10340 REAL(KIND(1D0)),
DIMENSION(NSURF + 1, NSURF - 1),
INTENT(IN) :: WaterDist
10341 REAL(KIND(1D0)),
DIMENSION(:),
INTENT(IN) :: Ts5mindata_ir
10344 REAL(KIND(1D0)),
DIMENSION(0:23, 2),
INTENT(IN) :: AHProf_24hr
10345 REAL(KIND(1D0)),
DIMENSION(0:23, 2),
INTENT(IN) :: HumActivity_24hr
10346 REAL(KIND(1D0)),
DIMENSION(0:23, 2),
INTENT(IN) :: PopProf_24hr
10347 REAL(KIND(1D0)),
DIMENSION(0:23, 2),
INTENT(IN) :: TraffProf_24hr
10348 REAL(KIND(1D0)),
DIMENSION(0:23, 2),
INTENT(IN) :: WUProfA_24hr
10349 REAL(KIND(1D0)),
DIMENSION(0:23, 2),
INTENT(IN) :: WUProfM_24hr
10355 REAL(KIND(1D0)),
INTENT(INOUT) :: qn_av
10356 REAL(KIND(1D0)),
INTENT(INOUT) :: dqndt
10357 REAL(KIND(1D0)),
INTENT(INOUT) :: qn_s_av
10358 REAL(KIND(1D0)),
INTENT(INOUT) :: dqnsdt
10361 REAL(KIND(1D0)),
INTENT(INOUT) :: SnowfallCum
10362 REAL(KIND(1D0)),
INTENT(INOUT) :: SnowAlb
10363 REAL(KIND(1D0)),
DIMENSION(NSURF),
INTENT(INOUT) :: IceFrac
10364 REAL(KIND(1D0)),
DIMENSION(NSURF),
INTENT(INOUT) :: SnowWater
10365 REAL(KIND(1D0)),
DIMENSION(NSURF),
INTENT(INOUT) :: SnowDens
10366 REAL(KIND(1D0)),
DIMENSION(NSURF),
INTENT(INOUT) :: SnowFrac
10367 REAL(KIND(1D0)),
DIMENSION(NSURF),
INTENT(INOUT) :: SnowPack
10370 REAL(KIND(1D0)),
DIMENSION(NSURF),
INTENT(INOUT) :: soilstore_surf
10371 REAL(KIND(1D0)),
DIMENSION(NSURF),
INTENT(INOUT) :: state_surf
10372 REAL(KIND(1D0)),
DIMENSION(6, NSURF),
INTENT(INOUT) :: StoreDrainPrm
10375 REAL(KIND(1D0)),
DIMENSION(NSURF),
INTENT(INOUT) :: alb
10376 REAL(KIND(1D0)),
DIMENSION(nvegsurf),
INTENT(INOUT) :: GDD_id
10377 REAL(KIND(1D0)),
DIMENSION(nvegsurf),
INTENT(INOUT) :: SDD_id
10378 REAL(KIND(1D0)),
DIMENSION(nvegsurf),
INTENT(INOUT) :: LAI_id
10379 REAL(KIND(1D0)),
INTENT(INOUT) :: DecidCap_id
10380 REAL(KIND(1D0)),
INTENT(INOUT) :: albDecTr_id
10381 REAL(KIND(1D0)),
INTENT(INOUT) :: albEveTr_id
10382 REAL(KIND(1D0)),
INTENT(INOUT) :: albGrass_id
10383 REAL(KIND(1D0)),
INTENT(INOUT) :: porosity_id
10384 REAL(KIND(1D0)),
INTENT(INOUT) :: Tmin_id
10385 REAL(KIND(1D0)),
INTENT(INOUT) :: Tmax_id
10386 REAL(KIND(1D0)),
INTENT(INOUT) :: lenday_id
10389 REAL(KIND(1D0)),
DIMENSION(12),
INTENT(INOUT) :: HDD_id
10392 REAL(KIND(1D0)),
DIMENSION(9),
INTENT(INOUT) :: WUDay_id
10395 REAL(KIND(1D0)),
INTENT(INOUT) :: Tair_av
10400 REAL(KIND(1D0)),
DIMENSION(nlayer),
INTENT(INOUT) :: tsfc_roof
10401 REAL(KIND(1D0)),
DIMENSION(nlayer) :: sfr_roof
10402 REAL(KIND(1D0)),
DIMENSION(nlayer),
INTENT(in) :: tin_roof
10403 REAL(KIND(1D0)),
DIMENSION(nlayer, ndepth),
INTENT(inout) :: temp_roof
10404 REAL(KIND(1D0)),
DIMENSION(nlayer, ndepth),
INTENT(in) :: k_roof
10405 REAL(KIND(1D0)),
DIMENSION(nlayer, ndepth),
INTENT(in) :: cp_roof
10406 REAL(KIND(1D0)),
DIMENSION(nlayer, ndepth),
INTENT(in) :: dz_roof
10409 REAL(KIND(1D0)),
DIMENSION(nlayer),
INTENT(INOUT) :: tsfc_wall
10410 REAL(KIND(1D0)),
DIMENSION(nlayer) :: sfr_wall
10411 REAL(KIND(1D0)),
DIMENSION(nlayer),
INTENT(in) :: tin_wall
10412 REAL(KIND(1D0)),
DIMENSION(nlayer, ndepth),
INTENT(inout) :: temp_wall
10413 REAL(KIND(1D0)),
DIMENSION(nlayer, ndepth),
INTENT(in) :: k_wall
10414 REAL(KIND(1D0)),
DIMENSION(nlayer, ndepth),
INTENT(in) :: cp_wall
10415 REAL(KIND(1D0)),
DIMENSION(nlayer, ndepth),
INTENT(in) :: dz_wall
10418 REAL(KIND(1D0)),
DIMENSION(nsurf),
INTENT(INOUT) :: tsfc_surf
10419 REAL(KIND(1D0)),
DIMENSION(nsurf),
INTENT(in) :: tin_surf
10420 REAL(KIND(1D0)),
DIMENSION(nsurf, ndepth),
INTENT(inout) :: temp_surf
10421 REAL(KIND(1D0)),
DIMENSION(nsurf, ndepth),
INTENT(in) :: k_surf
10422 REAL(KIND(1D0)),
DIMENSION(nsurf, ndepth),
INTENT(in) :: cp_surf
10423 REAL(KIND(1D0)),
DIMENSION(nsurf, ndepth),
INTENT(in) :: dz_surf
10426 INTEGER,
INTENT(IN) :: n_vegetation_region_urban, &
10427 n_stream_sw_urban, n_stream_lw_urban
10428 REAL(KIND(1D0)),
INTENT(IN) :: sw_dn_direct_frac, air_ext_sw, air_ssa_sw, &
10429 veg_ssa_sw, air_ext_lw, air_ssa_lw, veg_ssa_lw, &
10430 veg_fsd_const, veg_contact_fraction_const, &
10431 ground_albedo_dir_mult_fact
10432 LOGICAL,
INTENT(IN) :: use_sw_direct_albedo
10433 REAL(KIND(1D0)),
DIMENSION(nlayer + 1),
INTENT(IN) :: height
10434 REAL(KIND(1D0)),
DIMENSION(nlayer),
INTENT(IN) :: building_frac
10435 REAL(KIND(1D0)),
DIMENSION(nlayer),
INTENT(IN) :: veg_frac
10436 REAL(KIND(1D0)),
DIMENSION(nlayer),
INTENT(IN) :: building_scale
10437 REAL(KIND(1D0)),
DIMENSION(nlayer),
INTENT(IN) :: veg_scale
10438 REAL(KIND(1D0)),
DIMENSION(nlayer),
INTENT(IN) :: alb_roof
10439 REAL(KIND(1D0)),
DIMENSION(nlayer),
INTENT(IN) :: emis_roof
10440 REAL(KIND(1D0)),
DIMENSION(nlayer),
INTENT(IN) :: alb_wall
10441 REAL(KIND(1D0)),
DIMENSION(nlayer),
INTENT(IN) :: emis_wall
10442 REAL(KIND(1D0)),
DIMENSION(nspec, nlayer),
INTENT(IN) :: roof_albedo_dir_mult_fact
10443 REAL(KIND(1D0)),
DIMENSION(nspec, nlayer),
INTENT(IN) :: wall_specular_frac
10448 REAL(KIND(1D0)),
DIMENSION(nlayer),
INTENT(IN) :: SoilStoreCap_roof
10449 REAL(KIND(1D0)),
DIMENSION(nlayer),
INTENT(IN) :: StateLimit_roof
10450 REAL(KIND(1D0)),
DIMENSION(nlayer),
INTENT(IN) :: wetthresh_roof
10451 REAL(KIND(1D0)),
DIMENSION(nlayer),
INTENT(INOUT) :: soilstore_roof
10452 REAL(KIND(1D0)),
DIMENSION(nlayer),
INTENT(INOUT) :: state_roof
10454 REAL(KIND(1D0)),
DIMENSION(nlayer),
INTENT(IN) :: SoilStoreCap_wall
10455 REAL(KIND(1D0)),
DIMENSION(nlayer),
INTENT(IN) :: StateLimit_wall
10456 REAL(KIND(1D0)),
DIMENSION(nlayer),
INTENT(IN) :: wetthresh_wall
10457 REAL(KIND(1D0)),
DIMENSION(nlayer),
INTENT(INOUT) :: soilstore_wall
10458 REAL(KIND(1D0)),
DIMENSION(nlayer),
INTENT(INOUT) :: state_wall
10464 REAL(KIND(1D0)),
DIMENSION(len_sim, ncolumnsDataOutSUEWS) :: dataOutBlockSUEWS
10465 REAL(KIND(1D0)),
DIMENSION(len_sim, ncolumnsDataOutSnow) :: dataOutBlockSnow
10466 REAL(KIND(1D0)),
DIMENSION(len_sim, ncolumnsDataOutESTM) :: dataOutBlockESTM
10467 REAL(KIND(1D0)),
DIMENSION(len_sim, ncolumnsDataOutEHC) :: dataOutBlockEHC
10468 REAL(KIND(1D0)),
DIMENSION(len_sim, ncolumnsDataOutRSL) :: dataOutBlockRSL
10469 REAL(KIND(1D0)),
DIMENSION(len_sim, ncolumnsdataOutBEERS) :: dataOutBlockBEERS
10470 REAL(KIND(1D0)),
DIMENSION(len_sim, ncolumnsDataOutDebug) :: dataOutBlockDebug
10471 REAL(KIND(1D0)),
DIMENSION(len_sim, ncolumnsDataOutSPARTACUS) :: dataOutBlockSPARTACUS
10472 REAL(KIND(1D0)),
DIMENSION(len_sim, ncolumnsDataOutDailyState) :: dataOutBlockDailyState
10489 INTEGER,
PARAMETER :: gridiv_x = 1
10490 REAL(KIND(1D0)) :: qn1_obs
10491 REAL(KIND(1D0)) :: qh_obs
10492 REAL(KIND(1D0)) :: qe_obs
10493 REAL(KIND(1D0)) :: qs_obs
10494 REAL(KIND(1D0)) :: qf_obs
10495 REAL(KIND(1D0)) :: avu1
10496 REAL(KIND(1D0)) :: avrh
10497 REAL(KIND(1D0)) :: Temp_C
10498 REAL(KIND(1D0)) :: Press_hPa
10499 REAL(KIND(1D0)) :: Precip
10500 REAL(KIND(1D0)) :: avkdn
10501 REAL(KIND(1D0)) :: snowFrac_obs
10502 REAL(KIND(1D0)) :: ldown_obs
10503 REAL(KIND(1D0)) :: fcld_obs
10504 REAL(KIND(1D0)) :: wu_m3
10505 REAL(KIND(1D0)) :: xsmd
10506 REAL(KIND(1D0)) :: LAI_obs
10507 REAL(KIND(1D0)) :: kdiff
10508 REAL(KIND(1D0)) :: kdir
10509 REAL(KIND(1D0)) :: wdir
10511 REAL(KIND(1D0)),
DIMENSION(5) :: datetimeLine
10512 REAL(KIND(1D0)),
DIMENSION(ncolumnsDataOutSUEWS - 5) :: dataOutLineSUEWS
10513 REAL(KIND(1D0)),
DIMENSION(ncolumnsDataOutSnow - 5) :: dataOutLineSnow
10514 REAL(KIND(1D0)),
DIMENSION(ncolumnsDataOutESTM - 5) :: dataOutLineESTM
10515 REAL(KIND(1D0)),
DIMENSION(ncolumnsDataOutEHC - 5) :: dataOutLineEHC
10516 REAL(KIND(1D0)),
DIMENSION(ncolumnsDataOutRSL - 5) :: dataOutLineRSL
10518 REAL(KIND(1D0)),
DIMENSION(ncolumnsDataOutBEERS - 5) :: dataOutLineBEERS
10519 REAL(KIND(1D0)),
DIMENSION(ncolumnsDataOutDebug - 5) :: dataOutLinedebug
10520 REAL(KIND(1D0)),
DIMENSION(ncolumnsDataOutSPARTACUS - 5) :: dataOutLineSPARTACUS
10521 REAL(KIND(1D0)),
DIMENSION(ncolumnsDataOutDailyState - 5) :: dataOutLineDailyState
10523 REAL(KIND(1D0)),
DIMENSION(len_sim, ncolumnsDataOutSUEWS, 1) :: dataOutBlockSUEWS_X
10524 REAL(KIND(1D0)),
DIMENSION(len_sim, ncolumnsDataOutSnow, 1) :: dataOutBlockSnow_X
10525 REAL(KIND(1D0)),
DIMENSION(len_sim, ncolumnsDataOutESTM, 1) :: dataOutBlockESTM_X
10526 REAL(KIND(1D0)),
DIMENSION(len_sim, ncolumnsDataOutEHC, 1) :: dataOutBlockEHC_X
10527 REAL(KIND(1D0)),
DIMENSION(len_sim, ncolumnsDataOutRSL, 1) :: dataOutBlockRSL_X
10528 REAL(KIND(1D0)),
DIMENSION(len_sim, ncolumnsdataOutBEERS, 1) :: dataOutBlockBEERS_X
10529 REAL(KIND(1D0)),
DIMENSION(len_sim, ncolumnsDataOutDebug, 1) :: dataOutBlockDebug_X
10530 REAL(KIND(1D0)),
DIMENSION(len_sim, ncolumnsDataOutSPARTACUS, 1) :: dataOutBlockSPARTACUS_X
10531 REAL(KIND(1D0)),
DIMENSION(len_sim, ncolumnsDataOutDailyState, 1) :: dataOutBlockDailyState_X
10541 DO ir = 1, len_sim, 1
10545 iy = int(metforcingblock(ir, 1))
10546 id = int(metforcingblock(ir, 2))
10547 it = int(metforcingblock(ir, 3))
10548 imin = int(metforcingblock(ir, 4))
10550 qn1_obs = metforcingblock(ir, 5)
10551 qh_obs = metforcingblock(ir, 6)
10552 qe_obs = metforcingblock(ir, 7)
10553 qs_obs = metforcingblock(ir, 8)
10554 qf_obs = metforcingblock(ir, 9)
10555 avu1 = metforcingblock(ir, 10)
10556 avrh = metforcingblock(ir, 11)
10557 temp_c = metforcingblock(ir, 12)
10558 press_hpa = metforcingblock(ir, 13)
10559 precip = metforcingblock(ir, 14)
10560 avkdn = metforcingblock(ir, 15)
10561 snowfrac_obs = metforcingblock(ir, 16)
10562 ldown_obs = metforcingblock(ir, 17)
10563 fcld_obs = metforcingblock(ir, 18)
10564 wu_m3 = metforcingblock(ir, 19)
10565 xsmd = metforcingblock(ir, 20)
10566 lai_obs = metforcingblock(ir, 21)
10567 kdiff = metforcingblock(ir, 22)
10568 kdir = metforcingblock(ir, 23)
10569 wdir = metforcingblock(ir, 24)
10804 ah_min, ahprof_24hr, ah_slope_cooling, &
10805 ah_slope_heating, &
10806 alb, albmax_dectr, albmax_evetr, albmax_grass, &
10807 albmin_dectr, albmin_evetr, albmin_grass, &
10808 alpha_bioco2, alpha_enh_bioco2, alt, avkdn, avrh, avu1, baset, basete, &
10809 beta_bioco2, beta_enh_bioco2, bldgh, capmax_dec, capmin_dec, &
10810 chanohm, co2pointsource, cpanohm, crwmax, crwmin, daywat, daywatper, &
10811 dectreeh, diagmethod, diagnose, drainrt, &
10812 dt_since_start, dqndt, qn_av, dqnsdt, qn_s_av, &
10813 ef_umolco2perj, emis, emissionsmethod, enef_v_jkm, enddls, evetreeh, faibldg, &
10814 faidectree, faievetree, faimethod, faut, fcef_v_kgkm, fcld_obs, flowchange, &
10815 frfossilfuel_heat, frfossilfuel_nonheat, g_max, g_k, g_q_base, g_q_shape, g_t, g_sm, gdd_id, &
10816 gddfull, gridiv, gsmodel, h_maintain, hdd_id, humactivity_24hr, &
10817 icefrac, id, ie_a, ie_end, ie_m, ie_start, imin, &
10818 internalwateruse_h, &
10819 irrfracpaved, irrfracbldgs, &
10820 irrfracevetr, irrfracdectr, irrfracgrass, &
10821 irrfracbsoil, irrfracwater, &
10823 iy, kkanohm, kmax, lai_id, laimax, laimin, lai_obs, &
10824 laipower, laitype, lat, lenday_id, ldown_obs, lng, maxconductance, maxfcmetab, maxqfmetab, &
10825 snowwater, minfcmetab, minqfmetab, min_res_bioco2, &
10826 narp_emis_snow, narp_trans_site, netradiationmethod, &
10828 n_vegetation_region_urban, &
10829 n_stream_sw_urban, n_stream_lw_urban, &
10830 sw_dn_direct_frac, air_ext_sw, air_ssa_sw, &
10831 veg_ssa_sw, air_ext_lw, air_ssa_lw, veg_ssa_lw, &
10832 veg_fsd_const, veg_contact_fraction_const, &
10833 ground_albedo_dir_mult_fact, use_sw_direct_albedo, &
10834 height, building_frac, veg_frac, building_scale, veg_scale, &
10835 alb_roof, emis_roof, alb_wall, emis_wall, &
10836 roof_albedo_dir_mult_fact, wall_specular_frac, &
10837 ohm_coef, ohmincqf, ohm_threshsw, &
10838 ohm_threshwd, pipecapacity, popdensdaytime, &
10839 popdensnighttime, popprof_24hr, pormax_dec, pormin_dec, &
10840 precip, preciplimit, preciplimitalb, press_hpa, &
10841 qf0_beu, qf_a, qf_b, qf_c, &
10842 qn1_obs, qs_obs, qf_obs, &
10843 radmeltfact, raincover, rainmaxres, resp_a, resp_b, &
10844 roughlenheatmethod, roughlenmommethod, runofftowater, s1, s2, &
10845 sathydraulicconduct, sddfull, sdd_id, smdmethod, snowalb, snowalbmax, &
10846 snowalbmin, snowpacklimit, snowdens, snowdensmax, snowdensmin, snowfallcum, snowfrac, &
10847 snowlimbldg, snowlimpaved, snowfrac_obs, snowpack, snowprof_24hr, snowuse, soildepth, &
10848 stabilitymethod, startdls, &
10849 soilstore_surf, soilstorecap_surf, state_surf, statelimit_surf, wetthresh_surf, &
10850 soilstore_roof, soilstorecap_roof, state_roof, statelimit_roof, wetthresh_roof, &
10851 soilstore_wall, soilstorecap_wall, state_wall, statelimit_wall, wetthresh_wall, &
10852 storageheatmethod, storedrainprm, surfacearea, tair_av, tau_a, tau_f, tau_r, &
10853 tmax_id, tmin_id, &
10854 baset_cooling, baset_heating, temp_c, tempmeltfact, th, &
10855 theta_bioco2, timezone, tl, trafficrate, trafficunits, &
10857 tsfc_roof, tsfc_wall, tsfc_surf, &
10858 temp_roof, temp_wall, temp_surf, &
10859 tin_roof, tin_wall, tin_surf, &
10860 k_roof, k_wall, k_surf, &
10861 cp_roof, cp_wall, cp_surf, &
10862 dz_roof, dz_wall, dz_surf, &
10863 traffprof_24hr, ts5mindata_ir, tstep, tstep_prev, veg_type, &
10864 waterdist, waterusemethod, wu_m3, &
10865 wuday_id, decidcap_id, albdectr_id, albevetr_id, albgrass_id, porosity_id, &
10866 wuprofa_24hr, wuprofm_24hr, xsmd, z, z0m_in, zdm_in, &
10875 dt_since_start = dt_since_start + tstep
10878 dataoutblockdailystate(ir, :) = [output_line_suews%dataOutLineDailyState]
10883 snowuse, storageheatmethod, &
10886 output_line_suews%dataOutLineSUEWS, &
10887 output_line_suews%dataOutLineSnow, &
10888 output_line_suews%dataOutLineESTM, &
10889 output_line_suews%dataoutLineRSL, &
10890 output_line_suews%dataOutLineBEERS, &
10891 output_line_suews%dataOutLinedebug, &
10892 output_line_suews%dataOutLineSPARTACUS, &
10893 output_line_suews%dataOutLineEHC, &
10894 dataoutblocksuews_x, dataoutblocksnow_x, dataoutblockestm_x, &
10895 dataoutblockrsl_x, dataoutblockbeers_x, dataoutblockdebug_x, dataoutblockspartacus_x, dataoutblockehc_x)
10899 dataoutblocksuews = dataoutblocksuews_x(:, :, 1)
10900 dataoutblocksnow = dataoutblocksnow_x(:, :, 1)
10901 dataoutblockestm = dataoutblockestm_x(:, :, 1)
10902 dataoutblockehc = dataoutblockehc_x(:, :, 1)
10903 dataoutblockrsl = dataoutblockrsl_x(:, :, 1)
10904 dataoutblockbeers = dataoutblockbeers_x(:, :, 1)
10905 dataoutblockdebug = dataoutblockdebug_x(:, :, 1)
10906 dataoutblockspartacus = dataoutblockspartacus_x(:, :, 1)
10913 output_block_suews%dataOutBlockSUEWS = dataoutblocksuews
10914 output_block_suews%dataOutBlockSnow = dataoutblocksnow
10915 output_block_suews%dataOutBlockESTM = dataoutblockestm
10916 output_block_suews%dataOutBlockEHC = dataoutblockehc
10917 output_block_suews%dataOutBlockRSL = dataoutblockrsl
10918 output_block_suews%dataOutBlockBEERS = dataoutblockbeers
10919 output_block_suews%dataOutBlockDebug = dataoutblockdebug
10920 output_block_suews%dataOutBlockSPARTACUS = dataoutblockspartacus
10921 output_block_suews%dataOutBlockDailyState = dataoutblockdailystate