# # QSUB -s /bin/csh # QSUB -lT 2500 # QSUB -q prem # QSUB -nr # QSUB -lM 4MW # QSUB -o flxproc.flxcm3.$$ # QSUB -eo # # 8/98: this must be run on paiute, as that is the only machine # that still has the ccm processor. # cd $TMPDIR # # RUNCCM2PROC = 1 -> execute ccm2 processor with below inputs to make # a "List Sorted Data" (lsd) save tape for use by flxproc below. # (not necessary to run ccm2proc *if* ccmlsd (flxproc input) mss path # pre-exists) # (comment out or set to 0 to not run ccm2proc) # set RUNCCM2PROC = 0 set job = flxproc #batchname $job.$$ # # Run ccm2 processor and save lsd file: # if (! $?RUNCCM2PROC || $RUNCCM2PROC == 0) goto runflxproc cat >! parms << 'END' C C Mendillo Case B: 9/27-10/7 (days 270-280) C /ROBLE/csm/flxcm3/ccm3/hist/h0138 (ha0277) - h0148 (ha0297) C file=h0138, date=9/27, hist file day=636, ccmproc day=342 C (636-342=294 diff between model and proc)) C CTAPESA = '/ROBLE/csm/flxcm3/ccm3/hist/h0138' CDAYSA = 342.0 CTITLEA = 'Mendillo Case B day1 9/27' C CTAPESA = '/ROBLE/csm/flxcm3/ccm3/hist/h0143' CDAYSA = 347.0 CTITLEA = 'Mendillo Case B day5 10/2' C CTAPESA = '/ROBLE/csm/flxcm3/ccm3/hist/h0148' CDAYSA = 352.0 CTITLEA = 'Mendillo Case B day10 10/7' C CTAPESA = '/ROBLE/csm/flxcm3/ccm3/hist/ha0276', C '/ROBLE/csm/flxcm3/ccm3/hist/ha0277', C '/ROBLE/csm/flxcm3/ccm3/hist/ha0278', C '/ROBLE/csm/flxcm3/ccm3/hist/ha0279', C '/ROBLE/csm/flxcm3/ccm3/hist/ha0280', C '/ROBLE/csm/flxcm3/ccm3/hist/ha0281', C '/ROBLE/csm/flxcm3/ccm3/hist/ha0282', C '/ROBLE/csm/flxcm3/ccm3/hist/ha0283', C '/ROBLE/csm/flxcm3/ccm3/hist/ha0284', C '/ROBLE/csm/flxcm3/ccm3/hist/ha0285', C '/ROBLE/csm/flxcm3/ccm3/hist/ha0286' CDAYSA = 342.0,347.0, .041666 CTITLEA = 'Mendillo Case B 9/27-10/2' C TAPESA = '/ROBLE/csm/flxcm3/ccm3/hist/ha0286', '/ROBLE/csm/flxcm3/ccm3/hist/ha0287', '/ROBLE/csm/flxcm3/ccm3/hist/ha0288', '/ROBLE/csm/flxcm3/ccm3/hist/ha0289', '/ROBLE/csm/flxcm3/ccm3/hist/ha0290', '/ROBLE/csm/flxcm3/ccm3/hist/ha0291', '/ROBLE/csm/flxcm3/ccm3/hist/ha0292', '/ROBLE/csm/flxcm3/ccm3/hist/ha0293', '/ROBLE/csm/flxcm3/ccm3/hist/ha0294', '/ROBLE/csm/flxcm3/ccm3/hist/ha0295', '/ROBLE/csm/flxcm3/ccm3/hist/ha0296', '/ROBLE/csm/flxcm3/ccm3/hist/ha0297' DAYSA = 347.0,352.0, .041666 TITLEA = 'Mendillo Case B 10/2-10/7' C C Get tgcm/t21 coupled fields from ccm2: C FIELDA1 = 'T','U','V','Q','Z2','OMEGA' C C 10 values of PRESSLE for T21: C (mb pressures corresponding to zp -20.5 to -16.0, every 0.5) C PRESSLE = 4.3292312622E+02, 2.6258111572E+02, 1.5926348877E+02, 9.6598197937E+01, 5.8589767456E+01, 3.5536491394E+01, 2.1553972244E+01, 1.3073143959E+01, 7.9292626381E+00, 4.8093414307E+00 C C Run ccm2 proc in LSD mode: C DRVRTYP = 'LSD1' C C RTLS1A = 2 -> ordering of horiz slices = field level, time, ensemble member C RTLS1A = 2 C C PKLS1A = 1 -> no packing C PKLS1A = 1 C C SAVLS1A = mss path at which to store lsd save tape C (should be same as ccmlsd in processor inputs below) C CSAVLS1A = '/ROBLE/csm/flxcm3/ccm3/lsd/h0138' CSAVLS1A = '/ROBLE/csm/flxcm3/ccm3/lsd/h0143' CSAVLS1A = '/ROBLE/csm/flxcm3/ccm3/lsd/h0148' CSAVLS1A = '/ROBLE/csm/flxcm3/ccm3/lsd/day270-275' SAVLS1A = '/ROBLE/csm/flxcm3/ccm3/lsd/day275-280' ENDOFDATA 'END' /ccm/proc/Processor || goto ccmproc_fail echo ' ' echo '-------------------------------------------------------------------' echo ' ' echo Completed ccm2 processor run at `date` runflxproc: ja # # Get flxproc executable from mss: # msread $job /FOSTER/bin/$job chmod u+x $job # # Run flx processor: # cat >! $job.inp << 'EOF' &input_sut E ; ;---------------------------------------------------------------------- ; ; tgcm/t21: ; ; histvols = time-gcm history (-16 to +5 coupled run) ; ; Mendillo Case B: 9/27-10/7 (days 270-280) ; /ROBLE/flxcm3/CCM3S49 (270,0,0) - CCM3S54 (280,0,0) prim ; /ROBLE/flxcm3/CCMS502 (270,0,0) - CCMS522 (280,0,0) sec ; ;histvols = '/ROBLE/flxcm3/CCM3S49', mtimes= 270,0,0 ;histvols = '/ROBLE/flxcm3/CCM3S52', mtimes= 275,0,0 ;histvols = '/ROBLE/flxcm3/CCM3S54', mtimes= 280,0,0 ; histvols = '/ROBLE/flxcm3/CCMS502','/ROBLE/flxcm3/CCMS503', '/ROBLE/flxcm3/CCMS504','/ROBLE/flxcm3/CCMS505', '/ROBLE/flxcm3/CCMS506','/ROBLE/flxcm3/CCMS507', '/ROBLE/flxcm3/CCMS508','/ROBLE/flxcm3/CCMS509', '/ROBLE/flxcm3/CCMS510','/ROBLE/flxcm3/CCMS511' mtimes = 270,0,0,'to',275,0,0,'by',60 ; ;histvols = '/ROBLE/flxcm3/CCMS511','/ROBLE/flxcm3/CCMS512', ; '/ROBLE/flxcm3/CCMS513','/ROBLE/flxcm3/CCMS514', ; '/ROBLE/flxcm3/CCMS515','/ROBLE/flxcm3/CCMS516', ; '/ROBLE/flxcm3/CCMS517','/ROBLE/flxcm3/CCMS518', ; '/ROBLE/flxcm3/CCMS519','/ROBLE/flxcm3/CCMS520', ; '/ROBLE/flxcm3/CCMS521','/ROBLE/flxcm3/CCMS522' ;mtimes = 275,0,0,'to',280,0,0,'by',60 ; ; ccmlsd = mss path of ccm2 lsd save tape written by ccm2 processor ; (should be same as SAVLS1A above if ccm2 processor was executed) ; ; ccmlsd = '/ROBLE/csm/flxcm3/ccm3/lsd/h0138' ; ccmlsd = '/ROBLE/csm/flxcm3/ccm3/lsd/h0143' ; ccmlsd = '/ROBLE/csm/flxcm3/ccm3/lsd/h0148' ccmlsd = '/ROBLE/csm/flxcm3/ccm3/lsd/day270-275' ; ccmlsd = '/ROBLE/csm/flxcm3/ccm3/lsd/day275-280' ; ;---------------------------------------------------------------------- ; ; Request fields: cfields is string array with the below valid values ; (spelling is important) ; Coupled fields are bracketed, e.g., [TN] (CH4 coupled only with maccm) ; Quotes around cfield field names are not necessary. ; Do not include brackets or parens in cfields. ; ; [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 UNVN UIVI FOF2 HMF2 ; cfields = UN ;cfields = FOF2,HMF2 ;cfields = TN,UN,VN,H2O,W,Z ;cfields = TN,UN,VN,O1,N2,Z,NE,FOF2,HMF2 ; ;fmnmxint = 'TN',1.,0.,5. fmnmxint = 'UN',1.,0.,10. ; ; iden = 0 -> mass mix ratios (for coupled H2O) ; iden = 0 ionvel = 0 modelhts = 1 ; ;---------------------------------------------------------------------- ; 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') ;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 = 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 ; ; iyaxright: single integer flag for extra right-hand y-axes. ; (applies to ipltlat and ipltlon vertical slices only) ; For now, assume left hand y-axis is always linear Zp, then: ; iyaxright = 0 -> add no extra right-hand y-axes ; iyaxright = 1 -> add extra (non-linear) right-hand y-axis in height (km) ; iyaxright = 2 -> add extra right-hand y-axes in height *and* pressure (mb) ; iyaxright = 2 ; ;---------------------------------------------------------------------- ; Maps: ; (contour over map projections if ipltmaps=1) ; (if ipltmaps=0, remaining map parameters are ignored) ; (as of 11/3/95, interpolation to height surface not yet available -- ; if fmap_zpht > 5., it will be reset to 5.) ; ipltmaps = 0 ; fmap_zpht = -20.,-18.,-16.,-14.,-12.,-7. ; 1 int: make map projection fmap_zpht = -16. map_continents = 0 ; 1 int: continent 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 = 0 ; 1 int: add uivi to epot ivec_label = 0 ; 1 int: label vec mags ; map_top_anno = 'Day 1 of Mendillo 10-day Case B' ; map_top_anno = 'Day 5 of Mendillo 10-day Case B' map_top_anno = 'Mendillo Case B days 270-275' ; ; 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 = -30., 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_eradii = 1.5 ; 1 float: earth radii dist ; ;---------------------------------------------------------------------- ; ; Latitude slices (contour longitude vs zp or ht at selected latitudes) ; (more than one pair of flat_zprange is allowed for vertical scales) ; (if flat_zprange(1) is < bottom of model(s), or flat_zprange(2) is above ; model(s), it will be adjusted to match bottom and top of model(s), ; e.g., flat_zprange = -99.,99. will be adjusted to full vertical range) ; flat_htscale = linear height scale(s) for interpolation on y-axis ; (triplets of bot,top,delta heights for y-axis) ; (more than one triplet for flat_htscale is allowed for height scales) ; ipltlat = 0 ; lats 51.5 for Slough, -51.7 for P. Stanley, -27.5 for Brisbane flats = 51.5, -51.7, -27.5 flat_zprange = -99.,5., -99.,-8. ; flat_htscale = 5.,500.,10., 5.,100.,1., 5.,50.,.5 ilat_log10 = 1 ; 1 int: plot log10 of species ; lat_top_anno = 'Day 1 of Mendillo 10-day Case B' ; lat_top_anno = 'Day 5 of Mendillo 10-day Case B' lat_top_anno = 'Day 10 of Mendillo 10-day Case B' ; ;---------------------------------------------------------------------- ; ; Longitude slices (contour latitude vs zp or ht at selected longitudes) ; (more than one pair of flon_zprange is allowed for vertical scales) ; (if flon_zprange(1) is < bottom of model(s), or flon_zprange(2) is above ; model(s), it will be adjusted to match bottom and top of model(s), ; e.g., flon_zprange = -99.,99. will be adjusted to full vertical range) ; flon_htscale = linear height scale(s) for interpolation on y-axis ; (triplets of bot,top,delta heights for y-axis) ; (more than one triplet for flon_htscale is allowed for height scales) ; ; (flons = selected longitudes. Use flons = 'zm' for zonal means, ; e.g. flons = -70.,'zm' will make lon slices at 70W, and zonal means) ; 1/flxproc.day168 ipltlon = 0 ; lons 0. for Slough, -58. for Stanley, 153. for Brisbane ; flons = 0., -58., 153. ; flon_zprange = -99.,5., -99.,-8. flons = -58. flon_zprange = -99., -8. ; flon_htscale = 5.,500.,10., 5.,100.,1., 5.,50.,.5 ilon_log10 = 1 ; 1 int: plot log10 of species ; lon_top_anno = 'Day 1 of Mendillo 10-day Case B' ; lon_top_anno = 'Day 5 of Mendillo 10-day Case B' ; lon_top_anno = 'Day 10 of Mendillo 10-day Case B' lon_top_anno = 'Mendillo Case B days 270-275' ; lon_top_anno = 'Mendillo Case B days 275-280' ; ;---------------------------------------------------------------------- ; ; 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_locname = 'Port Stanley' utvert_locs = -51.7,-58. ;utvert_locname ='Slough', 'Moscow', 'Wallops Island','Hobart', ; 'Brisbane', 'Port Stanley', 'Wakkanai' ;utvert_locs = 51.5,359.4, 55.5,37.3, 37.8,284.5, -42.9,147.3, ; -27.5,152.9, -51.7,-58., 45.,142. ;utvert_zprange = -99.,99.,-99.,-8. ; bot,top for y-axis zp utvert_zprange = -99.,-8. ;utvert_htscale = 30.,100.,1. ; bot,top,delta for y-axis ht iutvert_log10 = 1 ; log10 of species ; ;---------------------------------------------------------------------- ; Send Files: ; ;sendcgm = 'vishnu.hao:/d/foster/usr/mendillo/flxcm3.sut.d270.cgm' ;sendcgm = 'vishnu.hao:/d/foster/usr/mendillo/flxcm3.sut.d275.cgm' ;sendcgm = 'vishnu.hao:/d/foster/usr/mendillo/flxcm3.d270-275.cgm' ;sendcgm = 'vishnu.hao:/d/foster/usr/mendillo/flxcm3.d275-280.cgm' ;sendcgm = 'vishnu.hao:/d/foster/usr/mendillo/tn.lon.d270-275.cgm' ;sendcgm = 'vishnu.hao:/d/foster/usr/mendillo/tn.lon.d275-280.cgm' ;sendcgm = 'vishnu.hao:/d/foster/usr/mendillo/un.lon.d270-275.cgm' ;sendcgm = 'vishnu.hao:/d/foster/usr/mendillo/un.lon.d275-280.cgm' ;sendcgm = 'vishnu.hao:/d/foster/usr/mendillo/vn.lon.d270-275.cgm' ;sendcgm = 'vishnu.hao:/d/foster/usr/mendillo/vn.lon.d275-280.cgm' ;sendcgm = 'vishnu.hao:/d/foster/usr/mendillo/f2.lon.d270-275.cgm' sendcgm = 'vishnu.hao:/d/foster/usr/mendillo/un.utzp.d270-275.cgm' ; ;---------------------------------------------------------------------- &end_sut 'EOF' echo Executing $job.. $job < $job.inp || goto xfail ja -st exit # ccmproc_fail: echo \>\>\> ccm2 processor execution failed exit # mkproc_fail: echo \>\>\> make of flx processor failed exit # xfail: echo \>\>\> Execution of flxproc failed exit