Changeset 1be76d8 in flexpart.git


Ignore:
Timestamp:
Oct 1, 2015, 3:29:19 PM (9 years ago)
Author:
Ignacio Pisso <Ignacio.Pisso@…>
Branches:
master, 10.4.1_pesei, GFS_025, bugfixes+enhancements, dev, release-10, release-10.4.1, scaling-bug
Children:
ac5b22b
Parents:
1333256
Message:

add 7 basic tests

Location:
tests/NILU
Files:
40 added
31 deleted
3 edited
4 moved

Legend:

Unmodified
Added
Removed

  • tests/NILU/run_test

    r1333256 r1be76d8  
    22
    33#input FP exec
    4 test_dir=$1
    5 FP_exec=$2
    6 base_dir=$3
    7 
    8 echo FP_exec=$FP_exec
     4test_dir=$1  # name of dir with pathname and options to be used
     5FP_exec=$2   # path to executable to be run
     6base_dir=$3  # base for the output
     7options_default=$4
     8 
     9#echo FP_exec=$FP_exec
    910echo $FLEXOUT_TMP
    1011
    … …  
    1415
    1516current_dir=$(pwd)
    16 echo test_tmp=$test_tmp
     17#echo test_tmp=$test_tmp
    1718mkdir $test_tmp
    1819
    19 #copy default optios
    20 cp -r options_default $test_tmp/options
     20echo copy default options $options_default  to temp dir $test_tmp/options
     21cp -r $options_default $test_tmp/options
    2122#override default options
    22 cp -r ${test_dir}/options/* $test_tmp/options/
     23echo copy this test option files to temp dir
     24cp  ${test_dir}/options/* $test_tmp/options/
     25echo copy SPECIES from ${test_dir}/options/SPECIES/* to $test_tmp/options/SPECIES/
     26cp  ${test_dir}/options/SPECIES/* $test_tmp/options/SPECIES/
    2327# copy paths form repo source exec tmp dir
    2428cp ${test_dir}/pathnames $test_tmp
    … …  
    2832
    2933cd $test_tmp
    30 echo $(pwd)
    31 echo $FP_exec
     34#echo $(pwd)
     35#echo $FP_exec
    3236#  '&>' o is bash for both out and err
    33 echo run $FP_exec   to  out.log #&
    34 $FP_exec   &> out.log #&
     37#echo run $FP_exec   to  out.log #&
     38#$FP_exec   &> out.log #&
     39$FP_exec > stdout.log 2> stderr.log #&
    3540
    3641

  • tests/NILU/run_tests

    r1333256 r1be76d8  
    11#!/bin/bash
    22
    3 FLEXOUT_TMP=/home/ignacio/repos/flextest/
     3#input:
     4# FP_exec : path to flexpart executable 
     5# path_flextest: directory  path to output base directory
     6# test_names (bash array): list of tests to be run
     7
     8# defaults
     9FP_exec=/home/ignacio/repos/flexpart/src/FP_ecmwf_gfortran
     10path_flextest=/home/ignacio/repos/flextest/
     11options_default=/home/ignacio/repos/flexpart/options/
     12#declare -a test_names=('1' 'HelloWorld' 'Fwd1' 'Fwd2' 'Bwd1' 'Volc' '2')
     13#test_names[0]='1'
     14test_names[0]='HelloWorld'
     15
     16group_name=tests
     17
     18source ./$1
     19
    420
    521suffix=$(date "+%y-%m-%d_%H%M%S")
    6 #test_tmp=${FLEXOUT_TMP}${base_dir}${test_dir}_${suffix}
    7 base_dir=${FLEXOUT_TMP}/tests_${suffix}/
    8 mkdir $base_dir
    9 
    10 FP_exec=/home/ignacio/repos/flexpart/src/FP_ecmwf_gfortran
    11 
    12 #test_dir=$1
    13 #FP_exec=$2
    14 #base_dir=$3
    15 
    16 ./run_test test_1  /home/ignacio/repos/flexpart/src/FP_ecmwf_gfortran $base_dir
    17 ./run_test test_2  /home/ignacio/repos/flexpart/src/FP_ecmwf_gfortran $base_dir
    18 
    19 # 6 species
    20 ./run_test test_Volc  /home/ignacio/repos/flexpart/src/FP_ecmwf_gfortran $base_dir
     22#tests_dir=tests_${suffix}/
     23tests_dir=${group_name}_${suffix}/
     24path_tests=${path_flextest}/$tests_dir
     25mkdir $path_tests
    2126
    2227
    23 #exit
     28#declare Bash Array
     29#declare -a test_names=('1' '2' 'HelloWorld')
     30#declare -a test_names=('1' '2' 'HelloWorld' 'Fwd1' 'Fwd2' 'Bwd1' 'Volc')
     31#declare -a test_names=('2' 'Volc')
     32#declare -a test_names=('1' 'HelloWorld' 'Fwd1' 'Fwd2' 'Bwd1' 'Volc' '2')
    2433
    25 ./run_test test_HelloWorld  /home/ignacio/repos/flexpart/src/FP_ecmwf_gfortran $base_dir
    26 ./run_test test_Fwd1  /home/ignacio/repos/flexpart/src/FP_ecmwf_gfortran $base_dir
    27 ./run_test test_Fwd2  /home/ignacio/repos/flexpart/src/FP_ecmwf_gfortran $base_dir
     34#blabla
     35#Print the Whole Bash Array
     36s_test_names=${test_names[@]}
     37# Length of the Bash Array
     38n_tests=${#test_names[@]}
     39#echo n_tests=$n_tests
     40echo run $n_tests  tests:  $s_test_names
     41
     42#test
     43#tests_dir='tests_15-08-19_113139/'
     44#path_tests=${path_flextest}/$tests_dir
     45
     46
     47outlog=stdout.log
     48i=0
     49i_end=$n_tests
     50while [ $i -lt $i_end ] ;
     51do
     52  let j=i+1
     53  name=${test_names[i]};
     54  echo test no. $j name= $name
     55  echo in $path_tests/test_$name
     56  ./run_test test_$name  $FP_exec $path_tests $options_default
     57  echo terminated
     58  success="$(cat $path_tests/test_$name/$outlog  | grep CONGRATULATIONS  | wc -l)"
     59  echo success=$success
     60let i=i+1;
     61done
     62 
     63tot_successes="$(cat $path_tests/test_*/$outlog  | grep CONGRATULATIONS  | wc -l)"
     64tot_runs="$(ls -1 $path_tests | wc -l)"
     65
     66echo ' '
     67echo finished  $tot_runs tests of $n_tests
     68echo $tot_successes seccessful
     69exit
     70
     71
     72
     73
     74
     75
     76
     77
     78
     79
     80
     81
     82
     83./run_test test_1  $FP_exec $path_tests
     84./run_test test_2  $FP_exec $path_tests
     85# 6 species
     86./run_test test_Volc  $FP_exec $path_tests
     87./run_test test_HelloWorld  $FP_exec $path_tests
     88./run_test test_Fwd1  $FP_exec $path_tests
     89./run_test test_Fwd2  $FP_exec $path_tests
     90./run_test test_Bwd1  $FP_exec $path_tests
     91
     92

  • tests/NILU/test_1/pathnames

    rdc69f6a r1be76d8  
    22./output/
    33/
    4 /flex_wrk/WIND_FIELDS/AVAILABLE_ECMWF_OPER_fields_global
     4/xnilu_wrk/flex_wrk/WIND_FIELDS/AVAILABLE_ECMWF_OPER_fields_global
    55============================================
    66

  • tests/NILU/test_2/options/SPECIES/SPECIES_026

    • Property mode changed from 100755 to 100644
    rdc69f6a r1be76d8  
    55*                                                                          *
    66****************************************************************************
    7 Kr-85              Tracer name
     7CH4                Tracer name
    88-999.9             Species half life
    9 -9.9E-09           Below cloud scavenging - A
    10                    Below cloud scavenging - B
     9-9.9E-09           Wet deposition - A
     10                   Wet deposition - B
    1111-9.9E-09           In-cloud scavenging - Ai (cl=Ai*prec**Bi)
    1212-9.9               In-cloud scavenging - Bi (cl=Ai*prec**Bi)
    … …  
    2020                   Dry deposition (particles) - dsig
    2121-9.99              Alternative: dry deposition velocity
    22  -9.99             molweight
    23 -9.9E-09           OH Reaction rate at 25 deg, [cm^3/sec]
     2216.000             molweight
     23 2.5E-15           OH Reaction rate at 25 deg, [cm^3/sec]
    2424-9                 number of associated specias (neg. none)
    2525-99.99             KOA - organic matter air partitioning

  • tests/NILU/test_Bwd1/options/COMMAND

    • Property mode changed from 100755 to 100644
    rdc69f6a r1be76d8  
    66********************************************************************************
    77
    8  1              LDIRECT           1 FOR FORWARD SIMULATION, -1 FOR BACKWARD SIMULATION
    9 20040720 000000 YYYYMMDD HHMISS   BEGINNING DATE OF SIMULATION
    10 20040721 120000 YYYYMMDD HHMISS   ENDING DATE OF SIMULATION
    11 10800           SSSSS             OUTPUT EVERY SSSSS SECONDS
    12 10800           SSSSS             TIME AVERAGE OF OUTPUT (IN SSSSS SECONDS)
    13 900             SSSSS             SAMPLING RATE OF OUTPUT (IN SSSSS SECONDS)
    14 9999999         SSSSSSS           TIME CONSTANT FOR PARTICLE SPLITTING (IN SECONDS)
    15 900             SSSSS             SYNCHRONISATION INTERVAL OF FLEXPART (IN SECONDS)
    16 -5.0            CTL               FACTOR, BY WHICH TIME STEP MUST BE SMALLER THAN TL
    17 4               IFINE             DECREASE OF TIME STEP FOR VERTICAL MOTION BY FACTOR IFINE
    18 3               IOUT              1 CONC. (RESID. TIME FOR BACKWARD RUNS) OUTPUT,2 MIX. RATIO OUTPUT,3 BOTH,4 PLUME TRAJECT.,5=1+4
    19 0               IPOUT             PARTICLE DUMP: 0 NO, 1 EVERY OUTPUT INTERVAL, 2 ONLY AT END
    20 1               LSUBGRID          SUBGRID TERRAIN EFFECT PARAMETERIZATION: 1 YES, 0 NO
    21 1               LCONVECTION       CONVECTION: 1 YES, 0 NO
    22 0               LAGESPECTRA       AGE SPECTRA: 1 YES, 0 NO
    23 0               IPIN              CONTINUE SIMULATION WITH DUMPED PARTICLE DATA: 1 YES, 0 NO
    24 0               IOUTPUTFOREACHREL CREATE AN OUPUT FILE FOR EACH RELEASE LOCATION: 1 YES, 0 NO
    25 0               IFLUX             CALCULATE FLUXES: 1 YES, 0 NO
    26 0               MDOMAINFILL       DOMAIN-FILLING TRAJECTORY OPTION: 1 YES, 0 NO
    27 1               IND_SOURCE        1=MASS UNIT , 2=MASS MIXING RATIO UNIT
    28 1               IND_RECEPTOR      1=MASS UNIT , 2=MASS MIXING RATIO UNIT
    29 0               MQUASILAG         QUASILAGRANGIAN MODE TO TRACK INDIVIDUAL PARTICLES: 1 YES, 0 NO
    30 0               NESTED_OUTPUT     SHALL NESTED OUTPUT BE USED? YES, 0 NO
    31 2               LINIT_COND        INITIAL COND. FOR BW RUNS: 0=NO,1=MASS UNIT,2=MASS MIXING RATIO UNIT
    32 0               SURF_ONLY         IF THIS IS SET TO 1, OUTPUT IS WRITTEN ONLY OUT FOR LOWEST LAYER
     81. __                3X, I2
     9   -1       
     10   LDIRECT           1 FOR FORWARD SIMULATION, -1 FOR BACKWARD SIMULATION
     11
     122. ________ ______   3X, I8, 1X, I6
     13   20070121 090000
     14   YYYYMMDD HHMISS   BEGINNING DATE OF SIMULATION
     15
     163. ________ ______   3X, I8, 1X, I6
     17   20070122 180000
     18   YYYYMMDD HHMISS   ENDING DATE OF SIMULATION
     19
     204. _____             3X, I5
     21    3600
     22   SSSSS             OUTPUT EVERY SSSSS SECONDS
     23
     245. _____             3X, I5
     25    3600
     26   SSSSS             TIME AVERAGE OF OUTPUT (IN SSSSS SECONDS)
     27
     286. _____             3X, I5
     29     300
     30   SSSSS             SAMPLING RATE OF OUTPUT (IN SSSSS SECONDS)
     31
     327. _________         3X, I9
     33   999999999
     34   SSSSSSSSS         TIME CONSTANT FOR PARTICLE SPLITTING (IN SECONDS)
     35
     368. _____             3X, I5
     37     300
     38   SSSSS             SYNCHRONISATION INTERVAL OF FLEXPART (IN SECONDS)
     39
     409.  ---.--           4X, F6.4
     41      3.0
     42    CTL              FACTOR, BY WHICH TIME STEP MUST BE SMALLER THAN TL
     43
     4410. ---              4X, I3
     45      4
     46    IFINE            DECREASE OF TIME STEP FOR VERTICAL MOTION BY FACTOR IFINE
     47
     4811. -                4X, I1
     49    1 
     50    IOUT             1 CONCENTRATION (RESIDENCE TIME FOR BACKWARD RUNS) OUTPUT, 2 MIXING RATIO OUTPUT, 3 BOTH,4 PLUME TRAJECT., 5=1+4
     51
     5212. -                4X, I1
     53    0 
     54    IPOUT            PARTICLE DUMP: 0 NO, 1 EVERY OUTPUT INTERVAL, 2 ONLY AT END
     55
     5613. _                4X, I1
     57    1
     58    LSUBGRID         SUBGRID TERRAIN EFFECT PARAMETERIZATION: 1 YES, 0 NO
     59
     6014. _                4X, I1
     61    0
     62    LCONVECTION      CONVECTION: 1 YES, 0 NO
     63
     6415. _                4X, I1
     65    0
     66    LAGESPECTRA      AGE SPECTRA: 1 YES, 0 NO
     67
     6816. _                4X, I1
     69    0
     70    IPIN             CONTINUE SIMULATION WITH DUMPED PARTICLE DATA: 1 YES, 0 NO
     71
     7217. _               
     73    1                4X,I1
     74    IOFR             IOUTPUTFOREACHREL CREATE AN OUPUT FILE FOR EACH RELEASE LOCATION: 1 YES, 0 NO
     75
     7618. _                4X, I1
     77    0
     78    IFLUX            CALCULATE FLUXES: 1 YES, 0 NO
     79
     8019. _                4X, I1
     81    0
     82    MDOMAINFILL      DOMAIN-FILLING TRAJECTORY OPTION: 1 YES, 0 NO, 2 STRAT. O3 TRACER
     83
     8420. _                4X, I1
     85    1
     86    IND_SOURCE       1=MASS UNIT , 2=MASS MIXING RATIO UNIT
     87
     8821. _                4X, I1
     89    1
     90    IND_RECEPTOR     1=MASS UNIT , 2=MASS MIXING RATIO UNIT
     91
     9222. _                4X, I1
     93    0
     94    MQUASILAG        QUASILAGRANGIAN MODE TO TRACK INDIVIDUAL PARTICLES: 1 YES, 0 NO
     95
     9623. _                4X, I1
     97    0
     98    NESTED_OUTPUT    SHALL NESTED OUTPUT BE USED? 1 YES, 0 NO
     99
     10024. _                4X, I1
     101    0
     102    LINIT_COND       INITIAL COND. FOR BW RUNS: 0=NO,1=MASS UNIT,2=MASS MIXING RATIO UNIT
     103
     10425. _                4X, I1
     105    0
     106    SURF_ONLY        IF THIS IS SET TO 1, OUTPUT IS WRITTEN ONLY OUT FOR LOWEST LAYER
    33107
    34108
    35 1. Simulation direction, 1 for forward, -1 for backward in time
     1091. Simulation direction, 1 for forward, -1 for backward in time
     110        (consult Seibert and Frank, 2004 for backward runs)
    36111
    371122. Beginning date and time of simulation. Must be given in format
    … …  
    100175
    10117620. IND_SOURCE switches between different units for concentrations at the source
    102     NOTE that in backward simulations the release of computational particles
     177    NOTE that in backward simulations the release of computational particles 
    103178    takes place at the "receptor" and the sampling of particles at the "source".
    104179          1=mass units (for bwd-runs = concentration)
    105           2=mass mixing ratio units
     180          2=mass mixing ratio units 
    10618121. IND_RECEPTOR switches between different units for concentrations at the receptor
    107182          1=mass units (concentrations)
    108           2=mass mixing ratio units
     183          2=mass mixing ratio units 
    109184
    11018522. MQUASILAG indicates whether particles shall be numbered consecutively (1) or

  • tests/NILU/test_HelloWorld/options/COMMAND

    • Property mode changed from 100755 to 100644
    rdc69f6a r1be76d8  
    1111
    12122. ________ ______   3X, I8, 1X, I6
    13    20110310 000000
     13   20070121 090000
    1414   YYYYMMDD HHMISS   BEGINNING DATE OF SIMULATION
    1515
    16163. ________ ______   3X, I8, 1X, I6
    17    20110310 120000
     17   20070121 180000
    1818   YYYYMMDD HHMISS   ENDING DATE OF SIMULATION
    1919
    20204. _____             3X, I5
    21    10800
     21    3600
    2222   SSSSS             OUTPUT EVERY SSSSS SECONDS
    2323
    24245. _____             3X, I5
    25    10800
     25    3600
    2626   SSSSS             TIME AVERAGE OF OUTPUT (IN SSSSS SECONDS)
    2727
    28286. _____             3X, I5
    29      900
     29     300
    3030   SSSSS             SAMPLING RATE OF OUTPUT (IN SSSSS SECONDS)
    3131
    … …  
    3535
    36368. _____             3X, I5
    37      900
     37     300
    3838   SSSSS             SYNCHRONISATION INTERVAL OF FLEXPART (IN SECONDS)
    3939
    … …  
    4747
    484811. -                4X, I1
    49     3 
     49    1 
    5050    IOUT             1 CONCENTRATION (RESIDENCE TIME FOR BACKWARD RUNS) OUTPUT, 2 MIXING RATIO OUTPUT, 3 BOTH,4 PLUME TRAJECT., 5=1+4
    5151
    … …  
    5959
    606014. _                4X, I1
    61     1
     61    0
    6262    LCONVECTION      CONVECTION: 1 YES, 0 NO
    6363

  • tests/NILU/test_HelloWorld/options/RELEASES

    • Property mode changed from 100755 to 100644
    rdc69f6a r1be76d8  
    1313___                        i3    Total number of species emitted
    1414
    15  24
     15  1
    1616___                        i3    Index of species in file SPECIES
    1717
    1818=========================================================================
    19 20040720      0
     1920070121 090000
    2020________ ______            i8,1x,i6 Beginning date and time of release
    2121
    22 20040720 120000
     2220070121 100000
    2323________ ______            i8,1x,i6 Ending date and time of release
    2424
    25    9.4048
     25   2.116
    2626____.____                  f9.4  Longitude [DEG] of lower left corner
    2727
    28   48.5060
     28  41.384
    2929____.____                  f9.4  Latitude [DEG] of lower left corner
    3030
    31    9.5067
     31   2.116
    3232____.____                  f9.4  Longitude [DEG] of upper right corner
    3333
    34   48.5158
     34  41.384
    3535____.____                  f9.4  Latitude [DEG] of upper right corner
    3636
    37         2
     37        1
    3838_________                  i9    1 for m above ground, 2 for m above sea level
    3939
    40  6933.60
     40   60.0
    4141_____.___                  f10.3 Lower z-level (in m agl or m asl)
    4242 
    43  6950.40
     43   60.0
    4444_____.___                  f10.3 Upper z-level (in m agl or m asl)
    4545 
    46     20000               
     46    1000               
    4747_________                  i9    Total number of particles to be released
    4848
    … …  
    5353________________________________________   character*40 comment
    5454+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
    55 20040720      0
    56 ________ ______            i8,1x,i6 Beginning date and time of release
    57 
    58 20040720 120000
    59 ________ ______            i8,1x,i6 Ending date and time of release
    60 
    61    9.3038
    62 ____.____                  f9.4  Longitude [DEG] of lower left corner
    63 
    64   48.5158
    65 ____.____                  f9.4  Latitude [DEG] of lower left corner
    66 
    67    9.4048
    68 ____.____                  f9.4  Longitude [DEG] of upper right corner
    69 
    70   48.5906
    71 ____.____                  f9.4  Latitude [DEG] of upper right corner
    72 
    73         2
    74 _________                  i9    1 for m above ground, 2 for m above sea level
    75 
    76  6833.50
    77 _____.___                  f10.3 Lower z-level (in m agl or m asl)
    78  
    79  6950.40
    80 _____.___                  f10.3 Upper z-level (in m agl or m asl)
    81  
    82     20000               
    83 _________                  i9    Total number of particles to be released
    84 
    85 1.0000E00
    86 _.____E__                  e9.4  Total mass emitted
    87 
    88 RELEASE_TEST2
    89 ________________________________________   character*40 comment
    90 +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
Note: See TracChangeset for help on using the changeset viewer.
hosted by ZAMG