Commit 985d9d5e authored by ph290's avatar ph290
Browse files

update

parent c4add0c5
......@@ -2444,8 +2444,8 @@ spec_hum2=0.62*vap/(airP(idmet)-0.38*vap)
!hl=0.985*5.67d-8*(surf_temp**4.0)*(0.39-0.05*vap**0.5)*(1.0-0.6d-4*cloud(idmet)**2.0) !original codes Longwave
!Net longwave radiation flux=IR_surface_emissivity_1_for_black_body(Stefan_Boltzmann_constant * T**4.8 - downwelling_longwave) Note T should really be skin temp
!hl=0.985*5.67d-8*(surf_temp**4.0)-(lw_rad_down0) !Longwave
hl=0.99*5.67d-8*(surf_temp**4.0)-(lw_rad_down0*0.75) !Longwave The 0.95 is
!playing around to account for lw being absorbed and loast in teh skin rather
hl=0.98*5.67d-8*(surf_temp**4.0)-(lw_rad_down0*0.95) !Longwave The 0.95 is
!playing around to account for lw being absorbed and lost in the skin rather
!than top meters
sk=1.45d-3*1.3*1004.0*wind_speed(idmet)*(temp_old(N)-airT(idmet)) !Sensible
ske=1.5d-3*1.3*wind_speed(idmet)*(spec_hum1-spec_hum2)*(2.5d6-2.3d3*temp_old(N)) !Evaporation
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment