agwx_biophys 0.0.1

checksums.yaml

@@ -0,0 +1,7 @@
+---
+SHA1:
+  metadata.gz: f9f340d516a0d61d5d391d7f346540c6af7e38e3
+  data.tar.gz: 29dd22ef00cc3ec68b4a9aff760b7324fdb624a2
+SHA512:
+  metadata.gz: 991e2c2afbeb889135e0c71344c2fdd6de0a53b339e969d5984c7f89d223b5ce41f8177c97afec6354e0b84adf8ce5b95a3ae9a5e579bf86ad5773ad3dbf6efc
+  data.tar.gz: 8abb2837d359283f91762f624c31f63a1b94d7d559ee524b894cbda52db001b8a95ebe1199c29f353869d03029a0a82bf31b51490073bbe2550fa12999ceab85

data/.gitignore

@@ -0,0 +1,18 @@
+*.gem
+*.rbc
+.bundle
+.config
+.yardoc
+Gemfile.lock
+InstalledFiles
+_yardoc
+coverage
+doc/
+lib/bundler/man
+pkg
+rdoc
+spec/reports
+test/tmp
+test/version_tmp
+test/.~lock.TestValues.xls#
+tmp

data/Gemfile

@@ -0,0 +1,4 @@
+source 'https://rubygems.org'
+
+# Specify your gem's dependencies in agwx_biophys.gemspec
+gemspec

data/LICENSE.txt

@@ -0,0 +1,22 @@
+Copyright (c) 2013 RickWayne
+
+MIT License
+
+Permission is hereby granted, free of charge, to any person obtaining
+a copy of this software and associated documentation files (the
+"Software"), to deal in the Software without restriction, including
+without limitation the rights to use, copy, modify, merge, publish,
+distribute, sublicense, and/or sell copies of the Software, and to
+permit persons to whom the Software is furnished to do so, subject to
+the following conditions:
+
+The above copyright notice and this permission notice shall be
+included in all copies or substantial portions of the Software.
+
+THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND,
+EXPRESS OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF
+MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND
+NONINFRINGEMENT. IN NO EVENT SHALL THE AUTHORS OR COPYRIGHT HOLDERS BE
+LIABLE FOR ANY CLAIM, DAMAGES OR OTHER LIABILITY, WHETHER IN AN ACTION
+OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, OUT OF OR IN CONNECTION
+WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE SOFTWARE.

data/README.md

@@ -0,0 +1,29 @@
+# AgwxBiophys
+
+TODO: Write a gem description
+
+## Installation
+
+Add this line to your application's Gemfile:
+
+    gem 'agwx_biophys'
+
+And then execute:
+
+    $ bundle
+
+Or install it yourself as:
+
+    $ gem install agwx_biophys
+
+## Usage
+
+TODO: Write usage instructions here
+
+## Contributing
+
+1. Fork it
+2. Create your feature branch (`git checkout -b my-new-feature`)
+3. Commit your changes (`git commit -am 'Add some feature'`)
+4. Push to the branch (`git push origin my-new-feature`)
+5. Create new Pull Request

data/Rakefile

@@ -0,0 +1,9 @@
+require "bundler/gem_tasks"
+require 'rake/testtask'
+
+Rake::TestTask.new do |t|
+  t.libs << 'test'
+end
+
+desc "Run tests"
+task default: :test

data/agwx_biophys.gemspec

@@ -0,0 +1,23 @@
+# coding: utf-8
+lib = File.expand_path('../lib', __FILE__)
+$LOAD_PATH.unshift(lib) unless $LOAD_PATH.include?(lib)
+require 'agwx_biophys/version'
+
+Gem::Specification.new do |spec|
+  spec.name          = "agwx_biophys"
+  spec.version       = AgwxBiophys::VERSION
+  spec.authors       = ["RickWayne"]
+  spec.email         = ["fewayne@wisc.edu"]
+  spec.description   = %q{Biophysical calculators primarily useful in agriculture and land management}
+  spec.summary       = %q{So far, just an ET calculator (Priestly-Taylor)}
+  spec.homepage      = ""
+  spec.license       = "MIT"
+
+  spec.files         = `git ls-files`.split($/)
+  spec.executables   = spec.files.grep(%r{^bin/}) { |f| File.basename(f) }
+  spec.test_files    = spec.files.grep(%r{^(test|spec|features)/})
+  spec.require_paths = ["lib"]
+
+  spec.add_development_dependency "bundler", "~> 1.3"
+  spec.add_development_dependency "rake"
+end

