PM-JPL 1.2.1__py3-none-any.whl → 1.7.0__py3-none-any.whl

PMJPL/ECOv002-static-tower-PM-JPL-inputs.csv

@@ -0,0 +1,122 @@
+ID,name,elevation_km,IGBP,gl_sh,gl_e_wv,RBL_min,RBL_max,CL,Tmin_open,Tmin_closed,VPD_open,VPD_closed,geometry
+US-NC3,NC_Clearcut#3,0.005,12,0.02,0.02,20.0,50.0,0.007,12.02,-8.0,650.0,4500.0,POINT (-76.656 35.799)
+PE-QFR,Quistococha Forest Reserve,0.104,2,0.01,0.01,70.0,100.0,0.0025,9.09,-8.0,1000.0,4000.0,POINT (-73.319 -3.8344)
+US-Mi3,LTAR UCB (Upper Chesapeake Bay) Miscanthus 3,0.27,8,0.08,0.08,25.0,45.0,0.0065,11.39,-8.0,650.0,3500.0,POINT (-80.637 41.8222)
+US-NC4,NC_AlligatorRiver,0.001,5,0.04,0.04,65.0,95.0,0.0025,9.5,-7.0,650.0,2900.0,POINT (-75.9038 35.7879)
+CA-DB2,Delta Burns Bog 2,0.004,8,0.08,0.08,25.0,45.0,0.0065,11.39,-8.0,650.0,3500.0,POINT (-122.9951 49.119)
+US-Sne,Sherman Island Restored Wetland,-0.005,13,0.02,0.02,20.0,50.0,0.007,12.02,-8.0,650.0,4500.0,POINT (-121.7547 38.0369)
+US-Mi1,LTAR UCB (Upper Chesapeake Bay) Miscanthus 1,0.29,12,0.02,0.02,20.0,50.0,0.007,12.02,-8.0,650.0,4500.0,POINT (-80.6313 41.7727)
+US-PFe,NW4 Lake-1 CHEESEHEAD 2019,0.48,5,0.04,0.04,65.0,95.0,0.0025,9.5,-7.0,650.0,2900.0,POINT (-90.3004 45.9793)
+US-NR1,Niwot Ridge Forest (LTER NWT1),3.05,8,0.08,0.08,25.0,45.0,0.0065,11.39,-8.0,650.0,3500.0,POINT (-105.5464 40.0329)
+US-Vcp,Valles Caldera Ponderosa Pine,2.542,10,0.02,0.02,20.0,50.0,0.007,12.02,-8.0,650.0,4200.0,POINT (-106.5974 35.8624)
+US-xAB,NEON Abby Road (ABBY),0.363,1,0.04,0.04,65.0,95.0,0.0032,8.31,-8.0,650.0,3000.0,POINT (-122.3303 45.7624)
+US-HBK,Hubbard Brook Experimental Forest,0.367,4,0.01,0.01,65.0,100.0,0.0028,9.94,-6.0,650.0,2900.0,POINT (-71.7181 43.9397)
+US-EDN,Eden Landing Ecological Reserve,,13,0.02,0.02,20.0,50.0,0.007,12.02,-8.0,650.0,4500.0,POINT (-122.114 37.6156)
+US-PFh,NE2 Pine-3 CHEESEHEAD 2019,0.463,5,0.04,0.04,65.0,95.0,0.0025,9.5,-7.0,650.0,2900.0,POINT (-90.2406 45.9557)
+US-Me6,Metolius Young Pine Burn,0.998,8,0.08,0.08,25.0,45.0,0.0065,11.39,-8.0,650.0,3500.0,POINT (-121.6078 44.3233)
+US-NC2,NC_Loblolly Plantation,0.005,5,0.04,0.04,65.0,95.0,0.0025,9.5,-7.0,650.0,2900.0,POINT (-76.6685 35.803)
+US-Whs,Walnut Gulch Lucky Hills Shrub,1.37,10,0.02,0.02,20.0,50.0,0.007,12.02,-8.0,650.0,4200.0,POINT (-110.0522 31.7438)
+US-SRM,Santa Rita Mesquite,1.12,7,0.04,0.04,20.0,55.0,0.0065,8.8,-8.0,650.0,4400.0,POINT (-110.8661 31.8214)
+US-CS1,Central Sands Irrigated Agricultural Field,0.328,12,0.02,0.02,20.0,50.0,0.007,12.02,-8.0,650.0,4500.0,POINT (-89.5379 44.1031)
+US-PFs,SE5 Aspen-5 CHEESEHEAD 2019,0.489,5,0.04,0.04,65.0,95.0,0.0025,9.5,-7.0,650.0,2900.0,POINT (-90.2382 45.9381)
+US-PFg,NE1 Pine-2 CHEESEHEAD 2019,0.475,5,0.04,0.04,65.0,95.0,0.0025,9.5,-7.0,650.0,2900.0,POINT (-90.2723 45.9735)
+US-Ha1,Harvard Forest EMS Tower (HFR1),0.34,4,0.01,0.01,65.0,100.0,0.0028,9.94,-6.0,650.0,2900.0,POINT (-72.1715 42.5378)
+US-Ne2,Mead - irrigated maize-soybean rotation site,0.362,12,0.02,0.02,20.0,50.0,0.007,12.02,-8.0,650.0,4500.0,POINT (-96.4701 41.1649)
+US-PFm,SW3 Hardwood-2 CHEESEHEAD 2019,0.484,5,0.04,0.04,65.0,95.0,0.0025,9.5,-7.0,650.0,2900.0,POINT (-90.3099 45.9207)
+US-PFr,SE4 Tussock-2 CHEESEHEAD 2019,0.49,5,0.04,0.04,65.0,95.0,0.0025,9.5,-7.0,650.0,2900.0,POINT (-90.2475 45.9245)
+US-PFL,SW2 Aspen-3 CHEESEHEAD 2019,0.464,5,0.04,0.04,65.0,95.0,0.0025,9.5,-7.0,650.0,2900.0,POINT (-90.3177 45.9409)
+US-PFj,NE4 Maple-1 CHEESEHEAD 2019,0.484,5,0.04,0.04,65.0,95.0,0.0025,9.5,-7.0,650.0,2900.0,POINT (-90.227 45.9619)
+US-CC2,Coloma Corn 2,0.32,12,0.02,0.02,20.0,50.0,0.007,12.02,-8.0,650.0,4500.0,POINT (-89.6196 44.1039)
+US-NR3,Niwot Ridge Alpine (T-Van West),3.504,8,0.08,0.08,25.0,45.0,0.0065,11.39,-8.0,650.0,3500.0,POINT (-105.5864 40.052)
+US-PFt,SE6 Pine-4 CHEESEHEAD 2019,0.466,5,0.04,0.04,65.0,95.0,0.0025,9.5,-7.0,650.0,2900.0,POINT (-90.2288 45.9197)
+CA-Cbo,"Ontario - Mixed Deciduous, Borden Forest Site",0.12,9,0.08,0.08,25.0,45.0,0.0065,11.39,-8.0,650.0,3600.0,POINT (-79.9333 44.3167)
+PR-xLA,NEON Lajas Experimental Station (LAJA),0.024,12,0.02,0.02,20.0,50.0,0.007,12.02,-8.0,650.0,4500.0,POINT (-67.0769 18.0212)
+US-Ne3,Mead - rainfed maize-soybean rotation site,0.363,12,0.02,0.02,20.0,50.0,0.007,12.02,-8.0,650.0,4500.0,POINT (-96.4397 41.1797)
+US-PHM,Plum Island High Marsh,0.0014,8,0.08,0.08,25.0,45.0,0.0065,11.39,-8.0,650.0,3500.0,POINT (-70.8301 42.7423)
+US-Tw1,Twitchell Wetland West Pond,-0.005,12,0.02,0.02,20.0,50.0,0.007,12.02,-8.0,650.0,4500.0,POINT (-121.6469 38.1074)
+US-Hn3,Hanford 100H sagebrush,0.12090000000000001,10,0.02,0.02,20.0,50.0,0.007,12.02,-8.0,650.0,4200.0,POINT (-119.4614 46.6878)
+US-PFk,SW1 Aspen-2 CHEESEHEAD 2019,0.476,5,0.04,0.04,65.0,95.0,0.0025,9.5,-7.0,650.0,2900.0,POINT (-90.3425 45.9149)
+US-CS2,Tri county school Pine Forest,0.328,12,0.02,0.02,20.0,50.0,0.007,12.02,-8.0,650.0,4500.0,POINT (-89.5002 44.1467)
+US-ONA,Florida pine flatwoods,0.025,9,0.08,0.08,25.0,45.0,0.0065,11.39,-8.0,650.0,3600.0,POINT (-81.9509 27.3836)
+US-CC1,Coloma Corn 1,0.314,8,0.08,0.08,25.0,45.0,0.0065,11.39,-8.0,650.0,3500.0,POINT (-89.6787 44.0732)
+US-Rls,RCEW Low Sagebrush,1.608,10,0.02,0.02,20.0,50.0,0.007,12.02,-8.0,650.0,4200.0,POINT (-116.7356 43.1439)
+US-UMd,UMBS Disturbance,0.239,4,0.01,0.01,65.0,100.0,0.0028,9.94,-6.0,650.0,2900.0,POINT (-84.6975 45.5625)
+US-PFq,SE3 Aspen-4 CHEESEHEAD 2019,0.473,5,0.04,0.04,65.0,95.0,0.0025,9.5,-7.0,650.0,2900.0,POINT (-90.2475 45.9271)
+US-PFn,SW4 Hardwood-3 CHEESEHEAD 2019,0.478,5,0.04,0.04,65.0,95.0,0.0025,9.5,-7.0,650.0,2900.0,POINT (-90.2823 45.9392)
+US-HB3,Hobcaw Barony Longleaf Pine Restoration,0.0073,1,0.04,0.04,65.0,95.0,0.0032,8.31,-8.0,650.0,3000.0,POINT (-79.2322 33.3482)
+US-DFC,"US Dairy Forage Research Center, Prairie du Sac",0.26489999999999997,12,0.02,0.02,20.0,50.0,0.007,12.02,-8.0,650.0,4500.0,POINT (-89.7117 43.3448)
+US-Wkg,Walnut Gulch Kendall Grasslands,1.531,10,0.02,0.02,20.0,50.0,0.007,12.02,-8.0,650.0,4200.0,POINT (-109.9419 31.7365)
+US-CF1,CAF-LTAR Cook East,0.794,12,0.02,0.02,20.0,50.0,0.007,12.02,-8.0,650.0,4500.0,POINT (-117.0821 46.7815)
+US-Ho2,Howland Forest (west tower),0.061,5,0.04,0.04,65.0,95.0,0.0025,9.5,-7.0,650.0,2900.0,POINT (-68.747 45.2091)
+US-Snf,Sherman Barn,-0.004,0,0.02,0.02,20.0,50.0,0.007,12.02,-8.0,650.0,4500.0,POINT (-121.7272 38.0402)
+US-Me2,Metolius mature ponderosa pine,1.253,1,0.04,0.04,65.0,95.0,0.0032,8.31,-8.0,650.0,3000.0,POINT (-121.5574 44.4523)
+US-CF3,CAF-LTAR Boyd North,0.795,12,0.02,0.02,20.0,50.0,0.007,12.02,-8.0,650.0,4500.0,POINT (-117.1261 46.7551)
+US-Ne1,Mead - irrigated continuous maize site,0.361,12,0.02,0.02,20.0,50.0,0.007,12.02,-8.0,650.0,4500.0,POINT (-96.4766 41.1651)
+US-Ha2,Harvard Forest Hemlock Site,0.36,4,0.01,0.01,65.0,100.0,0.0028,9.94,-6.0,650.0,2900.0,POINT (-72.1779 42.5393)
+CA-Ca3,British Columbia - Pole sapling Douglas-fir stand,,1,0.04,0.04,65.0,95.0,0.0032,8.31,-8.0,650.0,3000.0,POINT (-124.9004 49.5346)
+US-SP1,Slashpine-Austin Cary- 65yrs nat regen,0.05,8,0.08,0.08,25.0,45.0,0.0065,11.39,-8.0,650.0,3500.0,POINT (-82.2188 29.7381)
+US-UMB,Univ. of Mich. Biological Station,0.234,4,0.01,0.01,65.0,100.0,0.0028,9.94,-6.0,650.0,2900.0,POINT (-84.7138 45.5598)
+US-Bi1,Bouldin Island Alfalfa,-0.0027,12,0.02,0.02,20.0,50.0,0.007,12.02,-8.0,650.0,4500.0,POINT (-121.4993 38.0992)
+US-PFi,NE3 Hardwood-1 CHEESEHEAD 2019,0.488,5,0.04,0.04,65.0,95.0,0.0025,9.5,-7.0,650.0,2900.0,POINT (-90.2327 45.9749)
+US-CF2,CAF-LTAR Cook West,0.807,12,0.02,0.02,20.0,50.0,0.007,12.02,-8.0,650.0,4500.0,POINT (-117.0908 46.784)
+US-KFS,Kansas Field Station,0.31,9,0.08,0.08,25.0,45.0,0.0065,11.39,-8.0,650.0,3600.0,POINT (-95.1907 39.0561)
+US-CMW,Charleston Mesquite Woodland,1.199,10,0.02,0.02,20.0,50.0,0.007,12.02,-8.0,650.0,4200.0,POINT (-110.1777 31.6637)
+US-Ro4,Rosemount Prairie,0.274,12,0.02,0.02,20.0,50.0,0.007,12.02,-8.0,650.0,4500.0,POINT (-93.0723 44.6781)
+US-PAS,"Florida, Paspalum notatum pasture",0.027100000000000003,9,0.08,0.08,25.0,45.0,0.0065,11.39,-8.0,650.0,3600.0,POINT (-81.951 27.3944)
+US-xYE,NEON Yellowstone Northern Range (Frog Rock) (YELL),2.116,10,0.02,0.02,20.0,50.0,0.007,12.02,-8.0,650.0,4200.0,POINT (-110.5391 44.9535)
+US-Los,Lost Creek,0.48,5,0.04,0.04,65.0,95.0,0.0025,9.5,-7.0,650.0,2900.0,POINT (-89.9792 46.0827)
+US-PFd,NW3 Tussock-1 CHEESEHEAD 2019,0.473,5,0.04,0.04,65.0,95.0,0.0025,9.5,-7.0,650.0,2900.0,POINT (-90.301 45.9689)
+US-KM4,KBS Marshall Farms Smooth Brome Grass (Ref),0.24630000000000002,9,0.08,0.08,25.0,45.0,0.0065,11.39,-8.0,650.0,3600.0,POINT (-85.3301 42.4423)
+US-Syv,Sylvania Wilderness Area,0.54,5,0.04,0.04,65.0,95.0,0.0025,9.5,-7.0,650.0,2900.0,POINT (-89.3477 46.242)
+US-WCr,Willow Creek,0.52,4,0.01,0.01,65.0,100.0,0.0028,9.94,-6.0,650.0,2900.0,POINT (-90.0799 45.8059)
+US-Rws,Reynolds Creek Wyoming big sagebrush,1.425,10,0.02,0.02,20.0,50.0,0.007,12.02,-8.0,650.0,4200.0,POINT (-116.7132 43.1675)
+US-xJE,NEON Jones Ecological Research Center (JERC),0.044,8,0.08,0.08,25.0,45.0,0.0065,11.39,-8.0,650.0,3500.0,POINT (-84.4686 31.1948)
+US-SRG,Santa Rita Grassland,1.291,7,0.04,0.04,20.0,55.0,0.0065,8.8,-8.0,650.0,4400.0,POINT (-110.8277 31.7894)
+US-xNW,NEON Niwot Ridge Mountain Research Station (NIWO),3.513,8,0.08,0.08,25.0,45.0,0.0065,11.39,-8.0,650.0,3500.0,POINT (-105.5824 40.0543)
+US-Rms,RCEW Mountain Big Sagebrush,2.111,10,0.02,0.02,20.0,50.0,0.007,12.02,-8.0,650.0,4200.0,POINT (-116.7486 43.0645)
+US-NR4,Niwot Ridge Alpine (T-Van East),3.502,8,0.08,0.08,25.0,45.0,0.0065,11.39,-8.0,650.0,3500.0,POINT (-105.5859 40.052)
+US-xDS,NEON Disney Wilderness Preserve (DSNY),0.015,9,0.08,0.08,25.0,45.0,0.0065,11.39,-8.0,650.0,3600.0,POINT (-81.4362 28.125)
+US-xJR,NEON Jornada LTER (JORN),1.329,7,0.04,0.04,20.0,55.0,0.0065,8.8,-8.0,650.0,4400.0,POINT (-106.8425 32.5907)
+US-Bar,Bartlett Experimental Forest,0.272,4,0.01,0.01,65.0,100.0,0.0028,9.94,-6.0,650.0,2900.0,POINT (-71.2881 44.0646)
+US-ALQ,Allequash Creek Site,,5,0.04,0.04,65.0,95.0,0.0025,9.5,-7.0,650.0,2900.0,POINT (-89.6067 46.0308)
+US-MMS,Morgan Monroe State Forest,0.275,4,0.01,0.01,65.0,100.0,0.0028,9.94,-6.0,650.0,2900.0,POINT (-86.4131 39.3232)
+US-UC1,LTAR UCB (Upper Chesapeake Bay) EC1,0.39,4,0.01,0.01,65.0,100.0,0.0028,9.94,-6.0,650.0,2900.0,POINT (-78.0056 40.7536)
+US-PFb,NW1 Pine-1 CHEESEHEAD 2019,0.474,5,0.04,0.04,65.0,95.0,0.0025,9.5,-7.0,650.0,2900.0,POINT (-90.3232 45.972)
+US-xDL,NEON Dead Lake (DELA),0.022,13,0.02,0.02,20.0,50.0,0.007,12.02,-8.0,650.0,4500.0,POINT (-87.8039 32.5417)
+US-xPU,NEON Pu'u Maka'ala Natural Area Reserve (PUUM),1.685,8,0.08,0.08,25.0,45.0,0.0065,11.39,-8.0,650.0,3500.0,POINT (-155.3173 19.5531)
+US-xBL,NEON Blandy Experimental Farm (BLAN),0.183,9,0.08,0.08,25.0,45.0,0.0065,11.39,-8.0,650.0,3600.0,POINT (-78.0716 39.0603)
+US-xST,NEON Steigerwaldt Land Services (STEI),0.481,4,0.01,0.01,65.0,100.0,0.0028,9.94,-6.0,650.0,2900.0,POINT (-89.5864 45.5089)
+US-xSB,NEON Ordway-Swisher Biological Station (OSBS),0.045,8,0.08,0.08,25.0,45.0,0.0065,11.39,-8.0,650.0,3500.0,POINT (-81.9934 29.6893)
+US-xRN,NEON Oak Ridge National Lab (ORNL),0.334,4,0.01,0.01,65.0,100.0,0.0028,9.94,-6.0,650.0,2900.0,POINT (-84.2826 35.9641)
+US-UC2,LTAR UCB (Upper Chesapeake Bay) EC2,0.396,4,0.01,0.01,65.0,100.0,0.0028,9.94,-6.0,650.0,2900.0,POINT (-77.9998 40.7559)
+US-xTR,NEON Treehaven (TREE),0.472,4,0.01,0.01,65.0,100.0,0.0028,9.94,-6.0,650.0,2900.0,POINT (-89.5857 45.4937)
+US-xRM,"NEON Rocky Mountain National Park, CASTNET (RMNP)",2.743,8,0.08,0.08,25.0,45.0,0.0065,11.39,-8.0,650.0,3500.0,POINT (-105.5459 40.2759)
+US-Tw4,Twitchell East End Wetland,-0.005,12,0.02,0.02,20.0,50.0,0.007,12.02,-8.0,650.0,4500.0,POINT (-121.6413 38.1027)
+US-Bi2,Bouldin Island corn,-0.005,12,0.02,0.02,20.0,50.0,0.007,12.02,-8.0,650.0,4500.0,POINT (-121.5351 38.1091)
+US-Jo2,Jornada Experimental Range Mixed Shrubland,1.469,7,0.04,0.04,20.0,55.0,0.0065,8.8,-8.0,650.0,4400.0,POINT (-106.6032 32.5849)
+US-CS5,Central Sands Irrigated Agricultural Field,0.328,12,0.02,0.02,20.0,50.0,0.007,12.02,-8.0,650.0,4500.0,POINT (-89.5377 44.1095)
+US-Mpj,Mountainair Pinyon-Juniper Woodland,2.138,10,0.02,0.02,20.0,50.0,0.007,12.02,-8.0,650.0,4200.0,POINT (-106.2377 34.4384)
+US-xWR,NEON Wind River Experimental Forest (WREF),0.407,1,0.04,0.04,65.0,95.0,0.0032,8.31,-8.0,650.0,3000.0,POINT (-121.9519 45.8205)
+US-xUK,NEON The University of Kansas Field Station (UKFS),0.335,10,0.02,0.02,20.0,50.0,0.007,12.02,-8.0,650.0,4200.0,POINT (-95.1921 39.0404)
+US-xUN,NEON University of Notre Dame Environmental Research Center (UNDE),0.518,5,0.04,0.04,65.0,95.0,0.0025,9.5,-7.0,650.0,2900.0,POINT (-89.5373 46.2339)
+US-xBR,NEON Bartlett Experimental Forest (BART),0.232,4,0.01,0.01,65.0,100.0,0.0028,9.94,-6.0,650.0,2900.0,POINT (-71.2873 44.0639)
+US-CF4,CAF-LTAR Boyd South,0.795,12,0.02,0.02,20.0,50.0,0.007,12.02,-8.0,650.0,4500.0,POINT (-117.1285 46.7518)
+US-xCL,NEON LBJ National Grassland (CLBJ),0.259,10,0.02,0.02,20.0,50.0,0.007,12.02,-8.0,650.0,4200.0,POINT (-97.57 33.4012)
+CA-DBB,Delta Burns Bog,0.004,8,0.08,0.08,25.0,45.0,0.0065,11.39,-8.0,650.0,3500.0,POINT (-122.9849 49.1293)
+US-Hn2,Hanford 100H grassland,0.1175,10,0.02,0.02,20.0,50.0,0.007,12.02,-8.0,650.0,4200.0,POINT (-119.4641 46.6889)
+US-PFp,SE2 Hardwood-4 CHEESEHEAD 2019,0.497,5,0.04,0.04,65.0,95.0,0.0025,9.5,-7.0,650.0,2900.0,POINT (-90.2641 45.9365)
+US-xKA,NEON Konza Prairie Biological Station - Relocatable (KONA),1.329,12,0.02,0.02,20.0,50.0,0.007,12.02,-8.0,650.0,4500.0,POINT (-96.6129 39.1104)
+US-CS3,Central Sands Irrigated Agricultural Field,0.328,12,0.02,0.02,20.0,50.0,0.007,12.02,-8.0,650.0,4500.0,POINT (-89.5727 44.1394)
+US-xTE,NEON Lower Teakettle (TEAK),2.147,9,0.08,0.08,25.0,45.0,0.0065,11.39,-8.0,650.0,3600.0,POINT (-119.006 37.0058)
+US-Rwf,RCEW Upper Sheep Prescibed Fire,1.878,10,0.02,0.02,20.0,50.0,0.007,12.02,-8.0,650.0,4200.0,POINT (-116.7231 43.1207)
+US-xML,NEON Mountain Lake Biological Station (MLBS),1.126,4,0.01,0.01,65.0,100.0,0.0028,9.94,-6.0,650.0,2900.0,POINT (-80.5248 37.3783)
+US-UiB,University of Illinois Miscanthus,0.224,12,0.02,0.02,20.0,50.0,0.007,12.02,-8.0,650.0,4500.0,POINT (-88.1984 40.0628)
+PR-xGU,NEON Guanica Forest (GUAN),0.143,8,0.08,0.08,25.0,45.0,0.0065,11.39,-8.0,650.0,3500.0,POINT (-66.8687 17.9696)
+US-ARM,ARM Southern Great Plains site- Lamont,0.314,12,0.02,0.02,20.0,50.0,0.007,12.02,-8.0,650.0,4500.0,POINT (-97.4888 36.6058)
+US-Kon,Konza Prairie LTER (KNZ),0.417,10,0.02,0.02,20.0,50.0,0.007,12.02,-8.0,650.0,4200.0,POINT (-96.5603 39.0824)
+US-HB2,Hobcaw Barony Mature Longleaf Pine,0.0047,1,0.04,0.04,65.0,95.0,0.0032,8.31,-8.0,650.0,3000.0,POINT (-79.244 33.3242)
+US-xSL,"NEON North Sterling, CO (STER)",1.364,10,0.02,0.02,20.0,50.0,0.007,12.02,-8.0,650.0,4200.0,POINT (-103.0293 40.4619)
+US-xWD,NEON Woodworth (WOOD),0.579,12,0.02,0.02,20.0,50.0,0.007,12.02,-8.0,650.0,4500.0,POINT (-99.2414 47.1282)
+US-CS4,Central Sands Irrigated Agricultural Field,0.328,12,0.02,0.02,20.0,50.0,0.007,12.02,-8.0,650.0,4500.0,POINT (-89.5475 44.1597)
+US-xAE,NEON Klemme Range Research Station (OAES),0.516,10,0.02,0.02,20.0,50.0,0.007,12.02,-8.0,650.0,4200.0,POINT (-99.0588 35.4106)
+US-GLE,GLEES,3.197,10,0.02,0.02,20.0,50.0,0.007,12.02,-8.0,650.0,4200.0,POINT (-106.2399 41.3665)

