完善蒸散量计算示例程序

liurunyu

2025-08-28 7d42f93be8077ec79b7569cb57a03d2a1f61d433

 pipIrr-platform/pipIrr-web/pipIrr-web-model/src/main/java/com/dy/pipIrrModel/modelCalculate/Hargreaves.java

@@ -120,30 +120,31 @@
        return result ;
    }


    public static void main(String[] args) {
        Double lat = 38.561976140977116 ;
        Integer dayIndex = 180 ;
        Double lat = 38.56064095143279 ;//纬度
        Integer dayIndex = 239 ;//年内日序数（比如1月1日为1，12月31日为365或366）
        Double kc = 0.41 ;//作物系数
        Double maxT = 40.1 ;//一日内最高温度
        Double minT = 40.1 ;//一日内最低温度
        Double minT = 28.1 ;//一日内最低温度

        Double fai = rad(lat);
        System.out.println(fai);
        System.out.println("弧度:" + fai);

        Double sunMagnetismAngular = sunMagnetismAngular(dayIndex);
        System.out.println(sunMagnetismAngular);

        Double sunEarthDistance = sunEarthDistance(dayIndex);
        System.out.println(sunEarthDistance);
        System.out.println("太阳磁偏角:" + sunMagnetismAngular);

        Double sunTimeAngular = sunTimeAngular(fai, sunMagnetismAngular);
        System.out.println(sunTimeAngular);
        System.out.println("太阳时角:" + sunTimeAngular);

        Double sunEarthDistance = sunEarthDistance(dayIndex);
        System.out.println("日地间相对距离的倒数:" + sunEarthDistance);

        Double zenithRadiation = zenithRadiation(sunEarthDistance, sunTimeAngular, fai, sunMagnetismAngular);
        System.out.println(zenithRadiation);
        System.out.println("天顶辐射:" + zenithRadiation);

        Double et0 = ET0(kc, maxT, minT, zenithRadiation);
        System.out.println(et0);
        System.out.println("蒸散量:" + et0);
    }

}