&proc_input ;tmpdir = '/hao/tgcm/tiegcm1.6_yr02/sy2r' ;tmpdir = '/ptmp/roble/tgcmproc' tmpdir = '/hao/roble/timegcm1.2.ncep' ; ; histvols = '/hao/tgcm/timegcm1.2_yr04/py4n/ROBLE.timegcm1.2.py4n001.nc','to' ; '/hao/tgcm/timegcm1.2_yr04/py4n/ROBLE.timegcm1.2.py4n037.nc','by','1' ; mtimes = 1,0,0,'to', 365,0,0,'by',1440 ; ; histvols = '/hao/tgcm/timegcm1.2_yr04/sy4n/ROBLE.timegcm1.2.sy4n001.nc','to' ; '/hao/tgcm/timegcm1.2_yr04/sy4n/ROBLE.timegcm1.2.sy4n194.nc','by','1' ; mtimes = 1,15,0,'to',364,15,0,'by',1440 ; mtimes = 1,7,0,'to',364,7,0,'by',1440 ; mtimes = 1,13,0,'to',364,13,0,'by',1440 ; mtimes = 1,5,0,'to',364,5,0,'by',1440 ; ; mtimes = 1,0,0,'to',5,0,0,'by',60 ; histvols = '/ROBLE/tgcm24rh/py3n001.nc','to', ; '/ROBLE/tgcm24rh/py3n037.nc','by','1' ; mtimes = 1,0,0,'to', 365,0,0,'by',1440 ; ; ; ; ; histvols = '/ROBLE/timegcm1.2/py2n001.nc','to', ; '/ROBLE/timegcm1.2/py2n037.nc','by','1' ; mtimes = 1,0,0,'to', 365,0,0,'by',1440 ; ; histvols = '/ROBLE/timegcm1.2/py3m030.nc','to', ; '/ROBLE/timegcm1.2/py3m037.nc','by','1' ; mtimes = 300,0,0,'to', 310,0,0,'by',1440 histvols = '/ROBLE/timegcm1.2/sy3m156.nc','to', '/ROBLE/timegcm1.2/sy3m159.nc','by','1' mtimes = 300,0,0,'to', 305,0,0,'by',60 ; histvols = '/ROBLE/timegcm1.2/py5n001.nc','to', ; '/ROBLE/timegcm1.2/py5n037.nc','by','1' ; mtimes = 1,0,0,'to', 100,0,0,'by',1440 ;---------------------------------------------------------------------- ;---------------------------------------------------------------------- ;------------------------------------------------------------------ ; ; 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) iden = 1 idifpercent = 0 ; ; ionvel=1 -> ExB ion velocities')") ; ionvel=2 -> ExB+unvn ion velocities')") ; ionvel=3 -> E0xB+unvn ion velocities')") ionvel = 2 modelhts = 1 ; ; 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,O1,O3,OH,H,E5577,EO200,OHB-62 cfields = TN,UN,VN,OHB-62,EO200,E5577,O2,OH,O1 ; ; fmnmxint (optional) = groups of 4: 'field', cmin, cmax, cint ; ('field' may be one of cfields as in above comments) ; (if field = UN+VN or UI+VI then cmin,max,int = vlow,vhigh,vscale) ; ;fmnmxint = 'UN+VN',4.,320.,250., 'UI+VI',4.,320.,250., ; 'TN',480.,1640.,40., 'POTEN',-18000.,14000.,2000., ; 'WI',-65.,35.,5., 'FOF2',1.,20.,1., 'HMF2',100.,675.,25., ; 'UN',-140.,220.,20., 'NE',3.4,6.6,.2 ;fmnmxint = 'UN+VN',.5,350.,250., 'TN',185.,800.,10. ; ;---------------------------------------------------------------------- ; Volume Emissions (airglow) fields: ; (note: these fields are calculated only from time-gcm histories) ; (note: iden must be 1 (cm3 species) for these calculations) ; (note: solred is called if ie6300=1, or ie5577(3,4, or 5)=1 ; (solred is cpu intensive over several histories) ; iemis_integ = ; height-integration of emission fields (ipltmaps>0) ; E6300 E5577 EO200 EOH83 ECO215u ENO53u 0, 1, 1, 0, 0, 0 ; ; (add SR63 via solred) ie6300 = 0 ; ; idoppler=1 -> generate doppler t,u,v for each emissions field except ; ECO215u and ENO53u. (ipltxyut type line plots only) idoppler = 1 ; ; 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, will use 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 ; ; iyd,f107a,f107d needed only if ie5577(3 or 4) > 0, i.e., photoe or airglow ; (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,1,0,0,0,1,0,0,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 = 62 ; ;---------------------------------------------------------------------- ; Plot options: ; icolor = 1 ; make color fill contours if icolor=1 ibox_clabs = 0 ; box contour line labels outplt = 'cgm','ps' ; plot output type(s) ('cgm' and/or 'ps') ;outplt = 'ps' ; plot output type(s) ('cgm' and/or 'ps') ;outplt = 'cgm' ; plot output type(s) ('cgm' and/or 'ps') ;psmode = 'port' ; for ps: 'port' for portrait or 'land' for landscape ;psmode = 'land' ; for ps: 'port' for portrait or 'land' for landscape ; ; multiplt = 1 -> put multiple plots on each frame ; ipltrowcol = nrows,ncols of plots on each frame (multiplt > 0 only) ; multiadvfr = 0 -> advance frame *only* when page is full according to ; ipltrowcol. ; multiadvfr = 1 -> advance frame with full page, *and* between plot ; types (maps, slices, xy plots) and ut's ; multiplt = 0 ; multiple plots per frame (0 or 1 toggle) ipltrowcol = 2,2 ; # rows,cols of plots per frame (multiplt > 0) multiadvfr = 1 ; frame adv flag for multiplt ; ;---------------------------------------------------------------------- ipltmaps=0 ; plot flag for maps map_global=0 ; CE projections (cyl equidistant) map_global_cenlon=90 ; center x-axis longitude map_global_censlt=12 ; center x-axis local time map_polar=0 ; ST projections (polar stereographic) fmap_polar_perimlat = -40.,40. map_satview=0 ; SV projections (sat view) ;---------------------------------------------------------------------- ; ; ipltutvert = 1 -> contour ut vs zp/ht at utvert_locs ; ('ZM' not allowed, but 'LTxx.xx' is allowed) ; ipltutvert = 1 ; ; utvert_locs = lat,lon pairs of locations (lon can be 'LTxx.xx') ; E.g.: utvert_locs = 0.,'LT12.' -> lat=0., local noon ; utvert_locs = 60.,-70. -> lat=60.,lon = -70. ; ; utvert_locs = -35.,'LT00.',-23.5,'LT00.', 18.3,'LT00.', ; 22.5,'LT00.', 35.,'LT00.', 42.5,'LT00.' ; utvert_locs = 42.87,278.62, 42.43,25.62, 37.2,356.7, 29.16,278.89 utvert_locname = 'DELAWARE', 'STARA ZAGORA', 'SIERRA NEVADA', 'DAYTONA BEACH' ; utvert_locs = -35.,138., -23.5,133., 18.3,-66.8, 22.5,-156., ; 35.,-106., 42.5,-105. ; utvert_locname = 'ADELAIDE', 'ALICE SPRINGS', 'ARECIBO', 'MAUI', ; 'ALBUQUERQUE', 'FORT COLLINS' ; ; utvert_locname = 'ADELAIDE', 'ALICE SPRINGS' ; utvert_locs = -35.,138., -23.5,133. ; ; utvert_locname = 'ALBUQUERQUE','FORT COLLINS' ; utvert_locs = 35.,-106., 42.5,-105. ; ; utvert_locname = 'MAUI' ; utvert_locs = 22.5,-156. ; ; utvert_locname = 'ARECIBO' ; utvert_locs = 18.3,-66.8 ; ; utvert_locs = 52.,-1., -52.,-58. ; utvert_locname = 'CHILTON', 'PORT STANLEY' ; ;utvert_zprange = -7.,+7. ; bot,top for y-axis zp utvert_htscale = 70.,105.,1. ; bot,top,delta for y-axis ht ;utvert_htscale = 100.,400.,5. ; bot,top,delta for y-axis ht iutvert_log10 = 0 ; log10 of species ; ;---------------------------------------------------------------------- ; ; if ipltutlat >= 1 -> contour ut vs latitude at selected longitudes ; ('ZM' not allowed, but 'LTxx.xx' is allowed) ; ipltutlat = 0 ; ; utlat_zphtlon = pairs of zp or ht and lon (lon can be 'LTxx.xx') ; (1st of pair may be zp or ht, 2nd of pair may be ; longitude or 'LTxx.xx' for xx.xx local time) ; utlat_zphtlon = 85.,0., 96., 0., 111.,0. ; utlat_zphtlon = 90.,0., 95.,0., 100.,0. , 111.,0. ; utlat_zphtlon = 85.,-156., 90.,-156., 95.,-156. ; utlat_zphtlon = 85.,-105., 90.,-105., 95.,-105. ; utlat_zphtlon = 85.,-111., 90.,-111., 95.,-111. ; utlat_zphtlon = 300.,0.,120.,0. ; utlat_zphtlon = 50.,'LT20.',70.,'LT20.',85.,'LT20.',95.,'LT20.', ; 110.,'LT20.',200.,'LT20.' ; utlat_zphtlon = 97.,'LT00.' ;iutlat_log10 = 0 ; log10 of species ; ; iutlat_sltmaps: Integer flag, 0 or 1. If iutlat_sltmaps = 1, then ; processor will make map projections at fixed local times given ; by utlat_zphtlon (e.g., utlat_zphtlon = -4.,'LT12.'). These maps ; are made only when a full day, no more no less, is being processed. ; Map projections are made based on map_global, map_polar, map_satview, ; and map_mollweide (regardless of value of ipltmaps flag). ; (default iutlat_sltmaps = 0) ;iutlat_sltmaps = 0 ; ;---------------------------------------------------------------------- ; ; if ipltxyut >= 1 -> xy line plot with ut on x-axis, field on y-axis ; at selected grid points ; ('ZM' not allowed, but 'LTxx.xx' is allowed) ipltxyut = 1 ; ; xyut_locs = pairs of glat, glon (or slt) pairs of locations ; ; xyut_locs = -35.,'LT00.',-23.5,'LT00.', 18.3,'LT00.', ; 22.5,'LT00.', 35.,'LT00.', 42.5,'LT00.' ; xyut_locs = -35.,138., -23.5,133., 18.3,-66.8, 22.5,-156., ; 35.,-106., 42.5,-105. ; xyut_locname = 'ADELAIDE', 'ALICE SPRINGS', 'ARECIBO', 'MAUI', ; 'ALBUQUERQUE', 'FORT COLLINS' ; ; xyut_locname = 'ADELAIDE', 'ALICE SPRINGS' ; xyut_locs = -35.,138., -23.5,133. ; ; xyut_locname = 'ALBUQUERQUE','FORT COLLINS' ; xyut_locs = 35.,-106., 42.5,-105. ; ; xyut_locname = 'MAUI' ; xyut_locs = 22.5,-156. ; ; xyut_locname = 'ARECIBO' ; xyut_locs = 18.3,-66.8 ; ; xyut_locs = 52.,-1., -52.,-58., 56.,37., 45.,142., ; 39.,-77., 30.5,-83.5, -43.,147. ; xyut_locname = 'CHILTON', 'PORT STANLEY', 'MOSCOW', 'WAKKANAI', ; 'WALLOPS', 'EGLIN', 'HOBART' ; xyut_locs = 42.87,278.62, 42.43,25.62, 37.2,356.7, 29.16,278.89 xyut_locname = 'DELAWARE', 'STARA ZAGORA', 'SIERRA NEVADA', 'DAYTONA BEACH' ; ; xyut_zpht = pressures or heights for each xyut line plot ; (if not provided, only ht-integ emissions (iemis_integ ; above) and doppler (idoppler=1 above) plots are made) ; ;xyut_zpht = 90. ;xyut_zpht = -17.,-15.,-13.,-9.,-7.,-4.,+2. ; ; ixyut_log10 = 0 -> linear axis, linear field ; ixyut_log10 = 1 -> linear axis, log10 field ; ixyut_log10 = 2 -> log10 axis, linear field ; ixyut_log10 = 0 ; ; ;---------------------------------------------------------------------- ipltxyloc = 0 xylocs = 'zm','zm' xyloc_locname = 'global means' ; xyloc_zprange = -999.,999. xyloc_zprange = -7.,+7. ; xyloc_htscale = 30.,105.,2. ixyloc_log10 = 1 ixyloc_yaxright = 1 ;------------------------------------------------------------------ ; Send Files: ; ; sendcgm = 'ouray.hao:/e/roble/marianna/yr03.sta.marianna.d300-305hr' sendps = 'ouray.hao:/e/roble/marianna/yr03.sta.marianna.d300-305hr.ps' ; sendps = 'ouray.hao:/e/roble/marianna/yr03.sta.marianna.d300-365d.ps' ;---------------------------------------------------------------------- ; /