C+ C NAME: C iReadProxyMap C PURPOSE: C Read a proxy map for velocity, density or g-level. C CATEGORY: C Data processing C CALLING SEQUENCE: function iReadProxyMap(iq,Rot,nLng,nLat,Z) C INPUTS: C iq integer Type of map C iq = 1 velocity map (vmap.dat) C iq = 2 density map (dmap.dat) C iq = 3 g-level map (gmap.dat) C nLng integer # longitudes points in GMAP C nLat integer # latitude points in GMAP C ROT real Modified Carrington variable C OUTPUTS: C AMAP(nLng,nLat) real Proxy map read from file C INCLUDE: include 'dirspec.h' C CALLS: C bOpenFile, iFilePath, iFltArr C PROCEDURE: C MODIFICATION HISTORY: C NOV, 1995 B. Jackson (STEL,UCSD) C- parameter (nq = 3) character cFileQ(nq)*10 /'vmap.dat','d????.map','g????.map'/ real ZBad(nq) /-1,0.,0./ real Z(nLng,nLat) character cFile*512 external fncWZ integer nPad(4) /4*0/ I = iFilePath('$dat',0,' ',cFileQ(iq),cFile) I = index(cFile,'????') if (I .ne. 0) write (cFile(I:I+3),'(I4)') nint(Rot) nX = 0 nY = 0 iReadProxyMap = iFltArr(cFile,0,0.,1.,fncWZ,BadR4(),nLng*nLat,nX,nY,nPad,Z) c call ArrR4Mask(nLng*nLat,Z,ZBad(iq),BadR4(),0.,0.,1.,Z) return end