# # 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_sut E ; ;---------------------------------------------------------------------- ; ; Inputs for single-ut (iproc_sut) "snapshot" processor: ; (see hao ~foster/tgcmvis/tgcmproc/doc/tgcmproc.doc for details on ; all valid input parameters) ; iproc_sut = 1 ; ;---------------------------------------------------------------------- ; ; History and fields: ; histvols = mss paths to history volume(s) ; mtimes = integer triplets specifying model day,hour,min ; histvols = '/ROBLE/RGR95/TE7ESN1','/ROBLE/RGR95/TE7ESX2' ; ; mtimes = integer triplets, where each triplet represents model ; day, hour, and minute. The processor will sequentially ; search histvols for mtimes histories. ; mtimes = 6,0,0, 8,0,0 ; ;---------------------------------------------------------------------- ; ; 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 = '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. ;fmnmxint = 'TN',190.,255.,5. ;fmnmxint = 'TN',0.,0.,5. ;fmnmxint = 'TN',230.,240.,2. ; ;---------------------------------------------------------------------- ; Volume Emissions (airglow) fields: ; (note: these fields are calculated only from time-gcm histories) ; (note: iden must be 1 (cm3 species) for these calculations) ; iemis_integ = ; height-integration of emission fields (ipltmaps>0) ; E6300 E5577 EO200 EOH83 ECO215u ENO53u 1, 1, 1, 1, 1, 1 ; ; (add SR63) ie6300 = 1 ; (o1 recomb) (o2+ recomb) (photoe) (airglow) (o2 ly-beta) ie5577 = 1, 1, 1, 1, 1 ; ; 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 type ipltxyloc only, zonal and global means not ; available. (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 ; ; 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: ; 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 ; ;---------------------------------------------------------------------- ; Maps: ; (contour over map projections if ipltmaps=1) ; (if ipltmaps=0, remaining map parameters are ignored) ; ipltmaps = 1 ; 1 int: make map projection1 ; fmap_zpht = -10.,-9.,-8.,-7.,-6. ; n floats: zp or ht surfaces ; fmap_zpht = 80.,85.,90.,95.,100. ; n floats: zp or ht surfaces ; fmap_zpht = -16.,-15.5,-11.,-4.,2. fmap_zpht = -7.,7. map_continents = 0 ; 1 int: continent outlines map_tn_unvn = 1 ; 1 int: add unvn to tn map_ht_unvn = 0 ; 1 int: add unvn to z map_ep_uivi = 1 ; 1 int: add uivi to epot ivec_label = 0 ; 1 int: label vec mags ; ; Map projections and related options: ; map_global = 1 ; 1 int: CE projection map_global_cenlon = 0 ; 1 int: center longitude ; map_global_censlt = 12 ; 1 int: center local time map_polar = 0 ; 1 int: ST projection fmap_polar_perimlat = 45. ; n floats: perimeter lat(s) map_satview = 0 ; 1 int: SATV projection fmap_satview_latlon = -55.,0. ; 2 floats: center of satv ; fmap_satview_eradii = 1.5 ; 1 float: earth radii dist ; map_top_anno = 'MAP ANNOTATION' ; top label on plots ; ;---------------------------------------------------------------------- ; Longitude slices (contour latitude vs zp or ht at selected longitudes) ; ipltlon = 1 ; 1 intzp/ht vs lat at flon(s) ; flons = 0. ; n floats: lons for ipltlon flons = 'zm' ; fslts = 0.,6.,12.,18. ; n floats: local times for ipltlon fslts = 0.,12. ; n floats: local times for ipltlon flon_zprange = -99.,99. ; 2 floats: zp range for y-axis ; flon_htscale = 80.,130.,1. ; 3 floats: ht range, delta for y-axis ilon_log10 = 0 ; 1 int: plot log10 of species ; lon_top_anno = 'LON SLICE ANNOTATION'; top label on plots ; ;---------------------------------------------------------------------- ; Latitude slices (contour longitude vs zp or ht at selected latitudes) ; ipltlat = 1 ; flats = -90.,-60.,0.,60.,90. flats = -60.,0.,60. ; flats = 0.,42.5,72.5 ; n floats: latitudes for ipltlat flat_zprange = -99.,99. ; 2 floats: zp range for y-axis ; flat_htscale = 80.,130.,1. ; 3 floats: ht range, delta for y-axis ilat_log10 = 0 ; 1 int: plot log10 of species ; lat_top_anno = 'TOP ANNOTATION FOR LAT SLICES' ; ;---------------------------------------------------------------------- ; XY vertical profile plots at specified locations (field vs zp): ; ipltxyloc = 0 ; xylocs = 0.,'ZM', 40.,-70., 'zm','zm' xylocs = 67.5,-150. ; xylocs = 0.,'LT0.', 15.,'LT0.', 30.,'LT0.', 45.,'LT0.', 60.,'LT0.', ; -15.,'LT0.', -30.,'LT0.', -45.,'LT0.', -60.,'LT0.' xyloc_zprange = -99.,99. ; 2 floats: zp range for y-axis xyloc_htscale = 75.,100.,1. ; 3 floats: ht range, delta for y-axis ; ; ixyloc_log10 = 0 -> linear axis, linear field ; ixyloc_log10 = 1 -> linear axis, log10 field ; ixyloc_log10 = 2 -> log10 axis, linear field ; ixyloc_log10 = 2 ; log10 flag ; xyloc_top_anno = 'TOP ANNOTATION FOR XY LOC' ; ;---------------------------------------------------------------------- ; Send output data 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_sut.cgm' ; metafile ;sendps = 'cedar.hao:/d/login/tgcmproc/tgcmproc_sut.ps' ; postscript senddat = 'cedar.hao:/d/login/tgcmproc/tgcmproc_sut.dat' ; ascii data senddatms = '/USER/tgcmproc/tgcmproc_sut.dat' ; ascii on mss sendxdr = 'cedar.hao:/d/login/tgcmproc/tgcmproc_sut.xdr' ; xdr data sendxdrms = '/USER/tgcmproc/tgcmproc_sut.xdr' ; xdr to mss ; ;---------------------------------------------------------------------- ; End inputs for single-ut processor ; &end_sut 'EOF' ja -st exit badx: echo 'Execution failed' rcp $job.cgm cedar.hao:/d/login/$job/$job.fail.cgm exit