PMJPL/ECOv002_calval_PMJPL_inputs.py

@@ -0,0 +1,19 @@
+import os
+import pandas as pd
+
+def load_ECOv002_calval_PMJPL_inputs() -> pd.DataFrame:
+    """
+    Load the input data for the PM-JPL model from the ECOSTRESS Collection 2 Cal-Val dataset.
+
+    Returns:
+        pd.DataFrame: A DataFrame containing the input data.
+    """
+
+    # Define the path to the input CSV file relative to this module's directory
+    module_dir = os.path.dirname(os.path.abspath(__file__))
+    input_file_path = os.path.join(module_dir, "ECOv002-cal-val-PM-JPL-inputs.csv")
+
+    # Load the input data into a DataFrame
+    inputs_df = pd.read_csv(input_file_path)
+
+    return inputs_df

PMJPL/ECOv002_static_tower_PMJPL_inputs.py

@@ -0,0 +1,19 @@
+import os
+import pandas as pd
+
+def load_ECOv002_static_tower_PMJPL_inputs() -> pd.DataFrame:
+    """
+    Load the input data for the PM-JPL model from the ECOSTRESS Collection 2 Cal-Val dataset.
+
+    Returns:
+        pd.DataFrame: A DataFrame containing the input data.
+    """
+
+    # Define the path to the input CSV file relative to this module's directory
+    module_dir = os.path.dirname(os.path.abspath(__file__))
+    input_file_path = os.path.join(module_dir, "ECOv002-static-tower-PM-JPL-inputs.csv")
+
+    # Load the input data into a DataFrame
+    inputs_df = pd.read_csv(input_file_path)
+
+    return inputs_df