data/lib/agwx_biophys/et.rb

@@ -0,0 +1,149 @@
+module AgwxBiophys
+  module ET
+    SOLCON    = 1367.0
+    STEFAN    = (0.0864*0.0000000567)
+    SFCEMISS  = 0.96
+    ALBEDO    = 0.25
+    KELVIN    = 273.15
+    
+    def rad_lat(lat); lat * Math::PI / 180.0;end
+
+    def daily_to_sol(total_sol); 0.0864 * total_sol; end
+
+    def declin(doy)
+      0.41 * 
+        Math::cos(
+          2 * Math::PI * ( doy - 172.0 ) / 365.0
+        )
+    end
+
+    def sunrise_angle(doy,lat)
+      Math.acos(
+        -1 *
+        Math.tan(declin(doy)) *
+        Math.tan(rad_lat(lat))
+      )
+    end
+
+    def sunrise_hour(doy,lat); 12 - ( 12 / Math::PI ) * sunrise_angle(doy,lat); end
+
+    def day_hours(doy,lat); 24 - 2 * sunrise_hour(doy,lat); end
+
+    def av_eir(doy)
+      SOLCON * ( 1 + 0.035 * Math.cos( 2 * Math::PI * doy / 365.0 ) )
+    end
+    
+    def to_eir(doy,lat)
+      ( 0.0864 / Math::PI ) * av_eir(doy) *
+
+      (
+        sunrise_angle(doy,lat) *
+         Math.sin( declin(doy) ) *
+         Math.sin( rad_lat(lat) ) +
+         
+        Math.cos( declin(doy) ) *
+         Math.cos( rad_lat(lat) ) *
+         Math.sin( sunrise_angle(doy,lat) )
+      )
+    end
+    
+    def to_clr(doy,lat)
+      to_eir(doy,lat) * (-0.7 + 0.86*day_hours(doy,lat)) / day_hours(doy,lat)
+    end
+
+    def c_to_k(temperature); temperature + KELVIN; end
+
+    # Temps are in Deg C
+    
+    def lwu(min_temp,max_temp)
+      avg = (min_temp+max_temp) / 2.0
+      SFCEMISS * STEFAN * (273.15+avg) ** 4
+    end
+    
+    def sfactor(min_temp,max_temp)
+      avg = (min_temp+max_temp) / 2.0
+      0.398 + 0.0171 * avg - 0.000142 * avg * avg
+    end
+    
+    def sky_emiss(avg_v_press,min_temp,max_temp,avg_temp)
+
+      if ( avg_v_press > 0.5)
+        0.7 + (5.95e-4) * avg_v_press * Math.exp( 1500/(273+avg_temp) )
+      else
+        avg_t = (min_temp+max_temp) / 2.0
+        (1 - 0.261*Math.exp(-0.000777*avg_t*avg_t))
+      end
+    end
+
+    def angstrom(avg_v_press,min_temp,max_temp,avg_temp)
+      1.0 - sky_emiss(avg_v_press,min_temp,max_temp,avg_temp) / SFCEMISS
+    end
+
+    def clr_ratio(d_to_sol,doy,lat)
+      dts = daily_to_sol(d_to_sol)
+      tc = to_clr(doy,lat)
+      # Never return higher than 1
+      [dts/tc,1.0].min
+    end
+    
+    def lwnet(avg_v_press,min_temp,max_temp,avg_temp,d_to_sol,doy,lat)
+      angstrom(avg_v_press,min_temp,max_temp,avg_temp) * lwu(min_temp,max_temp)  * clr_ratio(d_to_sol,doy,lat)
+    end
+    
+    def et(max_temp,min_temp,avg_temp,avg_v_press,d_to_sol,doy,lat)
+      lwnet = lwnet(avg_v_press,min_temp,max_temp,avg_temp,d_to_sol,doy,lat  )
+      lwd = lwu(min_temp,max_temp) - lwnet
+      r_net_1 = (1.0 - ALBEDO) * daily_to_sol(d_to_sol) - lwnet
+      ret1 = 1.28 * sfactor(min_temp,max_temp) * r_net_1
+      # puts "lwnet #{lwnet}, lwd #{lwd}, r_net_1 #{r_net_1}, ret1 #{ret1} refet #{ret1 / 62.3}"
+      [ret1 / 62.3,clr_ratio(d_to_sol,doy,lat) * 100.0]
+    end
+#     double LWD = LWU - LWnet;
+#     double Rnet1 = ( 1-ALBEDO ) * DailyToSol - LWnet; /*  ;(RESULTANT ET)=MJ/M2/DY
+#                                                ;(RESULTANT ET)=REFN ET IN/DY */
+#     double RET1 = 1.28 * Sfactor * Rnet1;
+#     DRefET = RET1 / 62.3;
+#     DClrRatio = ClrRatio*100;
+
+    # Here's the old Java code
+    # double rad_lat =       lat *Math.PI/180.0;
+    # double DailyToSol =    0.0864 * DToSol;
+    # double Declin =    0.41 *Math.cos( 2 * Math.PI * ( dayOfYear -172. ) / 365. );
+    # double Sunrise_Angle = Math.acos( -1*Math.tan(Declin) * Math.tan(rad_lat));
+    # double Sunrise_Hour = 12 - ( 12 / Math.PI ) * Sunrise_Angle;
+    # double Day_Hours = 24 - 2 * Sunrise_Hour;
+
+    # double AvEIR = SOLCON * ( 1 + 0.035 * Math.cos( 2 * Math.PI * dayOfYear / 365. ) );
+    # double ToEIR = ( 0.0864 / Math.PI ) * AvEIR *
+    #     ( Sunrise_Angle * Math.sin( Declin ) * Math.sin( rad_lat )
+    #     + Math.cos( Declin ) * Math.cos( rad_lat ) * Math.sin( Sunrise_Angle ) );
+    #     /* MJ/SQM/DY */
+    # 
+    # double ToClr = ToEIR * (-.7 + 0.86*Day_Hours) / Day_Hours ;
+    # 
+    # double T = ( DMxTAir + DMnTAir ) / 2;
+    # double LWU = SFCEMISS * STEFAN * Math.pow( 273.15+T,4  ); /*LONG WAVE RAD
+    #                          EMITTED OUT */
+    # double Sfactor = 0.398 + 0.0171 * T - 0.000142 * T * T;
+    # 
+    # double Angstrom = 1 - SkyEmiss / SFCEMISS;
+# 
+#     double ClrRatio;
+#     if (DailyToSol <= ToClr) {          /*CLRRATIO = PCTCLR, < 100%
+#                                               100 added to relieve problem */
+#         ClrRatio =  DailyToSol / ToClr;
+#     }
+#     else {
+#         ClrRatio = 1.0;
+#     }
+# 
+#     double LWnet = Angstrom * LWU  * ClrRatio;
+#     double LWD = LWU - LWnet;
+#     double Rnet1 = ( 1-ALBEDO ) * DailyToSol - LWnet; /*  ;(RESULTANT ET)=MJ/M2/DY
+#                                                ;(RESULTANT ET)=REFN ET IN/DY */
+#     double RET1 = 1.28 * Sfactor * Rnet1;
+#     DRefET = RET1 / 62.3;
+#     DClrRatio = ClrRatio*100;
+    
+  end # module ET
+end # AgwxBiophys

