&proc_input ;tmpdir = '/hao/roble/timegcm1.1/tmpdir' ;tmpdir = '/ptmp/roble/timegcm1.1/tmpdir' ;tmpdir = '/hao/roble/timegcm1.2/tmpdir' ;tmpdir = '/hao/roble/timegcm1/tmpdir' ;tmpdir = '/hao/tgcm/tgcm24rh/py2n' ;tmpdir = '/hao/roble/timegcm1' ;tmpdir = '/hao/tgcm/tgcm24rh/sy2v' ; ; histvols = '/ROBLE/timegcm1.3/seqst004.nc' mtimes = 4,0,0 ; ; ;histvols = '/ROBLE/timegcm1.3/seqsm007.nc' ;mtimes = 98,0,0 ; ; ;histvols_cntr = '/ROBLE/timegcm1.3/seqsm009.nc' ;mtimes_cntr = 50,0,0 ; ; ;histvols='/ROBLE/timegcm1.3/seqsm005.nc','to','/ROBLE/timegcm1.3/seqsm006.nc', ; 'by','1', ;mtimes = 85,0,0,'to',86,0,0,'by', 60 ; ; ;histvols = '/ROBLE/timegcm1.3/seqstn05.nc' ;mtimes = 12,0,0 ; ; ;histvols = '/ROBLE/timegcm1.3/sdssm009.nc' ;mtimes = 54,0,0 ; ; ; ;histvols='/ROBLE/timegcm1/sdsaz001.nc','to','/ROBLE/timegcm1/sdsaz002.nc', ; 'by','1', ;mtimes = 0,0,0,'to', 1,0,0,'by', 60 ; ; ; mtimes = integer triplets, where each triplet represents model ; day, hour, and minute. The processor will sequentially ; search histvols for mtimes histories. ; Note mtimes may specify: mtimes = begin_time, 'to', end_time, 'by', dmin ; where "begin_time" and "end_time" are starting and ending model times ; (3 ints), and "dmin" is delta time in minutes, e.g.: ; mtimes = 0,0,0,'to',1,0,0,'by',60 ; ; mtimes =35,0,0,'to',51,0,0,'by',1440 ; ;------------------------------------------------------------------ ;histvols_cntr = '/ROBLE/t24co2r/seqsm02.nc' ;mtimes_cntr = 70,0,0 ;------------------------------------------------------------------ maxdiskvols = 10 isearch = 0 ; ; Request fields: cfields is string array with the following valid values: ; TN UN VN O2 OX N4S NOZ CO ; CO2 H2O H2 HOX O+ CH4 AR HE ; O21D NO2 NO O3 O1 OH HO2 H ; N2D TI TE NE O2+ W Z POTEN ; UI VI WI N2 RHO UN+VN UI+VI FOF2 ; HMF2 O/O2 O/N2 N2/O O2/N2 O/O2+N2 E6300 E5577 ; EO200 EOH83 ECO215u ENO53u OHB-XX OHV-X ; ; 'NaS ','NaO ','NaO3 ','NaO2 ','NaOH ', ; 'NaCO3 ','NaHCO3 ','NaS+ ','NaN2+ ','NaCO2+ ', ; 'NaH2O+ ','NaO+ ','NaEMIS ' ; ;cfields = 'TN','UN','VN','UTNDCY','VTNDCY','TTNDCY','DIFKK','DIFKT' ; 'GWU','UMOM','RAYK' ; cfields = 'TN','UN','VN','Z','NE' ; ;cfields = 'UTNDCY' iprint_fldminmax = 0 ;modelhts = 1 modelhts = 0 idifpercent = 0 iden = 2 ;iden = 2 ; fscale = 'H2O',1.e+6,'CO2',1.e+6, 'CH4',1.e+6, 'O3',1.e+6, ; 'H2',1.e+6, 'HOX',1.e+6, 'NO',1.e+6 ; fmnmxint = 'TN',1.,0.,50., 'UN',1.,0.,10, ; 'VN',1.,0.,10., 'UTNDCY',1.,0.,10., ; 'DIFKK',1.,0.,10. ; fmnmxint = 'TN',1.,0.,10., 'H2',1.,0.,0.5, 'HOX',1.,0.,1.0, ; 'NO',1.,0.,10.,'CO2',1.,0.,25., 'H2O',1.,0.,0.5, ; 'CH4',1.,0.,0.2, 'O3',1.,0.,0.5 ; fmnmxint = 'TN',1.,0.,10., 'UN',1.,0.,10.0 ; fmnmxint = 'TN',1.,0.,5., 'UN',1.,0.,5.0, 'VN',1.,0.,1.0 ; fmnmxint = 'TN',1.,0.,5., 'H2O',1.,0.,10., ; 'CO2',1.,0.,5.,'CH4',1.,0.,10., ; 'O3',-100.,100.,10.,'H2',1.,0.,10., ; 'HOX',1.,0.,100.,'NO',1.,0.,100. ; fmnmxint = 'TN',1.,0.,10., 'H2O',1.,0.,0.5, ; 'UN',1.,0.,10., 'VN',1.,0.,10., ; 'CO2',1.,0.,25.,'CH4',1.,0.,0.1, ; 'O3',1.,0.,1. ; Species density conversion flag: ; iden=0 -> leave species as on history (most are mass mix ratios) ; iden=1 -> convert species to number densities (cm3) ; iden=2 -> convert species to number density mixing ratios ; iden=3 -> convert species to mass density (gm/cm3) ionvel = 1 ; ionvel = 1 -> ExB only ; ionvel = 2 -> ExB+unvn (total) ; ionvel = 3 -> ExB+unvn, with E=0 ; ionvel = 4 -> ExB+unvn, with unvnwn=0. ;------------------------------------------------------------------ ;---------------------------------------------------------------------- ; Volume Emissions (airglow) fields: ; (note: these fields are calculated only from time-gcm histories) ; (note: iden must be 1 (cm3 species) for these calculations) ; (fmnmxint and fscale will work for E-fields, but not ht-integrations) ; ; (add SR63) ie6300 = 0 ; ; The following fields are required on the history for e5577: ; if ie5577(1) > 0 -> need tn,o2,o,n2 ; if ie5577(2) > 0 -> need te,o2p,ne,o21d ; if ie5577(3) > 0 -> need tn,ht,o2,o,n2,o21d and call solred ; if ie5577(4) > 0 -> need tn,ht,o2,o,n2,o21d and call solred ; if ie5577(5) > 0 -> need tn,ht,o2,o,n2,o21d and call solred ; (if o21d is needed but is not available, it will be set to 1.e-20) ; (if te is needed but it not available, will use te = tn) ; ; (o1 recomb) (o2+ recomb) (photoe) (airglow) (o2 ly-beta) ie5577 = 1, 0, 0, 0, 0 ; ; idoppler=1 -> generate doppler t,u,v for each emissions field except ; ECO215u and ENO53u. (ipltxyut type line plots only) idoppler = 0 ; ; iyd,f107a,f107d are needed for solred calculation, i.e., if ; ie5577(3, 4, or 5) > 0, or ie6300 = 1. ; (if not provided here, will use values from history header) ;iyd = 94080 ;f107d = 67. ;f107a = 72. ; ; ieohv = plot flags for 10 oh vibrational states (ieohv(1) is ground) ; (each state plotted as separate field if iehov(i) > 0 ; for plot types ipltutvert and/or ipltxyut (ht-integrations)) ; (fmnmxint and fscale not implemented for ieohv) ; ;ieohv = 0,1,1,1,1,1,1,1,1,1 ;ieohv = 0,1,0,0,0,0,0,0,0,0 ;ieohv = 0,0,1,0,0,0,0,1,0,0 ;ieohv = 0,0,0,0,0,0,0,0,0,0 ; ; ibohv = plot flags for oh-v bands (9-8 to 9-3, 8-7 to 8-3, etc) ; (these will be in kilo-Rayleighs). The following are valid bands ; (ibohv = 2-digit number giving high and low range, e.g. 98 for 9-8) ; ; 1-0 2-1 3-2 4-3 5-4 6-5 7-6 8-7 9-8 ; 2-0 3-1 4-2 5-3 6-4 7-5 8-6 9-7 ; 3-0 4-1 5-2 6-3 7-4 8-5 9-6 ; 4-0 5-1 6-2 7-3 8-4 9-5 ; 5-0 6-1 7-2 8-3 9-4 ; 6-0 7-1 8-2 9-3 ; ibohv = 83,62 ; icolor = 1 ;outplt = 'cgm','ps' outplt = 'cgm' ;outplt = 'ps' multiplt = 0 ipltrowcol = 2,2 multiadvfr = 1 ; iwr_label = 1,1,1,0,0,0 ;------------------------------------------------------------------ ipltmaps = 1 ; plot flag for maps map_global = 1 ; CE projections (cyl equidistant) map_global_cenlon = 0 ; center x-axis longitude for CE ; map_global_censlt = 12 ; center x-axis local time for CE map_polar = 0 ; ST projections (polar stereographic) fmap_polar_perimlat = 45. map_satview = 0 ; SV projections (sat view) ; fmap_satview_latlon = 20.,-80. fmap_satview_latslt = 0.,12. fmap_satview_eradii = 10.0 map_mollweide = 1 ; MO projections (mollweide) fmap_mollweide_latlon = 0.,0. ; fmap_mollweide_latslt = -55.,12. map_continents = 1 ; continental outlines map_tn_unvn = 1 ; 1 int: add unvn to tn map_ht_unvn = 1 ; 1 int: add unvn to ht (z) map_ep_uivi = 1 ; 1 int: add uivi to epot map_log10 = 0 ; log10 of densities ;fmap_zpht = 50.,80.,100. ;fmap_zpht = -17.,-15.,-13.,-9.,-7.,-4.,+2. fmap_zpht = -17.,-15.,-7.,-4.,+2. ;fmap_zpht = -9.,-7.,-4. ;fmap_zpht = -17. ;------------------------------------------------------------------ ipltlat = 0 ;flats = -72.5,-62.5,-42.5,-17.4,-2.5,2.5,17.5,42.5,62.5,72.5 ;flats = -87.5,-42.5,-17.5,2.5,17.5,42.5,87.5 flats = 2.5,17.5,42.5,67.5 ;flats = 62.5,77.5 ;flat_zprange = -999.,999. ;flat_zprange = -17.,-1. ;flat_htscale = 30.,110.,2. ;flat_htscale = 30.,130.,2. flat_htscale = 50.,200.,2. ilat_yaxright = 1 ilat_log10 = 0 ;------------------------------------------------------------------ ipltlon = 0 flons = -180.,0.,'zm' ;flons = -180. ;flons = 'zm' ;flons = -90.,'zm' ;fslts = 0.,12. ;flon_zprange = -999.,999. ;flon_zprange = -17.,-1. ;flon_htscale = 30.,105.,2. flon_htscale = 50.,200.,2. ilon_yaxright = 1 ilon_log10 = 0 amphase = 0 istream = 0 ;------------------------------------------------------------------ ipltxyloc = 0 xylocs = 'zm','zm' xyloc_locname = 'global means' ; xyloc_zprange = -999.,999. ; xyloc_zprange = -17.,+5. xyloc_htscale = 50.,200.,2. ixyloc_log10 = 0 ixyloc_yaxright = 1 ;------------------------------------------------------------------ ipltutvert = 0 utvert_locs = 'zm','zm' ; utvert_zprange = -999.,999. utvert_htscale = 60.,120.,1. ; utvert_locname = 'Moscow' iutvert_log10 = 2 ;------------------------------------------------------------------ ipltutlat = 0 ; utlat_zphtlon = 100.,'LT0.', 100.,'LT6.', 100.,'LT12.', 100.,'LT18.' ; utlat_zphtlon = 100.,'LT0.' iutlat_log10 = 1 iutlat_sltmaps = 1 ;------------------------------------------------------------------ ipltxyut = 0 ; xyut_locname ='Slough', 'Moscow', 'Wallops Island','Christchurch', ; 'Hobart', 'Brisbane' ; xyut_locs = 51.5,359.4, 55.5,37.3, 37.8,284.5, -43.6,172.8 ; -42.9,147.3, -27.5,152.9 ; xyut_locname ='Moscow' ; xyut_locs = 55.5,37.3 xyut_locs = 'zm','zm' ; xyut_zpht = -4.,'hmf2',100.,200. xyut_zpht = 300.,'hmf2' ixyut_log10 = 0 ixyut_doppler = 0 ;------------------------------------------------------------------ sendcgm = 'ouray.hao:/e/roble/timegcm1.3/eqd.global.diag.d4.0' ;sendcgm = 'ouray.hao:/e/roble/timegcm1.3/eqd.uars.d98.0' ;sendcgm = 'ouray.hao:/e/roble/timegcm1.3/eqd.uars.utndcy.d85.0-86.0' ;sendcgm = 'ouray.hao:/e/roble/timegcm1.3/eqd.uars.diff.d60-50.0' ;sendcgm = 'ouray.hao:/e/roble/timegcm1.3/dsold.uars.d54.0' ; ; ;sendps = 'ouray.hao:/d/roble/tgcmproc/tgcm22.ps' ;sendcdf = 'ouray.hao:/d/roble/tgcmproc/tgcm22.nc' ;senddat = 'ouray.hao:/d/roble/tgcmproc/tgcm22.dat' ;sendxdr = 'ouray.hao:/d/roble/tgcmproc/tgcm22.xdr' ;sendv5d = 'ouray.hao:/d/roble/tgcmproc/tgcm22.v5d' /