PMJPL/PMJPL.py

@@ -1,3 +1,4 @@
+
 """
 MOD16 model of evapotranspiration
 
@@ -6,401 +7,67 @@ https://landweb.nascom.nasa.gov/QA_WWW/forPage/user_guide/MOD16UsersGuide2016.pd
 
 Developed by Gregory Halverson in the Jet Propulsion Laboratory Year-Round Internship Program (Columbus Technologies and Services), in coordination with the ECOSTRESS mission and master's thesis studies at California State University, Northridge.
 """
+
+# All imports moved to the top
 import logging
-from typing import Callable, Dict, Union
+from typing import Dict, Union
 from datetime import datetime
-from os.path import join, abspath, dirname, expanduser
 
 import numpy as np
-import pandas as pd
-from numpy import where, nan, exp, array, isnan, logical_and, clip, float32
-import warnings
+import rasters as rt
+from rasters import Raster, RasterGrid, RasterGeometry, VectorGeometry
 
-import rasters as rt 
+from check_distribution import check_distribution
+from GEOS5FP import GEOS5FP
+from NASADEM import NASADEM
+from verma_net_radiation import verma_net_radiation
+from SEBAL_soil_heat_flux import calculate_SEBAL_soil_heat_flux
+from MCD12C1_2019_v006 import load_MCD12C1_IGBP
 
