# # QSUB -s /bin/csh # QSUB -eo # QSUB -lT 900 # QSUB -q prem # QSUB -nr # QSUB -lM 6MW # # Run tgcmproc cray post-model processor mss executable: # ja cd $TMPDIR setenv NCPUS 1 set job = tgcmproc batchname $job.$$ # output file to be returned to ntwk dir msread $job /FOSTER/bin/$job # get executable from mss chmod u+x $job msread $job.res /FOSTER/src/$job.res # get ncarg resource file from mss $job << 'EOF' || goto badx # execute with below inputs ; &input_mut E ; ; Inputs for multi-ut (iproc_mut) "station/location" processor: ; (see hao ~foster/tgcmvis/tgcmproc/doc/tgcmproc.doc for details on ; all valid input parameters) ; iproc_mut = 1 ; ;---------------------------------------------------------------------- ; histvols = '/ROBLE/RGR95/TSEQDR1','/ROBLE/RGR95/TSEQDR2', '/ROBLE/RGR95/TSEQDR3','/ROBLE/RGR95/TSEQDR4', '/ROBLE/RGR95/TSEQDR5','/ROBLE/RGR95/TSEQDR6', '/ROBLE/RGR95/TSEQDR7' ; ; 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 = 0,0,0,'to',1,0,0,'by',60 ; ;---------------------------------------------------------------------- ; ; Control history volumes and model times (for difference fields) ; ;histvols_cntr = '/ROBLE/RGR95/SECCMJ4' ;mtimes_cntr = 29,15,0 ; ;---------------------------------------------------------------------- ; ; 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 ; ; 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 ; ; 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 ; cfields = TN,UN,VN,N2,O2,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 = '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_integ = ; plot height-integration of emission fields ; E6300 E5577 EO200 EOH83 ECO215u ENO53u 1, 1, 1, 1, 1, 1 ; ; (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, 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,62,42,31 ; ; if ibohv_watts=0, then oh-v bands are in photons/cm3-sec (vol emis) ; and integrations are in kilo-rayleighs (brightness) ; if ibohv_watts=1, then oh-v bands are in watts/cm3 (vol emis), ; and integrations are in watts/cm2-str (brightness) ; (default is ibohv_watts=0) ibohv_watts = 0 ; ;---------------------------------------------------------------------- ; Plot options (multi-ut): ; icolor = 0 ; 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') ;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 ; ;---------------------------------------------------------------------- ; ; 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 = 41.9,-111.4 utvert_locname = 'BEAR LAKE' ;utvert_zprange = -17.,5. ; bot,top for y-axis zp utvert_htscale = 35.,120.,1. ; 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 = 0 ; ; 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 ; ;---------------------------------------------------------------------- ; ; ipltutlon = 0 -> contour ut vs longitude or local time at selected ; latitudes and zp or ht: ; ipltutlon = 1 ; ; utlon_zphtlat = pairs of zp or ht and lats (1st of pair may be zp or ht) utlon_zphtlat = 2.,52.5, 350.,52.5 iutlon_yslt = 1 ; slt on y-axis rather than lon iutlon_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 = 0 ; ; 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, doppler (idoppler=1 above), ; and/or oh-vib (ieohv above) plots are made) ; xyut_zpht = 80. ; ; ixyut_log10 = 0 -> linear axis, linear field ; ixyut_log10 = 1 -> linear axis, log10 field ; ixyut_log10 = 2 -> log10 axis, linear field ; ixyut_log10 = 0 ; ;---------------------------------------------------------------------- ; Send Files: ; sendcgm = machine:dir or machine:dir/file to which metafile is sent (rcp) ; sendps = machine:dir or machine:dir/file to which postscript is sent (rcp) ; senddat = machine:dir or machine:dir/file to which ascii data is sent (rcp) ; senddatms = mass store path to which ascii data is to be disposed (mswrite) ; sendxdr = machine:dir or machine:dir/file to which xdr data is sent (rcp) ; sendxdrms = mass store path to which xdr data is to be disposed (mswrite) ; sendcgm = 'cedar.hao:/d/login/tgcmproc/tgcmproc_mut.cgm' ; metafile ;sendps = 'cedar.hao:/d/login/tgcmproc/tgcmproc_mut.ps' ; postscript senddat = 'cedar.hao:/d/login/tgcmproc/tgcmproc_mut.dat' ; ascii data senddatms = '/USER/tgcmproc/tgcmproc_mut.dat' ; ascii on mss sendxdr = 'cedar.hao:/d/login/tgcmproc/tgcmproc_mut.xdr' ; xdr data sendxdrms = '/USER/tgcmproc/tgcmproc_mut.xdr' ; xdr to mss ; ;---------------------------------------------------------------------- ; &end 'EOF' ja -st exit badx: echo 'Execution failed' rcp $job.cgm cedar.hao:/d/login/$job/$job.fail.cgm exit