&proc_input ;tmpdir = '/ptmp/roble/tgcmproc' ; ; ; ;**************************************************** ; ; ; histvols = 'ROBLE.timegcm1.42d.sy92ncepnmt.003.nc' mtimes = 4,0,0 ; ; ; ; histvols = '/ROBLE/timegcm1.42d/s96jandnmtb007.nc','to', ; '/ROBLE/timegcm1.42d/s96jandnmtb009.nc','by','1' ; ; histvols = '/ROBLE/timegcm1.42d/s96eqdnmt007.nc','to', ; '/ROBLE/timegcm1.42d/s96eqdnmt009.nc','by','1' ; mtimes = 39,0,0,'to', 40,0,0,'by',1440 ; mtimes = 29,0,0,'to', 30,0,0,'by',60 ; ; ; ;------------------------------------------------------------------ ;------------------------------------------------------------------ ; maxdiskvols = 50 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 ' ; ; ; RAW-MESOX ;cfields = TN,UN,VN,WN,Z,E5577,OHB-83,O1 ; cfields = TN,UN,VN,WN,O1,E5577,NE,E6300 ;cfields = Z cfields = E5577 ; ;cfields = TN,UN,VN,WN,O3,O1,NO,NO2, ; H2O,OH,HO2,H,CO2,CO,CH4,NE, ; QSOLAR,QJOULE,CHEMHEAT, ; RADCOOL,CO2_DEG,NO_DEG,O3_DEG, ; UTNDCY,VTNDCY,DIFKK,O2,H2 ; ; cfields = TN,UN,VN,O1,O2,NO,NO2,CO,CO2,H,OH,HO2,CH4,H2, ; E5577,E0200,H2O,O/O2+N2 ; ; PERCENT-MESOX ; cfields = TN,O1,O2,NO,NO2,CO,CO2,H,OH,HO2,CH4,H2, ; E5577,E0200,H2O,O/O2+N2,TNFP ; ; ;fscale = 'O1',1.E-11 ;fmnmxint = 'TN',1.,0.,10., 'UN',1.,0.,10.,'O1',1.,0.,0.5, ; 'E5577',1.,0.,50. ; idifpercent = 0 makes raw differnces ;idifpercent = 0 idifpercent = 1 iden = 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 = 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, 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, 1, 1, 1, 1 ; ; idoppler=1 -> generate doppler t,u,v for each emissions field except ; ECO215u and ENO53u. (ipltxyut type line plots only) idoppler = 1 ; ; 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 ; 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 = 0 ; 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 = 35. map_satview = 0 ; SV projections (sat view) fmap_satview_latlon = -60.,-90. ; fmap_satview_latslt = -55.,12. fmap_satview_eradii = 10.0 map_mollweide = 0 ; MO projections (mollweide) fmap_mollweide_latlon = 65.,0. ; fmap_mollweide_latslt = -55.,12. map_continents = 0 ; 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 = 60.,80.,95.,100.,110. ;fmap_zpht = -17.,-15.,-9.,-7. fmap_zpht = -17. ;------------------------------------------------------------------ ipltxyut = 0 ; xyut_locname ='Bear Lake', 'Haute Provence', 'Stockholm', ; 'Fritz Peak','Mt.Johns' ; xyut_locs = 41.9,-111.4, 41.0, 6.0, 60.0, 20.0, ; 40.5,-105.5, -44.0,170.4 ; xyut_locname = '2.5,slt12' xyut_locs = 2.5,'slt12', 40.,'slt12', 65.,'slt12' ixyut_log10 = 0 xyut_zpht = 95.,250. ;------------------------------------------------------------------ ;------------------------------------------------------------------ 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.,-5. flat_htscale = 30.,105.,2. ;flat_htscale = 30.,130.,2. ;flat_htscale = 30.,200.,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 = -9.,-5. flon_htscale = 70.,105.,1. ;flon_htscale = 50.,200.,2. ;flon_htscale = 30.,200.,2. ilon_yaxright = 1 ilon_log10 = 0 flon_xlatrange = -20.,20. ;amphase = 4 ;istream = 1 ;------------------------------------------------------------------ ipltutlat = 0 iutlat_log10=0 ; log10 flag ; utlat_zphtlon =90.,'zm',100.,'zm',110.,'zm',120.,'zm', ; 200.,'zm',300.,'zm' ; utlat_zphtlon = 35.,'zm',50.,'zm',60.,'zm',70.,'zm',80.,'zm', ; 90.,'zm',100.,'zm',120.,'zm',300.,'zm' ; utlat_zphtlon = -17.,'zm',-15.,'zm',-13.,'zm',-10.,'zm',-9.,'zm', ; -7.,'zm',-4.,'zm' ; utlat_zphtlon = -7.,'zm',-4.,'zm',-1.,'zm',+2.,'zm',+5.,'zm' ; utlat_zphtlon = 100.,-150., 120.,-150., 140.,-150., 180.,-150., ; 400.,-150. ; utlat_zphtlon = 110., -180.,350., -180., 110.,0., 350.,0. ; utlat_zphtlon = 110.,'zm', 150.,'zm', 350.,'zm' ,500.,'zm' utlat_zphtlon = -9.,'zm', -7.,'zm', -4.,'zm', -1.,'zm' , +2.,'zm',+5.,'zm' ; ;------------------------------------------------------------------ ipltutvert = 0 utvert_locs = 'zm','zm' ; utvert_zprange = -9.,+5. utvert_zprange = -17.,-5. ; utvert_htscale = 80.,400.,1. utvert_htscale = 30.,105.,1. ; utvert_htscale = 30.,100.,1. utvert_locname = 'GLOBAL MEAN' iutvert_log10 = 0 ;------------------------------------------------------------------ ipltxyloc = 1 ; xylocs = 'zm','zm' xylocs = 1.25,'SLT12.00.' ; xyloc_locname = 'global means' xyloc_locname = '1.0','slt00.' ; xyloc_zprange = -999.,999. ; xyloc_zprange = -17.,-5. xyloc_htscale = 70.,250.,1. ixyloc_log10 = 0 ixyloc_yaxright = 1 ;------------------------------------------------------------------ ;------------------------------------------------------------------ ;------------------------------------------------------------------ ; sendcgm = 'sofie.hao.ucar.edu:/sofie/d/roble/timegcm1.42d.gordonss/s96eqd.profile.airglow.d30' ; flnm_cgm = '/hao/sofie2/roble/timegcm1.42d.gordon.bright/jan92.profile.green.d4' flnm_ps = '/hao/sofie2/roble/timegcm1.42d.gordon.bright/jan92.profile.green.d4.ps' ; /