; &input_mut E ; ; Inputs for multi-ut (iproc_mut) "station/location" processor: ; iproc_mut = 1 ; ;---------------------------------------------------------------------- ; histvols = '/ROBLE/RGR94/TS10DS1','/ROBLE/RGR94/TS10DS2', '/ROBLE/RGR94/TS10DS3','/ROBLE/RGR94/TS10DS4', '/ROBLE/RGR94/TS10DS5','/ROBLE/RGR94/TS10DS6', '/ROBLE/RGR94/TS10DS7','/ROBLE/RGR94/TS10DS8', '/ROBLE/RGR94/TS10DS9','/ROBLE/RGR94/TS10DT1', '/ROBLE/RGR94/TS10DT2','/ROBLE/RGR94/TS10DT3', '/ROBLE/RGR94/TS10DT4','/ROBLE/RGR94/TS10DT5', '/ROBLE/RGR94/TS10DT6','/ROBLE/RGR94/TS10DT7', '/ROBLE/RGR94/TS10DT8','/ROBLE/RGR94/TS10DT9', '/ROBLE/RGR94/TS10DU8','/ROBLE/RGR94/TS10DU9' ; ; 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 = 20,0,0,'to',30,0,0,'by',60 ; ; 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 = 0 ; ; ionvel=1 -> ExB ion velocities')") ; ionvel=2 -> ExB+unvn ion velocities')") ; ionvel=3 -> E0xB+unvn ion velocities')") ; ionvel = 1 ; ; modelhts = 0 -> calculate heights at model grid from tn, mean mass, etc. ; modelhts = 1 -> use heights from model history ; modelhts = 1 ; ; Flag fields desired (ifields(i) = 1 -> plot corresponding field) ; ifields = ; TN UN VN O2 OX N4S NOZ CO CO2 H2O H2 HOX O+ CH4 O21D 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, ; NO2 NO O3 O1 OH HO2 H N2D TI TE NE O2+ W Z POTEN 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, ; UI VI WI N2 RHO UN+VN UI+VI FOF2 HMF2 0, 0, 0, 0, 0, 0, 0, 0, 0 ; ; fmnmxint (optional) = groups of 4: 'field', cmin, cmax, cint ; ('field' may be one of ifields as in above comments, or iemis as below) ; (if field = UN+VN or UI+VI then cmin,max,int = vlow,vhigh,vscale) ; ;fmnmxint = 'O1',0.,0.,.5, 'E5577',0.,0.,30. ; ; fscale (optional) = groups of 2: 'field', scalefac, where scalefac ; is to be scaled (multipled by) the field before contouring ; (note fmnmxint will be used *after* scaling by fscale) ; ;fscale = 'O1',1.e-11 ; ;---------------------------------------------------------------------- ; 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) ; iemis = ; emission field flags ; E6300 E5577 EO200 EOH83 ECO215u ENO53u 0, 0, 0, 0, 0, 0 ; iemis_integ = ; plot height-integration of emission fields 0, 0, 0, 0, 0, 0 ; ; (add SR63) ie6300 = 0 ; ; 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,0,0,0,0,0,0,0,0,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 = 0, 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. ; ;---------------------------------------------------------------------- ; Plot options (multi-ut): ; icolor = 0 ; make color fill contours if icolor=1 ibox_clabs = 0 ; box contour line labels ; ;---------------------------------------------------------------------- ; ; ipltutvert = 1 -> contour ut vs zp/ht at utvert_locs ; ('ZM' not allowed, but 'LTxx.x' is allowed) ipltutvert = 1 ; ; utvert_locs = lat,lon pairs of locations (lon can be 'LTxx.x') ; E.g.: utvert_locs = 0.,'LT12.' -> lat=0., local noon ; utvert_locs = 60.,-70. -> lat=60.,lon = -70. ; ;utvert_locs = 2.0,-157., 20.8,-156.5, 41.9,-111.4, 74.4,-94.5, ; 40.,-105.5, -44.,170.4, -36.,139. ;utvert_locname = 'CHRISTMAS', 'MAUI', 'BEAR LAKE', 'RESOLUTE BAY', ; 'FORT COLLINS','MTJOHNS','ADELAIDE' utvert_locs = 2.0,-157. utvert_locname = 'CHRISTMAS' ;utvert_zprange = -17.,5. ; bot,top for y-axis zp utvert_htscale = 35.,300.,5. ; bot,top,delta for y-axis ht iutvert_log10 = 0 ; log10 of species ; ;---------------------------------------------------------------------- ; ; ipltutlat = 1 -> contour ut vs latitude at selected longitudes ; ('ZM' not allowed, but 'LTxx.x' is allowed) ; ipltutlat = 1 ; ; utlat_zphtlon = pairs of zp or ht and lon (lon can be 'LTxx.x') ; (1st of pair may be zp or ht, 2nd of pair may be ; longitude or 'LTxx.x' for xx.x local time) utlat_zphtlon = 97.,'LT0.', -10.,180. iutlat_log10 = 0 ; log10 of species ; ;---------------------------------------------------------------------- ; ; ipltxyut = 1 -> xy line plot with ut on x-axis, field on y-axis ; at selected grid points and optionally zp/ht level. ; ('ZM' not allowed, but 'LTxx.x' is allowed) ipltxyut = 1 ; ; xyut_locs = pairs of floats: glat, glon (or slt) locations ; ;xyut_locs = 2.0,-157., 20.8,-156.5, 41.9,-111.4, 74.4,-94.5, ; 40.,-105.5, -44.,170.4, -36.,139. ;xyut_locs = 40.,-105.5, 0.,'LT12.' ;xyut_locs = 2.0,-157., 20.8,-156.5, 41.9,-111.4, 74.4,-94.5, ; 40.,-105.5, -44.,170.4, -36.,139. ; xyut_locs = 0.,'LT0.', 0.,'LT12.' ; ; xyut_zpht = pressures or heights for each xyut line plot ; (if not provided, only ht-integ emissions (iemis and iemis_integ ; above), doppler (idoppler=1 above), and/or oh-vib (ieohv above) ; plots are made) ; xyut_zpht = 80. ; ;---------------------------------------------------------------------- ; Send Files: ; sendcgm = machine:dir or machine:dir/file to which metafile is sent (rcp) ; senddat = machine:dir or machine:dir/file to which ascii data is sent (rcp) ; ;sendcgm = 'euchre.hao:/d/palo/tgcmproc/tgcmproc_mut.cgm' ; metafile ;senddat = 'euchre.hao:/d/palo/tgcmproc/tgcmproc_mut.dat' ; ascii data sendcgm = 'vishnu.hao:/d/foster/tgcmvis/tgcmproc/tgcmproc_mut.cgm' ;senddat = 'euchre.hao:/d/palo/tgcmproc/tgcmproc_mut.dat' ; ascii data ; ;---------------------------------------------------------------------- ; &end