-from rasters import Raster, RasterGrid, RasterGeometry
+from carlson_leaf_area_index import carlson_leaf_area_index
+from carlson_fractional_vegetation_cover import carlson_fractional_vegetation_cover
+from carlson_leaf_area_index import carlson_leaf_area_index
 
-from .MCD12C1.MCD12C1 import load_MCD12C1_IGBP
-from .parameters import MOD16_parameter_from_IGBP
-from .evapotranspiration_conversion.evapotranspiration_conversion import lambda_Jkg_from_Ta_C
-from .meteorology_conversion.meteorology_conversion import SVP_Pa_from_Ta_C, calculate_specific_heat, calculate_surface_pressure, celcius_to_kelvin
-from .penman_monteith.penman_monteith import calculate_gamma
-from .priestley_taylor.priestley_taylor import delta_Pa_from_Ta_C, delta_kPa_from_Ta_C
+from daylight_evapotranspiration import lambda_Jkg_from_Ta_C
 
-from .soil_heat_flux import calculate_soil_heat_flux
-from .evapotranspiration_conversion import daily_ET_from_daily_LE
-from .meteorology_conversion import kelvin_to_celsius, calculate_specific_humidity, calculate_air_density
-from .vegetation_conversion.vegetation_conversion import FVC_from_NDVI, LAI_from_NDVI
+from meteorology_conversion import SVP_Pa_from_Ta_C
+from meteorology_conversion import calculate_air_density
+from meteorology_conversion import calculate_specific_heat
+from meteorology_conversion import calculate_specific_humidity
+from meteorology_conversion import calculate_surface_pressure
+from meteorology_conversion import celcius_to_kelvin
 