data/lib/agwx_biophys/version.rb

@@ -0,0 +1,3 @@
+module AgwxBiophys
+  VERSION = "0.0.1"
+end

data/lib/agwx_biophys.rb

@@ -0,0 +1,6 @@
+require "agwx_biophys/version"
+require "agwx_biophys/et"
+
+
+module AgwxBiophys
+end

data/test/Et_testing.java

@@ -0,0 +1,117 @@
+public class Et_testing {
+
+static final double SOLCON     = 1367.0;
+static final double STEFAN     = (0.0864*0.0000000567);
+static final double SFCEMISS   = 0.96;
+static final double ALBEDO     = 0.25;
+
+  double DRefET, DClrRatio;
+  
+  public Et_testing() {
+   DRefET = -1.0;
+   DClrRatio = -1.0;
+  }
+  public Et_testing(double DToSol, int dayOfYear, double DMnTAir, double DMxTAir, 
+            double DAvVPre, double DAvTAir, double lat) {
+    this();
+    if (lat < 1.0 || lat > 90.0) /* Dont use this in Indonesia!!! */
+       return;
+    // wha***k?
+    /*
+    if ( ( dayOfYear >= 91 ) && ( dayOfYear <= 273 ) )
+         return;
+    */
+      
+    double rad_lat =       lat *Math.PI/180.0;
+    double DailyToSol =    0.0864 * DToSol;
+    double Declin =    0.41 *Math.cos( 2 * Math.PI * ( dayOfYear -172. ) / 365. );
+
+/*-----
+ * The Old sunrise angle was given by
+ *   Sunrise_Angle=acos(tan(Declin)*tan(rad_lat);formula used for ET days
+ *       + 0.0145/(cos(Declin)*cos(rad_lat)))   ;other than months ?? (4-9)
+ */
+
+    double Sunrise_Angle = Math.acos( -1*Math.tan(Declin) * Math.tan(rad_lat));
+    double Sunrise_Hour = 12 - ( 12 / Math.PI ) * Sunrise_Angle;
+    double Day_Hours = 24 - 2 * Sunrise_Hour;
+    double AvEIR = SOLCON * ( 1 + 0.035 * Math.cos( 2 * Math.PI * dayOfYear / 365. ) );
+    double ToEIR = ( 0.0864 / Math.PI ) * AvEIR *
+        ( Sunrise_Angle * Math.sin( Declin ) * Math.sin( rad_lat )
+        + Math.cos( Declin ) * Math.cos( rad_lat ) * Math.sin( Sunrise_Angle ) );
+        /* MJ/SQM/DY */
+
+    double ToClr = ToEIR * (-.7 + 0.86*Day_Hours) / Day_Hours ;
+
+    double T = ( DMxTAir + DMnTAir ) / 2;
+    double LWU = SFCEMISS * STEFAN * Math.pow( 273.15+T,4  ); /*LONG WAVE RAD
+                             EMITTED OUT */
+    double Sfactor = 0.398 + 0.0171 * T - 0.000142 * T * T;
+/*-------
+ * Old Sky Emiss equation
+ * SkyEmiss = ( 1 - 0.261 * Exp ( -0.000777 * T * T ) )
+ */ double SkyEmiss;
+    if ( DAvVPre > 0.5) {
+        SkyEmiss = 0.7 + (5.95e-4)*DAvVPre*
+            Math.exp(1500/(273+DAvTAir));
+    }
+    else {
+            SkyEmiss = (1-0.261*Math.exp(-0.000777*T*T));
+    }
+    double Angstrom = 1 - SkyEmiss / SFCEMISS;
+
+    double ClrRatio;
+    if (DailyToSol <= ToClr) {          /*CLRRATIO = PCTCLR, < 100%
+                                              100 added to relieve problem */
+        ClrRatio =  DailyToSol / ToClr;
+    }
+    else {
+        ClrRatio = 1.0;
+    }
+
+    double LWnet = Angstrom * LWU  * ClrRatio;
+    double LWD = LWU - LWnet;
+    double Rnet1 = ( 1-ALBEDO ) * DailyToSol - LWnet; /*  ;(RESULTANT ET)=MJ/M2/DY
+                                               ;(RESULTANT ET)=REFN ET IN/DY */
+    double RET1 = 1.28 * Sfactor * Rnet1;
+    System.out.print("ToEIR: "); System.out.println(ToEIR);
+    System.out.print("AvEIR: "); System.out.println(AvEIR);
+    System.out.print("SkyEmiss: "); System.out.println(SkyEmiss);
+    System.out.print("LWnet: "); System.out.println(LWnet);
+    System.out.print("Angstrom: "); System.out.println(Angstrom);
+    System.out.print("LWU: "); System.out.println(LWU);
+    System.out.print("ClrRatio: "); System.out.println(ClrRatio);
+    System.out.print("LWD: "); System.out.println(LWD);
+    System.out.print("Rnet1: "); System.out.println(Rnet1);
+    System.out.print("RET1: "); System.out.println(RET1);
+    DRefET = RET1 / 62.3;
+    DClrRatio = ClrRatio*100;
+  } // ET ctor
+  
+  public double getET() {
+   return DRefET;
+  }
+  public double getPctClr() {
+   return DClrRatio;
+  }
+
+  public static void main(String[] args) {
+    // public Et_testing(double DToSol, int dayOfYear, double DMnTAir, double DMxTAir, 
+    //           double DAvVPre, double DAvTAir, float lat) {
+    Et_testing eT = new Et_testing(336.0,172,16.0,24.0,1.0,(16.0+24.0)/2.0,45.0);
+    System.out.println(eT.getET());
+    System.out.println();
+    System.out.println();
+    System.out.println();
+    // double[] et_obs1 = {244.40,10.80,24.88,1.49,152,44.12,0.18,63.05};
+    eT = new Et_testing(285.40,153,5.23,15.91,0.86,10.60,44.12);
+    System.out.println(eT.getET());
+    System.out.println(eT.getET() - 0.18);
+    double[] et_obs2 = {285.40,5.23,15.91,0.86,153,44.12,0.17,73.49};
+    double[] et_obs3 = {319.60,3.63,21.95,0.92,154,44.12,0.20,82.15};
+    double[] et_obs4 = {45.76,9.95,13.56,1.25,155,44.12,0.03,11.74};
+    
+    eT.invokedStandalone = true;
+  }
+  private boolean invokedStandalone = false;
+} 

