#!/bin/sh #======================================================================= # CMAM dynamics and chemistry diagnostics --- mamdiag_strat_t31 --- # $Id: mamdiag_strat_t31_jobdef 618 2010-02-11 01:57:41Z acrnrls $ #======================================================================= # keyword :: mamdiag_strat_t31 # description :: CMAM stratospheric dynamics and chemistry diagnostics # hide :: yes # # NOTE: this jobdef contains 2 decks # # Standard CMAM dynamics diagnostic string # # Based on mamchemdiag-1.2.2 example.job set -a . betapath2 # * ........................... Parmsub Parameters ....................... crawork=${runid}_mamdiag; username="acrnxxx"; user="XXX"; days="MAM" runid=aaa; uxxx=uxxx; mamdiag_flabel_uxxx=$uxxx; mamdiag_model_uxxx=$uxxx jobname=mamdiag; mamdiag_flabel_prefix=${mamdiag_flabel_uxxx}_${runid} nqsprfx="x${runid}_"; nqsext=''; RUNID=`echo "${runid}"|tr '[a-z]' '[A-Z]'`; run=CMAM year="001"; mon="06"; monn="jan"; flabel="${mamdiag_flabel_prefix}_${year}_m${mon}_"; mamdiag_model_prefix=${mamdiag_model_uxxx}_${runid} model1="${mamdiag_model_prefix}_${year}_m${mon}_"; days="JUN"; obsday="JUN"; plpfn="${runid}_${year}_m${mon}_"; # Libraries lopgm="lopgmhr"; oldiag="diag4"; # Moist moist=" QHYB"; itrvar=" Q"; # Memory and cpu settings dtime="200"; gptime="200"; stime="600"; memory1="100mb"; memory0="16mb"; memory2="100mb"; memory3="400mb"; # User identification username="ACRNY04"; user="JONSSON"; # Time coordinate details t1="000000000"; t2="999999999"; t3="0030"; t4=" 1"; delt=" 720.0"; # Momentum, energy and water stats tied to ISBEG g1=" 0"; g2="999999999"; g3="0090"; # Energy stats tied to ISRAD r1=" 0"; r2="999999999"; r3="0030"; # Horizontal coordinate details lrt=" 31"; lmt=" 31"; typ=" 2"; lon=" 64"; lat=" 32"; npg=" 2"; d="B"; # Vertical coordinate details lml=" 995"; lay=" 2"; coord=" ET15"; topsig="-1.00"; plid=" .0575"; r=" 1.5"; r_1=" 0.5"; etatop=" 5.675E-7"; etalvs=on; pmin="-4810"; pmax=" 1000"; # Set of pressure levels for diagnostic output plv=" 63"; plv2=" 52"; p01="-4700"; p02="-3100"; p03="-3150"; p04="-3200"; p05="-3300"; p06="-3400"; p07="-3500"; p08="-3700"; p09="-2100"; p10="-2150"; p11="-2200"; p12="-2300"; p13="-2400"; p14="-2500"; p15="-2700"; p16="-1100"; p17="-1150"; p18="-1200"; p19="-1300"; p20="-1400"; p21="-1500"; p22="-1700"; p23="-0100"; p24="-0150"; p25="-0200"; p26="-0300"; p27="-0400"; p28="-0500"; p29="-0700"; p30=" 10"; p31=" 15"; p32=" 20"; p33=" 25"; p34=" 30"; p35=" 35"; p36=" 40"; p37=" 50"; p38=" 60"; p39=" 70"; p40=" 80"; p41=" 90"; p42=" 100"; p43=" 115"; p44=" 130"; p45=" 150"; p46=" 170"; p47=" 200"; p48=" 250"; p49=" 285"; p50=" 300"; p51=" 350"; p52=" 400"; p53=" 450"; p54=" 500"; p55=" 600"; p56=" 650"; p57=" 700"; p58=" 750"; p59=" 780"; p60=" 800"; p61=" 850"; p62=" 925"; p63=" 1000"; # Diagnostics datatype=specsig; # data type expected for diagnostic input # (secsig/specpr/gridpr) join=1; # the number of files file1, file2,...to be processed # EP-flux diagnostic options for decks epf_total and epf_stat # The following setup computes total and stationary EP fluxes # and EPFD for all waves, as well as for the sum of zonal wavenumbers # 1 to 3 (i.e., planetary waves): #nbands="0"; # to process all zonal wavenumbers (k) nbands="1"; # to process a single wavenumber band, in addition to all k mfirst_band1=" 1"; # starting k of first band mlast_band1=" 3"; # ending k of first band plotlab_band1="M=1-3" # label for plots maxf=" 32"; # maxf=lon/2 (needed for llafrb & used only if nbands > 0) # gwd diagnostics (for gwd_diag_prs.dk) tt1=" 30"; # tt1 = israd (in model job) tt2="999999999"; st3=" 30"; # st3= = israd (in model job) s3=" 1" pmaxl=" 1000"; # Plot options plot=off ; lxp="00000";kax="00001";kin="00001"; plunit="UQAM" # Chemistry diagnostics psclev=" 50"; # diagnostic pressure level (nearest to 50 hPa) on # which PSCs are output # # * ............................ Condef Parameters ....................... # nextjob=on auto=on noprint=off tnoprint=on splitfiles=on # debug=on # gcmtsav=on; # model saves spectral temperature specsig=on; # related to datatype gcm2plus=on; # ensure diagnostics are compatible with GCM2 and # gcm2plus on all later versions gerhosave=on; # density deck: save number density on model levels gepsave=on; # pressure deck: save pressure on model levels wxstats=on; # save vert eddy stats ztsave=on; # save zonal-time mean data to xp file isave=on; # save instantaneous (un-averaged) data to ie/ip files # # * ........................... Deck Definition ............................. # Chemistry standard post-processing diagnostic. # # ---------- compute beta/del and process DATA DESCRIPTION record . del.dk . modinfo.dk # ---------- diagnose U, V, q, T, PHI and w . gpint.dk # ---------- other stats. MSLP, TS, CLDS, ENERGY...etc . mslpr_pres.dk . relhum.dk . gpstats.dk . mom_stat.dk . eng_stat.dk . wat_stat.dk . misc_stat.dk . cloud8.dk . sfc_stat.dk # ----------- extra CCmVal (08) saves . tndsav_rad.dk # averages and saves TTPS/TTPL . pressure.dk # computes P on the eta surface . o3c1.dk # save time aeraged ozone column . gwd_diag_prs.dk # extracts gravity wave quantities # ---------- computes density on p and eta surfaces # (this chemistry version uses reduced packing density) # (calculate _gep field required by thstats3) . density_dyn_diag.dk # # ---------- potential temperature . thstats3.dk # # ---------- horizontal eddy flux of temperature . hrz_tflux.dk # # ---------- estimate PSC prevalence . diag_chem_PSC.dk # # ----------- special MA diagnostics . geqtz2.dk # computes q,T and Z on eta surfaces . resd_circ2.dk # evaluate VR and WR . epf_total.dk # do EP-flux diagnostics (total) . epf_stat.dk # do EP-flux diagnostics (stationary waves) . Ertelpv1.dk # computes PV . gpww.dk # computes w(m/s) # ---------- clean-up pressure-level fields from dynamics . gppurg.dk # ---------- clean-up density output files . gxpurg_dyn_diag.dk # ---------- clean-up remains in $CCRNTMP from splitfiles=on #. cleanup.cdk # #end_of_job # # Standard CMAM chemistry diagnostic string # # Based on mamchemdiag-1.2.2 example.job # #!/bin/sh set -a . betapath2 # * ........................... Parmsub Parameters ....................... crawork=${runid}_mamdiag; username="acrnxxx"; user="XXX"; days="MAM" runid=aaa; uxxx=uxxx; mamdiag_flabel_uxxx=$uxxx; mamdiag_model_uxxx=$uxxx jobname=mamdiag; mamdiag_flabel_prefix=${mamdiag_flabel_uxxx}_${runid} nqsprfx="${runid}_"; nqsext="_${runid}"; RUNID=`echo "${runid}"|tr '[a-z]' '[A-Z]'`; run=CMAM year="001"; mon="06"; monn="jan"; flabel="${mamdiag_flabel_prefix}_${year}_m${mon}_"; mamdiag_model_prefix=${mamdiag_model_uxxx}_${runid} model1="${mamdiag_model_prefix}_${year}_m${mon}_"; days="JUN"; obsday="JUN"; plpfn="${runid}_${year}_m${mon}_"; # Libraries lopgm="lopgmhr"; oldiag="diag4"; # Moist moist=" QHYB"; itrvar=" Q"; # Memory and cpu settings dtime="200"; gptime="200"; stime="600"; memory1="100mb"; memory0="16mb"; memory2="100mb"; memory3="400mb"; # User identification username="ACRNY04"; user="JONSSON"; # Time coordinate details t1="000000000"; t2="999999999"; t3="0090"; t4=" 1"; delt=" 720.0"; # Horizontal coordinate details lrt=" 31"; lmt=" 31"; typ=" 2"; lon=" 64"; lat=" 32"; npg=" 2"; d="B"; # Vertical coordinate details lml=" 995"; lay=" 2"; coord=" ET15"; topsig="-1.00"; plid=" .0575"; r=" 1.5"; etatop=" 5.675E-7"; pmin="-4810"; pmax=" 1000"; # Set of pressure levels for diagnostic output plv=" 63"; plv2=" 52"; p01="-4700"; p02="-3100"; p03="-3150"; p04="-3200"; p05="-3300"; p06="-3400"; p07="-3500"; p08="-3700"; p09="-2100"; p10="-2150"; p11="-2200"; p12="-2300"; p13="-2400"; p14="-2500"; p15="-2700"; p16="-1100"; p17="-1150"; p18="-1200"; p19="-1300"; p20="-1400"; p21="-1500"; p22="-1700"; p23="-0100"; p24="-0150"; p25="-0200"; p26="-0300"; p27="-0400"; p28="-0500"; p29="-0700"; p30=" 10"; p31=" 15"; p32=" 20"; p33=" 25"; p34=" 30"; p35=" 35"; p36=" 40"; p37=" 50"; p38=" 60"; p39=" 70"; p40=" 80"; p41=" 90"; p42=" 100"; p43=" 115"; p44=" 130"; p45=" 150"; p46=" 170"; p47=" 200"; p48=" 250"; p49=" 285"; p50=" 300"; p51=" 350"; p52=" 400"; p53=" 450"; p54=" 500"; p55=" 600"; p56=" 650"; p57=" 700"; p58=" 750"; p59=" 780"; p60=" 800"; p61=" 850"; p62=" 925"; p63=" 1000"; # Diagnostics datatype=specsig; # data type expected for diagnostic input # (secsig/specpr/gridpr) join=1; # the number of files file1, file2,...to be processed # Chemistry diagnostics # pmaxc=" 420"; # model level at lower limit of chemistry domain pmint=" 474"; # top level of "tropospheric" domain (=first model # level below $pmaxc) psclev=" 50"; # diagnostic pressure level (nearest to 50 hPa) on # which PSCs are output daylev=" -100"; # diagnostic pressure level (near stratopause) on # which daymask is derived # # * ............................ Condef Parameters ....................... # nextjob=on auto=on noprint=off tnoprint=on splitfiles=on # debug=on # gcmtsav=on; # model saves spectral temperature specsig=on; # related to datatype gcm2plus=on; # ensure diagnostics are compatible with GCM2 and # gcm2plus on all later versions gerhosave=on; # density deck: save number density on model levels # gepsave=on; # density deck: save pressure on model levels # wxstats=on; # save vert eddy stats ztsave=on; # save zonal-time mean data to xp file isave=on; # save instantaneous (un-averaged) data to ie/ip files # Chemistry diagnostics # chem=on; # chemistry diagnostics turned on dnsave=on; # save day/night zonal-time mean data to [c]xp file trchem=off; # is chemistry calculated on all model levels (only # implemented in some major chemistry decks: Hydrogen, # Nitrogen, Oxygen, Hydrocarbon); see also parameters # plv2 and pmaxc # # * ........................... Deck Definition ............................. # Chemistry standard post-processing diagnostic. # # ---------- initiate xp, gp and ie files if they were not already created by the preceding dynamics job #. initstdiag.dk # ---------- computes density on p and eta surfaces # (this chemistry version uses reduced packing density) . density_chem_diag.dk # ---------- numberdensity for ND to VMR conversion . numberdensity.dk # ---------- daymask for day/night time zonal averaging . daymask.dk # ---------- chemistry . diag_chem_nd.dk . diag_chem_Oxygen.dk #. diag_chem_Ox.dk #. diag_chem_Ozone.dk #. diag_chem_O3.dk #. diag_chem_O3col.dk . diag_chem_Hydrocarbons.dk #. diag_chem_CH4.dk #. diag_chem_CH4ox.dk #. diag_chem_CO.dk . diag_chem_Hydrogen.dk #. diag_chem_HOx.dk #. diag_chem_H2.dk #. diag_chem_H2O.dk #. diag_chem_H2Otot.dk . diag_chem_Nitrogen.dk #. diag_chem_NOy.dk #. diag_chem_N2O.dk . diag_chem_Chlorine.dk #. diag_chem_Cly.dk #. diag_chem_CFC.dk . diag_chem_Bromine.dk #. diag_chem_Bry.dk #. diag_chem_CH3Br.dk # ---------- age of air . diag_chem_age.dk # ---------- chemical tendencies and burdens . diag_chem_chtend.dk # ---------- clean-up daymask, density, and gptqz output files . gxpurg_chem_diag.dk # ---------- clean-up model history and leftover pressure level files ###. mclean.dk # ---------- clean-up remains in $CCRNTMP from splitfiles=on . cleanup.cdk # #end_of_job