diff MFLUX2.f MFLUX2.f.old > diff_MFLUX2.f.txt 133c133,135 < real :: windmks,znotm, znott --- > real :: yz,y1,y2,y3,y4,windmks,znotm > real :: g1, g2, g3, g4, g5 > real :: w1, w2, w3, w4, w5, w6, w7 138a141,168 > yz= 0.0001344 > y1= 3.015e-05 > y2= 1.517e-06 > y3= -3.567e-08 > y4= 2.046e-10 > > g1 = -5.0375e-10 > g2 = 1.2913e-07 > g3 = -1.2323e-05 > g4 = 4.9848e-04 > g5 = -0.0052 > > c w1 = 4.5142e-16 > c w2 = -6.7892e-13 > c w3 = 1.6040e-10 > c w4 = -1.5575e-08 > c w5 = 7.2416e-07 > c w6 = -1.5485e-05 > c w7 = 1.3757e-04 > > w1 = -4.1339e-15 > w2 = 1.2190e-12 > w3 = -1.3543e-10 > w4 = 6.9407e-09 > w5 = -1.5108e-07 > w6 = 1.0665e-06 > w7 = 1.7438e-05 > 166,169c196,232 < call znot_m(windmks,znotm) < zoc(i) = -100.*znotm < call znot_t(windmks,znott) < zot(i) = -100* znott --- > if(wind(i) .lt. 1250.0 ) then > zoc(i)=-(0.0185/9.8*(7.59e-8*wind(i)**2+ > & 2.46e-4*wind(i))**2)*100. > endif > if(wind(i) .ge. 1250.0 .and. wind(i) .lt. 3000.0) > & zoc(i)=-(.000739793 * wind(i) -0.58)/10 > > if(wind(i) .ge. 3000.0 .and. wind(i) .le. 4000.0) then > znotm=yz+windmks*y1+windmks**2*y2+windmks**3*y3+ > & windmks**4*y4 !powell 2003 > zoc(i) = -100.*znotm > endif > > if(wind(i) .gt. 4000.0 ) then > znotm = g5 + g4*windmks + g3*windmks**2 + g2*windmks**3 > & + g1*windmks**4 > zoc(i) = -100.*znotm > endif > > if(windmks .le. 7.0) then > znott = (0.0185/9.8*(7.59e-8*wind(i)**2+ > & 2.46e-4*wind(i))**2) > zot(i) = -100* znott > endif > > if(windmks .gt. 7.0 .and. windmks .lt. 20.0 ) then > znott=0.2375*exp(-0.5250*windmks)+ > & 0.0025*exp(-0.0211*windmks) > zot(i) = -znott > endif > > c > if(windmks .ge. 20.0 ) then > zot(i)= -(w1*windmks**6 + w2*windmks**5 + w3*windmks**4 + > & w4*windmks**3 + w5*windmks**2 + w6 *windmks + w7)*100. > endif > c