data/test/test_et.rb

@@ -0,0 +1,232 @@
+require 'test/unit'
+require 'agwx_biophys'
+
+class TestET < Test::Unit::TestCase
+  include AgwxBiophys::ET
+  def test_rad_lat
+    [
+      [0,0],
+      [30,0.5235983333],
+      [45,0.7853975],
+      [60,1.0471966667],
+      [90,1.570795],
+      [180,3.14159]
+    ].each do |lat,radians|
+      assert_in_delta(radians, rad_lat(lat), 0.00001)
+    end
+    assert_equal(0, rad_lat(0.0))
+    assert_equal(Math::PI, rad_lat(180.0))
+    assert_equal(Math::PI/2.0, rad_lat(90.0))
+  end
+  
+  def test_daily_to_sol
+    assert_equal(0.0864, daily_to_sol(1.0))
+  end
+  
+  def test_declin
+    [
+      [1,-0.4019921553],
+      [150,0.3809480272],
+      [172,0.41],
+      [200,0.3632890322],
+      [250,0.092707884]
+    ].each do |doy,declination|
+      assert_in_delta(declination, declin(doy), 2 ** -20)
+    end
+  end
+  
+  def test_sunrise_angle
+    [
+      [1,0,1.5707963268],
+      [150,30,1.8041439502],
+      [172,45,2.020424427],
+      [200,60,2.2895663737],
+      [1,60,0.7431090864],
+      [150,45,1.9828717405],
+      [172,30,1.8244415261],
+      [200,0,1.5707963268]
+      
+    ].each do |doy,latitude,sun_angle|
+      assert_in_delta(sun_angle, sunrise_angle(doy,latitude), 0.00001)
+    end
+  end
+  
+  def test_sunrise_hour
+    [
+      [1,0,5.999994932],
+      [150,30,5.1086719137],
+      [172,45,4.2825406486],
+      [200,60,3.254493271]
+    ].each do |doy,latitude,sun_hour|
+      assert_in_delta(sun_hour, sunrise_hour(doy,latitude), 0.00001)
+    end
+  end
+
+  def test_day_hours
+    [
+      [1,0,12.0],
+      [150,30,13.7826561727],
+      [172,45,15.4349187028],
+      [200,60,17.491013458]
+    ].each do |doy,latitude,day_hrs|
+      assert_in_delta(day_hrs, day_hours(doy,latitude), 0.00002)
+    end
+  end
+  def test_av_eir
+    [
+      [1,1414.8379112589],
+      [150,1326.4494599268],
+      [172,1319.9344503498],
+      [200,1321.3095743309]
+    ].each do |doy,aveir|
+      assert_in_delta(aveir, av_eir(doy), 0.00006)
+    end
+  end
+  
+  def test_to_eir
+    [
+      [1,0,35.8090258247],
+      [150,30,40.7680366801],
+      [172,45,41.8737911],
+      [200,60,38.3861488725]
+    ].each do |doy,latitude,total_eir|
+      assert_in_delta(total_eir, to_eir(doy,latitude), 0.00005)
+    end
+  end
+  
+  def test_to_clr
+    [
+      [1,0,28.7069041338],
+      [150,30,32.9899653948],
+      [172,45,34.1124418882],
+      [200,60,31.4758531933]
+    ].each { |doy,lat,clr| assert_in_delta(clr, to_clr(doy,lat), 0.0003) }
+  end
+  
+  def test_lwu
+    [
+      [5,15,30.2297320891],
+      [16,24,34.73182695],
+      [18,27,35.9318491504],
+      [16,24,34.7318269495]
+    ].each do |min,max,exp_lwu|
+      assert_in_delta(exp_lwu, lwu(min,max), 0.000001)
+    end
+  end
+  
+  def test_sfactor
+    [
+      [5,15,0.5548],
+      [16,24,0.6832],
+      [18,27,0.7108625]
+    ].each do |min,max,exp_sfactor|
+      assert_in_delta(exp_sfactor, sfactor(min,max), 2 ** -20)
+    end
+  end
+  
+  def test_exponentiation
+    assert_in_delta(2.71828182845904, Math.exp(1), 2 ** -20)
+    assert_in_delta(148.4131591026, Math.exp(5), 2 ** -20)
+    [
+      [5 ,15,10  ,283.0, 0.9252419575,200.4076128756],
+      [16,24,20  ,293.0, 0.7328604333,167.2440168899],
+      [18,27,22.5,295.5, 0.6747883114,160.1550057145]
+    ].each do |min,max,avg,avg_kelvin,lo,hi|
+      avg_temp = (min+max) / 2.0
+      assert_in_delta(avg, avg_temp, 2 ** -20)
+      assert_in_delta(avg_temp + 273.0, avg_kelvin, 2 ** -20)
+      avg_kelvin = avg_temp + 273
+      lo_exp = Math.exp(-0.000777*avg_temp*avg_temp)
+      hi_exp = Math.exp( 1500.0/(273.0+avg_temp))
+      assert_in_delta(hi, hi_exp, 2 ** -20)
+      assert_in_delta(lo, lo_exp, 2 ** -20)
+    end
+  end
+  
+  def test_sky_emiss
+    [
+      [0.1,5,15     ,0.7585118491],
+      [0.2,16,24    ,0.8087234269],
+      [0.499,18,27  ,0.8238802507],
+      [0.51,5,15    ,0.7608136901],
+      [1,16,24      ,0.79951019],
+      [2,18,27      ,0.8905844568],
+      [1.0,16.0,24.0,0.79951019]
+    ].each do |vp,min,max,sky|
+      assert_in_delta(sky, sky_emiss(vp,min,max,(min+max)/2.0), 0.000001,"Out of range for vp #{vp} min #{min} max #{max}")
+    end  
+  end
+  
+  def test_angstrom
+    [
+      [0.1,5,15     ,0.2098834905],
+      [0.2,16,24    ,0.1575797636],
+      [0.499,18,27  ,0.1417914055],
+      [0.51,5,15    ,0.2074857395],
+      [1,16,24      ,0.1671768854],
+      [2,18,27      ,0.0723078575],
+      [1.0,16.0,24.0,0.16718]
+    ].each { |vp,min,max,angs| assert_in_delta(angs, angstrom(vp,min,max,(min+max)/2.0), 0.000004) }
+    
+  end
+  def test_clr_ratio
+    [
+      [62,1,0,0.186603194],
+      [142,150,30,0.3718949036],
+      [336,172,45,0.8510208708],
+      [62,200,60,0.1701876028]
+    ].each { |dAvSol,doy,lat,clr_ratio| assert_in_delta(clr_ratio, clr_ratio(dAvSol,doy,lat), 2 ** -20) }
+  end
+  
+  def test_lwnet
+    [
+      [0.1,5,15,142,150,30         ,2.35957],
+      [0.2,16,24,336,172,45        ,4.65767],
+      [0.499,18,27,62,200,60       ,0.86708],
+      [1.0,16.0,24.0,336.0,172,45.0,4.94133]
+    ].each do |vp,min,max,sol,doy,lat,lwnet|
+      assert_in_delta(lwnet, lwnet(vp,min,max,(min+max)/2.0,sol,doy,lat), 0.000025,"lwnet for #{vp}")
+    end
+  end
+
+  def test_ets_with_obs
+    # DAvSol,DMnTAir,DMxTAir,DAvTAir,DAvVPre,theDate,latitude,DRefET,DPctClr
+  
+    et_obs = [ 
+      [336.0,16.0,24.0,20.0,1.0,'2013-06-21',45.0,0.24,85.10],
+      [244.40,10.80,24.88,18.32,1.49,'2013-06-01',44.12,0.18,63.05],
+      [285.40,5.23,15.91,10.60,0.86,'2013-06-02',44.12,0.17,73.49],
+      [319.60,3.63,21.95,13.82,0.92,'2013-06-03',44.12,0.20,82.15],
+      [45.76,9.95,13.56,11.41,1.25,'2013-06-04',44.12,0.03,11.74],
+      [97.60,9.88,15.20,11.28,1.27,'2013-06-05',44.12,0.07,25.01],
+      [62.95,9.55,13.01,11.00,1.26,'2013-06-06',44.12,0.04,16.10],
+      [176.90,8.48,19.95,13.43,1.28,'2013-06-07',44.12,0.12,45.20],
+      [131.80,9.21,19.31,14.57,1.35,'2013-06-08',44.12,0.09,33.63],
+      [158.40,11.08,23.55,17.22,1.51,'2013-06-09',44.12,0.12,40.37],
+      [234.00,14.20,26.73,20.02,1.78,'2013-06-10',44.12,0.19,59.58],
+      [220.00,12.74,26.06,20.13,1.84,'2013-06-11',44.12,0.18,55.96],
+      [142.00,13.74,25.72,18.78,1.86,'2013-06-12',44.12,0.12,36.09],
+      [329.10,13.61,24.73,19.16,1.46,'2013-06-13',44.12,0.25,83.57],
+      [336.40,11.01,24.85,18.69,1.34,'2013-06-14',44.12,0.24,85.37],
+      [190.20,14.12,26.06,18.89,1.93,'2013-06-15',44.12,0.16,48.24],
+      [349.00,13.27,27.64,20.92,1.52,'2013-06-16',44.12,0.27,88.48],
+      [297.40,13.74,27.97,20.09,1.55,'2013-06-17',44.12,0.23,75.37],
+      [345.90,9.55,22.27,15.83,1.14,'2013-06-18',44.12,0.24,87.64],
+      [348.00,6.63,25.58,17.61,1.24,'2013-06-19',44.12,0.24,88.16],
+      [331.50,14.34,29.81,22.59,1.52,'2013-06-20',44.12,0.26,83.98],
+      [130.30,16.84,25.53,20.66,2.06,'2013-06-21',44.12,0.11,33.01],
+      [134.30,18.37,25.59,20.98,2.23,'2013-06-22',44.12,0.12,34.03],
+      [297.30,17.58,29.69,23.44,2.24,'2013-06-23',44.12,0.26,75.35],
+      [198.10,19.17,28.10,22.58,2.27,'2013-06-24',44.12,0.18,50.23],
+      [195.10,17.85,25.98,21.67,2.30,'2013-06-25',44.12,0.17,49.49],
+      [235.70,20.29,28.08,23.26,2.37,'2013-06-26',44.12,0.22,59.82],
+    ].each { |arr| arr[5] = Date.parse(arr[5]).yday }
+      
+    et_obs.each do |sol,min,max,avg,vapr,yday,lat,et,pct|
+      # max_temp,min_temp,avg_v_press,d_to_sol,doy,lat
+      ref_et,pct_clr = et(max,min,avg,vapr,sol,yday,lat)
+      assert_in_delta(et, ref_et.round(2), 0.00001,"ET bludgered for #{Date.civil(2013,1,1)+yday-1}")
+      assert_in_delta(pct, pct_clr.round(2), 0.00001,"PctClr bludgered for #{Date.civil(2013,1,1)+yday-1}")
+    end
+  end
+end

