; Reading multiple files from OPeNDAP server: NY Galadanci
; This code loads MOPITT CO measurements from OPeNDAP server

 load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
 load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
 load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl" ; BRH

 begin

;--- Using the data files from url
       NASA_URL1 = "https://opendap.larc.nasa.gov/opendap/MOPITT/MOP03JM.008/2013.01.01/MOP03JM-201301-L3V95.6.3.he5"
	   NASA_URL2 = "https://opendap.larc.nasa.gov/opendap/MOPITT/MOP03JM.008/2013.02.01/MOP03JM-201302-L3V95.6.3.he5"
	   NASA_URL3 = "https://opendap.larc.nasa.gov/opendap/MOPITT/MOP03JM.008/2013.03.01/MOP03JM-201303-L3V95.6.3.he5"
	   NASA_URL4 = "https://opendap.larc.nasa.gov/opendap/MOPITT/MOP03JM.008/2013.04.01/MOP03JM-201304-L3V95.6.3.he5"
	   NASA_URL5 = "https://opendap.larc.nasa.gov/opendap/MOPITT/MOP03JM.008/2013.05.01/MOP03JM-201305-L3V95.6.3.he5"
	   NASA_URL6 = "https://opendap.larc.nasa.gov/opendap/MOPITT/MOP03JM.008/2013.06.01/MOP03JM-201306-L3V95.6.3.he5"
	   NASA_URL7 = "https://opendap.larc.nasa.gov/opendap/MOPITT/MOP03JM.008/2013.07.01/MOP03JM-201307-L3V95.6.3.he5"
	   NASA_URL8 = "https://opendap.larc.nasa.gov/opendap/MOPITT/MOP03JM.008/2013.08.01/MOP03JM-201308-L3V95.6.3.he5"
	   NASA_URL9 = "https://opendap.larc.nasa.gov/opendap/MOPITT/MOP03JM.008/2013.09.01/MOP03JM-201309-L3V95.6.3.he5"
	   NASA_URL10 = "https://opendap.larc.nasa.gov/opendap/MOPITT/MOP03JM.008/2013.10.01/MOP03JM-201310-L3V95.6.3.he5"
	   NASA_URL11 = "https://opendap.larc.nasa.gov/opendap/MOPITT/MOP03JM.008/2013.11.01/MOP03JM-201311-L3V95.6.3.he5"
	   NASA_URL12 = "https://opendap.larc.nasa.gov/opendap/MOPITT/MOP03JM.008/2013.12.01/MOP03JM-201312-L3V95.6.3.he5"
	   
	   NASA_FILP1 = "MOP03JM-201301-L3V95.6.3.he5"
	   NASA_FILP2 = "MOP03JM-201302-L3V95.6.3.he5"
	   NASA_FILP3 = "MOP03JM-201303-L3V95.6.3.he5"
	   NASA_FILP4 = "MOP03JM-201304-L3V95.6.3.he5"
	   NASA_FILP5 = "MOP03JM-201305-L3V95.6.3.he5"
	   NASA_FILP6 = "MOP03JM-201306-L3V95.6.3.he5"
	   NASA_FILP7 = "MOP03JM-201307-L3V95.6.3.he5"
	   NASA_FILP8 = "MOP03JM-201308-L3V95.6.3.he5"
	   NASA_FILP9 = "MOP03JM-201309-L3V95.6.3.he5"
	   NASA_FILP10 = "MOP03JM-201310-L3V95.6.3.he5"
	   NASA_FILP11 = "MOP03JM-201311-L3V95.6.3.he5"
	   NASA_FILP12 = "MOP03JM-201312-L3V95.6.3.he5"
	   
	   f_1 = NASA_URL1 + NASA_FILP1
	   f_2 = NASA_URL2 + NASA_FILP2
       f_3 = NASA_URL3 + NASA_FILP3
	   f_4 = NASA_URL4 + NASA_FILP4
	   f_5 = NASA_URL5 + NASA_FILP5
	   f_6 = NASA_URL6 + NASA_FILP6
	   f_7 = NASA_URL7 + NASA_FILP7
	   f_8 = NASA_URL8 + NASA_FILP8
	   f_9 = NASA_URL9 + NASA_FILP9
	   f_10 = NASA_URL10 + NASA_FILP10
	   f_11 = NASA_URL11 + NASA_FILP11
	   f_12 = NASA_URL12 + NASA_FILP12
	   files = (/f_1, f_2, f_3, f_4, f_5, f_6, f_7, f_8, f_9, f_10, f_11, f_12/)
    f_all = addfiles(files,"r")
  exists = isfilepresent(files)
  if(.not.exists) then
    print("OPeNDAP test unsuccessful.")
    print("Either the file doesn't exist, or NCL does")
    print("not have OPeNDAP cabilities on this system.")
  else
    f1 = addfile(NASA_URL1 + NASA_FILP1,"r")
	f2 = addfile(NASA_URL2 + NASA_FILP2,"r")
	f3 = addfile(NASA_URL3 + NASA_FILP3,"r")
	f4 = addfile(NASA_URL4 + NASA_FILP4,"r")
	f5 = addfile(NASA_URL5 + NASA_FILP5,"r")
	f6 = addfile(NASA_URL6 + NASA_FILP6,"r")
	f7 = addfile(NASA_URL7 + NASA_FILP7,"r")
	f8 = addfile(NASA_URL8 + NASA_FILP8,"r")
	f9 = addfile(NASA_URL9 + NASA_FILP9,"r")
	f10 = addfile(NASA_URL10 + NASA_FILP10,"r") 
	f11 = addfile(NASA_URL11 + NASA_FILP11,"r")
	f12 = addfile(NASA_URL12 + NASA_FILP12,"r")
	  vnames = getfilevarnames(f1)
	   vnames = getfilevarnames(f2)
	    vnames = getfilevarnames(f3)
		 vnames = getfilevarnames(f4)
		  vnames = getfilevarnames(f5)
		   vnames = getfilevarnames(f6)
		    vnames = getfilevarnames(f7)
			 vnames = getfilevarnames(f8)
			  vnames = getfilevarnames(f9)
			   vnames = getfilevarnames(f10)
			    vnames = getfilevarnames(f11)
				 vnames = getfilevarnames(f12)
    print(vnames)
  end if
 end