Displaying 2 results from an estimated 2 matches for "latitcr".
Did you mean:
latitcrs
2005 Apr 23
0
reading fortran binary file
...e field
& ,IB,JB,KB,LB ! starting indices
& ,IE,JE,KE,LE ! ending indices
& ,XTIME ! the ingegration or forcast time
& ,CRDT,CORDER,CHRDATE
& ,CNAME,CUNIT ! field name and unit description
if(cname.eq.'LATITCRS ')then
read(11)out2d
else
read(11)
endif
goto 10
elseif(iflag.eq.2)then
goto 20
endif
20 close(11)
end
2.compile the subroutine as share object file with -byteswapio optin
3."dyn.load" the subroutine in R
4...
2008 Jul 09
0
"Rotated Lat-Lon" projection in mapproj
...;t identified a projection in library mapproj that can handle this.
Any ideas on how to plot this field on a map?
A typical R script to plot other outputs in a lambert projection would be:
##Read lat/lon data
lat<-as.numeric(system(paste(pycmd,'control',"200108150000","LATITCRS",1,0,sep="
"),intern=TRUE))
lon<-as.numeric(system(paste(pycmd,'control',"200108150000","LONGICRS",1,0,sep="
"),intern=TRUE))
## Projected coordinates
coords_lamb<-mapproject(lon,lat,projection="lambert",parameters=c(60.,30.)...