--- TC1_eckman/PROG/elcirc5_01_02g.f90 2004-06-20 01:13:56.000000000 -0400 +++ TC1_eckman_nhdf/PROG/elcirc5_01_02g.f90 2004-09-10 13:27:43.000000000 -0400 @@ -1749,8 +1749,10 @@ if(nws.eq.2) then wtime1=0 wtime2=wtiminc +#ifdef USEHDF call get_wind(wtime1,windx1,windy1,pr1,airt1,shum1) call get_wind(wtime2,windx2,windy2,pr2,airt2,shum2) +#endif endif !nws=2 !------------------------------------------------------------------ endif !ihot=0 @@ -1831,8 +1833,10 @@ ninv=time/wtiminc wtime1=ninv*wtiminc wtime2=(ninv+1)*wtiminc +#ifdef USEHDF call get_wind(wtime1,windx1,windy1,pr1,airt1,shum1) call get_wind(wtime2,windx2,windy2,pr2,airt2,shum2) +#endif endif !nws=2 !... Find positions in t.h. files fort.5[0-3] @@ -1944,8 +1948,10 @@ !... Initialize heat budget model if(nws==2.and.ihconsv/=0) then +#ifdef USEHDF call surf_fluxes(wtime1,windx1,windy1,pr1,airt1,shum1,& &srad,fluxsu,fluxlu,hradu,hradd,tauxz,tauyz) +#endif do i=1,np sflux(i)=-fluxsu(i)-fluxlu(i)-(hradu(i)-hradd(i)) enddo @@ -2053,8 +2059,10 @@ if(time>=wtime2) then !... Heat budget & wind stresses if(ihconsv.ne.0) then +#ifdef USEHDF call surf_fluxes(wtime2,windx2,windy2,pr2,airt2,shum2,& &srad,fluxsu,fluxlu,hradu,hradd,tauxz,tauyz) +#endif do i=1,np sflux(i)=-fluxsu(i)-fluxlu(i)-(hradu(i)-hradd(i)) enddo @@ -2071,8 +2079,9 @@ airt1(i)=airt2(i) shum1(i)=shum2(i) enddo +#ifdef USEHDF call get_wind(wtime2,windx2,windy2,pr2,airt2,shum2) - +#endif endif !time.ge.wtime2 wtratio=(time-wtime1)/wtiminc