#!/bin/sh ################################################################################ #### UNIX Script Documentation Block # . . # Script name: exrtma_gsianl.sh.sms # Script description: Runs regional GSI variational analysis # # Author: Manuel Pondeca Org: NP22 Date: 2006-06-23 # # Script history log: # 2006-06-23 Manuel Pondeca # ################################################################################ set -x cd $DATA msg="JOB $job HAS BEGUN" postmsg "$jlogfile" "$msg" . $COMIN/${RUN}.t${cyc}z.envir.sh # Set resoltion and other dependent parameters export JCAP=${JCAP:-62} export LEVS=${LEVS:-60} export DELTIM=${DELTIM:-$((3600/($JCAP/20)))} cp $COMIN/${RUN}.t${cyc}z.gridname_input gridname_input # Copy in the prepbufr file echo "is the prepbufr file good?" echo `pwd` echo `ls -l $COMRTMA/rtma.t${cyc}z.prepbufr.tm00` cp $COMRTMA/rtma.t${cyc}z.prepbufr.tm00 prepbufr #this is the observation input file. Must know #where it's stored. Comes from Dennis Keyser's job echo $? echo `ls -l prepbufr` # Define the Analysis Date adate=$PDY$cyc # for plotting correlations, set an_test_output = .true. # variable names expected for var_plotcor are # st -- stream function # vp -- velocity potential # ps -- surface pressure # tv -- virtual temperature # q -- specific humidity # oz -- ozone # sst -- sea surface temperature # stl -- skin temp over land # sti -- skin temp over ice # cw -- cloud water # # coordinates of desired correlation test point are # # i_plotcor, j_plotcor, k_plotcor # # Make gsi namelist nhr_assimilation=12 cat << EOF > gsiparm.anl &SETUP miter=2,niter(1)=50,niter(2)=50 write_diag(1)=.true.,write_diag(2)=.true.,write_diag(3)=.true., gencode=78,qoption=1,tsensible=.true.,lgschmidt=.false., factqmin=0.0,factqmax=0.0,deltim=$DELTIM, ndat=5,iguess=-1, oneobtest=.false.,retrieval=.false., diag_rad=.false.,diag_pcp=.false.,diag_ozone=.false., nhr_assimilation=${nhr_assimilation},perturb_obs=.false.,perturb_fact=0.1 / &GRIDOPTS JCAP=$JCAP,NLAT=$NLAT,NLON=$LONA,nsig=$LEVS,hybrid=.true., wrf_nmm_regional=.false.,wrf_mass_regional=.false.,twodvar_regional=.true., diagnostic_reg=.false., filled_grid=.false.,half_grid=.true.,netcdf=.false., / &BKGERR as=0.35,0.35,0.35,0.50,0.50,0.80,1.00,1.00, hzscl=1.414,1.000,0.707, vs=0.5,bw=0.0, / &ANBKGERR anisotropic=.true.,an_vs=0.5,ngauss=1, an_flen_u=-5.,an_flen_t=3.,an_flen_z=-200., ifilt_ord=2,npass=3,normal=-200,grid_ratio=1.,nord_f2a=4, rtma_subdomain_option=.false.,triad4=.true.,nsmooth=0,nsmooth_shapiro=0,lreadnorm=.true. / &JCOPTS / &STRONGOPTS jcstrong=.false.,jcstrong_option=3,nstrong=1,nvmodes_keep=20,period_max=3., baldiag_full=.true.,baldiag_inc=.true., / &OBSQC dfact=0.75,dfact1=3.0,noiqc=.false.,oberrflg=.false., c_varqc=0.02,vadfile='prepbufr', / &OBS_INPUT dmesh(1)=60.0,dmesh(2)=60.0,dmesh(3)=60.0,dmesh(4)=60.0,time_window_max=6.0, dfile(01)='prepbufr', dtype(01)='ps', dplat(01)=' ', dsis(01)='ps', dval(01)=1.0, dthin(01)=0, dfile(02)='prepbufr' dtype(02)='t', dplat(02)=' ', dsis(02)='t', dval(02)=1.0, dthin(02)=0, dfile(03)='prepbufr', dtype(03)='q', dplat(03)=' ', dsis(03)='q', dval(03)=1.0, dthin(03)=0, dfile(04)='prepbufr', dtype(04)='uv', dplat(04)=' ', dsis(04)='uv', dval(04)=1.0, dthin(04)=0, dfile(05)='prepbufr', dtype(05)='spd', dplat(05)=' ', dsis(05)='spd', dval(05)=1.0, dthin(05)=0, / &SUPEROB_RADAR / &SINGLEOB_TEST maginnov=0.1,magoberr=0.1,oneob_type='v', oblat=36.,oblon=264.,obpres=1000.,obdattim=${adate}, obhourset=0., / &LAG_DATA / &HYBRID_ENSEMBLE / EOF cat << EOF > parmcard_input &parmcardanisof latdepend=.false., scalex1=1., scalex2=1., scalex3=1., afact0=1., hsteep=750., lsmoothterrain=.true., hsmooth_len=1.0, volpreserve=.false., lwater_scaleinfl=.false., water_scalefactpsi=1., water_scalefactchi=1., water_scalefacttemp=2., water_scalefactq=2., water_scalefactpsfc=1., nhscale_pass=1, rltop_wind=20000., rltop_temp=500., rltop_q=400., rltop_psfc=300., svpsi=0.35, svchi=0.72205, svpsfc=1.4, svtemp=2.5, svshum=1.5, sclpsi=0.1504, sclchi=0.1504, sclpsfc=0.20, scltemp=0.25, sclhum=0.25 / &parmcardreadprepb cgrid='guam' ladjusterr=.false., oberrinflfact=5.0, ineighbour=3, jneighbour=3 lshoreline=.true. slmland=0.9 / EOF # Set fixed files # berror = forecast model background error statistics # specoef = OPTRAN spectral coefficients # trncoef = OPTRAN transmittance coefficients # emiscoef = coefficients for IR sea surface emissivity model # satinfo = text file with information about assimilation of brightness temperatures # satangl = angle dependent bias correction file (fixed in time) # pcpinfo = text file with information about assimilation of prepcipitation rates # ozinfo = text file with information about assimilation of ozone data # errtable = text file with obs error for conventional data (regional only) # bufrtable= text file ONLY needed for single obs test (oneobstest=.true.) # bftab_sst= bufr table for sst ONLY needed for sst retrieval (retrieval=.true.) cp $FIXrtma/gurtma_regional_nmm_berror.f77 berror_stats cp $FIXrtma/gurtma_nam_errtable.r3dv errtable cp $FIXrtma/gurtma_regional_convinfo.txt convinfo cp $FIXrtma/rtma_prepobs_prep.bufrtable prepobs_prep.bufrtable #note: available in /nwprod/fix cp $COMIN/${RUN}.t${cyc}z.t_rejectlist t_rejectlist cp $COMIN/${RUN}.t${cyc}z.p_rejectlist p_rejectlist cp $COMIN/${RUN}.t${cyc}z.q_rejectlist q_rejectlist cp $COMIN/${RUN}.t${cyc}z.w_rejectlist w_rejectlist cp $FIXrtma/gurtma_slmask.dat rtma_slmask.dat cp $FIXrtma/gurtma_terrain.dat rtma_terrain.dat # Copy bias correction, sigma, and surface files # # *** NOTE: The regional gsi analysis is written to (over) # the input guess field file (wrf_inout) # cp $COMIN/${RUN}.t${cyc}z.2dvar_input ./wrf_inout if [[ "$usefgat" = ".true." ]] ; then cp $COMIN/${RUN}.t${cyc}z.2dvar_input_bi2 ./wrf_inou2 cp $COMIN/${RUN}.t${cyc}z.2dvar_input_bi3 ./wrf_inou3 cp $COMIN/${RUN}.t${cyc}z.2dvar_input_bi4 ./wrf_inou4 cp $COMIN/${RUN}.t${cyc}z.2dvar_input_bi5 ./wrf_inou5 cp $COMIN/${RUN}.t${cyc}z.2dvar_input_bi6 ./wrf_inou6 cp $COMIN/${RUN}.t${cyc}z.2dvar_input_bi7 ./wrf_inou7 cp $COMIN/${RUN}.t${cyc}z.2dvar_input_bi8 ./wrf_inou8 fi cp $FIXrtma/gurtma_fltnorm.dat_psi $DATA/fltnorm.dat_psi cp $FIXrtma/gurtma_fltnorm.dat_chi $DATA/fltnorm.dat_chi cp $FIXrtma/gurtma_fltnorm.dat_ps $DATA/fltnorm.dat_ps cp $FIXrtma/gurtma_fltnorm.dat_t $DATA/fltnorm.dat_t cp $FIXrtma/gurtma_fltnorm.dat_pseudorh $DATA/fltnorm.dat_pseudorh cp $FIXrtma/gurtma_fltnorm.dat_ist $DATA/fltnorm.dat_ist cp $FIXrtma/gurtma_fltnorm.dat_lst $DATA/fltnorm.dat_lst cp $FIXrtma/gurtma_fltnorm.dat_oz $DATA/fltnorm.dat_oz cp $FIXrtma/gurtma_fltnorm.dat_qw $DATA/fltnorm.dat_qw cp $FIXrtma/gurtma_fltnorm.dat_sst $DATA/fltnorm.dat_sst cp /nwprod/fix/rtma_random_flips $DATA/random_flips export pgm=gurtma_gsianl . prep_step poe $EXECrtma/gurtma_gsianl >$pmgout 2>errfile export err=$?; err_chk ls -lrt gradx.dat_01_* | grep -c gradx.dat > used_iterations.dat cat shoreline_obrelocation.dat* > all_shoreline_relocatedobs.dat if [ $SENDCOM = YES ] then typeset -2Z nhr_assimilation cp wrf_inout $COMOUT/${RUN}.t${cyc}z.wrfanl cp siganl $COMOUT/${RUN}.t${cyc}z.siganl cp sigf12 $COMOUT/${RUN}.t${cyc}z.sigf12 cp sigf13 $COMOUT/${RUN}.t${cyc}z.sigf13 cp sigf14 $COMOUT/${RUN}.t${cyc}z.sigf14 cp sigf15 $COMOUT/${RUN}.t${cyc}z.sigf15 cp sigf16 $COMOUT/${RUN}.t${cyc}z.sigf16 cp sigf17 $COMOUT/${RUN}.t${cyc}z.sigf17 cp sigf18 $COMOUT/${RUN}.t${cyc}z.sigf18 cp sigf19 $COMOUT/${RUN}.t${cyc}z.sigf19 cp sigf${nhr_assimilation} $COMOUT/${RUN}.t${cyc}z.sigges cp bckg_dxdy.dat $COMOUT/${RUN}.t${cyc}z.bckg_dxdy.dat cp bckg_qsat.dat $COMOUT/${RUN}.t${cyc}z.bckg_qsat.dat cp bckg_psfc.dat $COMOUT/${RUN}.t${cyc}z.bckg_psfc.dat cp bckgvar.dat_psi $COMOUT/${RUN}.t${cyc}z.bckgvar_psi.dat cp bckgvar.dat_chi $COMOUT/${RUN}.t${cyc}z.bckgvar_chi.dat cp bckgvar.dat_ps $COMOUT/${RUN}.t${cyc}z.bckgvar_ps.dat cp bckgvar.dat_t $COMOUT/${RUN}.t${cyc}z.bckgvar_t0.dat cp bckgvar.dat_pseudorh $COMOUT/${RUN}.t${cyc}z.bckgvar_pseudorh.dat cp bckg_z.dat $COMOUT/${RUN}.t${cyc}z.bckg_z.dat cp used_iterations.dat $COMOUT/${RUN}.t${cyc}z.used_iterations.dat cp all_shoreline_relocatedobs.dat $COMOUT/${RUN}.t${cyc}z.all_shoreline_relocatedobs.dat nfiles=`cat used_iterations.dat` ic=0 while [ $ic -le $(($nfiles-1)) ] ; do typeset -3Z ic cp gradx.dat_01_${ic} $COMOUT/${RUN}.t${cyc}z.gradx.dat_01_${ic} cp grady.dat_01_${ic} $COMOUT/${RUN}.t${cyc}z.grady.dat_01_${ic} let "ic=ic+1" done cp prepbufr $COMOUT/${RUN}.t${cyc}z.prepbufr cp fort.201 $COMOUT/${RUN}.t${cyc}z.pressfcfits cp fort.202 $COMOUT/${RUN}.t${cyc}z.windfits cp fort.203 $COMOUT/${RUN}.t${cyc}z.tempfits cp fort.204 $COMOUT/${RUN}.t${cyc}z.shumidfits cp fort.220 $COMOUT/${RUN}.t${cyc}z.penaltyinfo #cp p_rejectlist $COMOUT/${RUN}.t${cyc}z.p_rejectlist #cp q_rejectlist $COMOUT/${RUN}.t${cyc}z.q_rejectlist #cp t_rejectlist $COMOUT/${RUN}.t${cyc}z.t_rejectlist #cp w_rejectlist $COMOUT/${RUN}.t${cyc}z.w_rejectlist cp gsiparm.anl $COMOUT/${RUN}.t${cyc}z.gsiparm.anl cp parmcard_input $COMOUT/${RUN}.t${cyc}z.parmcard_input fi # Loop over first and last outer loops to generate innovation # diagnostic files for indicated observation types (groups) # # NOTE: Since we set miter=2 in GSI namelist SETUP, outer # loop 03 will contain innovations with respect to # the analysis. Creation of o-a innovation files # is triggered by write_diag(3)=.true. The setting # write_diag(1)=.true. turns on creation of o-g # innovation files. # loops="01 02 03" for loop in $loops do # Collect diagnostic files for obs types (groups) below list="conv" for type in $list do count=`ls ${DATA}/dir.*/${type}_${loop}* | wc -l` if [[ count -gt 0 ]]; then cat ${DATA}/dir.*/${type}_${loop}* > ${DATA}/diag_${type}_${loop} compress ${DATA}/diag_${type}_${loop} if [ $SENDCOM = YES ] then if [ "$loop" = "01" ] then cp ${DATA}/diag_${type}_${loop}.Z $COMOUT/${RUN}.t${cyc}z.diag_${type}_ges.Z elif [ "$loop" = "03" ] then cp ${DATA}/diag_${type}_${loop}.Z $COMOUT/${RUN}.t${cyc}z.diag_${type}_anl.Z else cp ${DATA}/diag_${type}_${loop}.Z $COMOUT/${RUN}.t${cyc}z.diag_${type}_${loop}.Z fi fi fi done done ######################################################## msg='$job ENDED NORMALLY.' postmsg "$jlogfile" "$msg" ################## END OF SCRIPT #######################