source: trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/tools/regtest.csh @ 1759

Last change on this file since 1759 was 11, checked in by aslmd, 14 years ago

spiga@svn-planeto:ajoute le modele meso-echelle martien

File size: 117.1 KB
Line 
1#!/bin/csh
2# @ job_type            = parallel
3# @ environment         = COPY_ALL;MP_EUILIB=us
4# @ job_name            = regtest.$(jobid)
5# @ output              = regtest_out
6# @ error               = regtest_err
7# @ network.MPI         = csss,shared,us
8# @ node_usage          = shared
9# @ checkpoint          = no
10# @ wall_clock_limit    = 21600
11# @ node                = 1
12# @ total_tasks         = 4
13###################  NCAR  ########################
14# @ ja_report           = yes
15# @ class               = share
16###################  NCAR  ########################
17###################  NCEP  ########################
18## @ class              = dev
19## @ group              = devqpri
20## @ preferences        = Feature == "dev"
21###################  NCEP  ########################
22# @ queue
23
24# #BSUB -x                                # exlusive use of node (not_shared)
25# #BSUB -a mpich_gm                       # at NCAR: lightning
26# #BSUB -R "span[ptile=2]"                # how many tasks per node (1 or 2)
27#BSUB -a poe                            # at NCAR: bluevista
28#BSUB -R "span[ptile=4]"                # how many tasks per node (up to 8)
29#BSUB -n 4                              # number of total tasks
30#BSUB -o reg.out                        # output filename (%J to add job id)
31#BSUB -e reg.err                        # error filename
32#BSUB -J reg.test                       # job name
33#BSUB -q premium                        # queue
34#BSUB -W 6:00                           # wallclock time
35#BSUB -P 64000400
36
37# QSUB -q ded_4             # submit to 4 proc
38# QSUB -l mpp_p=4           # request 4 processors
39# QSUB -lT  21600           # max. job time limit is 6 h
40# QSUB -lF 250Mw            # max. job file size limit is 250 Megawords
41# QSUB -eo                  # merge error and output into one file
42# QSUB -o  reg.out          # output file name
43# QSUB                      # there are no further QSUB commands
44
45#       This is a script to test the bit-for-bit reproducibility of
46#       the WRF model, when comparing single processor serial runs to
47#       OpenMP and MPI parallel runs.  There are several regression tests
48#       that are performed.  Failed comparisons get reported, but don't
49#       stop the script.  Failed builds or forecasts force an exit from
50#       the script.
51
52#       Approximate time for completion of full test suite
53#               Compaq 733 MHz   ev67 :  2.5 hours (empty)
54#               Intel  1.2 GHz (4-pe) :  3.0 hours (empty)
55#               IBM            P4     :  2.0 hours (empty)
56
57#       These need to be changed for your particular set of runs.  This is
58#       where email gets sent.
59
60if ( ( `uname` == AIX ) && ( ( `hostname | cut -c 1-2` != bs ) && ( `hostname | cut -c 1-2` != bv ) ) ) then
61        set FAIL_MAIL = ( ${user}@noaa.gov )
62        set GOOD_MAIL = ( ${user}@noaa.gov )
63
64        setenv MP_EAGER_LIMIT 65536
65        setenv MP_SHARED_MEMORY yes
66        setenv MP_SINGLE_THREAD yes
67        setenv MP_LABELIO yes
68        setenv MP_STDOUTMODE ordered
69
70        setenv OMP_NUM_THREADS 4
71        setenv XLSMPOPTS "parthds=4:spins=0:yields=0:stack=128000000:schedule=static"
72        setenv AIXTHREAD_SCOPE S
73        setenv AIXTHREAD_MNRATIO 1:1
74        setenv SPINLOOPTIME 1000
75        setenv YIELDLOOPTIME 1000
76else
77        set FAIL_MAIL = ( ${user}@ucar.edu )
78        set GOOD_MAIL = ( ${user}@ucar.edu )
79endif
80
81unalias cd cp rm ls pushd popd mv
82if ( `uname` == Linux ) alias banner echo
83
84#       Get the command line input
85
86set thedate = -999
87set thefile = "null"
88set thedata = "null"
89set clrm = 0          # compile local run mmmtmp, for using clsroom cluster and local disk
90
91#       If this is a batch job (NCAR's IBMs or FSL's Intel and Alpha), we need to muck with the "input"
92#       parameters a bit.
93
94if      ( ( `uname` == AIX ) || ( `hostname` == tempest ) || ( `hostname | cut -c 1-2` == ln ) ) then
95        set argv = ( -here )
96        set argv = ( -ftp )
97        set argv = ( -D today )
98        set argv = ( -env )
99        set WRFREGFILE = /mmm/users/gill/wrf.tar
100        if ( ( `uname` == AIX ) && ( ( `hostname | cut -c 1-2` != bs ) && ( `hostname | cut -c 1-2` != bv ) ) ) then
101                set argv = ( -f /nbns/meso/wx22tb/regression_tests/wrf.tar )
102        else
103                set argv = ( -f wrf.tar )
104        endif
105else if ( ( `uname` == OSF1 ) && ( `hostname` == maple ) && ( $user == michalak ) ) then
106        set clrm=1
107endif
108
109#       Where is the input data located - for a few known NCAR/MMM machines.
110
111if      ( ( `hostname` == master ) || (`hostname | cut -c 1-4` == node ) ) then
112        set WRFREGDATAEM = /big/users/gill/WRF-data-EM
113        set WRFREGDATANMM = /big/users/gill/WRF-data-NMM
114else if   ( `hostname` == jacaranda ) then
115        set WRFREGDATAEM = /jacaranda/users/gill/WRF-data-EM
116        set WRFREGDATANMM = /jacaranda/users/gill/WRF-data-NMM
117else if   ( `hostname` == duku ) then
118        set WRFREGDATAEM = /duku/users/gill/WRF-data-EM
119        set WRFREGDATANMM = /duku/users/gill/WRF-data-NMM
120else if   ( `hostname` == cape ) then
121        set WRFREGDATAEM = /cape/users/michalak/WRF-data-EM
122        set WRFREGDATANMM = /cape/users/michalak/WRF-data-NMM
123else if ( (`hostname | cut -c 1-6` == joshua ) || \
124          ( `hostname` == maple ) || (`hostname | cut -c 1-7` == service ) ) then
125        set WRFREGDATAEM = /users/gill/WRF-data-EM
126        set WRFREGDATANMM = /users/gill/WRF-data-NMM
127else if ( ( `hostname | cut -c 1-2` == bs ) || ( `hostname` == tempest ) || ( `hostname | cut -c 1-2` == ln ) || \
128          ( `hostname | cut -c 1-2` == bv ) ) then
129        set WRFREGDATAEM = /mmm/users/gill/WRF-data-EM
130        set WRFREGDATANMM = /mmm/users/gill/WRF-data-NMM
131else if ( ( `uname` == AIX ) && ( ( `hostname | cut -c 1-2` != bs ) && ( `hostname | cut -c 1-2` != bv ) ) ) then
132        set WRFREGDATAEM = /nbns/meso/wx22tb/regression_tests/WRF-data-EM
133        set WRFREGDATANMM = /nbns/meso/wx22tb/regression_tests/WRF-data-NMM
134else
135        if      ( ( -d /users/gill/WRF-data-EM ) && ( -d /users/gill/WRF-data-NMM ) ) then
136                set WRFREGDATAEM = /users/gill/WRF-data-EM
137                set WRFREGDATANMM = /users/gill/WRF-data-NMM
138        else if ( ( -d /mmm/users/gill/WRF-data-EM ) && ( -d /mmm/users/gill/WRF-data-NMM ) ) then
139                set WRFREGDATAEM = /mmm/users/gill/WRF-data-EM
140                set WRFREGDATANMM = /mmm/users/gill/WRF-data-NMM
141        else
142                echo "stick the WRF em and nmm data somewhere, and then fill in the shell vars"
143                echo "inside this script, you NEED WRFREGDATAEM and WRFREGDATANMM set"
144                exit ( 1 ) 
145        endif
146endif
147#DAVE###################################################
148echo DAVE em data is located at $WRFREGDATAEM
149ls -ls $WRFREGDATAEM
150echo DAVE nmm data is located at $WRFREGDATANMM
151ls -ls $WRFREGDATANMM
152banner 1
153#set ans = "$<"
154#DAVE###################################################
155
156if ( $#argv == 0 ) then
157        echo "Please enter either a date for cvs checkout. ex regtest.csh -D date"
158        echo " or a file name containing WRF. ex regtest.csh -f tarfile"
159        echo " or the -ftp flag for the script to pick code off anon ftp"
160        exit ( 2 ) 
161endif
162
163set theargs = 0
164foreach a ( $argv )
165        if ( "$a" == "-D" ) then
166
167                rsh -n maple.mmm.ucar.edu w >& /dev/null
168                if ( $status ) then
169                        echo "Cannot execute a remote shell on maple.mmm.ucar.edu, where the"
170                        echo "WRF code resides."
171                        echo "Please check that it is up and that you have permission to rsh"
172                        echo "to this host. (Create a .rhosts file)."
173                        ping -c 1 maple.mmm.ucar.edu
174                        exit 2
175                endif
176                setenv CVSROOT maple.mmm.ucar.edu:/data3/mp/wrfhelp/WRF
177               
178                set acquire_from = "cvs"
179                set thedate = $argv[2]
180
181        endif
182
183        if ( "$a" == "-f" ) then
184
185                set thefile = $argv[2]
186                #       Check for absolute path, if not, make it absolute
187                echo $thefile | grep '^/' > /dev/null
188                if ( $status != 0 ) set thefile = `pwd`/$thefile
189                set acquire_from = "filearg"
190
191        endif
192
193        if ( "$a" == "-ftp" ) then
194                set acquire_from = "ftp"
195echo "anon ftp temporarily disabled"
196exit ( 3 )
197        endif
198
199        if ( "$a" == "-here" ) then
200                set acquire_from = "here"
201        endif
202
203        if ( "$a" == "-env" ) then
204                set acquire_from = "environment"
205                set thefile = $WRFREGFILE
206        endif
207end
208
209#       Start recording everything - for debug purposes.
210
211set echo
212set date
213
214#       And to tell us how long we've spent on this whole regression test,
215#       we should remember when we started.
216
217set start = ( `date` )
218
219#####################################################################
220
221#       Initial set up values
222
223#       Is this a single domain regression test or is this nested.  Well, a nested one
224#       is a bit special.  It can only run on machines that have the WRF RSL-but-no-MPI
225#       option available.
226
227set NESTED = TRUE
228set NESTED = FALSE
229
230if ( ( $NESTED == TRUE ) && ( ( `uname` == OSF1 ) || ( `uname` == Linux ) || ( `uname` == AIX ) ) ) then
231        echo DOING a NESTED TEST
232else if ( $NESTED == TRUE ) then
233        echo NESTED option is only valid on DEC, Linux, or AIX machines
234        exit ( 1 ) 
235endif
236
237#       We can choose to do grid and obs nudging tests.
238
239set FDDA = TRUE
240set FDDA = FALSE
241
242#       The default floating point precision is either 4 bytes or 8 bytes.
243#       We assume that it is 4 (or the default for the architecture) unless
244#       REAL8 is set to TRUE.
245
246set REAL8 = TRUE
247set REAL8 = FALSE
248
249#       Are we shooting for a bit-for-bit run (serial vs OpenMP, serial vs MPI), or not?
250#       If you want to do a performance-only run, the forecasts are still short, but you
251#       get to insure that the optimized code builds and runs.
252
253set REG_TYPE = OPTIMIZED
254set REG_TYPE = BIT4BIT
255
256#       For a Linux machine, we can use PGI or Intel compilers.
257
258if ( `uname` == Linux ) then
259        set LINUX_COMP = INTEL
260        set LINUX_COMP = PGI
261endif
262
263#       Intel requires /usr/local/intel-8.0/lib in LD_LIBRARY_PATH
264#       and it has to be in the .cshrc file for the rsh for mpirun
265#       intel 8.0 Apr 2005
266
267if ( `uname` == Linux ) then
268        if ( $LINUX_COMP == INTEL ) then
269       
270                grep LD_LIBRARY_PATH ~/.cshrc >& /dev/null
271                set ok1 = $status
272                echo $LD_LIBRARY_PATH | grep intel >& /dev/null
273                set ok2 = $status
274                if ( ( $ok1 != 0 ) || ( $ok2 != 0 ) ) then
275                        echo You need to stick the following line in your .cshrc file
276                        echo setenv LD_LIBRARY_PATH /usr/lib:/usr/local/lib:/usr/local/intel-8.0/lib
277                        echo Otherwise mpirun cannot find libcxa.so.5
278                        exit ( 1 )
279                endif
280        endif
281endif
282
283#       Is this a WRF chem test? 
284
285if ( $NESTED != TRUE ) then
286        set CHEM = TRUE
287        set CHEM = FALSE
288else if ( $NESTED == TRUE ) then
289        set CHEM = FALSE
290endif
291if ( $CHEM == TRUE ) then
292        setenv WRF_CHEM 1
293else if ( $CHEM == FALSE ) then
294        setenv WRF_CHEM 0
295endif
296
297#       For the real data case, we can run either one of two data cases.  If this is
298#       a chemistry run, we are forced to use that data.
299
300set dataset = jun01
301set dataset = jan00
302if ( $CHEM == TRUE ) then
303        set dataset = chem
304endif
305
306#       Yet another local variable to change the name of where the data is located.
307
308set thedataem = ${WRFREGDATAEM}/${dataset}
309set thedatanmm = $WRFREGDATANMM
310
311#       The distributed memory option RSL_LITE may be selected.
312
313set RSL_LITE = TRUE
314set RSL_LITE = FALSE
315
316set COMBO_NEST_RSL__LITE = TRUE
317set COMBO_NEST_RSL__LITE = FALSE
318
319if ( $COMBO_NEST_RSL__LITE == TRUE ) then
320        set NESTED = TRUE
321        set RSL_LITE = TRUE
322endif
323
324#       A separately installed version of the latest ESMF library (NOT the
325#       ESMF library included in the WRF tarfile) can be tested by setting
326#       "ESMF_LIB" to "TRUE" below.  This test is not supported on all
327#       machines. 
328
329set ESMF_LIB = TRUE
330set ESMF_LIB = FALSE
331
332# serial and OMP are not tested with ESMF so always start with env vars cleared
333unsetenv ESMFLIB
334unsetenv ESMFINC
335
336if ( $ESMF_LIB == TRUE ) then
337        if ( ( `uname` == AIX ) && ( `hostname | cut -c 1-2` == bs ) ) then
338                echo "A separately installed version of the latest ESMF library"
339                echo "(NOT the ESMF library included in the WRF tarfile) will"
340                echo "be used for MPI tests"
341                setenv OBJECT_MODE 32
342                set ESMFLIBSAVE = /home/bluesky/hender/esmf/lib/libO/AIX.default.32.default
343                set ESMFINCSAVE = /home/bluesky/hender/esmf/mod/modO/AIX.default.32.default
344                echo "Setting ESMFLIB = ${ESMFLIBSAVE}"
345                echo "Setting ESMFINC = ${ESMFINCSAVE}"
346        else
347                echo "Only the ESMF library included in the WRF tarfile is"
348                echo "tested on this machine"
349                exit ( 3 ) 
350        endif
351        if ( $NESTED == TRUE ) then
352                echo "The ESMF library does not work with nesting."
353                exit ( 3 ) 
354        endif
355        if ( $RSL_LITE == TRUE ) then
356                echo "The ESMF library does not work with RSL_LITE."
357                exit ( 3 ) 
358        endif
359endif
360
361#       A single WRF output "quilt" server can be tested by setting "QUILT"  to
362#       "TRUE" below.  At the moment, testing of I/O quilt servers is not supported
363#       on all machines. 
364
365set QUILT = TRUE
366set QUILT = FALSE
367
368if ( $QUILT == TRUE ) then
369        if ( `uname` == AIX ) then
370                echo "One WRF output quilt server will be used for some tests"
371        else if ( ( `uname` == OSF1 ) && \
372                  ( ( `hostname` == duku    ) || \
373                    ( `hostname` == joshua1 ) || \
374                    ( `hostname` == joshua3 ) ) ) then
375                echo "One WRF output quilt server will be used for some tests"
376        else
377                echo "WRF output quilt servers are not tested on this machine"
378                exit ( 3 ) 
379        endif
380endif
381
382#       Baseline data sets can be generated and archived or compared against. 
383#       - To generate and archive, set GENERATE_BASELINE to a pathname that can
384#         be created by this script via "mkdir -p $GENERATE_BASELINE".  This
385#         directory must not already exist. 
386#         Set GENERATE_BASELINE = FALSE to avoid baseline generation. 
387#       - To compare with a previously archived baseline, set COMPARE_BASELINE
388#         to an existing directory that contains an archived baseline. 
389#         Set COMPARE_BASELINE = FALSE to avoid baseline comparison. 
390set GENERATE_BASELINE = FALSE
391set COMPARE_BASELINE = FALSE
392
393#       Baseline generation and comparison are only done when BIT4BIT is set. 
394if ( $GENERATE_BASELINE != FALSE ) then
395        if ( $REG_TYPE != BIT4BIT ) then
396                echo "ERROR:  Baseline generation can only be done during BIT4BIT tests."
397                exit ( 3 ) 
398        endif
399        if ( -d $GENERATE_BASELINE ) then
400                echo "ERROR:  Baseline directory ${GENERATE_BASELINE} already exists."
401                exit ( 3 ) 
402        else
403                # Archive serial output file to baseline
404                mkdir -p $GENERATE_BASELINE || ( echo "ERROR:  cannot mkdir ${GENERATE_BASELINE}"; exit 3 )
405        endif
406endif
407if ( $COMPARE_BASELINE != FALSE ) then
408        if ( $REG_TYPE != BIT4BIT ) then
409                echo "Baseline comparison can only be done during BIT4BIT tests."
410                exit ( 3 ) 
411        endif
412        if ( ! -d $COMPARE_BASELINE ) then
413                echo "${0}: ERROR:  Baseline directory ${COMPARE_BASELINE} does not exist"
414                exit ( 3 ) 
415        endif
416endif
417
418#       Set the input/output format type (currently 1, 2 or 5 OK).
419#       Binary          NetCDF                          PHDF, IBM       GriB, history only
420#       1               2               3               4               5
421
422set IO_FORM = 2
423set IO_FORM_NAME = ( io_bin io_netcdf io_dummy io_phdf5 io_grib1 )
424set IO_FORM_WHICH =( IO     IO        IO       IO       O        )
425
426#       There is a breakdown of cores to test depending on the various
427#       options that the user is testing.
428#       nested: cannot test NMM
429#       rsl_lite: cannot test anything with y periodic bc
430#       chem: em_real only
431#       esmf_lib: cannot test NMM
432#       grib output: cannot test NMM
433
434if      ( ( $NESTED == TRUE ) && ( $RSL_LITE != TRUE ) ) then
435        set CORES = ( em_real em_b_wave em_quarter_ss          )
436else if ( ( $NESTED == TRUE ) && ( $RSL_LITE == TRUE ) ) then
437        set CORES = ( em_real )
438else if ( ( $NESTED != TRUE ) && ( $RSL_LITE != TRUE ) ) then
439        set CORES = ( em_real em_b_wave em_quarter_ss nmm_real )
440        if ( $CHEM == TRUE ) then
441                set CORES = ( em_real em_real )
442        endif
443        if ( $ESMF_LIB == TRUE ) then
444                set CORES = ( em_real em_b_wave em_quarter_ss )
445        endif
446        if ( $IO_FORM_NAME[$IO_FORM] == io_grib1 ) then
447                set CORES = ( em_real em_b_wave em_quarter_ss )
448        endif
449        if ( $FDDA == TRUE ) then
450                set CORES = ( em_real )
451        endif
452else if ( ( $NESTED != TRUE ) && ( $RSL_LITE == TRUE ) ) then
453        set CORES = ( em_real em_quarter_ss nmm_real )
454        if ( $CHEM == TRUE ) then
455                set CORES = ( em_real em_real )
456        endif
457        if ( $ESMF_LIB == TRUE ) then
458                set CORES = ( em_real )
459        endif
460        if ( $IO_FORM_NAME[$IO_FORM] == io_grib1 ) then
461                set CORES = ( em_real )
462        endif
463        if ( $FDDA == TRUE ) then
464                set CORES = ( em_real )
465        endif
466endif
467
468#       The b_wave case has binary input (4-byte only), the nmm
469#       core has raw MPI calls, skip them if we are doing real*8 floats.
470#       Bump up the OMP stack size for real*8 on OSF1 architectures.
471
472if      ( $REAL8 == TRUE ) then
473        if ( `uname` == OSF1 ) then
474                setenv MP_STACK_SIZE 64000000
475        endif
476        set CORES = ( em_real em_quarter_ss )
477endif
478
479if      ( ( $CHEM != TRUE ) && ( $FDDA != TRUE ) && ( $NESTED != TRUE ) ) then
480        set PHYSOPTS =  ( 1 2 3 4 )
481else if ( ( $CHEM != TRUE ) && ( $FDDA != TRUE ) && ( $NESTED == TRUE ) ) then
482        set PHYSOPTS =  ( 1 2 3 )
483else if ( ( $CHEM != TRUE ) && ( $FDDA == TRUE ) ) then
484        set PHYSOPTS =  ( 1 2 3 )
485else if ( $CHEM == TRUE ) then
486        set PHYSOPTS =  ( 1 2 3 4 5 6 )
487endif
488
489#       This is selecting the ideal physics options - mostly selecting BC options.
490#       With no nesting, run all three ideal physics options.  When we have
491#       RSL_LITE (we are only doing em_quarter_ss), choose the first option only
492#       since it uses open boundaries.
493
494if      ( ( $NESTED == TRUE ) && ( $RSL_LITE != TRUE ) ) then
495        set Max_Ideal_Physics_Options = 2
496else if   ( $RSL_LITE == TRUE )                          then
497        set Max_Ideal_Physics_Options = 1
498else if ( ( $NESTED != TRUE ) && ( $RSL_LITE != TRUE ) ) then
499        set Max_Ideal_Physics_Options = 3
500endif
501
502set CUR_DIR = `pwd`
503
504#       How many domains to run (nest tests).  Only em_real and ideals use this.
505#       The max is 3 due to the number of columns in the namelist that are
506#       currently filled in.
507
508if      ( $NESTED == TRUE ) then
509if      ( $dataset == jan00 ) then
510cat >! dom_real << EOF
511 time_step                           = 180,
512 time_step_fract_num                 = 0,
513 time_step_fract_den                 = 1,
514 max_dom                             = 2,
515 s_we                                = 1,     1,     1,
516 e_we                                = 74,    31,    31,
517 s_sn                                = 1,     1,     1,
518 e_sn                                = 61,    31,    31,
519 s_vert                              = 1,     1,     1,
520 e_vert                              = 28,    28,    28,
521 dx                                  = 30000, 10000,  3333.333333,
522 dy                                  = 30000, 10000,  3333.333333,
523 grid_id                             = 1,     2,     3,
524 parent_id                           = 0,     1,     2,
525 i_parent_start                      = 0,     31,    11,
526 j_parent_start                      = 0,     17,    11,
527 parent_grid_ratio                   = 1,     3,     3,
528 parent_time_step_ratio              = 1,     3,     3,
529 feedback                            = 1,
530 smooth_option                       = 0
531 num_moves                           = 3
532 move_id                             = 2 , 2 , 2
533 move_interval                       = 3 , 6 , 9
534 move_cd_x                           = 1 , 1 , 1
535 move_cd_y                           = 1 , 1 , 1
536EOF
537else if ( $dataset == jun01 ) then
538cat >! dom_real << EOF
539 time_step                           = 60,
540 time_step_fract_num                 = 0,
541 time_step_fract_den                 = 1,
542 max_dom                             = 2,
543 s_we                                = 1,     1,     1,
544 e_we                                = 91,    31,    31,
545 s_sn                                = 1,     1,     1,
546 e_sn                                = 82,    31,    31,
547 s_vert                              = 1,     1,     1,
548 e_vert                              = 28,    28,    28,
549 dx                                  = 10000,  3333.333333,  1111.111111,
550 dy                                  = 10000,  3333.333333,  1111.111111,
551 grid_id                             = 1,     2,     3,
552 parent_id                           = 0,     1,     2,
553 i_parent_start                      = 0,     30,    11,
554 j_parent_start                      = 0,     20,    11,
555 parent_grid_ratio                   = 1,     3,     3,
556 parent_time_step_ratio              = 1,     3,     3,
557 feedback                            = 1,
558 smooth_option                       = 0
559 num_moves                           = 3
560 move_id                             = 2 , 2 , 2
561 move_interval                       = 1 , 2 , 3
562 move_cd_x                           = 1 , 1 , 1
563 move_cd_y                           = 1 , 1 , 1
564EOF
565endif
566cat >! dom_ideal << EOF
567 max_dom                             = 2,
568EOF
569else if ( $NESTED != TRUE ) then
570if      ( $dataset == jan00 ) then
571cat >! dom_real << EOF
572 time_step                           = 180,
573 time_step_fract_num                 = 0,
574 time_step_fract_den                 = 1,
575 max_dom                             = 1,
576 s_we                                = 1,     1,     1,
577 e_we                                = 74,    31,    31,
578 s_sn                                = 1,     1,     1,
579 e_sn                                = 61,    31,    31,
580 s_vert                              = 1,     1,     1,
581 e_vert                              = 28,    28,    28,
582 dx                                  = 30000, 10000,  3333,
583 dy                                  = 30000, 10000,  3333,
584 grid_id                             = 1,     2,     3,
585 parent_id                           = 0,     1,     2,
586 i_parent_start                      = 0,     31,    30,
587 j_parent_start                      = 0,     17,    30,
588 parent_grid_ratio                   = 1,     3,     3,
589 parent_time_step_ratio              = 1,     3,     3,
590 feedback                            = 1,
591 smooth_option                       = 0
592EOF
593else if ( $dataset == jun01 ) then
594cat >! dom_real << EOF
595 time_step                           = 60,
596 time_step_fract_num                 = 0,
597 time_step_fract_den                 = 1,
598 max_dom                             = 1,
599 s_we                                = 1,     1,     1,
600 e_we                                = 91,    31,    31,
601 s_sn                                = 1,     1,     1,
602 e_sn                                = 82,    31,    31,
603 s_vert                              = 1,     1,     1,
604 e_vert                              = 28,    28,    28,
605 dx                                  = 10000,  3333.333333,  1111.111111,
606 dy                                  = 10000,  3333.333333,  1111.111111,
607 grid_id                             = 1,     2,     3,
608 parent_id                           = 0,     1,     2,
609 i_parent_start                      = 0,     30,    11,
610 j_parent_start                      = 0,     20,    11,
611 parent_grid_ratio                   = 1,     3,     3,
612 parent_time_step_ratio              = 1,     3,     3,
613 feedback                            = 1,
614 smooth_option                       = 0
615EOF
616endif
617cat >! dom_ideal << EOF
618 max_dom                             = 1,
619EOF
620endif
621
622#       The em_real entire physics namelist.  Change what you want.
623
624cat >! phys_real_1  << EOF
625 mp_physics                          = 3,     3,     3,
626 ra_lw_physics                       = 1,     1,     1,
627 ra_sw_physics                       = 1,     1,     1,
628 radt                                = 30,    30,    30,
629 sf_sfclay_physics                   = 1,     1,     1,
630 sf_surface_physics                  = 1,     1,     1,
631 bl_pbl_physics                      = 1,     1,     1,
632 bldt                                = 0,     0,     0,
633 cu_physics                          = 1,     1,     0,
634 cudt                                = 5,     5,     5,
635 isfflx                              = 1,
636 ifsnow                              = 0,
637 icloud                              = 1,
638 surface_input_source                = 1,
639 num_soil_layers                     = 5,
640 mp_zero_out                         = 0,
641 maxiens                             = 1,
642 maxens                              = 3,
643 maxens2                             = 3,
644 maxens3                             = 16,
645 ensdim                              = 144,
646EOF
647
648cat >! dyn_real_SAFE  << EOF
649 pd_moist                            = .false., .false., .false.,
650 pd_scalar                           = .false., .false., .false.,
651 pd_chem                             = .false., .false., .false.,
652 pd_tke                              = .false., .false., .false.,
653EOF
654
655cat >! dyn_real_1  << EOF
656 pd_moist                            = .true.,  .true.,  .false.,
657 pd_scalar                           = .false., .false., .false.,
658 pd_chem                             = .false., .false., .false.,
659 pd_tke                              = .false., .false., .false.,
660EOF
661
662cat >! time_real_1  << EOF
663 auxinput1_inname                    = "met_em.d<domain>.<date>"
664EOF
665
666cat >! phys_real_2 << EOF
667 mp_physics                          = 4,     4,     4,
668 ra_lw_physics                       = 1,     1,     1,
669 ra_sw_physics                       = 2,     2,     2,
670 radt                                = 30,    30,    30,
671 sf_sfclay_physics                   = 2,     2,     2,
672 sf_surface_physics                  = 2,     2,     2,
673 bl_pbl_physics                      = 2,     2,     2,
674 bldt                                = 0,     0,     0,
675 cu_physics                          = 2,     2,     0,
676 cudt                                = 5,     5,     5,
677 isfflx                              = 1,
678 ifsnow                              = 0,
679 icloud                              = 1,
680 surface_input_source                = 1,
681 num_soil_layers                     = 4,
682 mp_zero_out                         = 0,
683 maxiens                             = 1,
684 maxens                              = 3,
685 maxens2                             = 3,
686 maxens3                             = 16,
687 ensdim                              = 144,
688EOF
689
690cat >! dyn_real_2  << EOF
691 pd_moist                            = .true.,  .true.,  .false.,
692 pd_scalar                           = .false., .false., .false.,
693 pd_chem                             = .false., .false., .false.,
694 pd_tke                              = .false., .false., .false.,
695EOF
696
697cat >! time_real_2  << EOF
698 auxinput1_inname                    = "wrf_real_input_em.d<domain>.<date>"
699EOF
700
701cat >! phys_real_3 << EOF
702 mp_physics                          = 5,     5,     5,
703 ra_lw_physics                       = 1,     1,     1,
704 ra_sw_physics                       = 2,     2,     2,
705 radt                                = 30,    30,    30,
706 sf_sfclay_physics                   = 2,     2,     2,
707 sf_surface_physics                  = 3,     3,     3,
708 bl_pbl_physics                      = 2,     2,     2,
709 bldt                                = 0,     0,     0,
710 cu_physics                          = 3,     3,     0,
711 cudt                                = 5,     5,     5,
712 isfflx                              = 1,
713 ifsnow                              = 0,
714 icloud                              = 1,
715 surface_input_source                = 1,
716 num_soil_layers                     = 6,
717 mp_zero_out                         = 0,
718 maxiens                             = 1,
719 maxens                              = 3,
720 maxens2                             = 3,
721 maxens3                             = 16,
722 ensdim                              = 144,
723EOF
724
725cat >! dyn_real_3  << EOF
726 pd_moist                            = .false., .false., .false.,
727 pd_scalar                           = .false., .false., .false.,
728 pd_chem                             = .false., .false., .false.,
729 pd_tke                              = .false., .false., .false.,
730EOF
731
732cat >! time_real_3  << EOF
733 auxinput1_inname                    = "met_em.d<domain>.<date>"
734EOF
735
736cat >! phys_real_4 << EOF
737 mp_physics                          = 4,     4,     4,
738 ra_lw_physics                       = 1,     1,     1,
739 ra_sw_physics                       = 2,     2,     2,
740 radt                                = 30,    30,    30,
741 sf_sfclay_physics                   = 2,     2,     2,
742 sf_surface_physics                  = 2,     2,     2,
743 bl_pbl_physics                      = 2,     2,     2,
744 bldt                                = 0,     0,     0,
745 cu_physics                          = 2,     2,     0,
746 cudt                                = 5,     5,     5,
747 isfflx                              = 1,
748 ifsnow                              = 0,
749 icloud                              = 1,
750 ucmcall                             = 1,
751 surface_input_source                = 1,
752 num_soil_layers                     = 4,
753 mp_zero_out                         = 0,
754 maxiens                             = 1,
755 maxens                              = 3,
756 maxens2                             = 3,
757 maxens3                             = 16,
758 ensdim                              = 144,
759EOF
760
761cat >! dyn_real_4  << EOF
762 pd_moist                            = .false., .false., .false.,
763 pd_scalar                           = .false., .false., .false.,
764 pd_chem                             = .false., .false., .false.,
765 pd_tke                              = .false., .false., .false.,
766EOF
767
768cat >! time_real_4  << EOF
769 auxinput1_inname                    = "wrf_real_input_em.d<domain>.<date>"
770EOF
771
772cat >! phys_real_5 << EOF
773 mp_physics                          = 4,     4,     4,
774 ra_lw_physics                       = 3,     3,     3,
775 ra_sw_physics                       = 3,     3,     3,
776 radt                                = 30,    30,    30,
777 sf_sfclay_physics                   = 2,     2,     2,
778 sf_surface_physics                  = 2,     2,     2,
779 bl_pbl_physics                      = 2,     2,     2,
780 bldt                                = 0,     0,     0,
781 cu_physics                          = 2,     2,     0,
782 cudt                                = 5,     5,     5,
783 isfflx                              = 1,
784 ifsnow                              = 0,
785 icloud                              = 1,
786 ucmcall                             = 1,
787 surface_input_source                = 1,
788 num_soil_layers                     = 4,
789 mp_zero_out                         = 0,
790 maxiens                             = 1,
791 maxens                              = 3,
792 maxens2                             = 3,
793 maxens3                             = 16,
794 ensdim                              = 144,
795 levsiz                              = 59
796 paerlev                             = 29
797 cam_abs_freq_s                      = 21600
798 cam_abs_dim1                        = 4
799 cam_abs_dim2                        = 28
800EOF
801
802cat >! dyn_real_5  << EOF
803 pd_moist                            = .false., .false., .false.,
804 pd_scalar                           = .false., .false., .false.,
805 pd_chem                             = .false., .false., .false.,
806 pd_tke                              = .false., .false., .false.,
807EOF
808
809cat >! time_real_5  << EOF
810 auxinput1_inname                    = "met_em.d<domain>.<date>"
811EOF
812
813cat >! fdda_real_1 << EOF
814 grid_fdda                           = 1,     1,     1,
815 gfdda_inname                        = "wrffdda_d<domain>",
816 gfdda_end_h                         = 24,    24,    24,
817 gfdda_interval_m                    = 360,   360,   360,
818 fgdt                                = 0,     0,     0,
819 if_no_pbl_nudging_uv                = 0,     0,     1,
820 if_no_pbl_nudging_t                 = 0,     0,     1,
821 if_no_pbl_nudging_q                 = 0,     0,     1,
822 if_zfac_uv                          = 0,     0,     1,
823  k_zfac_uv                          = 10,   10,     1,
824 if_zfac_t                           = 0,     0,     1,
825  k_zfac_t                           = 10,   10,     1,
826 if_zfac_q                           = 0,     0,     1,
827  k_zfac_q                           = 10,   10,     1,
828 guv                                 = 0.0003,     0.0003,     0.0003,
829 gt                                  = 0.0003,     0.0003,     0.0003,
830 gq                                  = 0.0003,     0.0003,     0.0003,
831 if_ramping                          = 1,
832 dtramp_min                          = 360.0,
833 io_form_gfdda                       = 2,
834EOF
835
836cat >! fdda_real_time_1 << EOF
837EOF
838
839cat >! fdda_real_2 << EOF
840 obs_nudge_opt                       = 1,1,1,1,1
841 max_obs                             = 150000,
842 nobs_ndg_vars                       = 5,
843 nobs_err_flds                       = 9,
844 obs_nudge_wind                      = 1,1,1,1,1
845 obs_coef_wind                       = 6.E-4,6.E-4,6.E-4,6.E-4,6.E-4
846 obs_nudge_temp                      = 1,1,1,1,1
847 obs_coef_temp                       = 6.E-4,6.E-4,6.E-4,6.E-4,6.E-4
848 obs_nudge_mois                      = 1,1,1,1,1
849 obs_coef_mois                       = 6.E-4,6.E-4,6.E-4,6.E-4,6.E-4
850 obs_rinxy                           = 240.,240.,180.,180,180
851 obs_rinsig                          = 0.1,
852 obs_twindo                          = 40.
853 obs_npfi                            = 10,
854 obs_ionf                            = 2,
855 obs_idynin                          = 0,
856 obs_dtramp                          = 40.,
857 obs_ipf_errob                       = .true.
858 obs_ipf_nudob                       = .true.
859 obs_ipf_in4dob                      = .true.
860EOF
861
862cat >! fdda_real_time_2 << EOF
863 auxinput11_interval_s               = 180
864 auxinput11_end_h                    = 6
865EOF
866
867cat >! fdda_real_3 << EOF
868 grid_fdda                           = 1,     1,     1,
869 gfdda_inname                        = "wrffdda_d<domain>",
870 gfdda_end_h                         = 24,    24,    24,
871 gfdda_interval_m                    = 360,   360,   360,
872 fgdt                                = 0,     0,     0,
873 if_no_pbl_nudging_uv                = 0,     0,     1,
874 if_no_pbl_nudging_t                 = 0,     0,     1,
875 if_no_pbl_nudging_q                 = 0,     0,     1,
876 if_zfac_uv                          = 0,     0,     1,
877  k_zfac_uv                          = 10,   10,     1,
878 if_zfac_t                           = 0,     0,     1,
879  k_zfac_t                           = 10,   10,     1,
880 if_zfac_q                           = 0,     0,     1,
881  k_zfac_q                           = 10,   10,     1,
882 guv                                 = 0.0003,     0.0003,     0.0003,
883 gt                                  = 0.0003,     0.0003,     0.0003,
884 gq                                  = 0.0003,     0.0003,     0.0003,
885 if_ramping                          = 1,
886 dtramp_min                          = 360.0,
887 io_form_gfdda                       = 2,
888 obs_nudge_opt                       = 1,1,1,1,1
889 max_obs                             = 150000,
890 nobs_ndg_vars                       = 5,
891 nobs_err_flds                       = 9,
892 obs_nudge_wind                      = 1,1,1,1,1
893 obs_coef_wind                       = 6.E-4,6.E-4,6.E-4,6.E-4,6.E-4
894 obs_nudge_temp                      = 1,1,1,1,1
895 obs_coef_temp                       = 6.E-4,6.E-4,6.E-4,6.E-4,6.E-4
896 obs_nudge_mois                      = 1,1,1,1,1
897 obs_coef_mois                       = 6.E-4,6.E-4,6.E-4,6.E-4,6.E-4
898 obs_rinxy                           = 240.,240.,180.,180,180
899 obs_rinsig                          = 0.1,
900 obs_twindo                          = 40.
901 obs_npfi                            = 10,
902 obs_ionf                            = 2,
903 obs_idynin                          = 0,
904 obs_dtramp                          = 40.,
905 obs_ipf_errob                       = .true.
906 obs_ipf_nudob                       = .true.
907 obs_ipf_in4dob                      = .true.
908EOF
909
910cat >! fdda_real_time_3 << EOF
911 auxinput11_interval_s               = 180
912 auxinput11_end_h                    = 6
913EOF
914
915#       Tested options for ideal case em_b_wave.  Modifying these
916#       parameters is acceptable.  Adding to these requires changes
917#       to the ideal namelist build below.
918
919cat >! phys_b_wave_1a << EOF
920 diff_opt                            = 1,
921 km_opt                              = 1,
922 damp_opt                            = 0,
923EOF
924cat >! phys_b_wave_1b << EOF
925 mp_physics                          = 1,     1,     1,
926EOF
927cat >! phys_b_wave_1c << EOF
928 non_hydrostatic                     = .true., .true., .true.,
929EOF
930
931cat >! phys_b_wave_2a << EOF
932 diff_opt                            = 1,
933 km_opt                              = 1,
934 damp_opt                            = 0,
935EOF
936cat >! phys_b_wave_2b << EOF
937 mp_physics                          = 1,     1,     1,
938EOF
939cat >! phys_b_wave_2c << EOF
940 non_hydrostatic                     = .false., .false., .false.,
941EOF
942
943cat >! phys_b_wave_3a << EOF
944 diff_opt                            = 1,
945 km_opt                              = 1,
946 damp_opt                            = 0,
947EOF
948cat >! phys_b_wave_3b << EOF
949 mp_physics                          = 2,     2,     2,
950EOF
951cat >! phys_b_wave_3c << EOF
952 non_hydrostatic                     = .false., .false., .false.,
953EOF
954
955#       Tested options for ideal case em_quarter_ss.  Modifying these
956#       parameters is acceptable.  Adding to these requires changes
957#       to the ideal namelist build below.
958
959cat >! phys_quarter_ss_1a << EOF
960 diff_opt                            = 1,
961 km_opt                              = 1,
962 damp_opt                            = 0,
963EOF
964cat >! phys_quarter_ss_1b << EOF
965 mp_physics                          = 1,     1,     1,
966EOF
967cat >! phys_quarter_ss_1c << EOF
968 non_hydrostatic                     = .true., .true., .true.,
969EOF
970cat >! phys_quarter_ss_1d << EOF
971 periodic_x                          = .false.,.false.,.false.,
972 open_xs                             = .true., .false.,.false.,
973 open_xe                             = .true., .false.,.false.,
974 periodic_y                          = .false.,.false.,.false.,
975 open_ys                             = .true., .false.,.false.,
976 open_ye                             = .true., .false.,.false.,
977EOF
978
979cat >! phys_quarter_ss_2a << EOF
980 diff_opt                            = 2,
981 km_opt                              = 2,
982 damp_opt                            = 1,
983EOF
984cat >! phys_quarter_ss_2b << EOF
985 mp_physics                          = 1,     1,     1,
986EOF
987cat >! phys_quarter_ss_2c << EOF
988 non_hydrostatic                     = .true., .true., .true.,
989EOF
990cat >! phys_quarter_ss_2d << EOF
991 periodic_x                          = .true., .false.,.false.,
992 open_xs                             = .false.,.false.,.false.,
993 open_xe                             = .false.,.false.,.false.,
994 periodic_y                          = .true., .false.,.false.,
995 open_ys                             = .false.,.false.,.false.,
996 open_ye                             = .false.,.false.,.false.,
997EOF
998
999cat >! phys_quarter_ss_3a << EOF
1000 diff_opt                            = 2,
1001 km_opt                              = 3,
1002 damp_opt                            = 1,
1003EOF
1004cat >! phys_quarter_ss_3b << EOF
1005 mp_physics                          = 2,     2,     2,
1006EOF
1007cat >! phys_quarter_ss_3c << EOF
1008 non_hydrostatic                     = .false., .false., .false.,
1009EOF
1010cat >! phys_quarter_ss_3d << EOF
1011 periodic_x                          = .true., .false.,.false.,
1012 open_xs                             = .false.,.false.,.false.,
1013 open_xe                             = .false.,.false.,.false.,
1014 periodic_y                          = .true., .false.,.false.,
1015 open_ys                             = .false.,.false.,.false.,
1016 open_ye                             = .false.,.false.,.false.,
1017EOF
1018
1019if      ( $IO_FORM_WHICH[$IO_FORM] == IO ) then
1020cat >! io_format << EOF
1021 io_form_history                     = $IO_FORM
1022 io_form_restart                     = $IO_FORM
1023 io_form_input                       = $IO_FORM
1024 io_form_boundary                    = $IO_FORM
1025EOF
1026else if ( $IO_FORM_WHICH[$IO_FORM] == I  ) then
1027cat >! io_format << EOF
1028 io_form_history                     = 2
1029 io_form_restart                     = 2
1030 io_form_input                       = $IO_FORM
1031 io_form_boundary                    = $IO_FORM
1032EOF
1033else if ( $IO_FORM_WHICH[$IO_FORM] ==  O ) then
1034cat >! io_format << EOF
1035 io_form_history                     = $IO_FORM
1036 io_form_restart                     = 2
1037 io_form_input                       = 2
1038 io_form_boundary                    = 2
1039EOF
1040endif
1041
1042
1043if ( $dataset == jun01 ) then
1044        set filetag_real=2001-06-11_12:00:00
1045else if ( $dataset == jan00 ) then
1046        set filetag_real=2000-01-24_12:00:00
1047else if ( $dataset == chem ) then
1048        set filetag_real = ( 2006-04-06_00:00:00  2006-04-06_12:00:00 )
1049endif
1050
1051set filetag_ideal=0001-01-01_00:00:00
1052#DAVE###################################################
1053echo did phys, set date to $filetag_real
1054banner 2
1055#set ans = "$<"
1056#DAVE###################################################
1057
1058#####################################################################
1059
1060#       Set up info for particular architectures
1061
1062set ARCH = ( `uname` )
1063
1064set ZAP_SERIAL          = FALSE
1065set ZAP_OPENMP          = FALSE
1066set SERIALRUNCOMMAND    = 
1067set OMPRUNCOMMAND       = 
1068set MPIRUNCOMMANDPOST   = 
1069
1070touch version_info
1071if ( $ARCH[1] == AIX ) then
1072        set DEF_DIR             = $home
1073        set TMPDIR              = /ptmp/$user
1074        # keep stuff out of $HOME and /ptmp/$USER
1075        # this allows multiple regressions tests to run simultaneously
1076        # extend this to other machines later
1077        if      ( ( `hostname | cut -c 1-2` == bs ) && ( ! $?LOADL_JOB_NAME ) ) then
1078                echo "${0}: ERROR::  This batch script must be submitted via"
1079                echo "${0}:          LoadLeveler on an AIX machine\!"
1080                exit
1081        else if   ( `hostname | cut -c 1-2` == bs ) then
1082                set job_id              = `echo ${LOADL_JOB_NAME} | cut -f2 -d'.'`
1083                set DEF_DIR             = /ptmp/$user/wrf_regression.${job_id}
1084                set TMPDIR              = $DEF_DIR
1085                if ( -d $DEF_DIR ) then
1086                        echo "${0}: ERROR::  Directory ${DEF_DIR} exists, please remove it"
1087                        exit ( 1 ) 
1088                else
1089                        mkdir -p $DEF_DIR
1090                        echo "See ${DEF_DIR}/wrftest.output and other files in ${DEF_DIR} for test results"
1091                endif
1092                set CUR_DIR = ${LOADL_STEP_INITDIR}
1093        else if   ( `hostname | cut -c 1-2` == bv ) then
1094                set job_id              = $LSB_JOBID
1095                set DEF_DIR             = /ptmp/$user/wrf_regression.${job_id}
1096                set TMPDIR              = $DEF_DIR
1097                if ( -d $DEF_DIR ) then
1098                        echo "${0}: ERROR::  Directory ${DEF_DIR} exists, please remove it"
1099                        exit ( 1 ) 
1100                else
1101                        mkdir -p $DEF_DIR
1102                        echo "See ${DEF_DIR}/wrftest.output and other files in ${DEF_DIR} for test results"
1103                endif
1104        else if ( ( ( `hostname | cut -c 1-2` != bs ) && ( `hostname | cut -c 1-2` != bv ) ) && ( ! $?LOADL_JOB_NAME ) ) then
1105                echo "${0}: ERROR::  This batch script must be submitted via"
1106                echo "${0}:          LoadLeveler on an AIX machine\!"
1107                exit
1108        else if   ( ( `hostname | cut -c 1-2` != bs ) && ( `hostname | cut -c 1-2` != bv ) ) then
1109                set job_id              = `echo ${LOADL_JOB_NAME} | cut -f2 -d'.'`
1110                set DEF_DIR             = /ptmp/$user/wrf_regression.${job_id}
1111                set TMPDIR              = $DEF_DIR
1112                if ( -d $DEF_DIR ) then
1113                        echo "${0}: ERROR::  Directory ${DEF_DIR} exists, please remove it"
1114                        exit ( 1 ) 
1115                else
1116                        mkdir -p $DEF_DIR
1117                        echo "See ${DEF_DIR}/wrftest.output and other files in ${DEF_DIR} for test results"
1118                endif
1119                set CUR_DIR = ${LOADL_STEP_INITDIR}
1120        endif
1121        if ( ! -d $TMPDIR ) mkdir $TMPDIR
1122        set MAIL                = /usr/bin/mailx
1123        if      ( ( $NESTED == TRUE ) && ( $RSL_LITE != TRUE ) ) then
1124                set COMPOPTS    = ( 9 10 4 )
1125        else if ( ( $NESTED != TRUE ) && ( $RSL_LITE != TRUE ) ) then
1126                set COMPOPTS    = ( 1 2 4 )
1127        else if ( ( $NESTED == TRUE ) && ( $RSL_LITE == TRUE ) ) then
1128                set COMPOPTS    = ( 9 10 3 )
1129        else if ( ( $NESTED != TRUE ) && ( $RSL_LITE == TRUE ) ) then
1130                set COMPOPTS    = ( 1 2 3 )
1131        endif
1132        set Num_Procs           = 4
1133        set OPENMP              = $Num_Procs
1134        setenv MP_PROCS  $Num_Procs
1135        setenv MP_RMPOOL 1
1136        if      ( `hostname | cut -c 1-2` == bs ) then
1137                set MPIRUNCOMMAND       =  poe
1138        else if ( `hostname | cut -c 1-2` == bv ) then
1139                set MPIRUNCOMMAND       =  mpirun.lsf
1140        else if ( ( `hostname | cut -c 1-2` != bs ) && ( `hostname | cut -c 1-2` != bv ) ) then
1141                set MPIRUNCOMMAND       =  poe
1142        endif
1143        if ( $CHEM == TRUE ) then
1144                set ZAP_OPENMP          = TRUE
1145        else if ( $CHEM == FALSE ) then
1146                set ZAP_OPENMP          = FALSE
1147        endif
1148# check compiler version, JM
1149        lslpp -i | grep xlf | grep ' xlfcmp ' | head -1
1150        set xlfvers=`lslpp -i | grep xlf | grep ' xlfcmp ' | head -1 | awk '{print $2}' | sed 's/\...*$//'`
1151        if ( $xlfvers < 10 ) then
1152                set ZAP_OPENMP          = TRUE
1153        endif
1154# end of compiler check, JM
1155        echo "Compiler version info: " >! version_info
1156        echo "FORTRAN:        " `lslpp -l | grep xlfrte | head -1 | awk '{print $1 "   " $2}'` >>! version_info
1157        echo " " >>! version_info
1158        echo "OS version info: " >>! version_info
1159        echo "AIX:            " `lslpp -l | grep bos.mp | head -1 | awk '{print $1 "   " $2}'` >>! version_info
1160        echo " " >>! version_info
1161        setenv MP_SHARED_MEMORY yes
1162else if ( $ARCH[1] == OSF1 && $clrm == 0 ) then
1163        if      ( ( `hostname` == duku )                && ( -d /data1/$user ) ) then
1164                set DEF_DIR     = /data1/$user
1165        else if ( ( `hostname | cut -c 1-6` == joshua ) && ( -d /data3/mp/$user ) ) then
1166                set DEF_DIR     = /data3/mp/${user}/`hostname`
1167                if ( ! -d $DEF_DIR ) mkdir $DEF_DIR
1168        else if ( ( `hostname | cut -c 1-6` == joshua ) && ( -d /data6a/md/$user ) ) then
1169                set DEF_DIR     = /data6a/md/${user}/`hostname`
1170                if ( ! -d $DEF_DIR ) mkdir $DEF_DIR
1171        else
1172                set DEF_DIR     = /mmmtmp/${user}/`hostname`
1173                if ( ! -d $DEF_DIR ) mkdir $DEF_DIR
1174        endif
1175        set TMPDIR              = .
1176        set MAIL                = /usr/bin/mailx
1177        if      ( ( $NESTED == TRUE ) && ( $RSL_LITE != TRUE ) ) then
1178                set COMPOPTS    = ( 2 4 6 )
1179        else if ( ( $NESTED != TRUE ) && ( $RSL_LITE != TRUE ) ) then
1180                set COMPOPTS    = ( 1 3 6 )
1181        else if ( ( $NESTED == TRUE ) && ( $RSL_LITE == TRUE ) ) then
1182                set COMPOPTS    = ( 2 4 11 )
1183        else if ( ( $NESTED != TRUE ) && ( $RSL_LITE == TRUE ) ) then
1184                set COMPOPTS    = ( 1 3 5 )
1185        endif
1186        set Num_Procs           = 4
1187        set OPENMP              = $Num_Procs
1188        cat >! `pwd`/machfile << EOF
1189`hostname`
1190`hostname`
1191`hostname`
1192`hostname`
1193EOF
1194        set Mach = `pwd`/machfile
1195        set SERIALRUNCOMMAND    = 
1196        set OMPRUNCOMMAND       = 
1197        set MPIRUNCOMMAND       = ( /usr/local/mpich/bin/mpirun -np $Num_Procs -machinefile $Mach )
1198        if ( $CHEM == TRUE ) then
1199                set ZAP_OPENMP          = TRUE
1200        else if ( $CHEM == FALSE ) then
1201                set ZAP_OPENMP          = FALSE
1202        endif
1203        echo "Compiler version info: " >! version_info
1204        f90 -version >>&! version_info
1205        echo " " >>! version_info
1206        echo "OS version info: " >>! version_info
1207        uname -a >>&! version_info
1208        echo " " >>! version_info
1209else if ( $ARCH[1] == OSF1 && $clrm == 1 ) then
1210        set DEF_DIR             = /`hostname | cut -d. -f1`/$user
1211        set TMPDIR              = /mmmtmp/$user
1212        set MAIL                = /usr/bin/mailx
1213        if      ( ( $NESTED == TRUE ) && ( $RSL_LITE != TRUE ) ) then
1214                set COMPOPTS    = ( 2 4 6 )
1215        else if ( ( $NESTED != TRUE ) && ( $RSL_LITE != TRUE ) ) then
1216                set COMPOPTS    = ( 1 3 6 )
1217        else if ( ( $NESTED == TRUE ) && ( $RSL_LITE == TRUE ) ) then
1218                set COMPOPTS    = ( 2 4 5 )
1219        else if ( ( $NESTED != TRUE ) && ( $RSL_LITE == TRUE ) ) then
1220                set COMPOPTS    = ( 1 3 5 )
1221        endif
1222        set Num_Procs           = 4
1223        set OPENMP              = 0
1224        set ZAP_OPENMP          = TRUE
1225        cat >! $TMPDIR/machfile << EOF
1226service03
1227service04
1228service05
1229service06
1230EOF
1231        set Mach = $TMPDIR/machfile
1232        set MPIRUNCOMMAND       = ( mpirun -np $Num_Procs -machinefile $Mach )
1233        echo "Compiler version info: " >! version_info
1234        f90 -version >>&! version_info
1235        echo " " >>! version_info
1236        echo "OS version info: " >>! version_info
1237        uname -a >>&! version_info
1238        echo " " >>! version_info
1239else if ( ( $ARCH[1] == Linux ) && ( `hostname` == bay-mmm ) ) then
1240        set DEF_DIR     = /data3/mp/${user}/`hostname`
1241        if ( ! -d $DEF_DIR ) mkdir $DEF_DIR
1242        set TMPDIR              = .
1243        set MAIL                = /bin/mail
1244        if      ( $LINUX_COMP == PGI ) then
1245                if      ( ( $NESTED == TRUE ) && ( $RSL_LITE != TRUE ) ) then
1246                        set COMPOPTS    = ( 2 4 5 )
1247                else if ( ( $NESTED != TRUE ) && ( $RSL_LITE != TRUE ) ) then
1248                        set COMPOPTS    = ( 1 3 5 )
1249                else if ( ( $NESTED == TRUE ) && ( $RSL_LITE == TRUE ) ) then
1250                        set COMPOPTS    = ( 2 4 6 )
1251                else if ( ( $NESTED != TRUE ) && ( $RSL_LITE == TRUE ) ) then
1252                        set COMPOPTS    = ( 1 3 6 )
1253                endif
1254        else if ( $LINUX_COMP == INTEL ) then
1255                if        ( $NESTED == TRUE )                            then
1256                        set COMPOPTS    = ( 8 10 11 )
1257                else if   ( $NESTED != TRUE )                            then
1258                        set COMPOPTS    = ( 7  9 11 )
1259                endif
1260        endif
1261        set Num_Procs           = 2
1262        set OPENMP              = $Num_Procs
1263        cat >! machfile << EOF
1264`hostname`
1265`hostname`
1266`hostname`
1267`hostname`
1268EOF
1269        set Mach                = `pwd`/machfile
1270        if ( $CHEM == TRUE ) then
1271                set ZAP_OPENMP          = TRUE
1272        else if ( $CHEM == FALSE ) then
1273                set ZAP_OPENMP          = FALSE
1274        endif
1275        if ( $LINUX_COMP == INTEL ) then
1276                set ZAP_OPENMP          = TRUE
1277        endif
1278        set MPIRUNCOMMAND       = ( mpirun -np $Num_Procs -machinefile $Mach )
1279        echo "Compiler version info: " >! version_info
1280        if      ( $LINUX_COMP == PGI ) then
1281                pgf90 -V >>&! version_info
1282        else if ( $LINUX_COMP == INTEL ) then
1283                ifort -v >>&! version_info
1284        endif
1285        echo " " >>! version_info
1286        echo "OS version info: " >>! version_info
1287        uname -a >>&! version_info
1288        echo " " >>! version_info
1289else if ( ( $ARCH[1] == Linux ) && ( `hostname | cut -c 1-2` ==  ln ) ) then
1290        set DEF_DIR     = /ptmp/${user}/wrf_regtest
1291        if ( ! -d $DEF_DIR ) mkdir $DEF_DIR
1292        set TMPDIR              = .
1293        set MAIL                = /bin/mail
1294        if      ( $LINUX_COMP == PGI ) then
1295                if      ( ( $NESTED == TRUE ) && ( $RSL_LITE != TRUE ) ) then
1296                        set COMPOPTS    = ( 4 2 3 )
1297                else if ( ( $NESTED != TRUE ) && ( $RSL_LITE != TRUE ) ) then
1298                        set COMPOPTS    = ( 1 2 3 )
1299                endif
1300        endif
1301        set Num_Procs           = 4
1302        set OPENMP              = 2
1303        cat >! machfile << EOF
1304`hostname`
1305`hostname`
1306`hostname`
1307`hostname`
1308EOF
1309        set Mach                = `pwd`/machfile
1310        set ZAP_OPENMP          = TRUE
1311        set MPIRUNCOMMAND       =  mpirun.lsf
1312        echo "Compiler version info: " >! version_info
1313        if      ( $LINUX_COMP == PGI ) then
1314                pgf90 -V >>&! version_info
1315        else if ( $LINUX_COMP == INTEL ) then
1316                ifort -v >>&! version_info
1317        endif
1318        echo " " >>! version_info
1319        echo "OS version info: " >>! version_info
1320        uname -a >>&! version_info
1321        echo " " >>! version_info
1322else if ( ( $ARCH[1] == Linux ) && ( `hostname` == loquat ) ) then
1323        set job_id              = $$
1324        set DEF_DIR             = /loquat2/$user/wrf_regression.${job_id}
1325        set TMPDIR              = $DEF_DIR
1326        if ( -d $DEF_DIR ) then
1327                echo "${0}: ERROR::  Directory ${DEF_DIR} exists, please remove it"
1328                exit ( 1 ) 
1329        else
1330                mkdir -p $DEF_DIR
1331                echo "See directory ${DEF_DIR}/ for wrftest.output and other test results"
1332        endif
1333        set MAIL                = /bin/mail
1334        if      ( $LINUX_COMP == PGI ) then
1335                if      ( ( $NESTED == TRUE ) && ( $RSL_LITE != TRUE ) ) then
1336                        set COMPOPTS    = ( 2 4 5 )
1337                else if ( ( $NESTED != TRUE ) && ( $RSL_LITE != TRUE ) ) then
1338                        set COMPOPTS    = ( 1 3 5 )
1339                else if ( ( $NESTED == TRUE ) && ( $RSL_LITE == TRUE ) ) then
1340                        set COMPOPTS    = ( 2 4 6 )
1341                else if ( ( $NESTED != TRUE ) && ( $RSL_LITE == TRUE ) ) then
1342                        set COMPOPTS    = ( 1 3 6 )
1343                endif
1344        else if ( $LINUX_COMP == INTEL ) then
1345                if        ( $NESTED == TRUE )                            then
1346                        set COMPOPTS    = ( 8 10 11 )
1347                else if   ( $NESTED != TRUE )                            then
1348                        set COMPOPTS    = ( 7  9 11 )
1349                endif
1350        endif
1351        set Num_Procs           = 2
1352        set OPENMP              = $Num_Procs
1353        cat >! machfile << EOF
1354`hostname`
1355`hostname`
1356`hostname`
1357`hostname`
1358EOF
1359        set Mach                = `pwd`/machfile
1360        if ( $CHEM == TRUE ) then
1361                set ZAP_OPENMP          = TRUE
1362        else if ( $CHEM == FALSE ) then
1363                set ZAP_OPENMP          = FALSE
1364        endif
1365        if ( $LINUX_COMP == INTEL ) then
1366                set ZAP_OPENMP          = TRUE
1367        endif
1368        set MPIRUNCOMMAND       = ( mpirun -np $Num_Procs -machinefile $Mach )
1369        echo "Compiler version info: " >! version_info
1370        if      ( $LINUX_COMP == PGI ) then
1371                pgf90 -V >>&! version_info
1372        else if ( $LINUX_COMP == INTEL ) then
1373                ifort -v >>&! version_info
1374        endif
1375        echo " " >>! version_info
1376        echo "OS version info: " >>! version_info
1377        uname -a >>&! version_info
1378        echo " " >>! version_info
1379else if ( `hostname` == tempest ) then
1380        set DEF_DIR     = /ptmp/${user}/wrf_regtest.${QSUB_REQID}
1381        if ( ! -d $DEF_DIR ) mkdir $DEF_DIR
1382        set TMPDIR              = .
1383        set MAIL                = /usr/sbin/Mail
1384        set COMPOPTS            = ( 1 2 3 )
1385        set Num_Procs           = 2
1386        set OPENMP              = $Num_Procs
1387        set Mach                = `pwd`/machfile
1388        set ZAP_OPENMP          = TRUE
1389        set MPIRUNCOMMAND       = ( mpirun -np $Num_Procs )
1390        echo "Compiler version info: " >! version_info
1391        f90 -version >>&! version_info
1392        echo " " >>! version_info
1393        echo "OS version info: " >>! version_info
1394        uname -a >>&! version_info
1395        echo " " >>! version_info
1396else if ( ( $ARCH[1] == Linux ) && ( `hostname` == master ) ) then
1397        set DEF_DIR             = /big6/gill/DO_NOT_REMOVE_DIR
1398        set TMPDIR              = .
1399        set MAIL                = /bin/mail
1400        if      ( $LINUX_COMP == PGI ) then
1401                if        ( $NESTED == TRUE )                            then
1402                        set COMPOPTS    = ( 2 4 5 )
1403                else if ( ( $NESTED != TRUE ) && ( $RSL_LITE == TRUE ) ) then
1404                        set COMPOPTS    = ( 1 3 6 )
1405                else if ( ( $NESTED != TRUE ) && ( $RSL_LITE != TRUE ) ) then
1406                        set COMPOPTS    = ( 1 3 5 )
1407                endif
1408        else if ( $LINUX_COMP == INTEL ) then
1409                if        ( $NESTED == TRUE )                            then
1410                        set COMPOPTS    = ( 8 10 11 )
1411                else if   ( $NESTED != TRUE )                            then
1412                        set COMPOPTS    = ( 7  9 11 )
1413                endif
1414        endif
1415        set Num_Procs           = 4
1416        set OPENMP              = 2
1417        cat >! machfile << EOF
1418node3
1419node3
1420node4
1421node4
1422EOF
1423        set Mach                = `pwd`/machfile
1424        if ( $CHEM == TRUE ) then
1425                set ZAP_OPENMP          = TRUE
1426        else if ( $CHEM == FALSE ) then
1427                set ZAP_OPENMP          = FALSE
1428        endif
1429        if ( $LINUX_COMP == INTEL ) then
1430                set ZAP_OPENMP          = TRUE
1431        endif
1432        set MPIRUNCOMMAND       = ( mpirun -v -np $Num_Procs -machinefile $Mach -nolocal )
1433        set MPIRUNCOMMANDPOST   = "< /dev/null"
1434        echo "Compiler version info: " >! version_info
1435        if      ( $LINUX_COMP == PGI ) then
1436                pgf90 -V >>&! version_info
1437        else if ( $LINUX_COMP == INTEL ) then
1438                ifort -v >>&! version_info
1439        endif
1440        echo " " >>! version_info
1441        echo "OS version info: " >>! version_info
1442        uname -a >>&! version_info
1443        echo " " >>! version_info
1444else
1445        echo "Unrecognized architecture for regression test"  >! error_message
1446        echo `uname`                                          >> error_message
1447        echo `hostname`                                       >> error_message
1448        $MAIL -s "Unknown architecture $ARCH[1] " $FAIL_MAIL   < error_message
1449        exit ( 1 )
1450endif
1451
1452#####################################################################
1453#DAVE###################################################
1454echo did the arch specific stuff
1455banner 3
1456#set ans = "$<"
1457#DAVE###################################################
1458
1459#       First of all, in which particular directory do we start.
1460
1461cd $DEF_DIR
1462
1463#       We want to keep the old regression stuff around
1464
1465if ( -d regression_test ) then
1466        if ( -d regression_test.old ) then
1467                /bin/rm -fr regression_test.old
1468        endif
1469        /bin/mv regression_test regression_test.old
1470endif
1471
1472#       Go to the regression test directory
1473
1474mkdir regression_test
1475set ok = $status
1476if ( $ok != 0 ) then
1477        echo "Gee, I cannot make a directory in $DEF_DIR"  >! error_message
1478        echo `pwd`                                         >> error_message
1479        echo `\ls -ls`                                     >> error_message
1480        $MAIL -s "$DEF_DIR not writable $ARCH[1] " $FAIL_MAIL < error_message
1481        exit ( 1 )
1482else
1483        pushd regression_test
1484endif
1485
1486if      ( $acquire_from == "cvs" ) then
1487
1488        #       Checkout the most recent version of WRF from the NCAR cvs repository,
1489        #       and pick up the required input data from the anonymous ftp site.
1490
1491        cvs checkout -D $thedate WRFV2
1492        find ./WRFV2 -exec touch \{\} \;
1493        ftp -n ftp.ucar.edu < ftp_script_data
1494
1495
1496else if ( $acquire_from == "filearg" ) then
1497
1498        #       A tar file of the WRF source was provided, so that is used, along with
1499        #       the required input data files from the ftp site.
1500
1501        tar xvf $thefile
1502        cd WRFV2
1503        clean -a
1504        cd ..
1505        ftp -n ftp.ucar.edu < ftp_script_data
1506
1507else if ( $acquire_from == "environment" ) then
1508
1509        #       A tar file of WRF is assumed to be available.
1510
1511        tar xvf $thefile
1512
1513endif
1514
1515#       And we can stick the input data where we want, the WRFV2 directory has been created.
1516
1517( cd WRFV2/test/em_real  ; ln -sf $thedataem/* . ) 
1518( cd WRFV2/test/nmm_real ; ln -s $thedatanmm/wrf_real* . ; cp $thedatanmm/namelist.input.regtest . )
1519#DAVE###################################################
1520( cd WRFV2/test/em_real ; ls -ls )
1521( cd WRFV2/test/nmm_real ; ls -ls )
1522banner 4
1523#set ans = "$<"
1524#DAVE###################################################
1525
1526#       John-specific stuff for maple is the else; part of the "using service machines".
1527
1528if ( ! $clrm ) then
1529        pushd WRFV2
1530else
1531        if ( ! -d $TMPDIR ) then
1532                echo something wrong 1
1533        endif
1534        if ( ! -d $TMPDIR/RUN ) then
1535                mkdir $TMPDIR/RUN
1536                /bin/rm -fr $TMPDIR/RUN/*
1537        endif
1538        if ( -d $TMPDIR/RUN ) then
1539                tar cf - ./WRFV2/test ./WRFV2/main | ( cd $TMPDIR/RUN ; tar xvf - )
1540                pushd WRFV2
1541        else
1542                echo something wrong 2
1543                exit
1544        endif
1545endif
1546
1547#       Here we initialize our output message.
1548
1549if ( -e ${DEF_DIR}/wrftest.output ) rm ${DEF_DIR}/wrftest.output
1550echo "Architecute: $ARCH[1]      machine: `hostname`" >>! ${DEF_DIR}/wrftest.output
1551echo "WRFV2 source from: $acquire_from " >>! ${DEF_DIR}/wrftest.output
1552echo "Number of OpenMP processes to use: $OPENMP" >>! ${DEF_DIR}/wrftest.output
1553echo "Number of MPI    processes to use: $Num_Procs" >>! ${DEF_DIR}/wrftest.output
1554set name = ( `grep ^${user}: /etc/passwd | cut -d: -f5` ) 
1555echo "Test conducted by $name" >>! ${DEF_DIR}/wrftest.output
1556echo " " >>! ${DEF_DIR}/wrftest.output
1557echo "Test run from directory ${CUR_DIR}" >>! ${DEF_DIR}/wrftest.output
1558echo " " >>! ${DEF_DIR}/wrftest.output
1559if ( $?LOADL_JOB_NAME ) then
1560        echo "Loadlever job name = ${LOADL_JOB_NAME}" >>! ${DEF_DIR}/wrftest.output
1561        echo " " >>! ${DEF_DIR}/wrftest.output
1562endif
1563echo "Real data case for EM is from $dataset " >>! ${DEF_DIR}/wrftest.output
1564echo " " >>! ${DEF_DIR}/wrftest.output
1565echo "The em real and ideal forecasts will be nested: $NESTED " >>! ${DEF_DIR}/wrftest.output
1566echo " " >>! ${DEF_DIR}/wrftest.output
1567if ( $REG_TYPE == BIT4BIT ) then
1568        echo "This is a bit-wise (traditional) regression test. " >>! ${DEF_DIR}/wrftest.output
1569        echo " " >>! ${DEF_DIR}/wrftest.output
1570else if ( $REG_TYPE == OPTIMIZED ) then
1571        echo "This is a fully optimized regression test. " >>! ${DEF_DIR}/wrftest.output
1572        echo "No inter-comparisons are made. " >>! ${DEF_DIR}/wrftest.output
1573        echo " " >>! ${DEF_DIR}/wrftest.output
1574endif
1575if ( $REAL8 == TRUE ) then
1576        echo "Floating point precision is 8-bytes" >>! ${DEF_DIR}/wrftest.output
1577        echo " " >>! ${DEF_DIR}/wrftest.output
1578endif
1579if ( $CHEM == TRUE ) then
1580        echo "WRF_CHEM tests run for em_real core only, no other cores run" >>! ${DEF_DIR}/wrftest.output
1581        echo " " >>! ${DEF_DIR}/wrftest.output
1582endif
1583if ( $RSL_LITE == TRUE ) then
1584        echo "Parallel DM portion using RSL_LITE build option" >>! ${DEF_DIR}/wrftest.output
1585        echo " " >>! ${DEF_DIR}/wrftest.output
1586else if ( $RSL_LITE != TRUE ) then
1587        echo "Parallel DM portion using RSL build option" >>! ${DEF_DIR}/wrftest.output
1588        echo " " >>! ${DEF_DIR}/wrftest.output
1589endif
1590if ( $ESMF_LIB == TRUE ) then
1591        echo "A separately installed version of the latest ESMF library" >>! ${DEF_DIR}/wrftest.output
1592        echo "(NOT the ESMF library included in the WRF tarfile) will" >>! ${DEF_DIR}/wrftest.output
1593        echo "be used for MPI tests" >>! ${DEF_DIR}/wrftest.output
1594        echo "Setting ESMFLIB = ${ESMFLIBSAVE}" >>! ${DEF_DIR}/wrftest.output
1595        echo "Setting ESMFINC = ${ESMFINCSAVE}" >>! ${DEF_DIR}/wrftest.output
1596        echo " " >>! ${DEF_DIR}/wrftest.output
1597endif
1598if ( $QUILT == TRUE ) then
1599        echo "One WRF output quilt server will be used for some tests" >>! ${DEF_DIR}/wrftest.output
1600        echo " " >>! ${DEF_DIR}/wrftest.output
1601endif
1602if ( $FDDA == TRUE ) then
1603        echo "Running FDDA tests" >>! ${DEF_DIR}/wrftest.output
1604        echo " " >>! ${DEF_DIR}/wrftest.output
1605endif
1606if ( $GENERATE_BASELINE != FALSE ) then
1607        echo "WRF output will be archived in baseline directory ${GENERATE_BASELINE} for some tests" >>! \
1608             ${DEF_DIR}/wrftest.output
1609        echo " " >>! ${DEF_DIR}/wrftest.output
1610endif
1611if ( $COMPARE_BASELINE != FALSE ) then
1612        echo "WRF output will be compared with files in baseline directory ${COMPARE_BASELINE} for some tests" >>! \
1613             ${DEF_DIR}/wrftest.output
1614        echo " " >>! ${DEF_DIR}/wrftest.output
1615endif
1616echo "The selected I/O option is $IO_FORM ($IO_FORM_NAME[$IO_FORM])" >>! ${DEF_DIR}/wrftest.output
1617if      ( $IO_FORM_WHICH[$IO_FORM] == IO ) then
1618        echo "This option is for both input and history files" >>! ${DEF_DIR}/wrftest.output
1619        echo " " >>! ${DEF_DIR}/wrftest.output
1620else if ( $IO_FORM_WHICH[$IO_FORM] == I  ) then
1621        echo "This option is for input files only" >>! ${DEF_DIR}/wrftest.output
1622        echo " " >>! ${DEF_DIR}/wrftest.output
1623else if ( $IO_FORM_WHICH[$IO_FORM] ==  O ) then
1624        echo "This option is for history files only" >>! ${DEF_DIR}/wrftest.output
1625        echo " " >>! ${DEF_DIR}/wrftest.output
1626endif
1627
1628cat ${CUR_DIR}/version_info >>! ${DEF_DIR}/wrftest.output
1629
1630#       There are three WRF em executables to be considered that can run in threaded and
1631#       distributed memory.  The 2d hills and 2d squall lines cannot be parallelized with
1632#       MPI, and are therefore not considered in this shell.  The nmm is only run with
1633#       distributed memory (1 vs 4 procs).
1634
1635set first_time_in = TRUE
1636foreach core ( $CORES )
1637#DAVE###################################################
1638echo doing core $core
1639banner 5
1640#set ans = "$<"
1641#DAVE###################################################
1642
1643        #       Some sleight of hand is required for the chemistry tests because we need to
1644        #       build it twice.  But normally, we only build with different real vs ideal, or em vs nmm.
1645        #       What to do, what to do?  Well, we ask for em_real TWICE.  The first time we build without chemistry
1646        #       activated, the second time with it activated.  The first time, we do a single
1647        #       test, the second time through, we do the other 5 tests.
1648
1649        if      ( ( $CHEM == TRUE ) && ( ${#CORES} == 2 ) && ( $first_time_in == TRUE ) ) then
1650                setenv WRF_CHEM 0
1651                set PHYSOPTS =  ( 1 )
1652                set first_time_in = FALSE
1653        else if ( ( $CHEM == TRUE ) && ( ${#CORES} == 2 ) && ( $first_time_in != TRUE ) ) then
1654                setenv WRF_CHEM 1
1655                set PHYSOPTS =  ( 2 3 4 5 6 )
1656        endif
1657
1658        #       Cores to test.
1659
1660        set ZAP_SERIAL_FOR_THIS_CORE          = FALSE
1661        set ZAP_OPENMP_FOR_THIS_CORE          = FALSE
1662        if      ( `echo $core | cut -c 1-2` == em ) then
1663                setenv WRF_EM_CORE      1
1664                setenv WRF_NMM_CORE     0
1665                setenv WRF_COAMPS_CORE  0
1666                setenv WRF_EXP_CORE     0
1667                set ZAP_SERIAL_FOR_THIS_CORE          = FALSE
1668                set ZAP_OPENMP_FOR_THIS_CORE          = FALSE
1669        else if ( `echo $core | cut -c 1-3` == nmm ) then
1670                setenv WRF_EM_CORE      0
1671                setenv WRF_NMM_CORE     1
1672                setenv WRF_COAMPS_CORE  0
1673                setenv WRF_EXP_CORE     0
1674                set ZAP_SERIAL_FOR_THIS_CORE          = TRUE
1675                set ZAP_OPENMP_FOR_THIS_CORE          = TRUE
1676        endif
1677
1678        #       Here we are looping over all of the various compilation configurations,
1679        #       such as serial only, OpenMP only, MPI only, etc.  Each architecture
1680        #       has its own list of these options.  We build each of the executables for
1681        #       this particular ${core}.
1682       
1683        foreach compopt ( $COMPOPTS )
1684#DAVE###################################################
1685echo doing compile option $compopt
1686banner 6
1687#set ans = "$<"
1688#DAVE###################################################
1689
1690                #       We sometimes are interested in bypassing the OpenMP option.
1691
1692                if ( $compopt == $COMPOPTS[2] ) then
1693                        if ( $ZAP_OPENMP == TRUE || $ZAP_OPENMP_FOR_THIS_CORE == TRUE ) then
1694                                goto GOT_THIS_EXEC
1695                        endif
1696                endif
1697
1698                #       NMM only runs parallel
1699                if ( $compopt == $COMPOPTS[1] ) then
1700                        if ( $ZAP_SERIAL == TRUE || $ZAP_SERIAL_FOR_THIS_CORE == TRUE ) then
1701                                goto GOT_THIS_EXEC
1702                        endif
1703                endif
1704
1705                if ( `uname` == AIX ) goto BUILD_REGARDLESS
1706       
1707                #       Did we already build this one?
1708
1709                if  ( $core == em_real )  then
1710                        if      ( ( $compopt == $COMPOPTS[1] ) && \
1711                                  ( -e main/wrf_${core}.exe.$compopt ) && \
1712                                  ( -e main/real_${core}.exe.1 ) && \
1713                                  ( -e ${DEF_DIR}/regression_test/WRFV2/external/$IO_FORM_NAME[$IO_FORM]/diffwrf ) ) then
1714                                goto GOT_THIS_EXEC
1715                        else if ( ( $compopt != $COMPOPTS[1] ) && \
1716                                  ( -e main/wrf_${core}.exe.$compopt ) && \
1717                                  ( -e ${DEF_DIR}/regression_test/WRFV2/external/$IO_FORM_NAME[$IO_FORM]/diffwrf ) ) then
1718                                goto GOT_THIS_EXEC
1719                        endif
1720                else
1721                        if      ( ( $compopt == $COMPOPTS[1] ) && \
1722                                  ( -e main/wrf_${core}.exe.$compopt ) && \
1723                                  ( -e main/ideal_${core}.exe.1 ) && \
1724                                  ( -e ${DEF_DIR}/regression_test/WRFV2/external/$IO_FORM_NAME[$IO_FORM]/diffwrf ) ) then
1725                                goto GOT_THIS_EXEC
1726                        else if ( ( $compopt != $COMPOPTS[1] ) && \
1727                                  ( -e main/wrf_${core}.exe.$compopt ) && \
1728                                  ( -e ${DEF_DIR}/regression_test/WRFV2/external/$IO_FORM_NAME[$IO_FORM]/diffwrf ) ) then
1729                                goto GOT_THIS_EXEC
1730                        endif
1731                endif
1732
1733                BUILD_REGARDLESS:
1734       
1735                #       The WRF configuration file works with a single integer
1736                #       input, which is the compiler option.  By convention, option $COMPOPTS[1] is
1737                #       serial, $COMPOPTS[2] is OMP, and $COMPOPTS[3] is MPI.
1738
1739                #       Print info about use of separately installed ESMF library. 
1740                set esmf_lib_str = " - - - - - - - - - - - - - "
1741                if ( $ESMF_LIB == TRUE ) then
1742                        # only test ESMF with MPI
1743                        if ( $compopt == $COMPOPTS[3] ) then
1744                                echo "A separately installed version of the latest ESMF library" >>! ${DEF_DIR}/wrftest.output
1745                                echo "(NOT the ESMF library included in the WRF tarfile) is" >>! ${DEF_DIR}/wrftest.output
1746                                echo "being used for this test of $core parallel $compopt..." >>! ${DEF_DIR}/wrftest.output
1747                                set esmf_lib_str = "using separate ESMF library"
1748                                setenv ESMFLIB $ESMFLIBSAVE
1749                                setenv ESMFINC $ESMFINCSAVE
1750                        else
1751                                unsetenv ESMFLIB
1752                                unsetenv ESMFINC
1753                        endif
1754                endif
1755       
1756#DAVE###################################################
1757echo start build mechanism
1758banner 7
1759#set ans = "$<"
1760#DAVE###################################################
1761                ./clean -a
1762                echo $compopt | ./configure
1763       
1764                #       Decide whether this a bit-for-bit run or an fully optimized run.  We are just
1765                #       tinkering with the configure.wrf file optimization here.
1766
1767                if ( $REG_TYPE == BIT4BIT ) then
1768                        if ( `uname` == AIX ) then
1769                                if ( ( $compopt == $COMPOPTS[1] ) || ( $compopt == $COMPOPTS[3] ) ) then
1770                                        sed -e '/^OMP/d' -e '/^FCOPTIM/d' -e '/^FCDEBUG/s/#/-O0 /g' ./configure.wrf >! foo ; /bin/mv foo configure.wrf
1771                                else
1772                                        sed -e '/^OMP/s/noauto/noauto:noopt/' \
1773                                            -e '/^FCOPTIM/d' -e '/^FCDEBUG/s/#/-O0 /g' ./configure.wrf >! foo ; /bin/mv foo configure.wrf
1774                                endif
1775                                sed -e 's/-lmassv//g'  -e 's/-lmass//g'  -e 's/-DNATIVE_MASSV//g' -e '/^FCBASEOPTS/s/#//g' ./configure.wrf >! foo ; /bin/mv foo configure.wrf
1776                        else if ( `uname` == Linux ) then
1777                                if ( ( $compopt == $COMPOPTS[1] ) || ( $compopt == $COMPOPTS[3] ) ) then
1778                                        sed -e '/^OMP/d' -e '/^FCOPTIM/d' -e '/^FCDEBUG/s/#-g/-g/g' ./configure.wrf >! foo ; /bin/mv foo configure.wrf
1779                                else
1780                                        sed              -e '/^FCOPTIM/d' -e '/^FCDEBUG/s/#-g/-g/g' ./configure.wrf >! foo ; /bin/mv foo configure.wrf
1781                                endif
1782                                sed -e '/^#PGI  /s/#PGI /       /g' ./configure.wrf >! foo ; /bin/mv foo configure.wrf
1783                        else if ( `uname` == OSF1 ) then
1784                                if ( ( $compopt == $COMPOPTS[1] ) || ( $compopt == $COMPOPTS[3] ) ) then
1785                                        sed -e '/^OMP/d' -e '/^FCOPTIM/d' -e '/^FCDEBUG/s/# -g/-O0/g' -e '/^FCDEBUG/s/# -O0/-O0/g' ./configure.wrf >! foo ; /bin/mv foo configure.wrf
1786                                else
1787                                        sed              -e '/^FCOPTIM/d' -e '/^FCDEBUG/s/# -g/-O0/g' -e '/^FCDEBUG/s/# -O0/-O0/g' ./configure.wrf >! foo ; /bin/mv foo configure.wrf
1788                                endif
1789                        else
1790                                if ( ( $compopt == $COMPOPTS[1] ) || ( $compopt == $COMPOPTS[3] ) ) then
1791                                        sed -e '/^OMP/d' -e '/^FCOPTIM/d' -e '/^FCDEBUG/s/#//g' ./configure.wrf >! foo ; /bin/mv foo configure.wrf
1792                                else
1793                                        sed              -e '/^FCOPTIM/d' -e '/^FCDEBUG/s/#//g' ./configure.wrf >! foo ; /bin/mv foo configure.wrf
1794                                endif
1795                        endif
1796                endif
1797
1798                #       We also need to modify the configure.wrf file based on whether or not there is
1799                #       going to be nesting.  All regtest em_real nesting runs utilize a moving nest.  All of the
1800                #       ARCHFLAGS macros need to have -DMOVE_NESTS added to the list of options.
1801
1802                if ( ( $NESTED == TRUE ) && ( $core == em_real ) ) then
1803                        sed -e '/^ARCHFLAGS/s/=/=       -DMOVE_NESTS/' ./configure.wrf >! foo ; /bin/mv foo configure.wrf
1804                endif
1805
1806                #       The configure.wrf file needs to be adjusted as to whether we are requesting real*4 or real*8
1807                #       as the default floating precision.
1808
1809                if ( $REAL8 == TRUE ) then
1810                        sed -e '/^RWORDSIZE/s/\$(NATIVE_RWORDSIZE)/8/'  configure.wrf > ! foo ; /bin/mv foo configure.wrf
1811                endif
1812
1813#DAVE###################################################
1814echo configure built with optim mods removed, ready to compile
1815banner 8
1816#set ans = "$<"
1817#DAVE###################################################
1818       
1819                # The WRF_SRC_ROOT_DIR hack is only used by the OSF1 build. 
1820                # It works around the annoying fact that in OSF1 $(PWD) does
1821                # not change during execution of regtest.csh, despite the "cd"
1822                # and "pushd" commands. 
1823                setenv WRF_SRC_ROOT_DIR "${DEF_DIR}/regression_test/WRFV2"
1824                #       Build this executable
1825               
1826                ./compile $core
1827#DAVE###################################################
1828echo compile done
1829banner 9
1830#set ans = "$<"
1831#DAVE###################################################
1832       
1833                #       Did the compile work?  Check the expected executable names and locations.
1834
1835                set ok = $status
1836                if ( ! -x main/wrf.exe ) set ok = 1
1837
1838                if      ( ( $core ==  em_real ) && ( $compopt == $COMPOPTS[1] ) ) then
1839                        if ( ! -e main/real.exe ) set ok = 1
1840                else if ( ( $core == nmm_real ) && ( $compopt == $COMPOPTS[3] ) ) then
1841                        if ( ! -e main/real_nmm.exe ) set ok = 1
1842                else if ( $compopt == $COMPOPTS[1] ) then
1843                        if ( ! -e main/ideal.exe ) set ok = 1
1844                endif
1845
1846                if ( ! -x external/io_netcdf/diffwrf ) set ok = 1
1847                if ( ! -x external/io_int/diffwrf ) set ok = 1
1848
1849                if ( $ok != 0 ) then
1850                        echo "SUMMARY compilation    for $core           parallel $compopt $esmf_lib_str FAIL" >>! ${DEF_DIR}/wrftest.output
1851                        $MAIL -s "REGRESSION FAILURE $ARCH[1] " $FAIL_MAIL < ${DEF_DIR}/wrftest.output
1852                        exit ( 3 )
1853                else
1854                        echo "SUMMARY compilation    for $core           parallel $compopt $esmf_lib_str PASS" >>! ${DEF_DIR}/wrftest.output
1855                        echo "-------------------------------------------------------------" >> ${DEF_DIR}/wrftest.output
1856                        mv main/wrf.exe main/wrf_${core}.exe.$compopt
1857                        if      ( ( $core ==  em_real ) && ( $compopt == $COMPOPTS[1] ) ) then
1858                                mv main/real.exe main/real_${core}.exe.1
1859                        else if ( ( $core == nmm_real ) && ( $compopt == $COMPOPTS[3] ) ) then
1860                                mv main/real_nmm.exe main/real_${core}.exe.$COMPOPTS[3]
1861                        else if ( $compopt == $COMPOPTS[1] ) then
1862                                mv main/ideal.exe main/ideal_${core}.exe.1
1863                        endif
1864#DAVE###################################################
1865echo exec exists
1866ls -ls main/*.exe*
1867banner 10
1868#set ans = "$<"
1869#DAVE###################################################
1870                endif
1871
1872                GOT_THIS_EXEC:
1873
1874                if ( $clrm ) then
1875                        cp main/*exe* $TMPDIR/RUN/WRFV2/main
1876                endif
1877       
1878        end
1879
1880        if ( $clrm ) then
1881                pushd $TMPDIR/RUN/WRFV2
1882        endif
1883       
1884        #       We have all of the executables built, now we run'em.  This is a loop
1885        #       over all of the various physics options for this particular
1886        #       ${core}.  Inside the physics loop, we loop over the parallel options.
1887        #       This allows us to use the same WRF input files for each of the parallel
1888        #       choices for a single physics loop index.
1889
1890        foreach phys_option ( $PHYSOPTS )
1891#DAVE###################################################
1892echo which phys option $phys_option
1893banner 11
1894#set ans = "$<"
1895#DAVE###################################################
1896       
1897                #       For each of the executables, we need to run several physics
1898                #       options. 
1899
1900                if (  $core == em_real )  then
1901
1902                        if ( $CHEM != TRUE ) then
1903                                set filetag=$filetag_real
1904                        else if ( $CHEM == TRUE ) then
1905                                if ( $phys_option <= 3 ) then
1906                                        set filetag=$filetag_real[1]
1907                                else
1908                                        set filetag=$filetag_real[2]
1909                                endif
1910                        endif
1911
1912                        foreach compopt ( $COMPOPTS )
1913#DAVE###################################################
1914echo real if filetag is $filetag
1915echo compopt = $compopt
1916banner 12
1917#set ans = "$<"
1918#DAVE###################################################
1919
1920                                #       We sometimes are interested in bypassing the OpenMP option.
1921
1922                                if ( $compopt == $COMPOPTS[2] ) then
1923                                        if ( $ZAP_OPENMP == TRUE || $ZAP_OPENMP_FOR_THIS_CORE == TRUE ) then
1924                                                goto BYPASS_COMP_LOOP_REAL
1925                                        endif
1926                                endif
1927
1928                                if ( $compopt == $COMPOPTS[1] ) then
1929                                        if ( $ZAP_SERIAL == TRUE || $ZAP_SERIAL_FOR_THIS_CORE == TRUE ) then
1930                                                goto BYPASS_COMP_LOOP_REAL
1931                                        endif
1932                                endif
1933
1934                                pushd test/$core
1935                               
1936                                #
1937                                #       Create the correct namelist.input file for real data cases.
1938                                #
1939
1940                                if ( $CHEM != TRUE ) then
1941                                        cp ${CUR_DIR}/phys_real_${phys_option} phys_opt
1942                                        if ( ( $RSL_LITE == TRUE ) && ( $NESTED != TRUE ) ) then
1943                                                cp ${CUR_DIR}/dyn_real_${phys_option} dyn_opt
1944                                        else
1945                                                cp ${CUR_DIR}/dyn_real_SAFE dyn_opt
1946                                        endif
1947                                        cp ${CUR_DIR}/time_real_${phys_option} time_opt
1948                                        cp ${CUR_DIR}/dom_real dom_real
1949                                        if ( -e fdda_opt ) rm fdda_opt
1950                                        cat " grid_fdda=0" > fdda_opt
1951                                        if ( -e fdda_time ) rm fdda_time
1952                                        cat " auxinput4_interval=0" > fdda_time
1953
1954                                        if ( $FDDA == TRUE ) then
1955                                                cp ${CUR_DIR}/fdda_real_${phys_option} fdda_opt
1956                                                cp ${CUR_DIR}/fdda_real_time_${phys_option} fdda_time
1957                                        endif
1958       
1959                                        set time_step = `awk ' /^ time_step /{ print $3 } ' namelist.input.$dataset | cut -d, -f1`
1960       
1961                                        #       Wanna do more/less time steps on the real cases?  Easy. Those last two numbers
1962                                        #       in the eqns are all you need.  Their product must be 60.  So, instead of 3 and 20,
1963                                        #       (3 coarse grid timesteps), you could use 20 and 3 (20 coarse grid time steps).
1964       
1965                                        if      ( $NESTED == TRUE ) then
1966                                                @ run_seconds = $time_step * 3
1967                                                @ history_interval = $time_step / 20
1968                                        else if ( $NESTED != TRUE ) then
1969                                                @ run_seconds = $time_step * 10
1970                                                @ history_interval = $time_step / 6
1971                                        endif
1972                                        rm ed_in namelist.input.temp
1973                                        cat >! ed_in << EOF
1974g/run_seconds/s/[0-9]/$run_seconds
1975g/history_interval/s/[0-9][0-9][0-9]/$history_interval
1976w namelist.input.temp
1977q
1978EOF
1979                                        ed namelist.input.$dataset < ed_in
1980       
1981                                        cp ${CUR_DIR}/io_format io_format
1982                                        sed -e '/^ mp_physics/,/ensdim/d' -e '/^ &physics/r ./phys_opt' \
1983                                            -e '/^ pd_moist/,/pd_tke/d' -e '/^ v_sca_adv_order/r ./dyn_opt' \
1984                                            -e '/^ auxinput1_inname/d' -e '/^ debug_level/r ./time_opt' \
1985                                            -e '/^ time_step /,/^ smooth_option/d' -e '/^ &domains/r ./dom_real' \
1986                                            -e '/^ io_form_history /,/^ io_form_boundary/d' -e '/^ restart_interval/r ./io_format' \
1987                                            -e 's/ frames_per_outfile *= [0-9][0-9]*/ frames_per_outfile = 200/g' \
1988                                            -e 's/ run_days *= [0-9][0-9]*/ run_days = 0/g' \
1989                                            -e 's/ run_hours *= [0-9][0-9]*/ run_hours = 0/g' \
1990                                            -e 's/ run_minutes *= [0-9][0-9]*/ run_minutes = 0/g' \
1991                                            -e '/^ &fdda/r fdda_opt' \
1992                                            -e '/^ debug_level/r fdda_time' \
1993                                        namelist.input.temp >! namelist.input
1994                               
1995                                #       The chem run has its own namelist, due to special input files (io_form not tested for chem)
1996
1997                                else if ( $CHEM == TRUE ) then
1998                                        cp namelist.input.chem_test_${phys_option} namelist.input
1999                                        if ( -e wrf_real_input_em.d01.${filetag} ) then
2000                                                \rm wrf_real_input_em.d01.*
2001                                        endif
2002                                        if ( ${phys_option} <= 3 ) then
2003                                                ln -s 00z/wrf_real* .
2004                                        else
2005                                                ln -s 12z/wrf_real* .
2006                                        endif
2007
2008                                endif
2009
2010                                # WRF output quilt servers are only tested for MPI configuration. 
2011                                # Currently, only one WRF output quilt server is used. 
2012
2013                                if ( $QUILT == TRUE ) then
2014                                        if ( $compopt == $COMPOPTS[3] ) then
2015                                                #       For now, test only one group of one output quilt servers. 
2016                                                sed -e 's/ nio_tasks_per_group *= *[0-9][0-9]*/ nio_tasks_per_group = 1/g' \
2017                                                    -e 's/ nio_groups *= *[0-9][0-9]*/ nio_groups = 1/g' \
2018                                                namelist.input >! namelist.input.temp
2019                                                mv -f namelist.input.temp namelist.input
2020                                                echo "Building namelist.input.$core.${phys_option}.$compopt with one I/O quilt server enabled."
2021                                                echo "NOTE  one I/O quilt server enabled for $core physics $phys_option parallel $compopt..." >>! ${DEF_DIR}/wrftest.output
2022                                        endif
2023                                endif
2024
2025                                /bin/cp namelist.input $TMPDIR/namelist.input.$core.${phys_option}.$compopt
2026#DAVE###################################################
2027echo built namelist $TMPDIR/namelist.input.$core.${phys_option}.$compopt
2028cat  $TMPDIR/namelist.input.$core.${phys_option}.$compopt
2029echo need history interval to be 30
2030banner 13
2031#set ans = "$<"
2032#DAVE###################################################
2033#DAVE###################################################
2034echo skipped link of data files, we push them elsewhere
2035ls -ls wrf_real*
2036banner 14
2037#set ans = "$<"
2038#DAVE###################################################
2039
2040                                #       If this is the serial code, generate the IC and BC.  The real.exe program is not
2041                                #       parallelized, so the data is generated and saved for the rest of the parallel tests.
2042                                #       This data is necessarily updated for each of the physics tests.
2043
2044                                if ( $compopt == $COMPOPTS[1] ) then
2045
2046                                        #       Zap any old input data laying around.
2047
2048                                        rm wrfinput_d01 >& /dev/null
2049                                        rm wrfbdy_d01   >& /dev/null
2050
2051                                        if ( $NESTED == TRUE ) then
2052                                                setenv OMP_NUM_THREADS 1
2053                                                if ( `uname` == AIX ) then
2054                                                        setenv XLSMPOPTS "parthds=1"
2055                                                endif
2056                                                $SERIALRUNCOMMAND ../../main/real_${core}.exe.1 >! print.out.real_${core}_Phys=${phys_option}_Parallel=${compopt}
2057                                        else if ( $NESTED != TRUE ) then
2058                                                ../../main/real_${core}.exe.1 >! print.out.real_${core}_Phys=${phys_option}_Parallel=${compopt}
2059                                        endif
2060#DAVE###################################################
2061echo finished real
2062banner 15
2063#set ans = "$<"
2064#DAVE###################################################
2065
2066                                        grep "SUCCESS COMPLETE" print.out.real_${core}_Phys=${phys_option}_Parallel=${compopt} >& /dev/null
2067                                        set success = $status
2068
2069                                        #       Did making the IC BC files work?
2070
2071                                        if ( ( -e wrfinput_d01 ) && ( -e wrfbdy_d01 ) && ( $success == 0 ) ) then
2072                                                echo "SUMMARY generate IC/BC for $core physics $phys_option parallel $compopt $esmf_lib_str PASS" >>! ${DEF_DIR}/wrftest.output
2073                                                echo "-------------------------------------------------------------" >> ${DEF_DIR}/wrftest.output
2074                                        else
2075                                                echo "SUMMARY generate IC/BC for $core physics $phys_option parallel $compopt $esmf_lib_str FAIL" >>! ${DEF_DIR}/wrftest.output
2076                                                $MAIL -s "WRF FAIL making IC/BC $ARCH[1] " $FAIL_MAIL < ${DEF_DIR}/wrftest.output
2077                                                exit ( 4 )
2078                                        endif
2079#DAVE###################################################
2080echo IC BC must be OK
2081ls -ls wrfi* wrfb*
2082if ( $IO_FORM_NAME[$IO_FORM] == io_netcdf ) then
2083        ncdump -v Times wrfb* | tail -20
2084endif
2085banner 16
2086#set ans = "$<"
2087#DAVE###################################################
2088                                endif
2089               
2090                                #       Run the forecast for this core, physics package and parallel option
2091
2092                                rm $TMPDIR/wrfout_d01_${filetag}.${core}.${phys_option}.$compopt >& /dev/null
2093                               
2094                                #       The chem run has its own set of namelists, due to special input files.
2095
2096                                if ( $CHEM == TRUE ) then
2097
2098                                        # WRF output quilt servers are only tested for MPI configuration. 
2099                                        # Currently, only one WRF output quilt server is used. 
2100       
2101                                        if ( $QUILT == TRUE ) then
2102                                                if ( $compopt == $COMPOPTS[3] ) then
2103                                                        #       For now, test only one group of one output quilt servers. 
2104                                                        sed -e 's/ nio_tasks_per_group *= *[0-9][0-9]*/ nio_tasks_per_group = 1/g' \
2105                                                            -e 's/ nio_groups *= *[0-9][0-9]*/ nio_groups = 1/g' \
2106                                                        namelist.input >! namelist.input.temp
2107                                                        mv -f namelist.input.temp namelist.input
2108                                                        echo "Building namelist.input.$core.${phys_option}.$compopt with one I/O quilt server enabled."
2109                                                        echo "NOTE  one I/O quilt server enabled for $core physics $phys_option parallel $compopt..." >>! ${DEF_DIR}/wrftest.output
2110                                                endif
2111                                        endif
2112                                endif
2113
2114                                if      ( $compopt == $COMPOPTS[1] ) then
2115                                        setenv OMP_NUM_THREADS 1
2116                                        if ( `uname` == AIX ) then
2117                                                setenv XLSMPOPTS "parthds=1"
2118                                        endif
2119                                        if ( $NESTED == TRUE ) then
2120                                                $SERIALRUNCOMMAND ../../main/wrf_${core}.exe.$compopt >! print.out.wrf_${core}_Phys=${phys_option}_Parallel=${compopt}.exe
2121                                        else if ( $NESTED != TRUE ) then
2122                                                ../../main/wrf_${core}.exe.$compopt >! print.out.wrf_${core}_Phys=${phys_option}_Parallel=${compopt}.exe
2123                                        endif
2124                                else if ( $compopt == $COMPOPTS[2] ) then
2125                                        setenv OMP_NUM_THREADS $OPENMP
2126                                        if ( `uname` == AIX ) then
2127                                                setenv XLSMPOPTS "parthds=${OPENMP}"
2128                                        endif
2129                                        if ( $NESTED == TRUE ) then
2130                                                $OMPRUNCOMMAND ../../main/wrf_${core}.exe.$compopt >! print.out.wrf_${core}_Phys=${phys_option}_Parallel=${compopt}.exe
2131                                        else if ( $NESTED != TRUE ) then
2132                                                ../../main/wrf_${core}.exe.$compopt >! print.out.wrf_${core}_Phys=${phys_option}_Parallel=${compopt}.exe
2133                                        endif
2134                                else if ( $compopt == $COMPOPTS[3] ) then
2135                                        setenv OMP_NUM_THREADS 1
2136                                        if ( `uname` == AIX ) then
2137                                                setenv XLSMPOPTS "parthds=1"
2138                                        endif
2139                                        $MPIRUNCOMMAND ../../main/wrf_${core}.exe.$compopt $MPIRUNCOMMANDPOST
2140                                        mv rsl.error.0000 print.out.wrf_${core}_Phys=${phys_option}_Parallel=${compopt}.exe
2141                                endif
2142#DAVE###################################################
2143echo ran wrf fcst compopt = $compopt
2144banner 17
2145#set ans = "$<"
2146#DAVE###################################################
2147
2148                                grep "SUCCESS COMPLETE" print.out.wrf_${core}_Phys=${phys_option}_Parallel=${compopt}.exe
2149                                set success = $status
2150
2151                                #       Did making the forecast work, by that, we mean "is there an output file created", and "are there 2 times periods".
2152
2153                                if ( ( -e wrfout_d01_${filetag} ) && ( $success == 0 ) ) then
2154                                        if      ( $IO_FORM_NAME[$IO_FORM] == io_netcdf ) then
2155                                                ncdump -h wrfout_d01_${filetag} | grep Time | grep UNLIMITED | grep currently | grep -q 2
2156                                                set ok = $status
2157                                                set found_nans = 1
2158                                                if      ( `uname` == AIX   ) then
2159                                                        ncdump wrfout_d01_${filetag} | grep NaN >& /dev/null
2160                                                        set found_nans = $status
2161                                                else if ( `uname` == OSF1  ) then
2162                                                        ncdump wrfout_d01_${filetag} | grep nan >& /dev/null
2163                                        #               set found_nans = $status
2164                                                else if ( `uname` == Linux ) then
2165                                                        ncdump wrfout_d01_${filetag} | grep nan >& /dev/null
2166                                                        set found_nans = $status
2167                                                endif
2168                                                if ( $found_nans == 0 ) then
2169                                                        echo found nans
2170                                                        set ok = 1
2171                                                endif
2172
2173                                        else if ( $IO_FORM_NAME[$IO_FORM] == io_grib1  ) then
2174#                                               set joe_times = `../../external/io_grib1/wgrib -s -4yr wrfout_d01_${filetag} |  grep -v ":anl:" | wc -l`
2175#                                               if ( $joe_times >= 100 ) then
2176                                                set joe_times = `../../external/io_grib1/wgrib -s -4yr wrfout_d01_${filetag} |  grep "UGRD:10 m" | wc -l`
2177                                                if ( $joe_times == 2 ) then
2178                                                        set ok = 0
2179                                                else
2180                                                        set ok = 1
2181                                                endif
2182                                        endif
2183                                        if ( $ok == 0 ) then
2184                                                echo "SUMMARY generate FCST  for $core physics $phys_option parallel $compopt $esmf_lib_str PASS" >>! ${DEF_DIR}/wrftest.output
2185                                                echo "-------------------------------------------------------------" >> ${DEF_DIR}/wrftest.output
2186                                        else
2187                                                echo "SUMMARY generate FCST  for $core physics $phys_option parallel $compopt $esmf_lib_str FAIL" >>! ${DEF_DIR}/wrftest.output
2188                                                $MAIL -s "WRF FAIL FCST $ARCH[1] " $FAIL_MAIL < ${DEF_DIR}/wrftest.output
2189                                                exit ( 5 )
2190                                        endif
2191                                else
2192                                        echo "SUMMARY generate FCST  for $core physics $phys_option parallel $compopt $esmf_lib_str FAIL" >>! ${DEF_DIR}/wrftest.output
2193                                        $MAIL -s "WRF FAIL FCST $ARCH[1] " $FAIL_MAIL < ${DEF_DIR}/wrftest.output
2194                                        exit ( 6 )
2195                                endif
2196#DAVE###################################################
2197echo success or failure of fcst
2198if ( $IO_FORM_NAME[$IO_FORM] == io_netcdf ) then
2199        ncdump -v Times wrfout_d01_${filetag} | tail -20
2200endif
2201banner 18
2202#set ans = "$<"
2203#DAVE###################################################
2204
2205                                #       We have to save this output file for our biggy comparison after all of the
2206                                #       parallel options have been considered.
2207
2208                                mv wrfout_d01_${filetag} $TMPDIR/wrfout_d01_${filetag}.${core}.${phys_option}.$compopt
2209
2210                                #       To save space, we move the executables after we are finished with them.
2211
2212                                if ( $phys_option == $PHYSOPTS[${#PHYSOPTS}] ) then
2213                                        mv ../../main/wrf_${core}.exe.$compopt $TMPDIR/wrf_${core}.exe.$compopt
2214                                endif
2215
2216                                popd
2217
2218                                BYPASS_COMP_LOOP_REAL:
2219
2220                        end
2221
2222                else if ( $core == nmm_real )  then
2223#DAVE###################################################
2224echo doing nmm pre
2225banner 19
2226#set ans = "$<"
2227#DAVE###################################################
2228
2229                        set compopt = $COMPOPTS[3]   # ! parallel only
2230                        set filetag = 2005-01-23_00:00:00
2231                        set phys_option=1
2232                        pushd test/$core
2233
2234#DAVE###################################################
2235echo did rms
2236echo $filetag $phys_option
2237banner 19a
2238#set ans = "$<"
2239#DAVE###################################################
2240
2241                        #       Build NMM namelist
2242
2243                        cp ${CUR_DIR}/io_format io_format
2244                        sed -e '/^ io_form_history /,/^ io_form_boundary/d' -e '/^ restart_interval/r ./io_format' \
2245                            namelist.input.regtest >! namelist.input.temp
2246       
2247                        #       A fairly short forecast, 10 time steps
2248
2249                        sed -e 's/^ run_days *= *[0-9]*/ run_days = 0 /' \
2250                            -e 's/^ run_seconds *= *[0-9]*/ run_seconds = 900 /' \
2251                            -e 's/^ history_interval *= *[0-9][0-9]*/ history_interval = 15 /' \
2252                            -e 's/^ frames_per_outfile *= [0-9]*/ frames_per_outfile = 200/g' \
2253                            namelist.input.temp >! namelist.input
2254
2255#DAVE###################################################
2256echo did cp of namelist
2257ls -ls namelist.input
2258cat namelist.input
2259banner 19b
2260#set ans = "$<"
2261#DAVE###################################################
2262
2263                        #       Generate IC/BC for NMM run
2264
2265                        set RUNCOMMAND = `echo $MPIRUNCOMMAND | sed "s/$Num_Procs/1/"`
2266
2267                        #       Zap any old input data laying around.
2268
2269                        rm wrfinput_d01 >& /dev/null
2270                        rm wrfbdy_d01   >& /dev/null
2271
2272                        $RUNCOMMAND ../../main/real_${core}.exe.$COMPOPTS[3] >! print.out.real_${core}_Phys=${phys_option}_Parallel=$COMPOPTS[3]
2273#DAVE###################################################
2274echo finished real
2275banner 19c
2276#set ans = "$<"
2277#DAVE###################################################
2278
2279                        mv rsl.out.0000 print.out.real_${core}_Phys=${phys_option}_Parallel=${compopt}
2280                        grep "SUCCESS COMPLETE" print.out.real_${core}_Phys=${phys_option}_Parallel=${compopt} >& /dev/null
2281                        set success = $status
2282
2283                        #       Did making the IC BC files work?
2284
2285                        if ( ( -e wrfinput_d01 ) && ( -e wrfbdy_d01 ) && ( $success == 0 ) ) then
2286                                echo "SUMMARY generate IC/BC for $core physics $phys_option parallel $compopt $esmf_lib_str PASS" >>! ${DEF_DIR}/wrftest.output
2287                                echo "-------------------------------------------------------------" >> ${DEF_DIR}/wrftest.output
2288                        else
2289                                echo "SUMMARY generate IC/BC for $core physics $phys_option parallel $compopt $esmf_lib_str FAIL" >>! ${DEF_DIR}/wrftest.output
2290                                $MAIL -s "WRF FAIL making IC/BC $ARCH[1] " $FAIL_MAIL < ${DEF_DIR}/wrftest.output
2291                                exit ( 4 )
2292                        endif
2293#DAVE###################################################
2294echo IC BC must be OK
2295ls -lsL wrfinput* wrfb*
2296if ( $IO_FORM_NAME[$IO_FORM] == io_netcdf ) then
2297        ncdump -v Times wrfb* | tail -20
2298endif
2299banner 20
2300#set ans = "$<"
2301#DAVE###################################################
2302
2303                        #       Run on 1 and then on Num_Procs processors
2304
2305                        foreach n ( 1 $Num_Procs )
2306#DAVE###################################################
2307echo running nmm on $n procs
2308banner 21
2309#set ans = "$<"
2310#DAVE###################################################
2311
2312                                if      ( ( $n != 1 ) && ( $QUILT != TRUE ) ) then
2313                                        @ nmm_proc = $Num_Procs
2314                                        cat >! nproc_xy << EOF
2315 nproc_x = $nmm_proc
2316 nproc_y = 1
2317EOF
2318                                        sed -e '/^ numtiles/r nproc_xy' namelist.input >! file.foo
2319                                        mv file.foo namelist.input
2320                                else if ( ( $n != 1 ) && ( $QUILT == TRUE ) ) then
2321                                        @ nmm_proc = $Num_Procs - 1
2322                                        cat >! nproc_xy << EOF
2323 nproc_x = $nmm_proc
2324 nproc_y = 1
2325EOF
2326                                        sed -e '/^ numtiles/r nproc_xy' namelist.input >! file.foo
2327                                        mv file.foo namelist.input
2328                                endif
2329
2330                                if ( `uname` == AIX ) then
2331                                        set RUNCOMMAND = $MPIRUNCOMMAND
2332                                else
2333                                        set RUNCOMMAND = `echo $MPIRUNCOMMAND | sed "s/$Num_Procs/$n/"`
2334                                endif
2335
2336                                #       WRF output quilt servers are only tested for MPI configuration. 
2337                                #       Currently, only one WRF output quilt server is used. 
2338                       
2339                                if ( ( $QUILT == TRUE ) && ( $n == $Num_Procs ) ) then
2340                                        if ( $compopt == $COMPOPTS[3] ) then
2341                                                #       For now, test only one group of one output quilt servers. 
2342                                                sed -e 's/ nio_tasks_per_group *= *[0-9]*/ nio_tasks_per_group = 1/g' \
2343                                                    -e 's/ nio_groups *= *[0-9]*/ nio_groups = 1/g' \
2344                                                namelist.input >! namelist.input.temp
2345                                                mv -f namelist.input.temp namelist.input
2346                                                echo "Building namelist.input.$core.${phys_option}.$compopt with one I/O quilt server enabled."
2347                                                echo "NOTE  one I/O quilt server enabled for $core physics $phys_option parallel $compopt..." >>! ${DEF_DIR}/wrftest.output
2348                                        endif
2349                                endif
2350
2351                                #       NMM can fail on spurious fp exceptions that don't affect soln. Retry if necessary.
2352
2353                                set tries=0
2354                                while ( $tries < 2 )
2355#DAVE###################################################
2356echo try attempt $tries allowed to be less than 2
2357banner 22
2358#set ans = "$<"
2359#DAVE###################################################
2360                                        @ tries = $tries + 1
2361                                        $RUNCOMMAND ../../main/wrf_${core}.exe.$compopt $MPIRUNCOMMANDPOST
2362                                        mv rsl.error.0000 print.out.wrf_${core}_Phys=${phys_option}_Parallel=${compopt}.exe_${n}p
2363                                        grep "SUCCESS COMPLETE" print.out.wrf_${core}_Phys=${phys_option}_Parallel=${compopt}.exe_${n}p
2364                                        set success = $status
2365                                        set ok = $status
2366                                        if ( ( -e wrfout_d01_${filetag} ) && ( $success == 0 ) ) then
2367                                                if      ( $IO_FORM_NAME[$IO_FORM] == io_netcdf ) then
2368                                                        ncdump -h wrfout_d01_${filetag} | grep Time | grep UNLIMITED | grep currently | grep -q 2
2369                                                        set ok = $status
2370                                                        set found_nans = 1
2371                                                        if      ( `uname` == AIX   ) then
2372                                                                ncdump wrfout_d01_${filetag} | grep NaN >& /dev/null
2373                                                                set found_nans = $status
2374                                                        else if ( `uname` == OSF1  ) then
2375                                                                ncdump wrfout_d01_${filetag} | grep nan >& /dev/null
2376                                                #               set found_nans = $status
2377                                                        else if ( `uname` == Linux ) then
2378                                                                ncdump wrfout_d01_${filetag} | grep nan >& /dev/null
2379                                                                set found_nans = $status
2380                                                        endif
2381                                                        if ( $found_nans == 0 ) then
2382                                                                echo found nans
2383                                                                set ok = 1
2384                                                        endif
2385
2386                                                else if ( $IO_FORM_NAME[$IO_FORM] == io_grib1  ) then
2387                                                        ../../external/io_grib1/wgrib -s -4yr wrfout_d01_${filetag} |  grep "UGRD:10 m above gnd:3600 sec fcst"
2388                                                        set ok = $status
2389                                                endif
2390                                                if ( $ok == 0 ) then
2391                                                        echo "SUMMARY generate FCST  for $core physics $phys_option parallel $compopt $esmf_lib_str PASS" >>! ${DEF_DIR}/wrftest.output
2392                                                        echo "-------------------------------------------------------------" >> ${DEF_DIR}/wrftest.output
2393                                                        set tries=# success, bail from loop
2394                                                else
2395                                                        echo "SUMMARY generate FCST  for $core physics $phys_option parallel $compopt $esmf_lib_str FAIL" >>! ${DEF_DIR}/wrftest.output
2396                                                        $MAIL -s "WRF FAIL FCST $ARCH[1] " $FAIL_MAIL < ${DEF_DIR}/wrftest.output
2397                                                        if ( $tries == 2 ) exit ( 5 )
2398                                                endif
2399                                        else
2400                                                echo "SUMMARY generate FCST  for $core physics $phys_option parallel $compopt $esmf_lib_str FAIL" >>! ${DEF_DIR}/wrftest.output
2401                                                $MAIL -s "WRF FAIL FCST $ARCH[1] " $FAIL_MAIL < ${DEF_DIR}/wrftest.output
2402                                                if ( $tries == 2 ) exit ( 6 )
2403                                        endif
2404                                        mv wrfout_d01_${filetag} $TMPDIR/wrfout_d01_${filetag}.${core}.${phys_option}.${compopt}_${n}p
2405#DAVE###################################################
2406echo did nmm fcst
2407ls -ls $TMPDIR/wrfout_d01_${filetag}.${core}.${phys_option}.${compopt}_${n}p
2408banner 23
2409#set ans = "$<"
2410#DAVE###################################################
2411                                end
2412                        end
2413
2414                        popd
2415
2416                else
2417
2418#DAVE###################################################
2419echo doing ideal runs
2420banner 24
2421#set ans = "$<"
2422#DAVE###################################################
2423                        #       The ideal cases have different physics tests than the real cases.  If this is
2424                        #       more that the total number of ideal physics experiments that we were led to
2425                        #       believe would exist, jump to the end of the physics loop.
2426
2427                        if ( $phys_option > $Max_Ideal_Physics_Options ) then
2428                                goto BOTTOM_OF_PHYSICS_LOOP
2429                        endif
2430
2431                        set filetag=$filetag_ideal
2432
2433                        foreach compopt ( $COMPOPTS )
2434#DAVE###################################################
2435echo doing compopt = $compopt
2436echo filetag = $filetag
2437banner 25
2438#set ans = "$<"
2439#DAVE###################################################
2440
2441
2442                                #       We sometimes are interested in bypassing the OpenMP option.
2443
2444                                if ( $compopt == $COMPOPTS[2] ) then
2445                                        if ( $ZAP_OPENMP == TRUE ) then
2446                                                goto BYPASS_COMP_LOOP_IDEAL
2447                                        endif
2448                                endif
2449
2450                                pushd test/$core
2451
2452                                if ( ! -e namelist.input.template ) cp namelist.input namelist.input.template
2453
2454                                #       Create the correct namelist.input file.
2455
2456                                cp ${CUR_DIR}/dom_ideal dom_ideal
2457                                if      ( $core == em_quarter_ss )  then
2458                                        cp ${CUR_DIR}/phys_quarter_ss_${phys_option}a phys_tke
2459                                        cp ${CUR_DIR}/phys_quarter_ss_${phys_option}b phys_mp
2460                                        cp ${CUR_DIR}/phys_quarter_ss_${phys_option}c phys_nh
2461                                        cp ${CUR_DIR}/phys_quarter_ss_${phys_option}d phys_bc
2462                                else if ( $core == em_b_wave     ) then
2463                                        cp ${CUR_DIR}/phys_b_wave_${phys_option}a     phys_tke
2464                                        cp ${CUR_DIR}/phys_b_wave_${phys_option}b     phys_mp
2465                                        cp ${CUR_DIR}/phys_b_wave_${phys_option}c     phys_nh
2466                                endif
2467
2468                                cp ${CUR_DIR}/io_format io_format
2469                                if      ( $NESTED == TRUE ) then
2470                                        if      ( $core == em_quarter_ss )  then
2471                                                sed -e 's/ run_days *= *[0-9][0-9]*/ run_days = 00/g'                                           \
2472                                                    -e 's/ run_minutes *= *[0-9][0-9]*/ run_minutes = 1/g'                                      \
2473                                                    -e 's/ run_seconds *= *[0-9][0-9]*/ run_seconds = 0/g'                                      \
2474                                                    -e 's/ history_interval *= [0-9][0-9]*/ history_interval = 1/g'                             \
2475                                                    -e 's/ frames_per_outfile *= [0-9][0-9]*/ frames_per_outfile = 100/g'                       \
2476                                                    -e '/^ diff_opt/d' -e '/^ km_opt/d' -e '/^ damp_opt/d' -e '/^ rk_ord/r ./phys_tke'          \
2477                                                    -e '/^ io_form_history /,/^ io_form_boundary/d' -e '/^ restart_interval/r ./io_format'      \
2478                                                    -e '/^ mp_physics/d' -e '/^ &physics/r ./phys_mp'                                           \
2479                                                    -e '/^ non_hydrostatic/d' -e '/^ epssm/r ./phys_nh'                                         \
2480                                                    -e '/^ periodic_x/d' -e '/^ open_xs/d' -e '/^ open_xe/d'                                    \
2481                                                    -e '/^ periodic_y/d' -e '/^ open_ys/d' -e '/^ open_ye/d'                                    \
2482                                                    -e '/^ &bdy_control/r ./phys_bc'                                                            \
2483                                                    -e '/^ max_dom/d' -e '/^ time_step_fract_den/r ./dom_ideal'                                 \
2484                                                    ./namelist.input.template >! namelist.input
2485                                        else if ( $core == em_b_wave     ) then
2486                                                sed -e 's/ run_days *= *[0-9][0-9]*/ run_days = 00/g'                                           \
2487                                                    -e 's/ run_minutes *= *[0-9][0-9]*/ run_minutes = 20/g'                                     \
2488                                                    -e 's/ history_interval *= [0-9][0-9]*/ history_interval = 20/g'                            \
2489                                                    -e 's/ frames_per_outfile *= [0-9][0-9]*/ frames_per_outfile = 100/g'                       \
2490                                                    -e '/^ diff_opt/d' -e '/^ km_opt/d' -e '/^ damp_opt/d' -e '/^ rk_ord/r ./phys_tke'          \
2491                                                    -e '/^ io_form_history /,/^ io_form_boundary/d' -e '/^ restart_interval/r ./io_format'      \
2492                                                    -e '/^ mp_physics/d' -e '/^ &physics/r ./phys_mp'                                           \
2493                                                    -e '/^ non_hydrostatic/d' -e '/^ epssm/r ./phys_nh'                                         \
2494                                                    -e '/^ max_dom/d' -e '/^ time_step_fract_den/r ./dom_ideal'                                 \
2495                                                ./namelist.input.template >! namelist.input
2496                                        endif
2497                                else if ( $NESTED != TRUE ) then
2498                                        if      ( $core == em_quarter_ss )  then
2499                                                sed -e 's/ run_days *= *[0-9][0-9]*/ run_days = 00/g'                                           \
2500                                                    -e 's/ run_minutes *= *[0-9][0-9]*/ run_minutes = 2/g'                                      \
2501                                                    -e 's/ history_interval *= [0-9][0-9]*/ history_interval = 2/g'                             \
2502                                                    -e 's/ frames_per_outfile *= [0-9][0-9]*/ frames_per_outfile = 100/g'                       \
2503                                                    -e '/^ diff_opt/d' -e '/^ km_opt/d' -e '/^ damp_opt/d' -e '/^ rk_ord/r ./phys_tke'          \
2504                                                    -e '/^ io_form_history /,/^ io_form_boundary/d' -e '/^ restart_interval/r ./io_format'      \
2505                                                    -e '/^ mp_physics/d' -e '/^ &physics/r ./phys_mp'                                           \
2506                                                    -e '/^ non_hydrostatic/d' -e '/^ epssm/r ./phys_nh'                                         \
2507                                                    -e '/^ periodic_x/d' -e '/^ open_xs/d' -e '/^ open_xe/d'                                    \
2508                                                    -e '/^ periodic_y/d' -e '/^ open_ys/d' -e '/^ open_ye/d'                                    \
2509                                                    -e '/^ &bdy_control/r ./phys_bc'                                                            \
2510                                                    -e '/^ max_dom/d' -e '/^ time_step_fract_den/r ./dom_ideal'                                 \
2511                                                    ./namelist.input.template >! namelist.input
2512                                        else if ( $core == em_b_wave     ) then
2513                                                sed -e 's/ run_days *= *[0-9][0-9]*/ run_days = 00/g'                                           \
2514                                                    -e 's/ run_minutes *= *[0-9][0-9]*/ run_minutes = 100/g'                                    \
2515                                                    -e 's/ history_interval *= [0-9][0-9]*/ history_interval = 100/g'                           \
2516                                                    -e 's/ frames_per_outfile *= [0-9][0-9]*/ frames_per_outfile = 100/g'                       \
2517                                                    -e '/^ diff_opt/d' -e '/^ km_opt/d' -e '/^ damp_opt/d' -e '/^ rk_ord/r ./phys_tke'          \
2518                                                    -e '/^ io_form_history /,/^ io_form_boundary/d' -e '/^ restart_interval/r ./io_format'      \
2519                                                    -e '/^ mp_physics/d' -e '/^ &physics/r ./phys_mp'                                           \
2520                                                    -e '/^ non_hydrostatic/d' -e '/^ epssm/r ./phys_nh'                                         \
2521                                                    -e '/^ max_dom/d' -e '/^ time_step_fract_den/r ./dom_ideal'                                 \
2522                                                ./namelist.input.template >! namelist.input
2523                                        endif
2524                                endif
2525
2526                                # WRF output quilt servers are only tested for MPI configuration. 
2527                                # Currently, only one WRF output quilt server is used. 
2528
2529                                if ( $QUILT == TRUE ) then
2530                                        if ( $compopt == $COMPOPTS[3] ) then
2531                                                #       For now, test only one group of one output quilt servers. 
2532                                                sed -e 's/ nio_tasks_per_group *= *[0-9][0-9]*/ nio_tasks_per_group = 1/g' \
2533                                                    -e 's/ nio_groups *= *[0-9][0-9]*/ nio_groups = 1/g' \
2534                                                namelist.input >! namelist.input.temp
2535                                                mv -f namelist.input.temp namelist.input
2536                                                echo "Building namelist.input.$core.${phys_option}.$compopt with one I/O quilt server enabled."
2537                                                echo "NOTE  one I/O quilt server enabled for $core physics $phys_option parallel $compopt..." >>! ${DEF_DIR}/wrftest.output
2538                                        endif
2539                                endif
2540#DAVE###################################################
2541echo built namelist
2542ls -ls namelist.input
2543banner 26
2544#set ans = "$<"
2545#DAVE###################################################
2546
2547                                #       If this is the serial code, generate the IC and BC.  The ideal.exe program is not
2548                                #       parallelized, so the data is generated and saved for the rest of the parallel tests.
2549
2550                                if ( $compopt == $COMPOPTS[1] ) then
2551
2552                                        #       Zap any old input data laying around.
2553
2554                                        rm wrfinput_d01 >& /dev/null
2555                                        rm wrfbdy_d01   >& /dev/null
2556
2557                                        ../../main/ideal_${core}.exe.1 >! print.out.ideal_${core}_Parallel=${compopt}
2558#DAVE###################################################
2559echo ran ideal
2560ls -ls wrfiput*
2561banner 27
2562#set ans = "$<"
2563#DAVE###################################################
2564
2565                                        grep "SUCCESS COMPLETE" print.out.ideal_${core}_Parallel=${compopt} >& /dev/null
2566                                        set success = $status
2567
2568                                        #       Did making the IC BC files work?
2569
2570                                        if ( ( -e wrfinput_d01 ) && ( $success == 0 ) ) then
2571                                                echo "SUMMARY generate IC/BC for $core physics $phys_option parallel $compopt $esmf_lib_str PASS" >>! ${DEF_DIR}/wrftest.output
2572                                                echo "-------------------------------------------------------------" >> ${DEF_DIR}/wrftest.output
2573                                        else
2574                                                echo "SUMMARY generate IC/BC for $core physics $phys_option parallel $compopt $esmf_lib_str FAIL" >>! ${DEF_DIR}/wrftest.output
2575                                                $MAIL -s "WRF FAIL making IC/BC $ARCH[1] " $FAIL_MAIL < ${DEF_DIR}/wrftest.output
2576                                                exit ( 7 )
2577                                        endif
2578                                endif
2579               
2580                                #       Run the forecast for this core and parallel option
2581
2582                                rm $TMPDIR/wrfout_d01_${filetag}.${core}.${phys_option}.$compopt >& /dev/null
2583
2584                                if      ( $compopt == $COMPOPTS[1] ) then
2585                                        setenv OMP_NUM_THREADS 1
2586                                        if ( `uname` == AIX ) then
2587                                                setenv XLSMPOPTS "parthds=1"
2588                                        endif
2589                                        if ( $NESTED == TRUE ) then
2590                                                $SERIALRUNCOMMAND ../../main/wrf_${core}.exe.$compopt >! print.out.wrf_${core}_Parallel=${compopt}.exe
2591                                        else if ( $NESTED != TRUE ) then
2592                                                ../../main/wrf_${core}.exe.$compopt >! print.out.wrf_${core}_Parallel=${compopt}.exe
2593                                        endif
2594                                else if ( $compopt == $COMPOPTS[2] ) then
2595                                        setenv OMP_NUM_THREADS $OPENMP
2596                                        if ( `uname` == AIX ) then
2597                                                setenv XLSMPOPTS "parthds=${OPENMP}"
2598                                        endif
2599                                        if ( $NESTED == TRUE ) then
2600                                                $OMPRUNCOMMAND ../../main/wrf_${core}.exe.$compopt >! print.out.wrf_${core}_Parallel=${compopt}.exe
2601                                        else if ( $NESTED != TRUE ) then
2602                                                ../../main/wrf_${core}.exe.$compopt >! print.out.wrf_${core}_Parallel=${compopt}.exe
2603                                        endif
2604                                else if ( $compopt == $COMPOPTS[3] ) then
2605                                        setenv OMP_NUM_THREADS 1
2606                                        if ( `uname` == AIX ) then
2607                                                setenv XLSMPOPTS "parthds=1"
2608                                        endif
2609                                        $MPIRUNCOMMAND ../../main/wrf_${core}.exe.$compopt $MPIRUNCOMMANDPOST
2610                                        mv rsl.error.0000 print.out.wrf_${core}_Parallel=${compopt}.exe
2611                                endif
2612#DAVE###################################################
2613echo ran ideal fcst
2614banner 28
2615#set ans = "$<"
2616#DAVE###################################################
2617
2618                                grep "SUCCESS COMPLETE" print.out.wrf_${core}_Parallel=${compopt}.exe
2619                                set success = $status
2620
2621                                #       Did making the forecast work, by that, we mean "is there an output file created?"
2622
2623                                if ( ( -e wrfout_d01_${filetag} ) && ( $success == 0 ) ) then
2624                                        if      ( $IO_FORM_NAME[$IO_FORM] == io_netcdf ) then
2625                                                ncdump -h wrfout_d01_${filetag} | grep Time | grep UNLIMITED | grep currently | grep -q 2
2626                                                set ok = $status
2627                                                set found_nans = 1
2628                                                if      ( `uname` == AIX   ) then
2629                                                        ncdump wrfout_d01_${filetag} | grep NaN >& /dev/null
2630                                                        set found_nans = $status
2631                                                else if ( `uname` == OSF1  ) then
2632                                                        ncdump wrfout_d01_${filetag} | grep nan >& /dev/null
2633                                        #               set found_nans = $status
2634                                                else if ( `uname` == Linux ) then
2635                                                        ncdump wrfout_d01_${filetag} | grep nan >& /dev/null
2636                                                        set found_nans = $status
2637                                                endif
2638                                                if ( $found_nans == 0 ) then
2639                                                        echo found nans
2640                                                        set ok = 1
2641                                                endif
2642
2643                                        else if ( $IO_FORM_NAME[$IO_FORM] == io_grib1  ) then
2644#                                               set joe_times = `../../external/io_grib1/wgrib -s -4yr wrfout_d01_${filetag} |  grep -v ":anl:" | wc -l`
2645#                                               if ( $joe_times >= 100 ) then
2646                                                set joe_times = `../../external/io_grib1/wgrib -s -4yr wrfout_d01_${filetag} |  grep "UGRD:10 m" | wc -l`
2647                                                if ( $joe_times == 2 ) then
2648                                                        set ok = 0
2649                                                else
2650                                                        set ok = 1
2651                                                endif
2652                                        endif
2653                                        if ( $ok == 0 ) then
2654                                                echo "SUMMARY generate FCST  for $core physics $phys_option parallel $compopt $esmf_lib_str PASS" >>! ${DEF_DIR}/wrftest.output
2655                                                echo "-------------------------------------------------------------" >> ${DEF_DIR}/wrftest.output
2656                                        else
2657                                                echo "SUMMARY generate FCST  for $core physics $phys_option parallel $compopt $esmf_lib_str FAIL" >>! ${DEF_DIR}/wrftest.output
2658                                                $MAIL -s "WRF FAIL FCST $ARCH[1] " $FAIL_MAIL < ${DEF_DIR}/wrftest.output
2659                                                exit ( 8 )
2660                                        endif
2661                                else
2662                                        echo "SUMMARY generate FCST  for $core physics $phys_option parallel $compopt $esmf_lib_str FAIL" >>! ${DEF_DIR}/wrftest.output
2663                                        $MAIL -s "WRF FAIL FCST $ARCH[1] " $FAIL_MAIL < ${DEF_DIR}/wrftest.output
2664                                        exit ( 9 )
2665                                endif
2666
2667                                #       We have to save this output file for our biggy comparison after all of the
2668                                #       parallel options have been considered.
2669
2670                                mv wrfout_d01_${filetag} $TMPDIR/wrfout_d01_${filetag}.${core}.${phys_option}.$compopt
2671
2672                                #       To save space, we move the executables after we are finished with them.
2673
2674                                if ( $phys_option == $Max_Ideal_Physics_Options ) then
2675                                        mv ../../main/wrf_${core}.exe.$compopt $TMPDIR/wrf_${core}.exe.$compopt
2676                                endif
2677#DAVE###################################################
2678echo fcst was a success
2679ls -ls $TMPDIR/wrfout_d01_${filetag}.${core}.${phys_option}.$compopt
2680if ( $IO_FORM_NAME[$IO_FORM] == io_netcdf ) then
2681        ncdump -v Times $TMPDIR/wrfout_d01_${filetag}.${core}.${phys_option}.$compopt | tail -20
2682endif
2683banner 29
2684#set ans = "$<"
2685#DAVE###################################################
2686
2687                                popd
2688
2689                                BYPASS_COMP_LOOP_IDEAL:
2690
2691                        end
2692
2693                endif
2694
2695                #       OK, once more, we gotta check if this is a BIT4BIT run.  If so then there
2696                #       are a number of comparisons to do.  If this is a an OPTIMIZED run, then the
2697                #       comparisons will fail the bit-wise comparisons.
2698
2699                if ( $REG_TYPE == BIT4BIT) then
2700
2701                        if ( $core == nmm_real ) then
2702       
2703                                pushd ${DEF_DIR}/regression_test/WRFV2/test/$core
2704                                set DIFFWRF = ${DEF_DIR}/regression_test/WRFV2/external/$IO_FORM_NAME[$IO_FORM]/diffwrf
2705       
2706                                if ( ( -e $TMPDIR/wrfout_d01_${filetag}.${core}.${phys_option}.$COMPOPTS[3]_1p ) && \
2707                                     ( -e $TMPDIR/wrfout_d01_${filetag}.${core}.${phys_option}.$COMPOPTS[3]_${Num_Procs}p ) ) then
2708                                        set foo1 = ( ` \ls -ls $TMPDIR/wrfout_d01_${filetag}.${core}.${phys_option}.$COMPOPTS[3]_1p `)
2709                                        set foo2 = ( ` \ls -ls $TMPDIR/wrfout_d01_${filetag}.${core}.${phys_option}.$COMPOPTS[3]_${Num_Procs}p `)
2710                                        set size1 = $foo1[6]
2711                                        set size2 = $foo2[6]
2712                                        if ( $size1 == $size2 ) then
2713                                                set RIGHT_SIZE = TRUE
2714                                        else
2715                                                set RIGHT_SIZE = FALSE
2716                                        endif
2717                                else
2718                                        set RIGHT_SIZE = FALSE
2719                                endif
2720       
2721                                #       1p vs Num_Procs MPI
2722       
2723                                rm fort.88 fort.98 >& /dev/null
2724                                if ( ( -e $TMPDIR/wrfout_d01_${filetag}.${core}.${phys_option}.$COMPOPTS[3]_1p ) && \
2725                                     ( -e $TMPDIR/wrfout_d01_${filetag}.${core}.${phys_option}.$COMPOPTS[3]_${Num_Procs}p ) && \
2726                                     ( $RIGHT_SIZE == TRUE ) ) then
2727                                        $DIFFWRF $TMPDIR/wrfout_d01_${filetag}.${core}.${phys_option}.$COMPOPTS[3]_1p \
2728                                                 $TMPDIR/wrfout_d01_${filetag}.${core}.${phys_option}.$COMPOPTS[3]_${Num_Procs}p >& /dev/null
2729                                else
2730                                        touch fort.88 fort.98
2731                                endif
2732                                if ( ! -e fort.88 ) then
2733                                        echo "SUMMARY 1 vs $Num_Procs MPI  for $core physics $phys_option $esmf_lib_str            PASS" >>! ${DEF_DIR}/wrftest.output
2734                                        echo "-------------------------------------------------------------" >> ${DEF_DIR}/wrftest.output
2735                                else
2736                                        echo "SUMMARY 1 vs $Num_Procs MPI  for $core physics $phys_option $esmf_lib_str            FAIL" >>! ${DEF_DIR}/wrftest.output
2737                                        echo "-------------------------------------------------------------" >> ${DEF_DIR}/wrftest.output
2738                                endif
2739       
2740                                popd
2741       
2742                        else if ( ${#COMPOPTS} != 1 ) then
2743       
2744                        #       If there is only a single parallel option, then we are just trying to
2745                        #       build the serial code.  That implies no comparisons are needed.
2746       
2747                                #       All of the forecasts for this set of physics and core have been
2748                                #       generated.  We now compare the WRF model output files to see
2749                                #       if they are S^2D^2.
2750               
2751                                pushd ${DEF_DIR}/regression_test/WRFV2/test/$core
2752                                set DIFFWRF = ${DEF_DIR}/regression_test/WRFV2/external/$IO_FORM_NAME[$IO_FORM]/diffwrf
2753       
2754                                #       Are we skipping the OpenMP runs?
2755       
2756                                if ( $ZAP_OPENMP == TRUE ) then
2757                                        goto BYPASS_OPENMP_SUMMARY1
2758                                endif
2759       
2760                                #       Are the files the same size?  If not, then only the initial times
2761                                #       will be compared.  That means, on a failure to run a forecast, the
2762                                #       diffwrf will give a pass.  We need to root out this evil.
2763       
2764                                if ( ( -e $TMPDIR/wrfout_d01_${filetag}.${core}.${phys_option}.$COMPOPTS[1] ) && \
2765                                     ( -e $TMPDIR/wrfout_d01_${filetag}.${core}.${phys_option}.$COMPOPTS[2] ) ) then
2766                                        set foo1 = ( ` \ls -ls $TMPDIR/wrfout_d01_${filetag}.${core}.${phys_option}.$COMPOPTS[1] `)
2767                                        set foo2 = ( ` \ls -ls $TMPDIR/wrfout_d01_${filetag}.${core}.${phys_option}.$COMPOPTS[2] `)
2768                                        set size1 = $foo1[6]
2769                                        set size2 = $foo2[6]
2770                                        if ( $size1 == $size2 ) then
2771                                                set RIGHT_SIZE_OMP = TRUE
2772                                        else
2773                                                set RIGHT_SIZE_OMP = FALSE
2774                                        endif
2775                                else
2776                                        set RIGHT_SIZE_OMP = FALSE
2777                                endif
2778       
2779                                BYPASS_OPENMP_SUMMARY1:
2780       
2781                                if ( ( -e $TMPDIR/wrfout_d01_${filetag}.${core}.${phys_option}.$COMPOPTS[1] ) && \
2782                                     ( -e $TMPDIR/wrfout_d01_${filetag}.${core}.${phys_option}.$COMPOPTS[3] ) ) then
2783                                        set foo1 = ( ` \ls -ls $TMPDIR/wrfout_d01_${filetag}.${core}.${phys_option}.$COMPOPTS[1] `)
2784                                        set foo3 = ( ` \ls -ls $TMPDIR/wrfout_d01_${filetag}.${core}.${phys_option}.$COMPOPTS[3] `)
2785                                        set size1 = $foo1[6]
2786                                        set size3 = $foo3[6]
2787                                        if ( $size1 == $size3 ) then
2788                                                set RIGHT_SIZE_MPI = TRUE
2789                                        else
2790                                                set RIGHT_SIZE_MPI = FALSE
2791                                        endif
2792                                else
2793                                        set RIGHT_SIZE_MPI = FALSE
2794                                endif
2795       
2796                                #       Are we skipping the OpenMP runs?
2797       
2798                                if ( $ZAP_OPENMP == TRUE ) then
2799                                        goto BYPASS_OPENMP_SUMMARY2
2800                                endif
2801               
2802                                #       Serial vs OpenMP
2803               
2804                                rm fort.88 fort.98 >& /dev/null
2805                                if ( ( -e $TMPDIR/wrfout_d01_${filetag}.${core}.${phys_option}.$COMPOPTS[1] ) && \
2806                                     ( -e $TMPDIR/wrfout_d01_${filetag}.${core}.${phys_option}.$COMPOPTS[2] ) && \
2807                                     ( $RIGHT_SIZE_OMP == TRUE ) ) then
2808                                        $DIFFWRF $TMPDIR/wrfout_d01_${filetag}.${core}.${phys_option}.$COMPOPTS[1] \
2809                                                 $TMPDIR/wrfout_d01_${filetag}.${core}.${phys_option}.$COMPOPTS[2] >& /dev/null
2810                                else
2811                                        touch fort.88 fort.98
2812                                endif
2813                                if ( ! -e fort.88 ) then
2814                                        echo "SUMMARY serial vs OMP  for $core physics $phys_option $esmf_lib_str            PASS" >>! ${DEF_DIR}/wrftest.output
2815                                        echo "-------------------------------------------------------------" >> ${DEF_DIR}/wrftest.output
2816                                else
2817                                        echo "SUMMARY serial vs OMP  for $core physics $phys_option $esmf_lib_str            FAIL" >>! ${DEF_DIR}/wrftest.output
2818                                        echo "-------------------------------------------------------------" >> ${DEF_DIR}/wrftest.output
2819                                endif
2820       
2821                                BYPASS_OPENMP_SUMMARY2:
2822               
2823                                #       Serial vs MPI
2824               
2825                                rm fort.88 fort.98 >& /dev/null
2826                                if ( ( -e $TMPDIR/wrfout_d01_${filetag}.${core}.${phys_option}.$COMPOPTS[1] ) && \
2827                                     ( -e $TMPDIR/wrfout_d01_${filetag}.${core}.${phys_option}.$COMPOPTS[3] ) && \
2828                                     ( $RIGHT_SIZE_MPI == TRUE ) ) then
2829                                        $DIFFWRF $TMPDIR/wrfout_d01_${filetag}.${core}.${phys_option}.$COMPOPTS[1] \
2830                                                 $TMPDIR/wrfout_d01_${filetag}.${core}.${phys_option}.$COMPOPTS[3] >& /dev/null
2831                                else
2832                                        touch fort.88 fort.98
2833                                endif
2834                                if ( ! -e fort.88 ) then
2835                                        echo "SUMMARY serial vs MPI  for $core physics $phys_option $esmf_lib_str            PASS" >>! ${DEF_DIR}/wrftest.output
2836                                        echo "-------------------------------------------------------------" >> ${DEF_DIR}/wrftest.output
2837                                else
2838                                        echo "SUMMARY serial vs MPI  for $core physics $phys_option $esmf_lib_str            FAIL" >>! ${DEF_DIR}/wrftest.output
2839                                        echo "-------------------------------------------------------------" >> ${DEF_DIR}/wrftest.output
2840                                endif
2841               
2842                                popd
2843       
2844                        endif
2845
2846                        # Generate and archive baseline or compare against baseline
2847
2848                        if ( $core != nmm_real ) then
2849                                if ( $GENERATE_BASELINE != FALSE ) then
2850                                        if ( ! -d $GENERATE_BASELINE ) then
2851                                                echo "ERROR:  Baseline directory ${GENERATE_BASELINE} does not exist" >>! \
2852                                                     ${DEF_DIR}/wrftest.output
2853                                                exit ( 10 ) 
2854                                        else
2855                                                # Archive serial output file to baseline
2856                                                pushd ${DEF_DIR}/regression_test/WRFV2/test/$core
2857                                                set basefilenm = wrfout_d01_${filetag}.${core}.${phys_option}
2858                                                set basefile = ${GENERATE_BASELINE}/${basefilenm}
2859                                                set outfile = $TMPDIR/${basefilenm}.$COMPOPTS[1]
2860                                                if ( -e $outfile ) then
2861                                                        cp $outfile $basefile || \
2862                                                          ( echo "ERROR:  cannot copy ${outfile} to ${basefile}" >>! \
2863                                                           ${DEF_DIR}/wrftest.output; exit 10 )
2864                                                else
2865                                                        echo "ERROR:  Cannot archive baseline, file $outfile does not exist" >>! \
2866                                                             ${DEF_DIR}/wrftest.output
2867                                        echo "SUMMARY baseline archived in ${GENERATE_BASELINE} for $core physics $phys_option  FAIL" >>! \
2868                                                             ${DEF_DIR}/wrftest.output
2869                                                        echo "-------------------------------------------------------------" >> \
2870                                                             ${DEF_DIR}/wrftest.output
2871                                                        exit ( 10 ) 
2872                                                endif
2873                                        echo "SUMMARY baseline archived in ${GENERATE_BASELINE} for $core physics $phys_option  PASS" >>! \
2874                                                     ${DEF_DIR}/wrftest.output
2875                                                echo "-------------------------------------------------------------" >> \
2876                                                     ${DEF_DIR}/wrftest.output
2877                                                popd
2878                                        endif
2879                                endif
2880                                if ( $COMPARE_BASELINE != FALSE ) then
2881                                        if ( ! -d $COMPARE_BASELINE ) then
2882                                                echo "${0}: ERROR::  Baseline directory ${COMPARE_BASELINE} does not exist" >>! \
2883                                                     ${DEF_DIR}/wrftest.output
2884                                                exit ( 10 ) 
2885                                        else
2886                                                # Compare against baseline output file
2887                                                set basefilenm = wrfout_d01_${filetag}.${core}.${phys_option}
2888                                                set basefile = ${COMPARE_BASELINE}/${basefilenm}
2889                                                set DIFFWRF = ${DEF_DIR}/regression_test/WRFV2/external/$IO_FORM_NAME[$IO_FORM]/diffwrf
2890                                                set testdir = ${DEF_DIR}/regression_test/WRFV2/test/$core
2891                                                pushd ${testdir}
2892                                                foreach compopt ( $COMPOPTS )
2893                                                        set cmpfile = $TMPDIR/${basefilenm}.$compopt
2894                                                        rm fort.88 fort.98 >& /dev/null
2895                                                        if ( ( -e ${basefile} ) && ( -e ${cmpfile} ) ) then
2896                                                                #       Are the files the same size?  If not, then only the initial times
2897                                                                #       will be compared.  That means, on a failure to run a forecast, the
2898                                                                #       diffwrf will give a pass.  We need to root out this evil.
2899                                                                set foob = ( ` \ls -ls ${basefile} `)
2900                                                                set fooc = ( ` \ls -ls ${cmpfile} `)
2901                                                                set sizeb = $foob[6]
2902                                                                set sizec = $fooc[6]
2903                                                                if ( $sizeb == $sizec ) then
2904                                                                        $DIFFWRF ${basefile} ${cmpfile} >& /dev/null
2905                                                                        if ( -e fort.88 ) then
2906                        echo "FAIL:  Baseline comparison, file ${cmpfile} did not match ${basefile} according to diffwrf" >>! \
2907                                                                                ${DEF_DIR}/wrftest.output
2908                                                                        endif
2909                                                                else
2910                                                                        touch fort.88 fort.98
2911                        echo "FAIL:  Baseline comparison, files ${cmpfile} and ${basefile} have different sizes" >>! \
2912                                                                        ${DEF_DIR}/wrftest.output
2913                                                                endif
2914                                                        else
2915                        echo "FAIL:  Baseline comparison, either ${cmpfile} or ${basefile} does not exist" >>! \
2916                                                                        ${DEF_DIR}/wrftest.output
2917                                                                touch fort.88 fort.98
2918                                                        endif
2919                                                        if ( ! -e fort.88 ) then
2920                        echo "SUMMARY compare vs baseline ${COMPARE_BASELINE} for $core physics $phys_option compopt $compopt $esmf_lib_str  PASS" >>! \
2921                                                                     ${DEF_DIR}/wrftest.output
2922                                                                echo "-------------------------------------------------------------" >> \
2923                                                                     ${DEF_DIR}/wrftest.output
2924                                                        else
2925                        echo "SUMMARY compare vs baseline ${COMPARE_BASELINE} for $core physics $phys_option compopt $compopt $esmf_lib_str  FAIL" >>! \
2926                                                                     ${DEF_DIR}/wrftest.output
2927                                                                echo "-------------------------------------------------------------" >> \
2928                                                                     ${DEF_DIR}/wrftest.output
2929                                                        endif
2930                                                end
2931                                                popd
2932                                        endif
2933                                endif
2934                        else if ( $core == nmm_real ) then
2935                                if ( $GENERATE_BASELINE != FALSE ) then
2936                                        if ( ! -d $GENERATE_BASELINE ) then
2937                                                echo "ERROR:  Baseline directory ${GENERATE_BASELINE} does not exist" >>! \
2938                                                     ${DEF_DIR}/wrftest.output
2939                                                exit ( 10 ) 
2940                                        else
2941                                                # Archive serial output file to baseline
2942                                                pushd ${DEF_DIR}/regression_test/WRFV2/test/$core
2943                                                set basefilenm = wrfout_d01_${filetag}.${core}.${phys_option}
2944                                                set basefile = ${GENERATE_BASELINE}/${basefilenm}
2945                                                set outfile = $TMPDIR/${basefilenm}.$COMPOPTS[3]_1p
2946                                                if ( -e $outfile ) then
2947                                                        cp $outfile $basefile || \
2948                                                          ( echo "ERROR:  cannot copy ${outfile} to ${basefile}" >>! \
2949                                                           ${DEF_DIR}/wrftest.output; exit 10 )
2950                                                else
2951                                                        echo "ERROR:  Cannot archive baseline, file $outfile does not exist" >>! \
2952                                                             ${DEF_DIR}/wrftest.output
2953                                        echo "SUMMARY baseline archived in ${GENERATE_BASELINE} for $core physics $phys_option  FAIL" >>! \
2954                                                             ${DEF_DIR}/wrftest.output
2955                                                        echo "-------------------------------------------------------------" >> \
2956                                                             ${DEF_DIR}/wrftest.output
2957                                                        exit ( 10 ) 
2958                                                endif
2959                                        echo "SUMMARY baseline archived in ${GENERATE_BASELINE} for $core physics $phys_option  PASS" >>! \
2960                                                     ${DEF_DIR}/wrftest.output
2961                                                echo "-------------------------------------------------------------" >> \
2962                                                     ${DEF_DIR}/wrftest.output
2963                                                popd
2964                                        endif
2965                                endif
2966                                if ( $COMPARE_BASELINE != FALSE ) then
2967                                        if ( ! -d $COMPARE_BASELINE ) then
2968                                                echo "${0}: ERROR::  Baseline directory ${COMPARE_BASELINE} does not exist" >>! \
2969                                                     ${DEF_DIR}/wrftest.output
2970                                                exit ( 10 ) 
2971                                        else
2972                                                # Compare against baseline output file
2973                                                set basefilenm = wrfout_d01_${filetag}.${core}.${phys_option}
2974                                                set basefile = ${COMPARE_BASELINE}/${basefilenm}
2975                                                set DIFFWRF = ${DEF_DIR}/regression_test/WRFV2/external/$IO_FORM_NAME[$IO_FORM]/diffwrf
2976                                                set testdir = ${DEF_DIR}/regression_test/WRFV2/test/$core
2977                                                pushd ${testdir}
2978                                                set compopt = $COMPOPTS[3]
2979                                                foreach proc ( 1p 4p )
2980                                                        set cmpfile = $TMPDIR/${basefilenm}.${compopt}_${proc}
2981                                                        rm fort.88 fort.98 >& /dev/null
2982                                                        if ( ( -e ${basefile} ) && ( -e ${cmpfile} ) ) then
2983                                                                #       Are the files the same size?  If not, then only the initial times
2984                                                                #       will be compared.  That means, on a failure to run a forecast, the
2985                                                                #       diffwrf will give a pass.  We need to root out this evil.
2986                                                                set foob = ( ` \ls -ls ${basefile} `)
2987                                                                set fooc = ( ` \ls -ls ${cmpfile} `)
2988                                                                set sizeb = $foob[6]
2989                                                                set sizec = $fooc[6]
2990                                                                if ( $sizeb == $sizec ) then
2991                                                                        $DIFFWRF ${basefile} ${cmpfile} >& /dev/null
2992                                                                        if ( -e fort.88 ) then
2993                        echo "FAIL:  Baseline comparison, file ${cmpfile} did not match ${basefile} according to diffwrf" >>! \
2994                                                                                ${DEF_DIR}/wrftest.output
2995                                                                        endif
2996                                                                else
2997                                                                        touch fort.88 fort.98
2998                        echo "FAIL:  Baseline comparison, files ${cmpfile} and ${basefile} have different sizes" >>! \
2999                                                                        ${DEF_DIR}/wrftest.output
3000                                                                endif
3001                                                        else
3002                        echo "FAIL:  Baseline comparison, either ${cmpfile} or ${basefile} does not exist" >>! \
3003                                                                        ${DEF_DIR}/wrftest.output
3004                                                                touch fort.88 fort.98
3005                                                        endif
3006                                                        if ( ! -e fort.88 ) then
3007                        echo "SUMMARY compare vs baseline ${COMPARE_BASELINE} for $core physics $phys_option compopt $compopt $esmf_lib_str  PASS" >>! \
3008                                                                     ${DEF_DIR}/wrftest.output
3009                                                                echo "-------------------------------------------------------------" >> \
3010                                                                     ${DEF_DIR}/wrftest.output
3011                                                        else
3012                        echo "SUMMARY compare vs baseline ${COMPARE_BASELINE} for $core physics $phys_option compopt $compopt $esmf_lib_str  FAIL" >>! \
3013                                                                     ${DEF_DIR}/wrftest.output
3014                                                                echo "-------------------------------------------------------------" >> \
3015                                                                     ${DEF_DIR}/wrftest.output
3016                                                        endif
3017                                                end
3018                                                popd
3019                                        endif
3020                                endif
3021                                goto ALL_SHE_WROTE_FOR_NMM
3022                        endif
3023                        # End of generate and archive baseline or compare against baseline
3024
3025                endif
3026
3027                BOTTOM_OF_PHYSICS_LOOP:
3028               
3029        end
3030
3031ALL_SHE_WROTE_FOR_NMM:
3032
3033        echo "-------------------------------------------------------------" >> ${DEF_DIR}/wrftest.output
3034
3035        if ( $clrm ) then
3036          popd
3037        endif
3038               
3039end
3040
3041
3042#       How long did this take.
3043
3044set end = ( `date` )
3045echo "Start WRF Regression: $start " >> ${DEF_DIR}/wrftest.output
3046echo "End   WRF Regression: $end   " >> ${DEF_DIR}/wrftest.output
3047
3048#       We have done all of the tests, and placed the PASS FAIL labels in the
3049#       output file.  If there are any FAIL messages, we are in trouble.
3050
3051grep FAIL ${DEF_DIR}/wrftest.output
3052set ok = $status
3053
3054#       Send email of the status.
3055
3056if ( $ok == 0 ) then
3057        $MAIL -s "REGRESSION FAILURE $ARCH[1] " $FAIL_MAIL < ${DEF_DIR}/wrftest.output
3058else
3059        $MAIL -s "REGRESSION SUCCESS $ARCH[1] " $GOOD_MAIL < ${DEF_DIR}/wrftest.output
3060endif
3061
3062#       Clean left up detritus
3063
3064cd $CUR_DIR
3065
3066rm -rf dom_*eal >& /dev/null
3067rm -rf phys_real_* >& /dev/null
3068rm -rf phys_quarter_* >& /dev/null
3069rm -rf phys_b_wave_* >& /dev/null
3070rm -rf version_info >& /dev/null
3071rm -rf machfile >& /dev/null
3072rm -rf fdda_real* >& /dev/null
3073rm -rf time_real_* >& /dev/null
3074rm -rf dyn_real_* >& /dev/null
Note: See TracBrowser for help on using the repository browser.