program
main if(ind.eq.2)then nlh=nl/2+1 global_umax=0. do 50 k=nlh,nl do 50 lgns=1,lat2 do 50 mg=1,lon global_umax=max(global_umax, _ sqrt(ureal(mg,lgns,k)**2+vreal(mg,lgns,k)**2) ) 50 continue end if end program |
program
main real :: velocity if(ind.eq.2)then nlh=nl/2+1 global_umax=0. do 50 k=nlh,nl do 50 lgns=1,lat2 do 50 mg=1,lon global_umax=max(global_umax, _ velocity = vreal(mg,lgns,k) sqrt(ureal(mg,lgns,k)**2+velocity**2) ) 50 continue end if end program |