-from .constants import *
+from priestley_taylor import delta_Pa_from_Ta_C
+from PTJPL import calculate_relative_surface_wetness
+from PTJPL import RH_THRESHOLD, MIN_FWET
 
-from .fwet import calculate_fwet
+from .constants import *
+from .PMJPL_parameter_from_IGBP import PMJPL_parameter_from_IGBP
+from .calculate_gamma import calculate_gamma
 from .soil_moisture_constraint import calculate_fSM
-from .tmin_factor import calculate_tmin_factor
-from .correctance_factor import calculate_rcorr
+from .minimum_temperature_factor import Tmin_factor
+from .correctance_factor import calculate_correctance_factor
 from .VPD_factor import calculate_VPD_factor
-
 from .canopy_conductance import calculate_canopy_conductance
-
 from .wet_canopy_resistance import calculate_wet_canopy_resistance
-from .canopy_aerodynamic_resistance import calculate_rtotc
-
+from .canopy_aerodynamic_resistance import calculate_canopy_aerodynamic_resistance
 from .wet_soil_evaporation import calculate_wet_soil_evaporation
 from .potential_soil_evaporation import calculate_potential_soil_evaporation
 from .interception import calculate_interception
 from .transpiration import calculate_transpiration
+from .leaf_conductance_to_sensible_heat import leaf_conductance_to_sensible_heat
+from .potential_stomatal_conductance import potential_stomatal_conductance
+from .open_minimum_temperature import open_minimum_temperature
+from .closed_minimum_temperature import closed_minimum_temperature
+from .open_vapor_pressure_deficit import open_vapor_pressure_deficit
+from .closed_vapor_pressure_deficit import closed_vapor_pressure_deficit
+from .leaf_conductance_to_evaporated_water import leaf_conductance_to_evaporated_water
+from .maximum_boundary_layer_resistance import maximum_boundary_layer_resistance
+from .minimum_boundary_layer_resistance import minimum_boundary_layer_resistance
+from .model import PMJPL
+from .process_PMJPL_table import process_PMJPL_table
+from .ECOv002_static_tower_PMJPL_inputs import load_ECOv002_static_tower_PMJPL_inputs
+from .ECOv002_calval_PMJPL_inputs import load_ECOv002_calval_PMJPL_inputs
+from .verify import verify
 
 __author__ = 'Qiaozhen Mu, Maosheng Zhao, Steven W. Running, Gregory Halverson'
