Huang: These are the mss paths to the history volumes for the cirris tgcm run: ; Primary volumes Mday:ut Day Date ; ; RGR97/TSCIR01 27:00->28:12 by 12hrs 91116 4/26/91 ; RGR97/TSCIR02 29:00->30:12 by 12hrs ; RGR97/TSCIR03 31:00->31:12 by 12hrs ; ; ; Secondary volumes Mday:ut Day Date ; ; RGR97/SECIR01 27:00->27:07 by 1 hrs 91116 4/26/91 ; RGR97/SECIR02 27:08->27:15 by 1 hrs ; RGR97/SECIR03 27:16->27:23 by 1 hrs ; RGR97/SECIR04 28:00->28:07 by 1 hrs ; RGR97/SECIR05 28:08->28:15 by 1 hrs ; RGR97/SECIR06 28:16->28:23 by 1 hrs ; RGR97/SECIR07 29:00->29:07 by 1 hrs ; RGR97/SECIR08 29:08->29:15 by 1 hrs ; RGR97/SECIR09 29:16->29:23 by 1 hrs ; RGR97/SECIR10 30:00->30:07 by 1 hrs ; RGR97/SECIR11 30:08->30:15 by 1 hrs ; RGR97/SECIR12 30:16->30:23 by 1 hrs ; RGR97/SECIR13 31:00->31:07 by 1 hrs ; RGR97/SECIR14 31:08->31:12 by 1 hrs I made a short run of tgcmproc for history 27,3,0 (volume /ROBLE/RGR97/SECIR01, 4/26). There are vertical profiles on hao.hao.ucar.edu anon ftp area in the dir /pub/foster, file names huang.ps (postscript line plots), and huang.dat (ascii data corresponding to the plots). These are for O1,O2,N2,NO,O2+,NE,TN,TE. NO+, N+, and N2+ are not saved on the histories at all (they are calculated in the model, but not saved on the histories). N4S, O+ and N2D are on the primary histories, but NOT on this series of secondary histories (as you can see above, the primaries were saved only every 12 hours, whereas secondary histories are hourly). I will copy the tgcmproc input script below. These are cray jobs (tgcmproc runs only on the crays). The above profiles are at the full log pressure range of the model (zp -17 to +5), but you can change that to a height scale by setting xyloc_htscale, e.g.: xyloc_htscale = 75.,100.,1. ; 3 floats: ht range, delta for y-axis BTW, I wasnt sure whether the 169 deg lon you mentioned was east or west. We use -180 (west) to +180 (east). You can find out more about tgcmproc at: http://www.hao.ucar.edu/public/research/tiso/tgcm/proc/tgcmproc.html --Ben # # 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 iproc_sut = 1 ; ; Timegcm (tgcm20) run for CIRRIS comparisons (made Nov-Dec, 1996): ; (all mss paths are under root dir /ROBLE). The model was not ; advanced by calendar day, but time varying inputs were used ; for f107, hem power, byimf, and ctpoten, as below. ; ; Primary volumes Mday:ut Day Date ; ; RGR97/TSCIR01 27:00->28:12 by 12hrs 91116 4/26/91 ; RGR97/TSCIR02 29:00->30:12 by 12hrs ; RGR97/TSCIR03 31:00->31:12 by 12hrs ; ; ; Secondary volumes Mday:ut Day Date ; ; RGR97/SECIR01 27:00->27:07 by 1 hrs 91116 4/26/91 ; RGR97/SECIR02 27:08->27:15 by 1 hrs ; RGR97/SECIR03 27:16->27:23 by 1 hrs ; RGR97/SECIR04 28:00->28:07 by 1 hrs ; RGR97/SECIR05 28:08->28:15 by 1 hrs ; RGR97/SECIR06 28:16->28:23 by 1 hrs ; RGR97/SECIR07 29:00->29:07 by 1 hrs ; RGR97/SECIR08 29:08->29:15 by 1 hrs ; RGR97/SECIR09 29:16->29:23 by 1 hrs ; RGR97/SECIR10 30:00->30:07 by 1 hrs ; RGR97/SECIR11 30:08->30:15 by 1 hrs ; RGR97/SECIR12 30:16->30:23 by 1 hrs ; RGR97/SECIR13 31:00->31:07 by 1 hrs ; RGR97/SECIR14 31:08->31:12 by 1 hrs ; histvols = '/ROBLE/RGR97/SECIR01', mtimes = 27,3,0 ; ; 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) ; iden = 1 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 ; NAT 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 ; NaS NaO NaO3 NaO2 NaOH NaCO3 NaHCO3 NaS+ ; NaN2+ NaCO2+ NaH2O+ NaO+ NaEMIS ; E6300 E5577 EO200 EOH83 ECO215u ENO53u ; cfields = 'O1','O2','N2','NO','N4S','N2D','O2+','O+','NE','TN','TE' ; iemis_integ = ; height-integration of emission fields (ipltmaps>0) 1, 1, 1, 1, 1, 1 ; (add SR63) ie6300 = 0 ; (o1 recomb) (o2+ recomb) (photoe) (airglow) (o2 ly-beta) ie5577 = 1, 0, 0, 0, 0 ;iyd = 94080 ;f107d = 67. ;f107a = 72. ;ieohv = 0,0,0,1,1,1,1,1,1,1 ;ieohv = 0,1,0,0,0,0,0,0,0,0 ieohv = 0,0,0,0,0,0,0,0,0,0 ; icolor = 0 ; make color fill contours if icolor=1 ishadeneg = 1 ; softfill (shade) negative contour areas (icolor=0) ibox_clabs = 0 ; box contour line labels ;iboxplt = 1 ; draw perimeter box around each frame 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 = 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 = 0 ; 1 int: make map projection1 ; fmap_zpht = -16.,-13.,-10.,-7.,-4. ; fmap_zpht = -10.,-8.,-6.,-4. ; fmap_zpht = 90. fmap_zpht = -4. map_continents = 1 ; 1 int: continent outlines map_tn_unvn = 1 ; 1 int: add unvn to tn map_ht_unvn = 0 ; 1 int: add unvn to ht (z) map_ep_uivi = 0 ; 1 int: add uivi to epot ivec_label = 0 ; 1 int: label vec mags ; map_top_anno = 'TGCM/T21 COUPLED RUN CASE flxT21dy' ; ; Map projections and related options: ; map_global = 0 ; 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 = 30. ; n floats: perimeter lat(s) map_satview = 0 ; 1 int: SATV projection ; fmap_satview_latlon = -55.,0. ; 2 floats: center of satv fmap_satview_latslt = -55.,12. ; 2 floats: center of satv ; fmap_satview_eradii = 1.5 ; 1 float: earth radii dist map_mollweide = 1 fmap_mollweide_latlon = 0.,-90. ; 2 floats: center of mollw ; fmap_mollweide_latslt = 0.,12. ; 2 floats: center of mollw ; ;---------------------------------------------------------------------- ; Longitude slices (contour latitude vs zp or ht at selected longitudes) ; ipltlon = 0 ; 1 intzp/ht vs lat at flon(s) ; flons = -180.,'zm' ; flons = -180. flons = 'zm' ; fslts = 0.,6.,12.,18. ; fslts = 0.,3.,6.,9.,12.,15. ; n floats: local times for ipltlon ; fslts = 0.,6.,12.,18. ; n floats: local times for ipltlon flon_zprange = -10.,-2. ; 2 floats: zp range for y-axis ; flon_htscale = 60.,120.,2. ; 3 floats: ht range, delta for y-axis ; flon_xlatrange = -60.,60. ; 2 floats: range of lat on x-axis ilon_log10 = 0 ; 1 int: plot log10 of species ilon_yaxright = 2 ; extra right-hand y-axes ; lon_top_anno = 'LON SLICE ANNOTATION' ; ; Set amphase > 0 (max 4) to plot longitude slices of amplitudes and ; phases. Wave numbers 1 to amphase are plotted for each field. ; (ipltlon need not be > 0, but flon_zprange and/or flon_htscale, ; and flon_xlatrange are used when plotting amphase) amphase = 0 ; ; Set istream = 1 to plot stream functions of zonal mean VN and/or W. ; (ipltlon need not be > 0, but VN and/or W must be in cfields) ; (flon_zprange and/or flon_htscale and flon_xlatrange are used ; when plotting stream functions) istream = 0 ; ;---------------------------------------------------------------------- ; Latitude slices (contour longitude vs zp or ht at selected latitudes) ; ipltlat = 0 ; flats = -90.,-60.,-40.,-20.,0.,20.,40.,60.,90. ; flats = -27.5,27.5,-42.5,42.5 flats = 2.5,17.5,42.5 ; n floats: latitudes for ipltlat ; flat_zprange = -99.,99. ; 2 floats: zp range for y-axis flat_htscale = 75.,100.,1. ; 3 floats: ht range, delta for y-axis ilat_log10 = 0 ; 1 int: plot log10 of species ilat_yaxright = 2 ; extra right-hand y-axes ; lat_top_anno = 'TOP ANNOTATION FOR LAT SLICES' ; ;---------------------------------------------------------------------- ; XY vertical profile plots at specified locations (field vs zp): ; ipltxyloc = 1 xylocs = 55.,169. ; ; xyloc_zprange = -99.,99. ; 2 floats: zp range for y-axis xyloc_zprange = -99.,5. ; 2 floats: zp range for y-axis ; xyloc_htscale = 75.,100.,1. ; 3 floats: ht range, delta for y-axis ; xyloc_htscale = 30.,500.,10. ; 3 floats: ht range, delta for y-axis ixyloc_yaxright = 2 ; extra right-hand y-axes ; ; 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' ; ;---------------------------------------------------------------------- ; ;sendcgm = 'vishnu.hao:/d/foster/tgcmvis/tgcmproc/cirris/huang.cgm' ;senddat = 'vishnu.hao:/d/foster/tgcmvis/tgcmproc/cirris/huang.dat' sendcgm = 'cedar.hao:/d/dothe/cirris.cgm' sendps = 'cedar.hao:/d/dothe/cirris.ps' senddat = 'cedar.hao:/d/dothe/cirris.dat' &end_sut 'EOF' ja -st exit badx: echo 'Execution failed' rcp $job.cgm cedar.hao:/d/dothe/$job.fail.cgm exit