C+ C NAME: C iWriteProxyMap N C PURPOSE: C Write a proxy map for velocity, density or g-level. C This program is not called in the main program currently, but C something like it is expected to be used when the SMEI data analysis C is running in order that the C CATEGORY: C Data processing C CALLING SEQUENCE: C I = iWriteProxyMapN(iq,nLng,nLat,ROT,AMAP) 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 FUNCTIONS/SUBROUTINES: C bOpenFile C PROCEDURE: C MODIFICATION HISTORY: C NOV, 1995 B. Jackson (STEL,UCSD) C- function iWriteProxyMapN(iq,Rot,nLng,nLat,Z) parameter (nq = 3) character cFileQ(nq)*10 /'v????.???','d????.???','g????.map'/, & cFmt(nq) *7 /'(F10.3)','(F12.1)','(F10.4)'/, & cStr(nq)*80 /'Velocity proxy array','Density proxy array','G-level proxy array'/, & cFile*80, & cDir(2)*80 logical bMirror /.FALSE./ real ZBad(nq), & Z(nLng,nLat) include 'dirspec.h' iStr = 1 ! One string is output cDir(1) = 'dat' cDir(2) = 'proxy' C ZBad(1) = -999.99 9 ZBad(2) = -999999.0 ZBad(3) = -99.9999 I = iFilePath(cHome,2,cDir,cFileQ(iq),cFile) I = index(cFile,'????') if (I .ne. 0 .and. iq. eq. 1) write (cFile(I:I+7),'(F8.3)') Rot if (I .ne. 0 .and. iq. eq. 2) write (cFile(I:I+7),'(F8.3)') Rot if (I .ne. 0 .and. iq. eq. 3) write (cFile(I:I+3),'(I4)') nint(Rot) C call ArrR4Mask(nLng*nLat,Z,BadR4(),ZBad(iq),0.,0.,1.,Z) call wr2darr(1,nLng,nLat,Z,cFile,cFmt(iq),bMirror,iStr,cStr(iq)) C return end