;
; Name:
;      PSIMAS
;
;Purpose:
;
; Predictive Science MAS model (PSIMAS) for HMI or MDI (depending on date)
; from web page:  http://www.predsci.com/hmi/data_access.php
;
;Purpose:
;
;       ESSENTIALLY SAME AS NUMCUBE, BUT WITH GETMAS FLAG SET
;
;               IF CUBENAME AND READPRAMS NOT SET - this then looks up 
;		PSI MAS model according to date from web page
;
;               If GETMAS not set, acts exactly like NUMCUBEPRAMS
;
;
;Keyword inputs:
;
;  CUBENAME: Filename (inluding extention and path) of the numerical datacube
;            you want to look at. File must be in the correct format,
;            use make_my_cube.pro in this directory to convert your
;            cube. 
;
;
;  VELIMPOSE - impose a velocity of magnitude VELIMPOSE directed along the field
;		will overwrite any velocity field already loaded into the cube
;		if nonzero
;			UNITS KM/SEC
;                       DEFAULT 0.d0
;
;  DTHRES    - DOUBLE : Identify POP2 coronal points that are dense and can be treated
;                       as different population
;                        If unset or set to 0, will be ignored; if set to non-zero
;                        will be interpreted as density above which points are
;                        designated population 2.  (If set to 1 will default to 5d8 cgs)
;                        This will result in creation of POP2DENS zeroed everywhere except
;                        D>THRESD - this will be set up in NUMCUBE.PRO
;                        This will overwrite any POP2DENS tag already in CUBE.
;                        Note if COLDPTS also set, DTHRES will be ignored.
;                        DEFAULT 0.d0
;
;  COLDPTS     - DOUBLE : Identify points that are cold as different population
;                        If unset or set to 0, will be ignored; if set to non-zero
;                        will be interpreted as temperature below which points are
;                        designated population 2.  (If set to 1 will default to 3d4 K)
;                        This will result in creation of POP2DENS zeroed everywhere except
;                        T<COLDPTS - this will be set up in NUMCUBE.PRO
;                        This will overwrite any POP2DENS tag already in CUBE.
;                        DEFAULT 0.d0
;
; P2FILL   - DOUBLE :  A single number between 0 and 1 to multiply population 2 density by
;                       so that POP2FILLINGFACTOR=P2FILL in POP2 voxels, and 0 elsewhere, and
;                       FILLINGFACTOR=1-P2FILL in POP2 voxels, and 1 elsewhere.
;                       Will overwrite any FILLINGFACTOR and/or POP2FILLINGFACTOR tag on CUBE.
;                       Note if there is no POP2DENS or POP2TEMP in cube and both COLDPTS and DTHRES
;                       are not set, P2FILL will be ignored if set to anything other than default
;                       DEFAULT 1.d0
;
;