private void checkParametersAndRunEnergyBalance ( double [ ] rain , double [ ] [ ] T , double [ ] [ ] V , double [ ] [ ] P , double [ ] [ ] RH , double month , double day , double hour , double [ ] Abasin , double [ ] [ ] [ ] A , double [ ] [ ] [ ] EI , double [ ] [ ] DTd , double [ ] [ ] DTm , double [ ] [ ] canopy ) { double Dt = ( ( double ) tTimestep / ( double ) pInternaltimestep ) * 60.0 ; /* * some hardcoded variables */ boolean hasNoStations = false ; double zmes_T = 2.0 ; // quota misura temperatura,pressione e umidita' double zmes_U = 2.0 ; // quota misura velocita' vento [m] double z0T = 0.005 ; // [m] roughness length della temperatura double z0U = 0.05 ; // [m] roughness length del vento double K = ka * ka / ( ( log ( zmes_U / z0U ) ) * ( log ( zmes_T / z0T ) ) ) ; double eps = 0.98 ; // emissivita' neve double Lc = 0.05 ; // ritenzione capillare double Ksat = 3.0 ; // 5.55; // conducibilita' idraulica della neve a saturazione double Ks = 5.55E-5 ; // conducibilita' termica superficiale della neve double aep = 50.0 ; // albedo extinction parameter (kg/m2==mm) double rho_g = 1600 ; // densita' del suolo [kg/m3] double De = 0.4 ; // suolo termicamente attivo double C_g = 890.0 ; // capacita' termica del suolo [J/(kg K)] double albedo_land = 0.2 ; double Ts_min = - 20.0 ; double Ts_max = 20.0 ; // TODO check parameters and add to the model parameter latitude = 46.6 * Math . PI / 180.0 ; // [rad] longitude = 10.88 * Math . PI / 180.0 ; // [rad] standard_time = - 1.0 ; // differenza rispetto a UMT [h] /* * start calculations */ sun ( hour , day ) ; for ( int i = 0 ; i < basinNum ; i ++ ) { calculateEnergyBalance ( i , month , hasNoStations , V [ i ] , canopy , T [ i ] , P [ i ] , RH [ i ] , rain , pTrain , pTsnow , Dt , A , Abasin , EI , DTd [ i ] , DTm [ i ] , K , eps , Lc , pRhosnow , Ksat , rho_g , De , C_g , aep , albedo_land , Ks , Ts_min , Ts_max ) ; } }
Method to check the input parameters .