REM ? MCV-TDIFF -- Creates an RGB image showing IR clouds over a Basemap
REM ?      MCV-TDIFF
REM ? ----------
    
rem
rem --- User can define the name of the temporary data set by changing the value of t$
rem 
     keyin "DATALOC ADD MCV-MCB sos.ssec.wisc.edu"
     let t$ = "LOCAL/AREAS"
     let y$ = "MCV-MCB/AREAS"

rem
rem --- User can define the directory name for storing temporary data required to create composite
rem --- Change the string r$ from temp_directory to your choice.  
rem --- The default value is $HOME/mcidas/data/RAD_TEMP
rem
     let r$ = "MCV-TEMP"

rem
rem --- Determine value of $HOME - needed for REDIRECT commands
rem
          keyin "OS DEV=T USER_HOME R {echo $HOME}"
          open "USER_HOME",10
          input @10,a$
          h$ = mid$(a$,1,len(a$))
          r$ = h$ + "/mcidas/data/MCV-TEMP"
          close 10

rem
rem --- Create a temporary data set to store temporary images for creating composite
rem
     keyin "DSSERVE ADD "; t$ ;" AREA 9001 9299 {Temporary areas for creating McV Demo}"

rem
rem --- Create a temporary directory for the temporary images for creating composite
rem
     keyin "OS {mkdir " ; r$ ;"}"
     keyin "OS {mkdir " ; r$ ;"/final}"

rem
rem --- Set up a file redirection for temporary images
rem
     keyin "REDIRECT ADD AREA90* {"; r$ ;"}"

rem
rem --- We don't have real-time data so we'll use a canned dataset contains a 
rem --- global composite of infrared satellite, GFS surface temperatures and SST's
rem ---
rem --- MCV-MCB/AREAS.1 is the global IR composite
rem --- MCV-MCB/AREAS.2 is the GFS temperatures
rem --- MCV-MCB/AREAS.3 is the SST temperatues
rem --- MCV-MCB/AREAS.4 is the land sea mask
rem --- MCV-MCB/AREAS.5 is the red values for the Big Blue Marble
rem --- MCV-MCB/AREAS.6 is the green values for the Big Blue Marble
rem --- MCV-MCB/AREAS.7 is the blue values for the Big Blue Marble
rem

rem 
rem --- Copy the SST and SFC areas 
rem 
     keyin "IMGFILT "; y$ ;".2 "; y$ ;".4 "; t$ ;".10 SIZE=ALL SCALE=0 255 0 255 FILT=DIS 0 255 4 4 0"
     keyin "IMGFILT "; y$ ;".3 "; y$ ;".4 "; t$ ;".11 SIZE=ALL SCALE=0 255 0 255 FILT=DIS 0 255 8 8 0"
     keyin "IMGOPER "; t$ ;".10 "; t$ ;".11 "; t$ ;".12 SIZE=ALL FORM=MAX ZERO=DATA"
     keyin "PRDUTIL ADD "; t$ ;".12 2 T 320 200 0 255"
     keyin "IMGOPER "; y$ ;".1 "; t$ ;".13 SIZE=ALL SCAL=320 200 0 255 PROD=T UNIT=TEMP"

rem 
rem --- Subract the temperature values
rem 

     keyin "IMGOPER "; t$ ;".12 "; t$ ;".13 "; t$ ;".32 COEF=1 -1 FORM=ADD UNIT=T T SIZE=ALL SCALE=-10 100 0 255"
     keyin "IMGRGB "; t$ ;".32 "; y$ ;".1 "; t$ ;".40 FILT=TDIFF BASE="; y$ ;".5 "; y$ ;".6 "; y$ ;".7 " 
     keyin "OS {cp "; r$ ;"/AREA9040 "; r$ ;"/final"
     stop