metadata

@@ -0,0 +1,87 @@
+--- !ruby/object:Gem::Specification
+name: agwx_biophys
+version: !ruby/object:Gem::Version
+  version: 0.0.1
+platform: ruby
+authors:
+- RickWayne
+autorequire: 
+bindir: bin
+cert_chain: []
+date: 2013-12-04 00:00:00.000000000 Z
+dependencies:
+- !ruby/object:Gem::Dependency
+  name: bundler
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - ~>
+      - !ruby/object:Gem::Version
+        version: '1.3'
+  type: :development
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - ~>
+      - !ruby/object:Gem::Version
+        version: '1.3'
+- !ruby/object:Gem::Dependency
+  name: rake
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - '>='
+      - !ruby/object:Gem::Version
+        version: '0'
+  type: :development
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - '>='
+      - !ruby/object:Gem::Version
+        version: '0'
+description: Biophysical calculators primarily useful in agriculture and land management
+email:
+- fewayne@wisc.edu
+executables: []
+extensions: []
+extra_rdoc_files: []
+files:
+- .gitignore
+- Gemfile
+- LICENSE.txt
+- README.md
+- Rakefile
+- agwx_biophys.gemspec
+- lib/agwx_biophys.rb
+- lib/agwx_biophys/et.rb
+- lib/agwx_biophys/version.rb
+- test/Et_testing.java
+- test/TestValues.xls
+- test/test_et.rb
+homepage: ''
+licenses:
+- MIT
+metadata: {}
+post_install_message: 
+rdoc_options: []
+require_paths:
+- lib
+required_ruby_version: !ruby/object:Gem::Requirement
+  requirements:
+  - - '>='
+    - !ruby/object:Gem::Version
+      version: '0'
+required_rubygems_version: !ruby/object:Gem::Requirement
+  requirements:
+  - - '>='
+    - !ruby/object:Gem::Version
+      version: '0'
+requirements: []
+rubyforge_project: 
+rubygems_version: 2.0.6
+signing_key: 
+specification_version: 4
+summary: So far, just an ET calculator (Priestly-Taylor)
+test_files:
+- test/Et_testing.java
+- test/TestValues.xls
+- test/test_et.rb