| |
| | [No title] (Site not responding. Last check: 2007-11-02) |
 | | end ; transformation for lambert conformal conic projection if(ip(4) eq 1) then begin r=(tan(!pi/4.-!dtor*rlat/2.))^rp(10) xp=rp(11)*r*sin(rp(10)*(rp(7)-!dtor*rlon)) yp=-rp(11)*r*cos(rp(10)*(rp(7)-!dtor*rlon)) xp=xp-rp(12) yp=yp-rp(13) end ; transform to 'unprimed' system ; first rotate yd=yp*cos(alpha)-xp*sin(alpha) ; now translate yd=yd+ip(3)*rp(2) return,yd end ;---------------------------------------------------------------------- function rlat,x,y,ip,rp ; find lat for x,y (meters) on bathymetric grid alpha=rp(6)*!pi/180. |
 | | dl=rp(19)*xp+rp(20)*yp+rp(21)*xp*yp+rp(22)*xp*xp rl=rp(0)+dl end ; transformation for lambert conformal conic projection if(ip(4) eq 1) then begin xp=xp+rp(12) yp=yp+rp(13) xx=xp/rp(11) yy=yp/rp(11) r=sqrt(xx^2+yy^2) rl=360.*(!pi/4.-atan(r^(1./rp(10))))/!pi end return,rl end ;---------------------------------------------------------------------- function rlon,x,y,ip,rp ; find lon for x,y (meters) on bathymetric grid alpha=rp(6)*!pi/180. |
 | | dl=rp(15)*xp+rp(16)*yp+rp(17)*xp*yp+rp(18)*xp*xp rl=rp(1)-dl end ; transformation for lambert conformal conic projection if(ip(4) eq 1) then begin xp=xp+rp(12) yp=yp+rp(13) xx=xp/rp(11) yy=yp/rp(11) rl=180.*(rp(7)-atan(xx,-yy)/rp(10))/!pi end return,rl end ;---------------------------------------------------------------------- pro rgrid,fname,gname,dgrid,iparm,rparm ; Purpose: ; To read a bathymetric grid data file ; and return grid parameters and depths. |
| www.glerl.noaa.gov /eegle/resources/model_group/rgrid.pro (229 words) |
|