-
-logger = logging.getLogger(__name__)
-
-DEFAULT_WORKING_DIRECTORY = "."
-DEFAULT_MOD16_INTERMEDIATE = "MOD16_intermediate"
-
-DEFAULT_OUTPUT_VARIABLES = [
-    'LEi',
-    'LEc',
-    'LEs',
-    'LE',
-    'LE_daily',
-    'ET_daily_kg'
-]
-
-def PMJPL(
-        Rn: Union[Raster, np.ndarray],
-        G: Union[Raster, np.ndarray],
-        NDVI: Union[Raster, np.ndarray],
-        Ta_C: Union[Raster, np.ndarray],
-        Tmin_C: Union[Raster, np.ndarray],
-        RH: Union[Raster, np.ndarray],
-        IGBP: Union[Raster, np.ndarray],
-        Ps_Pa: Union[Raster, np.ndarray] = None,
-        elevation_m: Union[Raster, np.ndarray] = None,
-        delta_Pa: Union[Raster, np.ndarray] = None,
-        gamma_Jkg: Union[Raster, np.ndarray, float] = None) -> Dict[str, Raster]:
-    results = {}
-
-    # calculate leaf area index if it's not given
-    if LAI is None:
-        # calculate leaf area index from NDVI
-        LAI = LAI_from_NDVI(NDVI)
-
-    # calculate fraction of vegetation cover if it's not given
-    if FVC is None:
-        # calculate fraction of vegetation cover from NDVI
-        FVC = FVC_from_NDVI(NDVI)
-
-    # calculate surface air pressure if it's not given
-    if Ps_Pa is None:
-        # calculate surface air pressure is Pascal
-        Ps_Pa = calculate_surface_pressure(elevation_m=elevation_m, Ta_C=Ta_C)
-
-    # calculate Penman-Monteith/Priestley-Taylor delta term if it's not given
-    if delta_Pa is None:
-        # calculate Penman-Monteith/Priestley-Taylor delta term in Pascal per degree Celsius
-        delta_Pa = delta_Pa_from_Ta_C(Ta_C)
-
-    # calculate latent heat of vaporization if it's not given
-    if lambda_Jkg is None:
-        # calculate latent heat of vaporization in Joules per kilogram
-        lambda_Jkg = lambda_Jkg_from_Ta_C(Ta_C)
-
-    logger.info("calculating PM-MOD meteorology")
-
-    # calculate air temperature in Kelvin
-    Ta_K = celcius_to_kelvin(Ta_C)
-
-    # calculate saturation vapor pressure in Pascal from air temperature in Celsius
-    SVP_Pa = SVP_Pa_from_Ta_C(Ta_C)
-
-    # calculate vapor pressure in Pascal from releative humidity and saturation vapor pressure
-    Ea_Pa = RH * SVP_Pa
-
-    # specific humidity of air
-    # as a ratio of kilograms of water to kilograms of air and water
-    # from surface pressure and actual water vapor pressure
-    specific_humidity = calculate_specific_humidity(Ea_Pa, Ps_Pa)
-    results['specific_humidity'] = specific_humidity
-
-    # calculate air density (rho) in kilograms per cubic meter
-    rho_kgm3 = calculate_air_density(Ps_Pa, Ta_K, specific_humidity)
-    results["rho_kgm3"] = rho_kgm3
-
-    # calculate specific heat capacity of the air (Cp)
-    # in joules per kilogram per kelvin
-    # from specific heat of water vapor (CPW)
-    # and specific heat of dry air (CPD)
-    Cp_Jkg = calculate_specific_heat(specific_humidity)
-    results["Cp"] = Cp_Jkg    
-
-    # calculate delta term if it's not given
-    if delta_Pa is None:
-        # slope of saturation vapor pressure curve in Pascal per degree
-        delta_Pa = delta_Pa_from_Ta_C(Ta_C)
-
-    results["delta_Pa"] = delta_Pa
-
-    # calculate gamma term if it's not given
-    if gamma_Jkg is None:
-        # calculate psychrometric gamma in Joules per kilogram
-        gamma_Jkg = calculate_gamma(
-            Ta_C=Ta_C,
-            Ps_Pa=Ps_Pa,
-            Cp_Jkg=Cp_Jkg
-        )
-
-    # vapor pressure deficit in Pascal
-    VPD_Pa = rt.clip(SVP_Pa - Ea_Pa, 0.0, None)
-
-    # calculate relative surface wetness (fwet)
-    # from relative humidity
-    fwet = calculate_fwet(RH)
-    results['fwet'] = fwet
-
-    logger.info("calculating PM-MOD resistances")
-
-    # query leaf conductance to sensible heat (gl_sh) in seconds per meter
-    gl_sh = MOD16_parameter_from_IGBP(
-        variable="gl_sh", 
-        IGBP=IGBP
-    )
-
-    results['gl_sh'] = gl_sh
-
-    # calculate wet canopy resistance to sensible heat (rhc) in seconds per meter
-    # from leaf conductance to sensible heat (gl_sh), LAI, and relative surface wetness (fwet)
-    rhc = calculate_wet_canopy_resistance(gl_sh, LAI, fwet)
-    results['rhc'] = rhc
-
-    # calculate resistance to radiative heat transfer through air (rrc)
-    rrc = float32(rho_kgm3 * Cp_Jkg / (4.0 * SIGMA * Ta_K ** 3.0))
-    results['rrc'] = rrc
-
-    # calculate aerodynamic resistance (rhrc)
-    # in seconds per meter
-    # from wet canopy resistance to sensible heat
-    # and resistance to radiative heat transfer through air
-    rhrc = float32((rhc * rrc) / (rhc + rrc))
-    results['rhrc'] = rhrc
-
-    # calculate leaf conductance to evaporated water vapor (gl_e_wv)
-    gl_e_wv = MOD16_parameter_from_IGBP(
-        variable="gl_e_wv", 
-        IGBP=IGBP
-    )
-
-    results['gl_e_wv'] = gl_e_wv
-
-    rvc = calculate_wet_canopy_resistance(gl_e_wv, LAI, fwet)
-    results['rvc'] = rvc
-
-    # caluclate available radiation to the canopy (Ac)
-    # in watts per square meter
-    # this is the same as net radiation to the canopy in PT-JPL
-    Ac = Rn * FVC
-    results['Ac'] = Ac
-
-    # calculate wet latent heat flux (LEi)
-    # in watts per square meter
-    LEi = calculate_interception(
-        delta_Pa=delta_Pa, 
-        Ac=Ac, 
-        rho=rho_kgm3, 
-        Cp=Cp_Jkg, 
-        VPD_Pa=VPD_Pa, 
-        FVC=FVC, 
-        rhrc=rhrc, 
-        fwet=fwet, 
-        rvc=rvc,
-        gamma_Jkg=gamma_Jkg
-    )
-    
-    results['LEi'] = LEi
-
-    # calculate correctance factor (rcorr)
-    # for stomatal and cuticular conductances
-    # from surface pressure and air temperature
-    rcorr = calculate_rcorr(Ps_Pa, Ta_K)
-    results['rcorr'] = rcorr
-
-    # query biome-specific mean potential stomatal conductance per unit leaf area
-    CL = MOD16_parameter_from_IGBP(
-        variable="CL", 
-        IGBP=IGBP
-    )
-
-    results['CL'] = CL
-
-    # query open minimum temperature by land-cover
-    tmin_open = MOD16_parameter_from_IGBP(
-        variable="tmin_open", 
-        IGBP=IGBP
-    )
-
-    results['tmin_open'] = tmin_open
-
-    # query closed minimum temperature by land-cover
-    tmin_close = MOD16_parameter_from_IGBP(
-        variable="tmin_close", 
-        IGBP=IGBP
-    )
-
-    results['tmin_close'] = tmin_close
-
-    # calculate minimum temperature factor for stomatal conductance
-    mTmin = calculate_tmin_factor(Tmin_C, tmin_open, tmin_close, IGBP)
-    results['mTmin'] = mTmin
-
-    # query open vapor pressure deficit by land-cover
-    VPD_open = MOD16_parameter_from_IGBP(
-        variable="VPD_open", 
-        IGBP=IGBP
-    )
-
-    results['vpd_open'] = VPD_open
-
-    # query closed vapor pressure deficit by land-cover
-    VPD_close = MOD16_parameter_from_IGBP(
-        variable="VPD_close", 
-        IGBP=IGBP
-    )
-
-    results['vpd_close'] = VPD_close
-
-    # calculate vapor pressure deficit factor for stomatal conductance
-    mVPD = calculate_VPD_factor(VPD_open, VPD_close, VPD_Pa)
-    results['mVPD'] = mVPD
-
-    # calculate stomatal conductance (gs1)
-    gs1 = CL * mTmin * mVPD * rcorr
-    results['gs1'] = gs1
-
-    # correct cuticular conductance constant to leaf cuticular conductance (Gcu) using correction factor (rcorr)
-    Gcu = CUTICULAR_CONDUCTANCE * rcorr
-    results['Gcu'] = Gcu
-
-    # calculate canopy conductance
-    # equivalent to g_canopy
-    Cc = calculate_canopy_conductance(LAI, fwet, gl_sh, gs1, Gcu)
-    results['Cc'] = Cc
-
-    # calculate surface resistance to evapotranspiration as inverse of canopy conductance (Cc)
-    rs = rt.clip(1.0 / Cc, 0.0, MAX_RESISTANCE)
-    results['rs'] = rs
-
-    # calculate convective heat transfer as inverse of leaf conductance to sensible heat (gl_sh)
-    rh = 1.0 / gl_sh
-    results['rh'] = rs
-
-    # calculate radiative heat transfer (rr)
-    rr = rho_kgm3 * Cp_Jkg / (4.0 * SIGMA * Ta_K ** 3)
-    results['rr'] = rr
-
-    # calculate parallel resistance (ra)
-    # MOD16 user guide is not clear about what to call this
-    ra = (rh * rr) / (rh + rr)
-    results["ra"] = ra
-
-    # calculate transpiration
-    LEc = calculate_transpiration(
-        delta_Pa=delta_Pa, 
-        Ac=Ac, 
-        rho=rho_kgm3, 
-        Cp_Jkg=Cp_Jkg, 
-        VPD_Pa=VPD_Pa, 
-        FVC=FVC, 
-        ra=ra, 
-        fwet=fwet, 
-        rs=rs,
-        gamma_Jkg=gamma_Jkg
-    )
-    
-    results['LEc'] = LEc
-
-    # soil evaporation
-
-    # query aerodynamic resistant constraints from land-cover
-    rbl_max = MOD16_parameter_from_IGBP(
-        variable="rbl_max", 
-        IGBP=IGBP
-    )
-
-    results['rbl_max'] = rbl_max
-
-    rbl_min = MOD16_parameter_from_IGBP(
-        variable="rbl_min", 
-        IGBP=IGBP
-    )
-
-    results['rbl_min'] = rbl_min
-
-    rtotc = calculate_rtotc(VPD_Pa, VPD_open, VPD_close, rbl_max, rbl_min)
-    results['rtotc'] = rtotc
-
-    # calculate total aerodynamic resistance
-    # by applying correction to total canopy resistance
-    rtot = rcorr * rtotc
-    results['rtot'] = rtot
-
-    # calculate resistance to radiative heat transfer through air
-    rrs = float32(rho_kgm3 * Cp_Jkg / (4.0 * SIGMA * Ta_K ** 3))
-    results['rrs'] = rrs
-
-    # calculate aerodynamic resistance at the soil surface
-    ras = (rtot * rrs) / (rtot + rrs)
-    results['ras'] = ras
-
-    # calculate available radiation at the soil
-    Asoil = rt.clip((1.0 - FVC) * Rn - G, 0.0, None)
-    results['Asoil'] = Asoil
-
-    # separate wet soil evaporation and potential soil evaporation
-
-    # calculate wet soil evaporation
-    LE_soil_wet = calculate_wet_soil_evaporation(
-        delta_Pa=delta_Pa, 
-        Asoil=Asoil, 
-        rho=rho_kgm3, 
-        Cp_Jkg=Cp_Jkg, 
-        FVC=FVC, 
-        VPD_Pa=VPD_Pa, 
-        ras=ras, 
-        fwet=fwet, 
-        rtot=rtot,
-        gamma_Jkg=gamma_Jkg
-    )
-    
-    results['LE_soil_wet'] = LE_soil_wet
-
-    # calculate potential soil evaporation
-    LE_soil_pot = calculate_potential_soil_evaporation(
-        delta_Pa=delta_Pa, 
-        Asoil=Asoil, 
-        rho=rho_kgm3, 
-        Cp_Jkg=Cp_Jkg, 
-        FVC=FVC, 
-        VPD_Pa=VPD_Pa, 
-        ras=ras, 
-        fwet=fwet, 
-        rtot=rtot,
-        gamma_Jkg=gamma_Jkg
-    )
-    
-    results['LE_soil_pot'] = LE_soil_pot
-
-    # calculate soil moisture constraint
-    fSM = calculate_fSM(RH, VPD_Pa)
-    results['fSM'] = fSM
-
-    # calculate soil evaporation
-    LEs = rt.clip(LE_soil_wet + LE_soil_pot * fSM, 0.0, None)
-
-    # fill soil evaporation with zero
-    LEs = rt.where(isnan(LEs), 0.0, LEs)
-    results['LEs'] = LEs
-
-    # sum partitions into total latent heat flux
-    LE = rt.clip(LEi + LEc + LEs, 0.0, Rn)
-    results['LE'] = LE
-
-    return results

