character *69 flnm,end,start,time,model
	character *69 variableb,modelb,startb,endb,flnmb
	character *69 time2,resolutionb,flnm2b
	character *4 timespot2,timespot,format
	character *69 variable,units,resolution,flnm2,barunits
	character *69 minlon, maxlon, minlat, maxlat
	character *69 region,manual,colors
	character *20 stmp,undefined
	real lev(11)
	character *20 variable2,variable2b
	character *2 filelen,dummyvar,filelenb
	character *8 filetime
	character *8 stmpnocolon
	character *15 correction,correction2,correctionb
	character *15 correction2b
	character *10 startdate,enddate1,enddate2,enddate3
        character *6 startdatea,enddate1a,enddate2a,enddate3a
	character *2 startdateb,enddate1b,enddate2b,enddate3b
	character *8 enddate1tot,enddate2tot,enddate3tot
	character *8 enddate1big,enddate2big,enddate3big
	character *13 totstring
	character *8 startdatetot,startdatebig
	integer timeint,timeint2,icount
	real conmin,conmax,cint
	character *20 conminc,conmaxc
        open (unit=25, file='time1.file',status='old')
        read (25,'(A10)') stmp 
	read (25,'(A8)') filetime
        open (unit=26, file='timenocolon1.file',status='old')
	read (26,'(A8)') stmpnocolon
	open (unit=10, file=
     &  'variables1.'//filetime,status='old')
        open (unit=11, file=
     &  'script1.'//filetime//'.dat',status='unknown')
	open (unit=15, file=
     &  'script11.'//filetime//'.dat',status='unknown')
        open (unit=12, file=
     &  'ctlfile1.'//filetime//'.dat',status='unknown')
        open (unit=14, file=
     &  'ctlfile1b.'//filetime//'.dat',status='unknown')
        open (unit=35, file=
     &  'variables1b.'//filetime,status='old')
	open (unit=30, file='/LDAS/PROC/STARTTIME',
     &  status='old')
        open (unit=31, file='/LDAS/PROC/ENDTIME',
     &  status='old')
        open (unit=32, file='/LDAS/PROC/FORENDTIME1',
     &  status='old')
        open (unit=33, file='/LDAS/PROC/FORENDTIME2',
     &  status='old')
	open (unit=34, file='/LDAS/PROC/flag12.dat',
     &  status='old')
	open (unit=49, file='tempp.'//filetime//'.dat',status='unknown')
	conmax=10
	conmin=1
	read (10,100) variable
        read (10,100) model
        read (10,100) start
        read (10,100) end
        read (10,100) flnm
	read (10,100) flnm2
        read (10,100) time
	read (10,100) resolution
	read (10,100) minlon
        read (10,100) maxlon
        read (10,100) minlat
        read (10,100) maxlat
	read (10,100) region
	read (10,100) manual
	read (10,'(A20)') variable2
	read (10,'(A20)') conmaxc
	read (10,'(A20)') conminc
	read (10,'(A4)') format
	if (conmaxc.ne.'') then
	if (conminc.ne.'') then
	write (49,'(A20)') conmaxc
	write (49,'(A20)') conminc
	close (unit=49)
        open (unit=49, file='tempp.'//filetime//'.dat',status='unknown')
	read (49,*) conmax
	read (49,*) conmin
	close (unit=49)
	endif
	endif
        read (30,'(A10)') startdate
	startdatea=startdate(1:6)
	startdateb=startdate(9:10)
	startdatebig=startdatea//startdateb
        icount=1
        do i=1,8
                if(startdatebig(i:i).ne.' ') then
                        if (icount.eq.1) then
                                startdatetot=startdatebig(i:i)
                        else
                        startdatetot(icount:icount)=startdatebig(i:i)
                        endif
                        icount=icount+1
                endif
        enddo


	read (31,'(A10)') enddate1
        enddate1a=enddate1(1:6)
        enddate1b=enddate1(9:10)
        enddate1big=enddate1a//enddate1b
	icount=1
	do i=1,8
		if(enddate1big(i:i).ne.' ') then
			if (icount.eq.1) then
				enddate1tot=enddate1big(i:i)
			else
			enddate1tot(icount:icount)=enddate1big(i:i)
			endif
			icount=icount+1
		endif
	enddo
        read (32,'(A10)') enddate2
        enddate2a=enddate2(1:6)
        enddate2b=enddate2(9:10)
        enddate2big=enddate2a//enddate2b
        icount=1
        do i=1,8
                if(enddate2big(i:i).ne.' ') then
                        if (icount.eq.1) then
                                enddate2tot=enddate2big(i:i)
                        else
                        enddate2tot(icount:icount)=enddate2big(i:i)
                        endif
                        icount=icount+1
                endif
        enddo
        read (33,'(A10)') enddate3
        enddate3a=enddate3(1:6)
        enddate3b=enddate3(9:10)
        enddate3big=enddate3a//enddate3b
        icount=1
        do i=1,8
                if(enddate3big(i:i).ne.' ') then
                        if (icount.eq.1) then
                                enddate3tot=enddate3big(i:i)
                        else
                        enddate3tot(icount:icount)=enddate3big(i:i)
                        endif
                        icount=icount+1
                endif
        enddo
	read (34,'(A2)') dummyvar
	read (34,'(A2)') filelen
c       !!!!!!!!!!!!!!!!!!!!!!!!!!!!
c	setting filelen=48 always...ignoring file read from above
c!!!!!!!!!!!!!!!!!!!!
	filelen='48'
        filelenb='72'
        flnm2b='/LDAS/DATA/LDAS8th.OUTPUT/var2.eta.timese'


        read (35,100) variableb
        read (35,100) modelb
        read (35,100) startb
        read (35,100) endb
        read (35,100) flnmb
        read (35,100) flnm2b
        read (35,100) time2
        read (35,100) resolutionb
        read (35,'(A20)') variable2b

        if ((modelb.eq.'"2nd--For Difference Plots"').or.(
     &  modelb.eq.'"---------------------------------"')) then
            flnm2b='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
	endif



	if (format.eq.'Sn') then
		format='"Sn"'
	endif
	if (format.eq.'DM') then
		format='"DM"'
	endif
	if (format.eq.'FM') then
                format='"FM"'
        endif


 100    format (A60)
	cint=(conmax-conmin)/10
	do i=1,11
		lev(i)=0
	enddo
	
	do i=1,11
		lev(i)=((i-1)*cint)+conmin
	enddo

	if (manual.ne.'"Specify Max/Min"') then
		cint=0
	endif
	timespot2='"1"'
	timespot='"1"'
	timeint2=1
	timeint=1
	start='"11/30/98"'
	correction2='"1"'
	correction='""'
	correctionb='""'
	correction2b='"(1)"'
	end='"12/01/98"'

cif (resolution.eq.'"125"') then
c	resolution='"800"'
c       elseif (resolution.eq.'"200"') then
c                resolution='"1200"'
c        elseif (resolution.eq.'"100"') then
c                resolution='"600"'
celseif (resolution.eq.'"75"') then
c		resolution='"400"'
c	else
c		resolution='"200"'
c	endif
	if (time(2:10).eq.'15Z'//startdatetot(1:6)) then
		time='"15 Z"'
		timespot='"1"'
		timeint=1
	endif
        if (time(2:10).eq.'18Z'//startdatetot(1:6)) then
		time='"18 Z"'
		timespot='"2"'
		timeint=2
	endif
        if (time(2:10).eq.'21Z'//startdatetot(1:6)) then
		time='"21 Z"'
		timespot='"3"'
		timeint=3
	endif
        if (time(2:10).eq.'00Z'//enddate1tot(1:6)) then
		time='"00 Z"'
		timespot='"4"'
		timeint=4
	endif
        if (time(2:10).eq.'03Z'//enddate1tot(1:6)) then
		time='"03 Z"'
		timespot='"5"'
		timeint=5
	endif
        if (time(2:10).eq.'06Z'//enddate1tot(1:6)) then
		time='"06 Z"'
		timespot='"6"'
		timeint=6
	endif
        if (time(2:10).eq.'09Z'//enddate1tot(1:6)) then
		time='"09 Z"'
		timespot='"7"'
		timeint=7
	endif
        if (time(2:10).eq.'12Z'//enddate1tot(1:6)) then
		time='"12 Z"'
		timespot='"8"'
		timeint=8
	endif


        if (time2(2:10).eq.'15Z'//startdatetot(1:6)) then
                time2='"15 Z"'
                timespot2='"1"'
                timeint2=1
        endif
        if (time2(2:10).eq.'18Z'//startdatetot(1:6)) then
                time2='"18 Z"'
                timespot2='"2"'
                timeint2=2
        endif
        if (time2(2:10).eq.'21Z'//startdatetot(1:6)) then
                time2='"21 Z"'
                timespot2='"3"'
                timeint2=3
        endif
        if (time2(2:10).eq.'00Z'//enddate1tot(1:6)) then
                time2='"00 Z"'
                timespot2='"4"'
                timeint2=4
        endif
        if (time2(2:10).eq.'03Z'//enddate1tot(1:6)) then
                time2='"03 Z"'
                timespot2='"5"'
                timeint2=5
        endif
        if (time2(2:10).eq.'06Z'//enddate1tot(1:6)) then
                time2='"06 Z"'
                timespot2='"6"'
                timeint2=6
        endif
        if (time2(2:10).eq.'09Z'//enddate1tot(1:6)) then
                time2='"09 Z"'
                timespot2='"7"'
                timeint2=7
        endif
        if (time2(2:10).eq.'12Z'//enddate1tot(1:6)) then
                time2='"12 Z"'
                timespot2='"8"'
                timeint2=8
        endif







        if (time(2:10).eq.'15Z'//enddate1tot(1:6)) then
                time='"15 Z"'
                timespot='"9"'
                timeint=9
        endif
        if (time(2:10).eq.'18Z'//enddate1tot(1:6)) then
                time='"18 Z"'
                timespot='"10"'
                timeint=10
        endif
        if (time(2:10).eq.'21Z'//enddate1tot(1:6)) then
                time='"21 Z"'
                timespot='"11"'
                timeint=11
        endif
        if (time(2:10).eq.'00Z'//enddate2tot(1:6)) then
                time='"00 Z"'
                timespot='"12"'
                timeint=12
        endif
        if (time(2:10).eq.'03Z'//enddate2tot(1:6)) then
                time='"03 Z"'
                timespot='"13"'
                timeint=13
        endif
        if (time(2:10).eq.'06Z'//enddate2tot(1:6)) then
                time='"06 Z"'
                timespot='"14"'
                timeint=14
        endif
        if (time(2:10).eq.'09Z'//enddate2tot(1:6)) then
                time='"09 Z"'
                timespot='"15"'
                timeint=15
        endif
        if (time(2:10).eq.'12Z'//enddate2tot(1:6)) then
                time='"12 Z"'
                timespot='"16"'
                timeint=16
        endif
        if (time(2:10).eq.'15Z'//enddate2tot(1:6)) then
                time='"15 Z"'
                timespot='"17"'
                timeint=17
        endif
        if (time(2:10).eq.'18Z'//enddate2tot(1:6)) then
                time='"18 Z"'
                timespot='"18"'
                timeint=18
        endif
        if (time(2:10).eq.'21Z'//enddate2tot(1:6)) then
                time='"21 Z"'
                timespot='"19"'
                timeint=19
        endif
        if (time(2:10).eq.'00Z'//enddate3tot(1:6)) then
                time='"00 Z"'
                timespot='"20"'
                timeint=20
        endif
        if (time(2:10).eq.'03Z'//enddate3tot(1:6)) then
                time='"03 Z"'
                timespot='"21"'
                timeint=21
        endif
        if (time(2:10).eq.'06Z'//enddate3tot(1:6)) then
                time='"06 Z"'
                timespot='"22"'
                timeint=22
        endif
        if (time(2:10).eq.'09Z'//enddate3tot(1:6)) then
                time='"09 Z"'
                timespot='"23"'
                timeint=23
        endif
        if (time(2:10).eq.'12Z'//enddate3tot(1:6)) then
                time='"12 Z"'
                timespot='"24"'
                timeint=24
        endif




        if (time2(2:10).eq.'15Z'//enddate1tot(1:6)) then
                time2='"15 Z"'
                timespot2='"9"'
                timeint2=9
        endif
        if (time2(2:10).eq.'18Z'//enddate1tot(1:6)) then
                time2='"18 Z"'
                timespot2='"10"'
                timeint2=10
        endif
        if (time2(2:10).eq.'21Z'//enddate1tot(1:6)) then
                time2='"21 Z"'
                timespot2='"11"'
                timeint2=11
        endif
        if (time2(2:10).eq.'00Z'//enddate2tot(1:6)) then
                time2='"00 Z"'
                timespot2='"12"'
                timeint2=12
        endif
        if (time2(2:10).eq.'03Z'//enddate2tot(1:6)) then
                time2='"03 Z"'
                timespot2='"13"'
                timeint2=13
        endif
        if (time2(2:10).eq.'06Z'//enddate2tot(1:6)) then
                time2='"06 Z"'
                timespot2='"14"'
                timeint2=14
        endif
        if (time2(2:10).eq.'09Z'//enddate2tot(1:6)) then
                time2='"09 Z"'
                timespot2='"15"'
                timeint2=15
        endif
        if (time2(2:10).eq.'12Z'//enddate2tot(1:6)) then
                time2='"12 Z"'
                timespot2='"16"'
                timeint2=16
        endif
        if (time2(2:10).eq.'15Z'//enddate2tot(1:6)) then
                time2='"15 Z"'
                timespot2='"17"'
                timeint2=17
        endif
        if (time2(2:10).eq.'18Z'//enddate2tot(1:6)) then
                time2='"18 Z"'
                timespot2='"18"'
                timeint2=18
        endif
        if (time2(2:10).eq.'21Z'//enddate2tot(1:6)) then
                time2='"21 Z"'
                timespot2='"19"'
                timeint2=19
        endif
        if (time2(2:10).eq.'00Z'//enddate3tot(1:6)) then
                time2='"00 Z"'
                timespot2='"20"'
                timeint2=20
        endif
        if (time2(2:10).eq.'03Z'//enddate3tot(1:6)) then
                time2='"03 Z"'
                timespot2='"21"'
                timeint2=21
        endif
        if (time2(2:10).eq.'06Z'//enddate3tot(1:6)) then
                time2='"06 Z"'
                timespot2='"22"'
                timeint2=22
        endif
        if (time2(2:10).eq.'09Z'//enddate3tot(1:6)) then
                time2='"09 Z"'
                timespot2='"23"'
                timeint2=23
        endif
        if (time2(2:10).eq.'12Z'//enddate3tot(1:6)) then
                time2='"12 Z"'
                timespot2='"24"'
                timeint2=24
        endif



        undefined='-999.9'

        if (variable.eq.'"LW Radiation"') then
                units='"(W/m`a2`n)"'
        barunits='"100 125 150 175 200 225 250 275 300 325 350 375 400"'
                correction='"-0.0"'
                colors='"15 9 14 4 11 5 13 3 10 7 12 8 2 6"'
                if (model.eq.'"MFng"') then
                flnm2='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif

        endif

        if (variable.eq.'"IGBP Vegetation"') then
                units='""'
        barunits='"0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17"'
                correction='""'
	colors=
     &  '"38 28 32 25 24 43 23 36 20 51 22 34 26 52 40 41 39 53 38"'
                if (model.eq.'"MFng"') then
                flnm2='/usr/local/www/LDAS/LDASpredictions/igbp.grads'
                endif

        endif


        if (variable.eq.'"SW Radiation"') then
                units='"(W/m`a2`n)"'
        barunits='"0 60 120 180 240 300 360 420 480 540 600 660 720"'
                correction='"-0.0"'
                colors='"15 9 14 4 11 5 13 3 10 7 12 8 2 6"'
                if (model.eq.'"MFng"') then
                flnm2='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif

        endif


        if (variableb.eq.'"LW Radiation"') then
                correctionb='"-0.0"'
                if (modelb.eq.'"MFng"') then
                flnm2b='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif

        endif

        if (variable.eq.'"SW Radiation"') then
                correctionb='"-0.0"'
                if (modelb.eq.'"MFng"') then
                flnm2b='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif

        endif

        if (variable.eq.'"2 Meter Humidity"') then
                units='"(kg/kg)"'
		variable2='"MeterHumidity"'
        barunits='"0 2 4 6 8 10 12 14 16 18 20 22 24"'
                correction='"*1000.0"'
                colors='"15 9 14 4 11 5 13 3 10 7 12 8 2 6"'
                if (model.eq.'"MFng"') then
                flnm2='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif


        endif

        if (variableb.eq.'"2 Meter Humidity"') then
                correctionb='"*1000.0"'
		variable2b='"MeterHumidity"'
                if (modelb.eq.'"MFng"') then
                flnm2b='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif


        endif


        if (variable.eq.'"Canopy Humidity"') then
                units='"(kg/kg)"'
        barunits='"0 2 4 6 8 10 12 14 16 18 20 22 24"'
                correction='"*1000.0"'
                colors='"15 9 14 4 11 5 13 3 10 7 12 8 2 6"'
                if (model.eq.'"Mosaic"') then
                flnm2='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif
                if (model.eq.'"MFcst"') then
                flnm2='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif


        endif

        if (variableb.eq.'"Canopy Humidity"') then
                correctionb='"*1000.0"'
                if (modelb.eq.'"Mosaic"') then
                flnm2b='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif
                if (modelb.eq.'"MFcst"') then
                flnm2b='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif


        endif


	if (variable.eq.'"Canopy Temp"') then
		units='"(`3.`0K)"'
	barunits='"248 253 258 263 268 273 278 283 288 293 298 303 308"'
		correction='"-0.0"'
		colors='"15 9 14 4 11 5 13 3 10 7 12 8 2 6"'
		if (model.eq.'"Mosaic"') then
		flnm2='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
		endif
                if (model.eq.'"MFcst"') then
                flnm2='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif

	endif


	if (variable.eq.'"Soil Tmp 0-10 cm"') then
                units='"(`3.`0K)"'
	variable2='"SoilTmp1"'
        barunits='"248 253 258 263 268 273 278 283 288 293 298 303 308"'
                correction='"-0.0"'
                colors='"15 9 14 4 11 5 13 3 10 7 12 8 2 6"'
                if (model.eq.'"Eta"') then
                flnm2='/LDAS/DATA/LDAS8th.OUTPUT/var2.eta.timese'
                endif

        endif

        if (variable.eq.'"Soil Tmp 10-40 cm"') then
                units='"(`3.`0K)"'
        variable2='"SoilTmp2"'
        barunits='"248 253 258 263 268 273 278 283 288 293 298 303 308"'
                correction='"-0.0"'
                colors='"15 9 14 4 11 5 13 3 10 7 12 8 2 6"'
                if (model.eq.'"Eta"') then
                flnm2='/LDAS/DATA/LDAS8th.OUTPUT/var2.eta.timese'
                endif

        endif

        if (variable.eq.'"Soil Tmp 40-100 cm"') then
                units='"(`3.`0K)"'
        variable2='"SoilTmp3"'
        barunits='"248 253 258 263 268 273 278 283 288 293 298 303 308"'
                correction='"-0.0"'
                colors='"15 9 14 4 11 5 13 3 10 7 12 8 2 6"'
                if (model.eq.'"Eta"') then
                flnm2='/LDAS/DATA/LDAS8th.OUTPUT/var2.eta.timese'
                endif

        endif

        if (variable.eq.'"Soil Tmp 100-200 cm"') then
                units='"(`3.`0K)"'
        variable2='"SoilTmp4"'
        barunits='"248 253 258 263 268 273 278 283 288 293 298 303 308"'
                correction='"-0.0"'
                colors='"15 9 14 4 11 5 13 3 10 7 12 8 2 6"'
                if (model.eq.'"Eta"') then
                flnm2='/LDAS/DATA/LDAS8th.OUTPUT/var2.eta.timese'
                endif

        endif


        if (variableb.eq.'"Soil Tmp 0-10 cm"') then
        variable2b='"SoilTmp1"'
                correctionb='"-0.0"'
                if (modelb.eq.'"Eta"') then
                flnm2b='/LDAS/DATA/LDAS8th.OUTPUT/var2.eta.timese'
                endif

        endif

        if (variableb.eq.'"Soil Tmp 10-40 cm"') then
        variable2b='"SoilTmp2"'
                correctionb='"-0.0"'
                if (modelb.eq.'"Eta"') then
                flnm2b='/LDAS/DATA/LDAS8th.OUTPUT/var2.eta.timese'
                endif

        endif

        if (variableb.eq.'"Soil Tmp 40-100 cm"') then
        variable2b='"SoilTmp3"'
                correctionb='"-0.0"'
                if (modelb.eq.'"Eta"') then
                flnm2b='/LDAS/DATA/LDAS8th.OUTPUT/var2.eta.timese'
                endif

        endif

        if (variableb.eq.'"Soil Tmp 100-200 cm"') then
        variable2b='"SoilTmp4"'
                correctionb='"-0.0"'
                if (modelb.eq.'"Eta"') then
                flnm2b='/LDAS/DATA/LDAS8th.OUTPUT/var2.eta.timese'
                endif

        endif



        if (variable.eq.'"Soil Mst 0-10 cm"') then
                units='"(mm)"'
        variable2='"SoilMst1"'
        barunits='"0 6 12 18 24 30 36 42 48 54 60 66 72"'
                correction='"*100"'
                colors='"15 9 14 4 11 5 13 3 10 7 12 8 2 6"'
                if (model.eq.'"Eta"') then
                flnm2='/LDAS/DATA/LDAS8th.OUTPUT/var2.eta.timese'
                endif

        endif

        if (variable.eq.'"Soil Mst 10-40 cm"') then
                units='"(mm)"'
        variable2='"SoilMst2"'
        barunits='"0 15 30 45 60 75 90 105 120 135 150 165 180"'
                correction='"*300"'
                colors='"15 9 14 4 11 5 13 3 10 7 12 8 2 6"'
                if (model.eq.'"Eta"') then
                flnm2='/LDAS/DATA/LDAS8th.OUTPUT/var2.eta.timese'
                endif

        endif

        if (variable.eq.'"Soil Mst 40-100 cm"') then
                units='"(mm)"'
        variable2='"SoilMst3"'
        barunits='"0 30 60 90 120 150 180 210 240 270 300 330 360"'
                correction='"*600"'
                colors='"15 9 14 4 11 5 13 3 10 7 12 8 2 6"'
                if (model.eq.'"Eta"') then
                flnm2='/LDAS/DATA/LDAS8th.OUTPUT/var2.eta.timese'
                endif

        endif

        if (variable.eq.'"Soil Mst 100-200 cm"') then
                units='"(mm)"'
        variable2='"SoilMst4"'
        barunits='"0 60 120 180 240 300 360 420 480 540 600 660 720"'
                correction='"*1000"'
                colors='"15 9 14 4 11 5 13 3 10 7 12 8 2 6"'
                if (model.eq.'"Eta"') then
                flnm2='/LDAS/DATA/LDAS8th.OUTPUT/var2.eta.timese'
                endif

        endif




        if (variable.eq.'"Volumetric Soil Mst 0-10 cm"') then
                units='"(m`a3`n/m`a3`n)"'
        variable2='"SoilMst1"'
        barunits='"0.0 .05 .10 .15 .2 .25 .3 .35 .4 .45 .5 .55 .6"'
                correction='"*1"'
                colors='"15 9 14 4 11 5 13 3 10 7 12 8 2 6"'
                if (model.eq.'"Eta"') then
                flnm2='/LDAS/DATA/LDAS8th.OUTPUT/var2.eta.timese'
                endif

        endif

        if (variable.eq.'"Volumetric Soil Mst 10-40 cm"') then
                units='"(m`a3`n/m`a3`n)"'
        variable2='"SoilMst2"'
        barunits='"0.0 .05 .10 .15 .2 .25 .3 .35 .4 .45 .5 .55 .6"'
                correction='"*1"'
                colors='"15 9 14 4 11 5 13 3 10 7 12 8 2 6"'
                if (model.eq.'"Eta"') then
                flnm2='/LDAS/DATA/LDAS8th.OUTPUT/var2.eta.timese'
                endif

        endif

        if (variable.eq.'"Volumetric Soil Mst 40-100 cm"') then
                units='"(m`a3`n/m`a3`n)"'
        variable2='"SoilMst3"'
        barunits='"0.0 .05 .10 .15 .2 .25 .3 .35 .4 .45 .5 .55 .6"'
                correction='"*1"'
                colors='"15 9 14 4 11 5 13 3 10 7 12 8 2 6"'
                if (model.eq.'"Eta"') then
                flnm2='/LDAS/DATA/LDAS8th.OUTPUT/var2.eta.timese'
                endif

        endif

        if (variable.eq.'"Volumetric Soil Mst 100-200 cm"') then
                units='"(m`a3`n/m`a3`n)"'
        variable2='"SoilMst4"'
        barunits='"0.0 .05 .10 .15 .2 .25 .3 .35 .4 .45 .5 .55 .6"'
                correction='"*1"'
                colors='"15 9 14 4 11 5 13 3 10 7 12 8 2 6"'
                if (model.eq.'"Eta"') then
                flnm2='/LDAS/DATA/LDAS8th.OUTPUT/var2.eta.timese'
                endif

        endif





        if (variableb.eq.'"Soil Mst 0-10 cm"') then
        variable2b='"SoilMst1"'
                correctionb='"*100"'
                if (modelb.eq.'"Eta"') then
                flnm2b='/LDAS/DATA/LDAS8th.OUTPUT/var2.eta.timese'
                endif

        endif

        if (variableb.eq.'"Soil Mst 10-40 cm"') then
        variable2b='"SoilMst2"'
                correctionb='"*300"'
                if (modelb.eq.'"Eta"') then
                flnm2b='/LDAS/DATA/LDAS8th.OUTPUT/var2.eta.timese'
                endif

        endif

        if (variableb.eq.'"Soil Mst 40-100 cm"') then
        variable2b='"SoilMst3"'
                correctionb='"*600"'
                if (modelb.eq.'"Eta"') then
                flnm2b='/LDAS/DATA/LDAS8th.OUTPUT/var2.eta.timese'
                endif

        endif

        if (variableb.eq.'"Soil Mst 100-200 cm"') then
        variable2b='"SoilMst4"'
                correctionb='"*1000"'
                if (modelb.eq.'"Eta"') then
                flnm2b='/LDAS/DATA/LDAS8th.OUTPUT/var2.eta.timese'
                endif

        endif

        if (variableb.eq.'"Volumetric Soil Mst 0-10 cm"') then
        variable2b='"SoilMst1"'
                correctionb='"*1"'
                if (modelb.eq.'"Eta"') then
                flnm2b='/LDAS/DATA/LDAS8th.OUTPUT/var2.eta.timese'
                endif

        endif

        if (variableb.eq.'"Volumetric Soil Mst 10-40 cm"') then
        variable2b='"SoilMst2"'
                correctionb='"*1"'
                if (modelb.eq.'"Eta"') then
                flnm2b='/LDAS/DATA/LDAS8th.OUTPUT/var2.eta.timese'
                endif

        endif

        if (variableb.eq.'"Volumetric Soil Mst 40-100 cm"') then
        variable2b='"SoilMst3"'
                correctionb='"*1"'
                if (modelb.eq.'"Eta"') then
                flnm2b='/LDAS/DATA/LDAS8th.OUTPUT/var2.eta.timese'
                endif

        endif

        if (variableb.eq.'"Volumetric Soil Mst 100-200 cm"') then
        variable2b='"SoilMst4"'
                correctionb='"*1"'
                if (modelb.eq.'"Eta"') then
                flnm2b='/LDAS/DATA/LDAS8th.OUTPUT/var2.eta.timese'
                endif

        endif


        if (variableb.eq.'"Canopy Temp"') then
                correctionb='"-0.0"'
                if (modelb.eq.'"Mosaic"') then
                flnm2b='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif
                if (model.eq.'"MFcst"') then
                flnm2b='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif

        endif


        if (variable.eq.'"Deep Soil Temp"') then
                units='"(`3.`0K)"'
        barunits='"248 253 258 263 268 273 278 283 288 293 298 303 308"'
                correction='"-0.0"'
		colors='"15 9 14 4 11 5 13 3 10 7 12 8 2 6"'
                if (model.eq.'"Mosaic"') then
                      flnm2='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif
                if (model.eq.'"MFcst"') then
                flnm2='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif

        endif

        if (variableb.eq.'"Deep Soil Temp"') then
                correctionb='"-0.0"'
                if (modelb.eq.'"Mosaic"') then
                      flnm2b='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif
                if (modelb.eq.'"MFcst"') then
                flnm2b='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif

        endif


        if (variable.eq.'"Total Runoff"') then
                units='"(mm/3hr)"'
        barunits='".5 1 1.5 2 2.5 3 3.5 4 4.5 5 5.5 6 6.5"'
                colors='"15 9 14 4 11 5 13 3 10 7 12 8 2 6"'
                correction='"*10800"'
                if (model.eq.'"Mosaic"') then
                      flnm2='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif
                if (model.eq.'"MFcst"') then
                flnm2='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif

        endif

        if (variableb.eq.'"Total Runoff"') then
                correctionb='"*10800"'
                if (modelb.eq.'"Mosaic"') then
                      flnm2b='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif
                if (modelb.eq.'"MFcst"') then
                flnm2b='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif

        endif

        if (variable.eq.'"Surface Runoff"') then
                units='"(mm/3hr)"'
        barunits='".5 1 1.5 2 2.5 3 3.5 4 4.5 5 5.5 6 6.5"'
                colors='"15 9 14 4 11 5 13 3 10 7 12 8 2 6"'
                correction='"*10800"'
                if (model.eq.'"Mosaic"') then
                      flnm2='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif
                if (model.eq.'"MFcst"') then
                flnm2='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif

        endif

        if (variableb.eq.'"Surface Runoff"') then
                correctionb='"*10800"'
                if (modelb.eq.'"Mosaic"') then
                      flnm2b='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif
                if (modelb.eq.'"MFcst"') then
                flnm2b='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif

        endif


        if (variable.eq.'"Sfc Soil Mst"') then
                units='"(mm)"'
        barunits='"0.0 1 2 3 4 5 6 7 8 9 10 11 12"'
                correction='"*1"'
                colors='"15 9 14 4 11 5 13 3 10 7 12 8 2 6"'
                if (model.eq.'"Mosaic"') then
                      flnm2='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif
                if (model.eq.'"MFcst"') then
                flnm2='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif

        endif

        if (variableb.eq.'"Sfc Soil Mst"') then
                correctionb='"*1"'
                if (modelb.eq.'"Mosaic"') then
                      flnm2b='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif
                if (modelb.eq.'"MFcst"') then
                flnm2b='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif

        endif


        if (variable.eq.'"Volumetric Sfc Soil Mst"') then
		variable2='"SfcSoilWet"'
                units='"(m`a3`n/m`a3`n)"'
        barunits='"0.0 .05 .10 .15 .2 .25 .3 .35 .4 .45 .5 .55 .6"'
                correction='"*1"'
                colors='"15 9 14 4 11 5 13 3 10 7 12 8 2 6"'
                if (model.eq.'"Mosaic"') then
                      flnm2='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif
                if (model.eq.'"MFcst"') then
                flnm2='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif

        endif

        if (variableb.eq.'"Volumetric Sfc Soil Mst"') then
		variable2b='"SfcSoilWet"'
                correctionb='"*1"'
                if (modelb.eq.'"Mosaic"') then
                      flnm2b='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif
                if (modelb.eq.'"MFcst"') then
                flnm2b='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif

        endif


        if (variable.eq.'"Volumetric Mid Soil Mst"') then
                variable2='"MidSoilWet"'
                units='"(m`a3`n/m`a3`n)"'
        barunits='"0.0 .05 .10 .15 .2 .25 .3 .35 .4 .45 .5 .55 .6"'
                correction='"*1"'
                colors='"15 9 14 4 11 5 13 3 10 7 12 8 2 6"'
                if (model.eq.'"Mosaic"') then
                      flnm2='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif
                if (model.eq.'"MFcst"') then
                flnm2='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif

        endif

        if (variableb.eq.'"Volumetric Mid Soil Mst"') then
                variable2b='"MidSoilWet"'
                correctionb='"*1"'
                if (modelb.eq.'"Mosaic"') then
                      flnm2b='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif
                if (modelb.eq.'"MFcst"') then
                flnm2b='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif

        endif

        if (variable.eq.'"Volumetric Deep Soil Mst"') then
                variable2='"DeepSoilWet"'
                units='"(m`a3`n/m`a3`n)"'
        barunits='"0.0 .05 .10 .15 .2 .25 .3 .35 .4 .45 .5 .55 .6"'
                correction='"*1"'
                colors='"15 9 14 4 11 5 13 3 10 7 12 8 2 6"'
                if (model.eq.'"Mosaic"') then
                      flnm2='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif
                if (model.eq.'"MFcst"') then
                flnm2='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif

        endif

        if (variableb.eq.'"Volumetric Deep Soil Mst"') then
                variable2b='"DeepSoilWet"'
                correctionb='"*1"'
                if (modelb.eq.'"Mosaic"') then
                      flnm2b='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif
                if (modelb.eq.'"MFcst"') then
                flnm2b='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif

        endif




        if (variable.eq.'"Mid Soil Mst"') then
                units='"(mm)"'
        barunits='"0 50 100 150 200 250 300 350 400 450 500 550 600"'
                correction='"*1"'
                colors='"15 9 14 4 11 5 13 3 10 7 12 8 2 6"'
                if (model.eq.'"Mosaic"') then
                      flnm2='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif
                if (model.eq.'"MFcst"') then
                flnm2='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif

        endif

        if (variableb.eq.'"Mid Soil Mst"') then
                correctionb='"*1"'
                if (modelb.eq.'"Mosaic"') then
                      flnm2b='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif
                if (modelb.eq.'"MFcst"') then
                flnm2b='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif

        endif


        if (variable.eq.'"Deep Soil Mst"') then
                units='"(mm)"'
        barunits='"0 75 150 225 300 375 425 500 575 625 700 775 850"'
                correction='"*1"'
                colors='"15 9 14 4 11 5 13 3 10 7 12 8 2 6"'
                if (model.eq.'"Mosaic"') then
                      flnm2='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif
                if (model.eq.'"MFcst"') then
                flnm2='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif

        endif


        if (variableb.eq.'"Deep Soil Mst"') then
                correctionb='"*1"'
                if (modelb.eq.'"Mosaic"') then
                      flnm2b='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif
                if (modelb.eq.'"MFcst"') then
                flnm2b='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif

        endif


        if (variable.eq.'"Sensible Heat Flux"') then
                units='"(W/m`a2`n)"'
        barunits='"-240 -200 -160 -120 -80 -40 0 40 80 120 160 200 240"'
                correction='"*1"'
		variable2='"SensHeatFlux"'

                colors='"15 9 14 4 11 5 13 3 10 7 12 8 2 6"'
                if (model.eq.'"Mosaic"') then
                      flnm2='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif
                if (model.eq.'"MFcst"') then
                flnm2='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif
                if (model.eq.'"Eta"') then
                flnm2='/LDAS/DATA/LDAS8th.OUTPUT/var2.eta.timese'
                endif


        endif


        if (variableb.eq.'"Sensible Heat Flux"') then
                correctionb='"*1"'
                variable2b='"SensHeatFlux"'
                if (modelb.eq.'"Mosaic"') then
                      flnm2b='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif
                if (modelb.eq.'"MFcst"') then
                flnm2b='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif
                if (modelb.eq.'"Eta"') then
                flnm2b='/LDAS/DATA/LDAS8th.OUTPUT/var2.eta.timese'
                endif


        endif


        if (variable.eq.'"Latent Heat Flux"') then
                units='"(W/m`a2`n)"'
	barunits='"-180 -150 -120 -90 -60 -30 0 30 60 90 120 150 180"'
                correction='"*1"'
                colors='"15 9 14 4 11 5 13 3 10 7 12 8 2 6"'
                if (model.eq.'"Mosaic"') then
                      flnm2='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif
                if (model.eq.'"MFcst"') then
                flnm2='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif
                if (model.eq.'"Eta"') then
                flnm2='/LDAS/DATA/LDAS8th.OUTPUT/var2.eta.timese'
                endif


        endif

        if (variableb.eq.'"Latent Heat Flux"') then
                correctionb='"*1"'
                if (modelb.eq.'"Mosaic"') then
                      flnm2b='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif
                if (modelb.eq.'"MFcst"') then
                flnm2b='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif
                if (modelb.eq.'"Eta"') then
                flnm2b='/LDAS/DATA/LDAS8th.OUTPUT/var2.eta.timese'
                endif
        endif

       if (variable.eq.'"3hr Evaporation"') then
                variable2='LatentHeatFlux'
                units='"(mm)"'
        barunits='".05 .1 .15 .2 .25 .3 .35 .4 .45 .5 .55 .6 .65"'
                correction='"*-4.3619E-03"'
                colors='"15 9 14 4 11 5 13 3 10 7 12 8 2 6"'
                if (model.eq.'"Mosaic"') then
                      flnm2='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif
                if (model.eq.'"MFcst"') then
                flnm2='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif
                if (model.eq.'"Eta"') then
                flnm2='/LDAS/DATA/LDAS8th.OUTPUT/var2.eta.timese'
                endif


        endif

        if (variableb.eq.'"3hr Evaporation"') then
                variable2b='LatentHeatFlux'
                correctionb='"*4.3619E-03"'
                if (modelb.eq.'"Mosaic"') then
                      flnm2b='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif
                if (modelb.eq.'"MFcst"') then
                flnm2b='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif
                if (modelb.eq.'"Eta"') then
                flnm2b='/LDAS/DATA/LDAS8th.OUTPUT/var2.eta.timese'
                endif



        endif

       if (variable.eq.'"Water Balance Error"') then
                units='"(mm)"'
		variable2='WaterBalance'
        barunits='".25 .5 .75 1 1.25 1.5 1.75 2 2.25 2.5 2.75 3 3.25"'
                correction='"*1"'
                colors='"15 9 14 4 11 5 13 3 10 7 12 8 2 6"'
                if (model.eq.'"Mosaic"') then
                      flnm2='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif
                if (model.eq.'"MFcst"') then
                flnm2='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif
                if (model.eq.'"Eta"') then
                flnm2='/LDAS/DATA/LDAS8th.OUTPUT/var2.eta.timese'
                endif


        endif

        if (variableb.eq.'"Water Balance Error"') then
		variable2b='WaterBalance'
                correctionb='"*1"'
                if (modelb.eq.'"Mosaic"') then
                      flnm2b='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif
                if (modelb.eq.'"MFcst"') then
                flnm2b='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif
                if (modelb.eq.'"Eta"') then
                flnm2b='/LDAS/DATA/LDAS8th.OUTPUT/var2.eta.timese'
                endif



        endif

       if (variable.eq.'"Tot Colmn Soil Mst"') then
                units='"(mm)"'
                barunits=
     &  '"100 200 300 400 500 600 700 800 900 1000 1100 1200 1300"'
                correction='"*1"'
		variable2='SoilMstTot'
                colors='"15 9 14 4 11 5 13 3 10 7 12 8 2 6"'
                if (model.eq.'"Mosaic"') then
                      flnm2='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif
                if (model.eq.'"MFcst"') then
                flnm2='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif
                if (model.eq.'"Eta"') then
                flnm2='/LDAS/DATA/LDAS8th.OUTPUT/var2.eta.timese'
                endif


        endif

        if (variableb.eq.'"Tot Colmn Soil Mst"') then
                correctionb='"*1"'
		variable2b='SoilMstTot'
                if (modelb.eq.'"Mosaic"') then
                      flnm2b='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif
                if (modelb.eq.'"MFcst"') then
                flnm2b='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif
                if (modelb.eq.'"Eta"') then
                flnm2b='/LDAS/DATA/LDAS8th.OUTPUT/var2.eta.timese'
                endif
	endif




       if (variable.eq.'"Intercepted Water"') then
		variable2='InterWater'
                units='"(mm)"'
        barunits='".05 .1 .15 .2 .25 .3 .35 .4 .45 .5 .55 .6 .65"'
                correction='"*1"'
                colors='"15 9 14 4 11 5 13 3 10 7 12 8 2 6"'
                if (model.eq.'"Mosaic"') then
                      flnm2='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif
                if (model.eq.'"MFcst"') then
                flnm2='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif
                if (model.eq.'"Eta"') then
                flnm2='/LDAS/DATA/LDAS8th.OUTPUT/var2.eta.timese'
                endif


        endif

        if (variableb.eq.'"Intercepted Water"') then
		variable2b='InterWater'
                correctionb='"*1"'
                if (modelb.eq.'"Mosaic"') then
                      flnm2b='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif
                if (modelb.eq.'"MFcst"') then
                flnm2b='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif
                if (modelb.eq.'"Eta"') then
                flnm2b='/LDAS/DATA/LDAS8th.OUTPUT/var2.eta.timese'
                endif



        endif


       if (variable.eq.'"Grnd Water Runoff"') then
                units='"(mm)"'
        barunits='".5 1 1.5 2 2.5 3 3.5 4 4.5 5 5.5 6 6.5"'
                correction='"*1"'
                colors='"15 9 14 4 11 5 13 3 10 7 12 8 2 6"'
                if (model.eq.'"Mosaic"') then
                      flnm2='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif
                if (model.eq.'"MFcst"') then
                flnm2='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif
                if (model.eq.'"Eta"') then
                flnm2='/LDAS/DATA/LDAS8th.OUTPUT/var2.eta.timese'
                endif


        endif

        if (variableb.eq.'"Grnd Water Runoff"') then
                correctionb='"*1"'
                if (modelb.eq.'"Mosaic"') then
                      flnm2b='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif
                if (modelb.eq.'"MFcst"') then
                flnm2b='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif
                if (modelb.eq.'"Eta"') then
                flnm2b='/LDAS/DATA/LDAS8th.OUTPUT/var2.eta.timese'
                endif



        endif


       if (variable.eq.'"Storm Sfc Runoff"') then
                units='"(mm)"'
        barunits='".5 1 1.5 2 2.5 3 3.5 4 4.5 5 5.5 6 6.5"'
                correction='"*1"'
                colors='"15 9 14 4 11 5 13 3 10 7 12 8 2 6"'
                if (model.eq.'"Mosaic"') then
                      flnm2='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif
                if (model.eq.'"MFcst"') then
                flnm2='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif
                if (model.eq.'"Eta"') then
                flnm2='/LDAS/DATA/LDAS8th.OUTPUT/var2.eta.timese'
                endif


        endif

        if (variableb.eq.'"Storm Sfc Runoff"') then
                correctionb='"*1"'
                if (modelb.eq.'"Mosaic"') then
                      flnm2b='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif
                if (modelb.eq.'"MFcst"') then
                flnm2b='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif
                if (modelb.eq.'"Eta"') then
                flnm2b='/LDAS/DATA/LDAS8th.OUTPUT/var2.eta.timese'
                endif



        endif


        if (variable.eq.'"Ground Heat Flux"') then
                units='"(W/m`a2`n)"'
        barunits='"-120 -100 -80 -60 -40 -20 0 20 40 60 80 100 120"'
                correction='"*1"'
                colors='"15 9 14 4 11 5 13 3 10 7 12 8 2 6"'
                if (model.eq.'"Mosaic"') then
                      flnm2='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif
                if (model.eq.'"MFcst"') then
                flnm2='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif
                if (model.eq.'"Eta"') then
                flnm2='/LDAS/DATA/LDAS8th.OUTPUT/var2.eta.timese'
                endif
        endif

        if (variableb.eq.'"Ground Heat Flux"') then
                correctionb='"*1"'
                if (modelb.eq.'"Mosaic"') then
                      flnm2b='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif
                if (modelb.eq.'"MFcst"') then
                flnm2b='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif
                if (modelb.eq.'"Eta"') then
                flnm2b='/LDAS/DATA/LDAS8th.OUTPUT/var2.eta.timese'
                endif
        endif


        if (variable.eq.'"Snow Depth [liq]"') then
	variable2='SnowDepth'
                units='"(mm)"'
        barunits='"0 7 14 21 28 35 42 49 56 63 70 77 84"'
                correction='"*1"'
                colors='"15 9 14 4 11 5 13 3 10 7 12 8 2 6"'
                if (model.eq.'"Mosaic"') then
                      flnm2='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif
                if (model.eq.'"MFcst"') then
                flnm2='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif
                if (model.eq.'"Eta"') then
                flnm2='/LDAS/DATA/LDAS8th.OUTPUT/var2.eta.timese'
                endif


	endif

        if (variableb.eq.'"Snow Depth [liq]"') then
        	variable2b='SnowDepth'
                correctionb='"*1"'
                if (modelb.eq.'"Mosaic"') then
                      flnm2b='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif
                if (modelb.eq.'"MFcst"') then
                flnm2b='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif
                if (modelb.eq.'"Eta"') then
                flnm2b='/LDAS/DATA/LDAS8th.OUTPUT/var2.eta.timese'
                endif


        endif

        if (variable.eq.'"Snow Melt [liq]"') then
        variable2='SnowMelt'
                units='"(mm)"'
        barunits='"0 1 2 3 4 5 6 7 8 9 10 11 12"'
                correction='"*1"'
                colors='"15 9 14 4 11 5 13 3 10 7 12 8 2 6"'
                if (model.eq.'"Mosaic"') then
                      flnm2='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif
                if (model.eq.'"MFcst"') then
                flnm2='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif
                if (model.eq.'"Eta"') then
                flnm2='/LDAS/DATA/LDAS8th.OUTPUT/var2.eta.timese'
                endif


        endif

        if (variableb.eq.'"Snow Melt [liq]"') then
                variable2b='SnowMelt'
                correctionb='"*1"'
                if (modelb.eq.'"Mosaic"') then
                      flnm2b='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif
                if (modelb.eq.'"MFcst"') then
                flnm2b='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif
                if (modelb.eq.'"Eta"') then
                flnm2b='/LDAS/DATA/LDAS8th.OUTPUT/var2.eta.timese'
                endif


        endif

        if (variableb.eq.'"3hr Snow Fall [liq]"') then
                variable2b='SnowFall'
                correctionb='"*1"'
                if (modelb.eq.'"Mosaic"') then
                      flnm2b='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif
                if (modelb.eq.'"MFcst"') then
                flnm2b='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif
                if (modelb.eq.'"Eta"') then
                flnm2b='/LDAS/DATA/LDAS8th.OUTPUT/var2.eta.timese'
                endif


        endif

        if (variable.eq.'"3hr Snow Fall [liq]"') then
        variable2='SnowFall'
                units='"(mm)"'
        barunits='"0 3 6 9 12 15 18 21 24 27 30 33 36"'
                correction='"*1"'
                colors='"15 9 14 4 11 5 13 3 10 7 12 8 2 6"'
                if (model.eq.'"Mosaic"') then
                      flnm2='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif
                if (model.eq.'"MFcst"') then
                flnm2='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif
                if (model.eq.'"Eta"') then
                flnm2='/LDAS/DATA/LDAS8th.OUTPUT/var2.eta.timese'
                endif


        endif
        if (variableb.eq.'"Prec Type Snow"') then
                variable2b='Snow'
                correctionb='"*1"'
                if (modelb.eq.'"Mosaic"') then
                      flnm2b='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif
                if (modelb.eq.'"MFcst"') then
                flnm2b='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif
                if (modelb.eq.'"Eta"') then
                flnm2b='/LDAS/DATA/LDAS8th.OUTPUT/var2.eta.timese'
                endif


        endif

        if (variable.eq.'"Prec Type Snow"') then
        variable2='Snow'
                units='""'
        barunits='"0 1"'
                correction='"*1"'
                colors='"0 4 0"'
                if (model.eq.'"Mosaic"') then
                      flnm2='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif
                if (model.eq.'"MFcst"') then
                flnm2='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif
                if (model.eq.'"Eta"') then
                flnm2='/LDAS/DATA/LDAS8th.OUTPUT/var2.eta.timese'
                endif


        endif
        if (variableb.eq.'"Prec Type Rain"') then
                variable2b='Rain'
                correctionb='"*1"'
                if (modelb.eq.'"Mosaic"') then
                      flnm2b='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif
                if (modelb.eq.'"MFcst"') then
                flnm2b='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif
                if (modelb.eq.'"Eta"') then
                flnm2b='/LDAS/DATA/LDAS8th.OUTPUT/var2.eta.timese'
                endif


        endif

        if (variable.eq.'"Prec Type Rain"') then
        variable2='Rain'
                units='""'
        barunits='"0 1"'
                correction='"*1"'
                colors='"0 4 0"'
                if (model.eq.'"Mosaic"') then
                      flnm2='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif
                if (model.eq.'"MFcst"') then
                flnm2='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif
                if (model.eq.'"Eta"') then
                flnm2='/LDAS/DATA/LDAS8th.OUTPUT/var2.eta.timese'
                endif


        endif

        if (variableb.eq.'"Prec Type Frz Rain"') then
                variable2b='FreezingRain'
                correctionb='"*1"'
                if (modelb.eq.'"Mosaic"') then
                      flnm2b='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif
                if (modelb.eq.'"MFcst"') then
                flnm2b='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif
                if (modelb.eq.'"Eta"') then
                flnm2b='/LDAS/DATA/LDAS8th.OUTPUT/var2.eta.timese'
                endif


        endif

        if (variable.eq.'"Prec Type Frz Rain"') then
        variable2='FreezingRain'
                units='""'
        barunits='"0 1"'
                correction='"*1"'
                colors='"0 4 0"'
                if (model.eq.'"Mosaic"') then
                      flnm2='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif
                if (model.eq.'"MFcst"') then
                flnm2='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif
                if (model.eq.'"Eta"') then
                flnm2='/LDAS/DATA/LDAS8th.OUTPUT/var2.eta.timese'
                endif


        endif

        if (variableb.eq.'"Prec Type Sleet"') then
                variable2b='IcePellet'
                correctionb='"*1"'
                if (modelb.eq.'"Mosaic"') then
                      flnm2b='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif
                if (modelb.eq.'"MFcst"') then
                flnm2b='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif
                if (modelb.eq.'"Eta"') then
                flnm2b='/LDAS/DATA/LDAS8th.OUTPUT/var2.eta.timese'
                endif


        endif

        if (variable.eq.'"Prec Type Sleet"') then
        variable2='IcePellet'
                units='""'
        barunits='"0 1"'
                correction='"*1"'
                colors='"0 4 0"'
                if (model.eq.'"Mosaic"') then
                      flnm2='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif
                if (model.eq.'"MFcst"') then
                flnm2='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif
                if (model.eq.'"Eta"') then
                flnm2='/LDAS/DATA/LDAS8th.OUTPUT/var2.eta.timese'
                endif


        endif



        if (variable.eq.'"Snow Depth [sol]"') then
	        variable2='SnowDepth'
                units='"(Inches)"'
        barunits='"0 3 6 9 12 15 18 21 24 27 30 33 36"'
                correction='"*0.3937"'
                colors='"15 9 14 4 11 5 13 3 10 7 12 8 2 6"'
                if (model.eq.'"Mosaic"') then
                      flnm2='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif
                if (model.eq.'"MFcst"') then
                flnm2='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif

        endif

        if (variableb.eq.'"Snow Depth [sol]"') then
	        variable2b='SnowDepth'
                correctionb='"*0.3937"'
                if (modelb.eq.'"Mosaic"') then
                      flnm2b='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif
                if (modelb.eq.'"MFcst"') then
                flnm2b='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif

        endif

        if ((variable2b.eq.'"2MAirTemp[K]"').or.
     &  (variable2b.eq.'"2MAirTemp[F]"')) then
		variable2b='"AirTemp"'
	endif

        if (variable.eq.'"2 M Air Temp [K]"') then
                units='"(`3.`0K)"'
		variable='"2m Air Temp"'
		variable2='"AirTemp"'
       barunits='"238 245 252 259 266 273 280 287 294 301 308 315 322"'
                correction='""'
                colors='"15 9 14 4 11 5 13 3 10 7 12 8 2 6"'
                if (model.eq.'"MFng"') then
                      flnm2='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif
                if (model.eq.'"Eta"') then
                      flnm2='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif
        endif

        if (variableb.eq.'"2 M Air Temp [K]"') then
                variableb='"2m Air Temp"'
                variable2b='"AirTemp"'
                correctionb='""'
                if (modelb.eq.'"MFng"') then
                      flnm2b='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif
                if (model.eq.'"Eta"') then
                      flnm2='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif
        endif

        if (variable.eq.'"2 M Air Temp [C]"') then
                units='"(`3.`0C)"'
                variable='"2m Air Temp"'
                variable2='"AirTemp"'
       barunits='"-35 -28 -21 -14 -7 0 7 14 21 28 35 42 49"'
                correction='"-273.15"'
                colors='"15 9 14 4 11 5 13 3 10 7 12 8 2 6"'
                if (model.eq.'"MFng"') then
                      flnm2='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif
                if (model.eq.'"Eta"') then
                      flnm2='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif
        endif

        if (variableb.eq.'"2 M Air Temp [C]"') then
                variableb='"2m Air Temp"'
                variable2b='"AirTemp"'
                correctionb='"-273.15"'
                if (modelb.eq.'"MFng"') then
                      flnm2b='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif
                if (model.eq.'"Eta"') then
                      flnm2='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif
        endif



        if (variable.eq.'"2 M Air Temp [F]"') then
		variable='"2m Air Temp"'
                variable2='"AirTemp"'
                units='"(`3.`0F)"'
                barunits='"-18 -8 2 12 22 32 42 52 62 72 82 92 102"'
                correction='"-273.15"'
		correction2='"(9/5)+32"'
                colors='"15 9 14 4 11 5 13 3 10 7 12 8 2 6"'
                if (model.eq.'"MFng"') then
                flnm2='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif
                if (model.eq.'"Eta"') then
                      flnm2='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif
        endif


        if (variableb.eq.'"2 M Air Temp [F]"') then
                variableb='"2m Air Temp"'
                variable2b='"AirTemp"'
                correctionb='"-273.15"'
                correction2b='"(9/5)+32"'
                if (modelb.eq.'"MFng"') then
                      flnm2b='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif
                if (model.eq.'"Eta"') then
                      flnm2='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif
        endif


        if (variable.eq.'"Wind Speed"') then
                units='"(m/s)"'
                barunits='"0 1 2 3 4 5 6 7 8 9 10 11"'
                colors='"15 9 14 4 11 5 13 3 10 7 12 8 2 6"'
                if (model.eq.'"MFng"') then
                      flnm2='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif
        endif

        if (variableb.eq.'"Wind Speed"') then
                if (modelb.eq.'"MFng"') then
                      flnm2b='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif
        endif


        if (variable.eq.'"Surface Pressure"') then
                variable2='"pressure"'
                units='"(mb)"'
       barunits='"700 730 760 790 820 850 880 910 940 970 1000 1030"'
                correction='"/100"'
                colors='"15 9 4 11 5 13 3 10 7 12 8 2"'
                if (model.eq.'"MFng"') then
                      flnm2='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif
        endif

        if (variableb.eq.'"Surface Pressure"') then
                variable2b='"pressure"'
                correctionb='"/100"'
                if (modelb.eq.'"MFng"') then
                      flnm2b='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif
        endif

        if (variable.eq.'"Total Precip"') then
                variable2='"TotalPrecip"'
                units='"(mm/3hr)"'
       barunits='"0 1 2 3 4 5 6 7 8 9 10 11"'
                correction='"*10800"'
                colors='"15 9 4 11 5 13 3 10 7 12 8 2 6"'
                if (model.eq.'"MFng"') then
                      flnm2='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif
        endif


        if (variableb.eq.'"Total Precip"') then
                variable2b='"TotalPrecip"'
                correctionb='"*10800"'
                if (modelb.eq.'"MFng"') then
                      flnm2b='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif
        endif


        if (variable.eq.'"Conv Precip"') then
                variable2='"ConvPrecip"'
                units='"(mm/hour)"'
       barunits='"0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 5.5"'
                correction='"*3600"'
                colors='"15 9 4 11 5 13 3 10 7 12 8 2 6"'
                if (model.eq.'"MFng"') then
                      flnm2='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif
        endif


        if (variableb.eq.'"Conv Precip"') then
                variable2b='"ConvPrecip"'
                correctionb='"*3600"'
                if (modelb.eq.'"MFng"') then
                      flnm2b='/LDAS/DATA/LDAS8th.OUTPUT/var2.timeseries'
                endif
        endif


	
	
		if (region.eq.'"Entire LDAS Domain"') then
			minlon='"-124.9375"'
			maxlon='"-67.0625"'
			minlat='"25.0625"'
			maxlat='"52.9375"'
		endif
                if (region.eq.'"Eastern US"') then
                        minlon='"-96"'
                        maxlon='"-67.0625"'
                        minlat='"25.0625"'
                        maxlat='"50"'
                endif
                if (region.eq.'"Western US"') then
                        minlon='"-124.9375"'
                        maxlon='"-95"'
                        minlat='"25.0625"'
                        maxlat='"50"'
                endif
                if (region.eq.'"Northeast"') then
                        minlon='"-85"'
                        maxlon='"-67.0625"'
                        minlat='"36"'
                        maxlat='"50"'
                endif
                if (region.eq.'"Southeast"') then
                        minlon='"-96"'
                        maxlon='"-74"'
                        minlat='"25.0625"'
                        maxlat='"38"'
                endif
                if (region.eq.'"Southwest"') then
                        minlon='"-123"'
                        maxlon='"-94"'
                        minlat='"25.0625"'
                        maxlat='"38"'
                endif
                if (region.eq.'"Midwest"') then
                        minlon='"-103"'
                        maxlon='"-83"'
                        minlat='"36"'
                        maxlat='"50"'
                endif
                if (region.eq.'"Northwest"') then
                        minlon='"-124.9375"'
                        maxlon='"-101"'
                        minlat='"36"'
                        maxlat='"50"'
                endif
                if (region.eq.'"Washington DC Area"') then
                        minlon='"-78.3"'
                        maxlon='"-74.4"'
                        minlat='"37.8"'
                        maxlat='"40.0"'
                endif
                if (region.eq.'"Los Angeles Area"') then
                        minlon='"-120.7"'
                        maxlon='"-116.9"'
                        minlat='"32.7"'
                        maxlat='"35"'
                endif


	if ((manual.eq.'"Dynamic-Automatic"').or.
     &	(manual.eq.'"Specify Max/Min"')) then
		barunits='""'
	endif
	if (timeint.le.3) then
	 	write (11,*)'date=','"',startdate,'"'
	endif
	
	if ((timeint.ge.4).and.(timeint.le.11)) then
		write (11,*)'date=','"',enddate1,'"'
	endif
 
        if ((timeint.ge.12).and.(timeint.le.19)) then
                write (11,*)'date=','"',enddate2,'"'
        endif

        if ((timeint.ge.20).and.(timeint.le.24)) then
                write (11,*)'date=','"',enddate3,'"'
        endif

        if ((modelb.eq.'"2nd--For Difference Plots"').or.(
     &  modelb.eq.'"---------------------------------"')) then
                correctionb='"*0"'
		correction2b='"(1)"'
		variable2b='"CanopyTemp"'
        endif

	write (11,*)'startdate="',startdate,'"'
	write (11,*)'enddate1="',enddate1,'"'
	write (11,*)'enddate2="',enddate2,'"'
	write (11,*)'enddate3="',enddate3,'"'
	write (11,*)'flnm=',flnm
        write (11,*)'end=',end
        write (11,*)'start=',start
	write (11,*)'time=',time
	write (11,*)'time2=',time2
	write (11,*)'timespot=',timespot
	write (11,*)'timespot2=',timespot2
	write (11,*)'variable=',variable
	write (11,*)'units=',units
	write (11,*)'contours=',manual
cif (modelb.eq.'"2nd--For Difference Plots"') then
celse
	write (11,*)'barunits=',barunits
cendif
	write (11,*)'type=',format
	write (11,*)'minlon=',minlon
        write (11,*)'maxlon=',maxlon
        write (11,*)'minlat=',minlat
        write (11,*)'maxlat=',maxlat
	write (11,'(A5,A10)')'stmp=',stmp
	write (11,'(A13,A6,A1)')'stmpnocolon="',stmpnocolon,'"'
	write (11,*)'res=',resolution
	write (11,*)'colors=',colors
	write (11,*)'correction=',correction
	write (11,*)'correction2=',correction2
        write (11,*)'correctionb=',correctionb
        write (11,*)'correction2b=',correction2b
        if (variable.eq.'"IGBP Vegetation"') then
	timeint=1
	variable2b='"IGBPVegetation"'
	endif
	write (11,*)'var2=',variable2,'".1 (t=',timeint,')"'
	write (11,*)'varshort=',variable2,'".1"'
	write (11,*)'varshortb=',variable2b,'".2"'
        write (11,*)
     &  'var2b=',variable2b,'".2 (t=',timeint2,')"'
	write (11,*) 'modelb=',modelb
	write (11,*) 'v2=',variableb
        if (manual.eq.'"Specify Max/Min"') then
c                write (11,*)'conmax="',conmax,'"'
c                write (11,*)'conmin="',conmin,'"'
c                write (11,*)'cint="',cint,'"'
		write (11,*) 'lev1="',lev(1),'"'
                write (11,*) 'lev2="',lev(2),'"'
                write (11,*) 'lev3="',lev(3),'"'
                write (11,*) 'lev4="',lev(4),'"'
                write (11,*) 'lev5="',lev(5),'"'
                write (11,*) 'lev6="',lev(6),'"'
                write (11,*) 'lev7="',lev(7),'"'
                write (11,*) 'lev8="',lev(8),'"'
                write (11,*) 'lev9="',lev(9),'"'
                write (11,*) 'lev10="',lev(10),'"'
                write (11,*) 'lev11="',lev(11),'"'

	endif
	if (region.ne.'"Entire LDAS Domain"') then
c        write (11,*)'"set parea 2 9 0 8.5"'	
	endif
	write (11,'(A12,A8,A5)')'"open ldas1.',filetime,'.ctl"'
        write (11,'(A13,A8,A5)')'"open ldas1b.',filetime,'.ctl"'

        if (manual.eq.'"Specify Max/Min"') then
c                write (11,*)'"set cmax "conmax'
c                write (11,*)'"set cmin "conmin'
c                write (11,*)'"set cint "cint'
cwrite (11,*)'"set rbrange "conmin" "conmax'
	endif
	if (modelb.eq.'"2nd--For Difference Plots"') then
	if (manual.eq.'"Static For All Models"') then
	        write (11,*)'barunits=',barunits
		write (11,*)'"set clevs "barunits'
	endif
	endif
	write (11,*)'"set grads off"'	
        write (11,*)'"set mpdset hires"'
        write (11,*)'"set lat "minlat" "maxlat'
        write (11,*)'"set lon "minlon" "maxlon'
        write (11,*)'"set gxout grfill"'
        write (11,*)'"set ccols "colors'

        if (manual.eq.'"Specify Max/Min"') then
c                write (15,*)'"set cmax "conmax'
c                write (15,*)'"set cmin "conmin'
c                write (15,*)'"set cint "cint'
        endif

        if (manual.eq.'"Static For All Models"') then
                write (15,*)'"set clevs "barunits'
        endif



c	write (12,101)'DSET /LDAS/DATA/LDAS8th.OUTPUT/var1.timeseries'
	write (12,101)'DSET ',flnm2
	write (12,104)'UNDEF -999.9'
	write (12,103) 'TDEF    ',filelen,' LINEAR  00Z01jan87    3hr'
        if (variable.eq.'"IGBP Vegetation"') then
        flnm2b='/usr/local/www/LDAS/LDASpredictions/igbp.grads'
	endif
        write (14,101)'DSET ',flnm2b
        write (14,104)'UNDEF -999.9'
        write (14,103) 'TDEF    ',filelenb,' LINEAR  00Z01jan87    3hr'




 101   format (A5,A50)
c 101	format (A43)
 102   format (A6,A15)
 103   format (A8,A2,A26)
 104   format (A12)

	end