; &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' &end_sut