PMJPL/{parameters.py → PMJPL_parameter_from_IGBP.py}

@@ -3,13 +3,20 @@ from os.path import join, abspath, dirname
 import numpy as np
 import pandas as pd
 
-from rasters import Raster, RasterGeometry
+from rasters import Raster, RasterGeometry, SpatialGeometry
 
-from .MCD12C1.MCD12C1 import load_MCD12C1_IGBP
+from MCD12C1_2019_v006 import load_MCD12C1_IGBP
+
+from .constants import *
 
 LUT = pd.read_csv(join(abspath(dirname(__file__)), 'mod16.csv'))
 
-def MOD16_parameter_from_IGBP(variable: str, IGBP: Union[Raster, np.ndarray] = None, geometry: RasterGeometry = None) -> Union[Raster, np.ndarray]:
+def PMJPL_parameter_from_IGBP(
+        variable: str, 
+        IGBP: Union[Raster, np.ndarray] = None, 
+        geometry: SpatialGeometry = None,
+        IGBP_upsampling_resolution_meters: int = IGBP_UPSAMPLING_RESOLUTION_METERS,
+        resampling=IGBP_PARAMETER_RESAMPLING) -> Union[Raster, np.ndarray]:
     """
     Translates the IGBP (International Geosphere-Biosphere Programme) values to the corresponding values in the Look-Up Table (LUT) for a given variable.
 
@@ -28,14 +35,22 @@ def MOD16_parameter_from_IGBP(variable: str, IGBP: Union[Raster, np.ndarray] = N
 
     Returns:
         Union[np.ndarray, Raster]: The translated values.
-
     """
