#!/bin/sh #======================================================================= # Create AGCM model diagnostic files (*gp and *xp) --- mdiag_t47 --- # $Id: mdiag_t47_jobdef 617 2010-02-11 01:53:48Z acrnrls $ #======================================================================= # keyword :: mdiag_t47 # description :: agcm model diagnostics (T47) # hide :: yes set -a . betapath2 # * ........................... Parmsub Parameters ............................ crawork=${runid}_mdiag; username="acrnxxx"; user="XXX"; days="MAM" runid=aaa; uxxx=uxxx; mdiag_flabel_uxxx=$uxxx; mdiag_model_uxxx=$uxxx jobname=mdiag; mdiag_flabel_prefix=${mdiag_flabel_uxxx}_${runid} oldiag="diag4"; nqsprfx="x${runid}_"; nqsext=''; RUNID=`echo "${runid}"|tr '[a-z]' '[A-Z]'`; run=AGCM # use_month_names will flag the use of jan,feb,mar,... rather than m01,m02,m03,... use_month_names=0 year="001"; mon="06"; monn="jan"; flabel="${mdiag_flabel_prefix}_${year}_m${mon}_"; mdiag_model_prefix=${mdiag_model_uxxx}_${runid} model1="${mdiag_model_prefix}_${year}_m${mon}_"; days="JUN"; obsday="JUN"; plpfn="${runid}_${year}_m${mon}_"; # atmos_file and year_offset are only required for ts_save.dk # atmos_file="${mdiag_flabel_prefix}_${year}_m${mon}_ts"; year_offset="0"; resol="96_48"; # # General surface observations: # obsfile="obs_sfc_${resol}_v2"; obsfile2="pd_nmc_${resol}_1979_1988_m${mon}_"; if [ "$mon" = "12" ] ; then obsfile2="pd_nmc_${resol}_1979_1987_m${mon}_"; fi # # GPCP, ISCCP, NVAP, and Dasilva observational data: # # obsgpcp=pooled_88_96_gpcp_${resol}; obsgpcp=pooled_79_97_xie_arkin_${resol}; obsisccp=pooled_84_90_isccpc2_${resol}; obspw=pooled_88_92_nvap_${resol}; obsbeg=pooled_${resol}_dasilva_obs; # # # mask=land_mask_${resol}; tlabel="tph"; tendlbl=" "; t1=" 1"; t2="999999999"; t3=" 1"; s3=" 1"; r1=" 1"; r2="999999999"; r3=" 1"; g1=" 1"; g2="999999999"; g3=" 1"; a1=" 1"; a2="999999999"; a3=" 1"; nbl=" 4"; lt1=" 1"; lt2=" 12"; lt3=" 21"; lt4=" 30"; lt5=" "; lb1=" 12"; lb2=" 21"; lb3=" 30"; lb4=" 33"; lb5=" "; lml=" 995"; d="B"; lay=" 2"; coord=" ET15"; topsig="-1.00"; moist=" QHYB"; plid=" 50.0"; sref=" 6.71E-3"; spow=" 1."; m01=" 012"; itrvar="QHYB"; it_LWC="01"; xref_LWC=" 0."; xpow_LWC=" 1."; it_IWC="02"; xref_IWC=" 0."; xpow_IWC=" 1."; it_BCO="03"; xref_BCO=" 0."; xpow_BCO=" 1."; it_BCY="04"; xref_BCY=" 0."; xpow_BCY=" 1."; it_OCO="05"; xref_OCO=" 0."; xpow_OCO=" 1."; it_OCY="06"; xref_OCY=" 0."; xpow_OCY=" 1."; it_SSA="07"; xref_SSA=" 9.58E-9"; xpow_SSA=" 1."; it_SSC="08"; xref_SSC=" 6.54E-8"; xpow_SSC=" 1."; it_DUA="09"; xref_DUA=" 1.22E-7"; xpow_DUA=" 1."; it_DUC="10"; xref_DUC=" 1.84E-7"; xpow_DUC=" 1."; it_DMS="11"; xref_DMS=" 1.32E-10"; xpow_DMS=" 1."; it_SO2="12"; xref_SO2=" 1.24E-9"; xpow_SO2=" 1."; it_SO4="13"; xref_SO4=" 1.18E-9"; xpow_SO4=" 1."; it_HPO="14"; xref_HPO=" 0."; xpow_HPO=" 1."; it_H2O="15"; xref_H2O=" 0."; xpow_H2O=" 1."; plv=" 17"; p01=" 10"; p02=" 20"; p03=" 30"; p04=" 50"; p05=" 70"; p06=" 100"; p07=" 150"; p08=" 200"; p09=" 250"; p10=" 300"; p11=" 400"; p12=" 500"; p13=" 600"; p14=" 700"; p15=" 850"; p16=" 925"; p17=" 1000"; p18=" "; p19=" "; p20=" "; p21=" "; p22=" "; p23=" "; p24=" "; p25=" "; p26=" "; p27=" "; p28=" "; p29=" "; p30=" "; p31=" "; p32=" "; p33=" "; p34=" "; p35=" "; p36=" "; p37=" "; p38=" "; p39=" "; p40=" "; p41=" "; p42=" "; p43=" "; p44=" "; p45=" "; p46=" "; p47=" "; p48=" "; p49=" "; p50=" "; kax=" 1"; kin=" 1"; lxp=" 0"; map=" 0"; b=" +"; dtime=100; gptime=600; stime=600; memory1="800mb"; memory2="800mb"; memory3="800mb"; plunit=VIC; lopgm=lopgm; # It is sometimes more efficient to set CCRNTMP = $RUNPATH CCRNTMP=$CCRNTMP lrt=" 47"; lmt=" 47"; typ=" 2"; lon=" 96"; lat=" 48"; npg=" 2"; ncx=" 2"; ncy=" 2"; delt=" 1200.0"; # MAM stuff # r=" 1.5"; etatop=" 5.00E-4"; # # pmaxl=" 1000"; pmin="-0100"; pmax=" 1000"; # * ............................ Condef Parameters ............................ join=1 auto=on nextjob=on noprint=on tnoprint=on datatype=specsig gcmtsav=on gcm2plus=on wxstats=on xtracld=off obsdat=off xtradif=on colourplots=on shade=off begplot=off amip2=off tprhs=off qprhs=off uprhs=off vprhs=off rhsplot=off debug=off lmstat=off carbon=off ctem=off slab=off tmpsave=on splitfiles=on # switch on the splitting mode. # * ............................. Deck Definition ............................. # --------------------------------- Compute beta/del and process DATA DESCRIPTION records. . del.dk . modinfo.dk # --------------------------------- Interpolate model data on pressure levels. # . gpqtz.dk # . gpuv.dk # . gpw.dk # . gpvort.dk # gpint.dk is a merge of the above 4 decks . gpint.dk # --------------------------------- Basic model statistics. . mslpr_eta.dk . gpstats.dk # --------------------------------- Spectral statistics (optional) # . spvdqtz.dk # . spmodinfo.dk # . sphum.dk # . spvort.dk # . spdiv.dk # . sptemp.dk # . spintr.dk # . vdelphi.dk # . sppurg.dk # --------------------------------- Tracer statistics (optional). # #par_def # trac=" SSA SSC DUA DUC DMS SO2 SO4" # . xstats.dk # . constat.dk # . dust.dk # --------------------------------- Model physics statistics. . mom_stat.dk . eng_stat3.dk . wat_stat.dk . misc_stat.dk . cloud8.dk . sfc_stat.dk . lmlstat.dk # lowest model level (optional) # --------------------------------- Time series deck (optional) # . ts_save.dk # set lmstat=on to accumulate lowest # # model level quantities # --------------------------------- Delete pressure-level files. . gppurg.dk . sppurg.dk # --------------------------------- cleanup . cleanup.dk #end_of_job