----------------------------------------------------------- # Part 1: run_cedsyn (BEGIN) BASIC FRAMEWORK for READING/MASSAGING/GRAPHING SINGLE RADAR ----------------------------------------------------------- # # On the command line type "run_cedsyn" followed by HHMM of the analysis-time # For example, "run_cedsyn 1500" will execute with $1 replaced with 1500. # # Link input files to fortran units 71 & 72 which are read (READVOL) by Cedric # New DOW files are written (WRITVOL to 81 & 82) after some massaging of # reflectivity and velocity. The original DOW fields (DD and VE) are # retained by renaming them (RENAME) as DZreo and VEreo and creating # new fields named as DZ and VE. Final results of synthesis are written # (WRITVOL) to fortran unit 90. # ln -sf ./$1-DOW2.ced fort.71 ln -sf ./$1-DOW3.ced fort.72 # # Execute cedric_linux.e with # input (<<) from cedsyn.inp and output (>&) into cedsyn.log # Files are all moved into a 06Aug.$1 directory after execution # /h/eol/ljmill/bin/cedric_linux.e << cedsyn.inp >& cedsyn.log COMMENT ------------------------------------------------------- READ REORDER OUTPUT FILES, MASSAGE FIELDS, GRAPH, and WRITE NEW CEDRIC OUTPUT FILES FOR INPUT TO SYNTHESIS ------------------------------------------------------- COMEND *-------*-------*-------*-------*-------*-------*-------*-------*-------*------- VALUE SET RADAR to DOW2 VALUE SET FortInp to 71.0 VALUE SET FortOut to 81.0 EXPAND FIL_GRAF *-------*-------*-------*-------*-------*-------*-------*-------*-------*------- VALUE SET RADAR to DOW3 VALUE SET FortInp to 72.0 VALUE SET FortOut to 82.0 EXPAND FIL_GRAF *-------*-------*-------*-------*-------*-------*-------*-------*-------*------- ------------------------------------------------------------- # Part 2: run_cedsyn (FIL_GRAF) DEFINE BLOCK used for READING/MASSAGING/GRAPHING SINGLE RADAR ------------------------------------------------------------- DEFINE FIL_GRAF VALUE ON READVOL FortInp NEXT 000000. 240000. YES STATS PRINT Z 1.0 ALL BOTH FULL * Delete all unnecessary fields DELETE DCZ ZZ VVUF SW AZ EL CNT OCT DM * Fix bad field name for DOW3 RENAME DDDBZFIEas DD * * Rename the reflectivity (DD) and radial velocity (VE) fields * to retain those from REORDER without any massaging. New * reflectivity (DZ) and radial velocity (VE) fields will be * created with the PATCHER commands. * *ELETE DZ *ENAME DD as DZreo RENAME DZ as DZreo RENAME VE as VEreo * * Create new reflectivity (DZ) and radial velocity (VE) fields. * PATCHER DZ DZreo FILLCON 9.0 4.0 20.0 Z PATCHER DZ DZ FILLCON 5.0 4.0 20.0 Z PATCHER DZ DZ FILLCON 2.0 4.0 20.0 Z PATCHER VE VEreo DECILOC 3.0 4.0 5.0 Z PATCHER VE VE FILLCON 9.0 4.0 20.0 Z PATCHER VE VE FILLCON 5.0 4.0 20.0 Z PATCHER VE VE FILLCON 2.0 4.0 20.0 Z STATS PRINT Z 1.0 ALL BOTH FULL *INDOW DISTANCE-40.0 40.0 -30.0 50.0 WINDOW DISTANCE-35.0 35.0 -20.0 50.0 0.0 4.0 GRAPHICS01*01*01Z 1.0 REORDER Dorade sweep files WINDOW CONTOUR DZreo FALL 00011043 UNI -30.0 70.0 5.0 CONTOUR DZ FALL 00011043 UNI -30.0 70.0 5.0 CONTOUR VEreo FALL 00011043 UNI -25.0 25.0 2.0 CONTOUR VE FALL 00011043 UNI -25.0 25.0 2.0 END WRITVOL FortOut RADAR BEG VALUE OFF DEFEND FIL_GRAF ---------------------------------------------------- # Part 3: run_cedsyn (SYNUVW) DEFINE BLOCK used for SYNTHESIS, CREATING ADDITIONAL FIELDS, GRAPHING, AND WRITING THESE FINAL RESULTS ---------------------------------------------------- EXPAND SYNUVW DEFINE SYNUVW VALUE ON COMMENT -------------------------------------------------------------- SYNTHESIS OF RADIAL VELOCITIES WITH DIFFERENTIAL ADVECTION 1) NUMBER OF RADARS = 2 2) TWO EQUATION SOLUTION U3 = U2 + EWU*(W+Vt) V3 = V2 + EWV*(W+Vt) 3) (U,V,W) STD = (5.0,5.0,5.0) 4) DIFFERENTIAL ADVECTION WITH STORM DIR/SPD = (127,12) 5) TRANSFER REFL AND RADIAL VELOCITY DTEST1 - Both |EWU| and |EWV| .LT. DTEST1 Little or no vertical motion contribution DTEST2 - Both USTD and VSTD .LT. DTEST2 DTEST3 - WSTD .LT. DTEST3 Note: Large DTEST1-3 ==> UV winds close to the baseline -------------------------------------------------------------- P_______P_______P_______P_______P_______P_______P_______P_______P_______P_______ 1_______2_______3_______4_______5_______6_______7_______8_______9_______10______ COMEND *YNTHES NSNVOL NRADS DTEST1 DTEST2 DTEST3 ITWOEQ ITANAL STMDIR STMSPD SYNTHES SYN1250 2.0 2.0 2.0 2.0 YES 81.0 NEXT YES VE YES INPUT DZ VE OUTPUT DZ_2 VE_2 82.0 NEXT YES VE YES INPUT DZ VE OUTPUT DZ_3 VE_3 END * * HSTD - proportional to the angle between dual-Doppler radar beams * FUNCTIONNAME 1.0 FULL Usyn P LINEAR U 1.0 0.0 Vsyn P LINEAR V 1.0 0.0 HSTD P SQ+SQ USTD VSTD HSTD P SQRT HSTD END STATS PRINT Z 1.0 ALL BOTH FULL ----------------------------------------- # Part 4: run_cedsyn (PATCHER) COMMANDS for ADDITIONAL FILLING OF FIELDS ----------------------------------------- COMMENT ---------------------------------------------------------------------- PATCHER NAMOUF NAMINF DECIGLO SIGFAC IFIXAX IWINDO PATCHER NAMOUF NAMINF DECILOC MAXSTP MINPTS MAXDEV IFIXAX IWINDO PATCHER NAMOUF NAMINF FILLCON3MAXSTP MINQAD MINPTS DECFLD PATCHER NAMOUF NAMINF FILLCON MAXSTP MINQAD MINPTS DECFLD Z FULL FILTER NAMOUF NAMINF METHOD NSTEP IFIXAX IWINDO ---------------------------------------------------------------------- COMEND PATCHER U U FILLCON 15.0 4.0 20.0 Z PATCHER U U FILLCON 13.0 4.0 20.0 Z PATCHER U U FILLCON 11.0 4.0 20.0 Z PATCHER U U FILLCON 9.0 4.0 20.0 Z PATCHER U U FILLCON 5.0 4.0 20.0 Z PATCHER U U FILLCON 2.0 4.0 20.0 Z PATCHER V V FILLCON 15.0 4.0 20.0 Z PATCHER V V FILLCON 13.0 4.0 20.0 Z PATCHER V V FILLCON 11.0 4.0 20.0 Z PATCHER V V FILLCON 9.0 4.0 20.0 Z PATCHER V V FILLCON 5.0 4.0 20.0 Z PATCHER V V FILLCON 2.0 4.0 20.0 Z ---------------------------------------------- # Part 5: run_cedsyn (FILTER) COMMANDS for FILTERING and CREATING NEW FIELDS ---------------------------------------------- FILTER U U TWO FILTER V V TWO *ONVERGECONV U V 3.0 CONVERGECONV U V 5.0 FUNCTIONNAME 1.0 FULL UVSQ T SQ+SQ U V SPEED P SQRT UVSQ DZ_max P MAX DZ_2 DZ_3 DZ_max P ORELSE DZ_max DZ_2 DZ_max P ORELSE DZ_max DZ_3 VORT P DDI-DDJ V U 1.0 1.0 0.0 0.0 END COMMENT ------------------------------------------------------------------------ * Filter horizontal convergence and vertical vorticity to eliminate * two-delta noise that results from finite differencing. HAN (Hanning) * followed by INV (inverse Hanning) or TWO (Hanning-inverse) does this. * If more filtering is needed or desired, use the LEA (least-squares) * method going out two (2.0) or more grid points. ------------------------------------------------------------------------ COMEND FILTER CONV CONV HAN FILTER CONV CONV INV *ILTER CONV CONV TWO *ILTER CONV CONV LEA 2.0 FILTER VORT VORT HAN FILTER VORT VORT INV *ILTER VORT VORT TWO *ILTER VORT VORT LEA 2.0 ---------------------------------------------- # Part 6: run_cedsyn (COMPARE) COMMANDS for COMPARING INPUT and OUTPUT FIELDS ---------------------------------------------- VALUE ON VALUE SET X2 to 0.0 VALUE SET Y2 to 0.0 VALUE SET Z2 to 0.152 VALUE SET X3 to -0.28 VALUE SET Y3 to 30.58 VALUE SET Z3 to 0.173 VALUE SET StrSpc to 7.0 VALUE SET VecSpc to 5.0 VALUE SET VecScl to 2.5 SAMPLER VS_2 U V X2 Y2 Z2 SAMPLER VS_3 U V X3 Y3 Z3 STATS PRINT Z 1.0 PRI PLANE FULL *INDOW DISTANCE-40.0 40.0 -30.0 50.0 WINDOW DISTANCE-35.0 35.0 -20.0 50.0 0.0 4.0 GRAPHICS01*01*01Z 1.0 REORDER SYNTHESIZED WINDS WINDOW VECTOR Usyn Vsyn STREAM StrSpc WHITE VECTOR U V STREAM StrSpc WHITE VECTOR Usyn Vsyn ARROW VecSpc VecScl 1.0 PRO WHITE VECTOR U V ARROW VecSpc VecScl 1.0 PRO WHITE VECTOR DU V ARROW VecSpc VecScl 1.0 PRO BLACK CONTOUR DZ_max FALL 000110531.0 UNI -30.0 70.0 5.0 CONTOUR CONV FALL 001110531.0 UNI -10.5 10.5 1.0 CONTOUR VORT FALL 001110531.0 UNI -10.5 10.5 1.0 CONTOUR U FALL 000110531.0 UNI -25.0 25.0 2.0 CONTOUR V FALL 000110531.0 UNI -25.0 25.0 2.0 CONTOUR SPEED FALL 001110531.0 UNI 0.0 25.0 2.0 CONTOUR DZ_2 FALL 000110531.0 UNI -30.0 70.0 5.0 CONTOUR DZ_3 FALL 000110531.0 UNI -30.0 70.0 5.0 CONTOUR VE_2 FALL 000110531.0 UNI -25.0 25.0 2.0 CONTOUR VS_2 FALL 000110531.0 UNI -25.0 25.0 2.0 CONTOUR VE_3 FALL 000110531.0 UNI -25.0 25.0 2.0 CONTOUR VS_3 FALL 000110531.0 UNI -25.0 25.0 2.0 SCATTER DZ_2 -30.0 70.0 DZ_3 -30.0 70.0 1.0 SCATTER VE_2 -25.0 25.0 VS_2 -25.0 25.0 1.0 SCATTER VE_3 -25.0 25.0 VS_3 -25.0 25.0 1.0 END VALUE OFF --------------------------------------------------- # Part 7: run_cedsyn (FINISH) COMMANDS for FINISHING and WRITING SYNTHESIS OUTPUT --------------------------------------------------- DELETE CT EWU EWV PATCHER DZ_max DZ_max FILLCON 5 4 4 COMMENT ------------------------------------------------------------------------ FUNCTIONNAME 1.0 SUM T + DZDOW2 DZDOW3 1.0 1.0 TWO T CON 2.0 AVGDZ P / SUM TWO 1.0 TEMPDZ T ORELSE DZDOW2 DZDOW3 AVGDZ P ORELSE TEMPDZ AVGDZ END ------------------------------------------------------------------------ COMEND DELETE USTD VSTD *-------*-------*-------*-------*-------*-------*-------*-------*-------*------- *RITVOL LUNIT NAMTAP NAMVOL IPOS NFSKIP FEET FORMAT WRITVOL 90.0 SYNTHES synDOWs BEG DEFEND SYNUVW QUIT cedsyn.inp # # Cleanup after massaging input radial velocities # and doing the dual-Doppler wind synthesis. # Make a directory and move all files for this # synthesis time ($1) into that directory. # set syn_dir=06Aug.$1 mkdir ./$syn_dir /bin/mv gmeta ./$syn_dir/$1-cedsyn.plt /bin/mv cedsyn.log ./$syn_dir/$1-cedsyn.log /bin/mv fort.81 ./$syn_dir/$1-DOW2_cedsyn.ced /bin/mv fort.82 ./$syn_dir/$1-DOW3_cedsyn.ced /bin/mv fort.90 ./$syn_dir/$1-DUAL_cedsyn.ced /bin/mv fort.99 ./$syn_dir/$1-DUAL_cedsyn.log /bin/rm .cededit .cedremap .sync .async #