+    if geometry is None and isinstance(IGBP, Raster):
+        geometry = IGBP.geometry
+
+    if isinstance(geometry, RasterGeometry):
+        IGBP_geometry = geometry.UTM(IGBP_upsampling_resolution_meters)
+    else:
+        IGBP_geometry = geometry
+
     if IGBP is None:
-        IGBP = load_MCD12C1_IGBP(geometry=geometry)
+        IGBP = load_MCD12C1_IGBP(geometry=IGBP_geometry)
 
     result = np.float32(np.array(LUT[variable])[np.array(IGBP).astype(int)])
 
-    if isinstance(IGBP, Raster):
-        result = Raster(result, geometry=IGBP.geometry)
+    if isinstance(IGBP, Raster) or isinstance(geometry, RasterGeometry):
+        result = Raster(result, geometry=IGBP_geometry)
+        result = result.to_geometry(geometry, resampling=resampling)
     
     return result

PMJPL/VPD_factor.py

@@ -20,7 +20,7 @@ def calculate_VPD_factor(
     """
     # calculate VPD factor using queried open and closed VPD
     mVPD = rt.where(VPD <= VPD_open, 1.0, np.nan)
-    mVPD = rt.where(rt.logical_and(VPD_open < VPD, VPD < VPD_close), (VPD_close - VPD) / (VPD_close - VPD_open), mVPD)
+    mVPD = rt.where(np.logical_and(VPD_open < VPD, VPD < VPD_close), (VPD_close - VPD) / (VPD_close - VPD_open), mVPD)
     mVPD = rt.where(VPD >= VPD_close, 0.0, mVPD)
 
     return mVPD