Concentration matter and Contour display matter.

Topics about the HYSPLIT dispersion model.
Post Reply
sol
Posts: 6
Joined: February 5th, 2020, 9:24 pm
Registered HYSPLIT User: No

Concentration matter and Contour display matter.

Post by sol » February 5th, 2020, 10:04 pm

Dear all.
I ran Dispersion model.
I used 'gfs0p25 data' which from NCAR-UCAR homepage as a meteorological data.
I used 'api2arl_v1.exe' tools to convert gfs0p25 data to arl data.
But when I ran HYSPLIT Dispersion model, I got concentration data as 0 in every grid cell that i made. And also contour image is appear one dark line, not color area.
I attached contour image , message file and default_conc file.

here is message file code

Code: Select all

Simulation Date (CCYYMMDD): 20200205
 Simulation Time (HHMMSS.S): 101732.686
  NOTICE metini: turbulence undefined (kblt=0); setting to Kantha-Clayson
 
   KG   KT DATA  SIZE   NX   NY
    1    1 kwbc  27.8 1440  721
 
  NOTICE   main: number meteo grids and times -            1           1
  
 HYSPLIT4
 Last Changed Rev: 951
 Last Changed Date: 2018-04-21 18:15:01 -0400 (Sat, 21 Apr 2018)
 
 ------------- Start Namelist configuration -------------------
 Internal grid parameters (nlvl,aa,bb,cc):           23   20.00000    
  -15.00000       5.000000    
 &SETUP
 INITD   =           0,
 KHMAX   =        9999,
 NUMPAR  =        2500,
 MAXPAR  =       10000,
 MAXDIM  =           1,
 QCYCLE  =  0.0000000E+00,
 P10F    =   1.000000    ,
 KMIXD   =           0,
 KMIX0   =         250,
 FRME    =  0.1000000    ,
 FRMR    =  0.0000000E+00,
 KRND    =           6,
 DELT    =  0.0000000E+00,
 KDEF    =           0,
 TKERD   =  0.1800000    ,
 TKERN   =  0.1800000    ,
 NDUMP   =           0,
 NCYCL   =           0,
 TRATIO  =  0.7500000    ,
 ISOT    =         -99,
 MGMIN   =        1440,
 KMSL    =           0,
 NSTR    =           0,
 MHRS    =           0,
 NVER    =           0,
 CPACK   = 99*1,
 CMASS   =           0,
 ICHEM   =           0,
 EFILE   =                                                                                 ,
 VSCALE  =   200.0000    ,
 VSCALES =   200.0000    ,
 VSCALEU =   200.0000    ,
 K10M    =           1,
 VDIST   = VMSDIST                                                                         ,
 KHINP   =           0,
 PINBC   = PARINBC                                                                         ,
 TOUT    =           0,
 TM_PRES =           0,
 TM_TPOT =           0,
 TM_TAMB =           0,
 TM_RAIN =           0,
 TM_MIXD =           0,
 TM_TERR =           0,
 HSCALE  =   10800.00    ,
 TM_DSWF =           0,
 TM_RELH =           0,
 TM_SPHU =           0,
 TM_MIXR =           0,
 DXF     =   1.000000    ,
 DYF     =   1.000000    ,
 DZF     =  9.9999998E-03,
 PINPF   = PARINIT                                                                         ,
 POUTF   = PARDUMP                                                                         ,
 NINIT   =           1,
 KAGL    =           0,
 KPUFF   =           0,
 KBLT    =           0,
 KZMIX   =           0,
 TVMIX   =   1.000000    ,
 KBLS    =           1,
 CONAGE  =          48,
 KSPL    =           1,
 FRHS    =   1.000000    ,
 FRVS    =  9.9999998E-03,
 FRTS    =  0.1000000    ,
 FRHMAX  =   3.000000    ,
 SPLITF  =   1.000000    ,
 MESSG   = MESSAGE                                                                         ,
 NBPTYP  =           1,
 RHB     =          80,
 RHT     =          60,
 OUTDT   =           0,
 DECAY   =           1,
 KWET    =           1,
 KRAND   =           2,
 CAPEMIN =  -1.000000    ,
 TKEMIN  =  1.0000000E-03,
 WBWR    =  0.0000000E+00,
 WBWF    =  0.0000000E+00,
 WBBH    =  0.0000000E+00
 /
 ------------- End Namelist configuration -------------------
  
  NOTICE   main: pollutant initialization flags
      Dry deposition -  T
      Wet deposition -  T
        Resuspension -  T
         Radioactive -  T
  NOTICE   main: meteorological data flags
     Ten meter wind -  T
     Two meter temp -  T
     Exchange coeff -  F
          Heat flux -  T
      Momentum flux -  T
  Velocity variance -  F
       3D TKE field -  F
     Down shortwave -  T
  Friction velocity -  F
      Friction temp -  F
     Terrain height -  T
   Surface pressure -  T
  Mixed layer depth -  T
  Convective mixing -  F
  Profile averaging -            0
   Stability method -            1
  Horizontal mixing -            0
  PBL mixing scheme -            2
   Free Trop Mixing -    1.000000    
  Precip accumulate -          360
 ------------- Start computation messages -------------------
  NOTICE advrng: (kg ,xyr,xy1) -  1 1440  721    1    1
  NOTICE metpos: (mtime,ftime) -     61370100           0
  NOTICE metpos: (mtime,ftime) -     61370280           0
  NOTICE advpnt: (kg,nx,ny,nz) -     1 1440  721   23
  NOTICE sfcinp: reading ASCDATA.CFG
  NOTICE sfcinp: record terminator = CR/LF
  NOTICE metgrd: (kg, xyr,xy1) -     1 1440  721    1    1
  NOTICE metinp: kwbc 1    1 1440  721    1    1 61370100 16  9  7  3
  NOTICE metinp: kwbc 1  161 1440  721    1    1 61370280 16  9  7  6
  NOTICE advrng: (kg ,xyr,xy1) -  1 1440  721    1    1
  NOTICE   main: Initial time step (min)            1
  NOTICE emspnt: emissions started
  NOTICE   main:            1    61370101           1  1.7499998E+13
  NOTICE   main:            1    61370102           2  3.4961814E+13
  NOTICE   main:            1    61370103           4  3.4923713E+13
  NOTICE   main:            1    61370104           6  3.4923711E+13
  NOTICE   main:            1    61370105           8  3.4923709E+13
  NOTICE   main:            1    61370106          10  3.4923707E+13
  NOTICE   main:            1    61370107          12  3.4923704E+13
  NOTICE   main:            1    61370108          14  3.4923702E+13
  NOTICE   main:            1    61370109          16  3.4923700E+13
  NOTICE   main:            1    61370110          18  3.4923698E+13
  NOTICE   main:            1    61370111          20  3.4923696E+13
  NOTICE   main:            1    61370112          22  3.4923694E+13
  NOTICE   main:            1    61370113          24  3.4923692E+13
  NOTICE   main:            1    61370114          26  3.4923690E+13
  NOTICE   main:            1    61370115          28  3.4923688E+13
  NOTICE   main:            1    61370116          30  3.4923686E+13
  NOTICE   main:            1    61370117          32  3.4923684E+13
  NOTICE   main:            1    61370118          34  3.4923681E+13
  NOTICE   main:            1    61370119          36  3.4923679E+13
  NOTICE   main:            1    61370120          38  3.4923677E+13
  NOTICE   main:            1    61370121          40  3.4923675E+13
  NOTICE   main:            1    61370122          42  3.4923673E+13
  NOTICE   main:            1    61370123          44  3.4923671E+13
  NOTICE   main:            1    61370124          46  3.4923669E+13
  NOTICE   main:            1    61370125          48  3.4899017E+13
  NOTICE   main:            1    61370126          51  3.4861054E+13
  NOTICE   main:            1    61370127          54  3.4823173E+13
  NOTICE   main:            1    61370128          57  3.4823171E+13
  NOTICE   main:            1    61370129          60  3.4823169E+13
  NOTICE   main:            1    61370130          63  3.4823167E+13
  NOTICE   main:            1    61370131          66  3.4823165E+13
  NOTICE   main:            1    61370132          69  3.4823163E+13
  NOTICE   main:            1    61370133          72  3.4823161E+13
  NOTICE   main:            1    61370134          75  3.4823159E+13
  NOTICE   main:            1    61370135          78  3.4823157E+13
  NOTICE   main:            1    61370136          81  3.4823154E+13
  NOTICE   main:            1    61370137          84  3.4823152E+13
  NOTICE   main:            1    61370138          87  3.4823150E+13
  NOTICE   main:            1    61370139          90  3.4823148E+13
  NOTICE   main:            1    61370140          93  3.4798632E+13
  NOTICE   main:            1    61370141          97  3.4760888E+13
  NOTICE   main:            1    61370142         101  3.4723225E+13
  NOTICE   main:            1    61370143         105  3.4685644E+13
  NOTICE   main:            1    61370144         109  3.4685644E+13
  NOTICE   main:            1    61370145         113  3.4685640E+13
  NOTICE   main:            1    61370146         117  3.4685640E+13
  NOTICE   main:            1    61370147         121  3.4685636E+13
  NOTICE   main:            1    61370148         125  3.4685636E+13
  NOTICE   main:            1    61370149         129  3.4685632E+13
  NOTICE   main:            1    61370150         133  3.4685632E+13
  NOTICE   main:            1    61370151         137  3.4685627E+13
  NOTICE   main:            1    61370152         141  3.4685627E+13
  NOTICE   main:            1    61370153         145  3.4685623E+13
  NOTICE   main:            1    61370154         149  3.4685623E+13
  NOTICE   main:            1    61370155         154  3.4685619E+13
  NOTICE   main:            1    61370156         159  3.4685619E+13
  NOTICE   main:            1    61370157         164  3.4685615E+13
  NOTICE   main:            1    61370158         169  3.4685615E+13
  NOTICE   main:            1    61370159         174  3.4685611E+13
  NOTICE   main:            1    61370160         180  3.4685611E+13
  NOTICE advrng: (kg ,xyr,xy1) -  1 1440  721    1    1
  NOTICE   main:            2    61370166          71  8.7185596E+13
  NOTICE   main:            2    61370172          80  1.3968557E+14
  NOTICE   main:            2    61370178          91  1.9218553E+14
  NOTICE   main:            2    61370184         102  2.4468546E+14
  NOTICE   main:            2    61370190         114  2.9688233E+14
  NOTICE   main:            2    61370196         127  3.4907931E+14
  NOTICE   main:            2    61370202         143  4.0157917E+14
  NOTICE   main:            2    61370208         162  4.5407904E+14
  NOTICE   main:            2    61370214         181  5.0455917E+14
  NOTICE   main:            2    61370220         202  5.5569931E+14
  NOTICE advrng: (kg ,xyr,xy1) -  1 1440  721    1    1
  NOTICE   main:            3    61370226         210  6.0819931E+14
  NOTICE   main:            3    61370232         236  6.6025734E+14
  NOTICE   main:            3    61370238         266  7.1275727E+14
  NOTICE   main:            3    61370244         296  7.6495804E+14
  NOTICE   main:            3    61370250         329  8.1634215E+14
  NOTICE   main:            3    61370256         366  8.6884209E+14
  NOTICE   main:            3    61370262         409  9.2001434E+14
  NOTICE   main:            3    61370268         457  9.7140517E+14
  NOTICE   main:            3    61370274         509  1.0229329E+15
  NOTICE condsk:  1 output at    9   7   6   0
  NOTICE   main:            3    61370280         564  1.0748864E+15
  NOTICE advrng: (kg ,xyr,xy1) -  1 1440  721    1    1
  NOTICE metpos: (mtime,ftime) -     61370460    61370100
  NOTICE metinp: kwbc 1  321 1440  721    1    1 61370460 16  9  7  9
  NOTICE   main:            4    61370286         546  0.0000000E+00
  NOTICE   main:            4    61370292         622  0.0000000E+00
  NOTICE   main:            4    61370298         700  0.0000000E+00
  NOTICE   main:            4    61370304         780  0.0000000E+00
  NOTICE   main:            4    61370310         860  0.0000000E+00
  NOTICE   main:            4    61370316         940  0.0000000E+00
  NOTICE   main:            4    61370322        1020  0.0000000E+00
  NOTICE   main:            4    61370328        1100  0.0000000E+00
  NOTICE   main:            4    61370334        1180  0.0000000E+00
  NOTICE   main:            4    61370340        1260  0.0000000E+00
  NOTICE advrng: (kg ,xyr,xy1) -  1 1440  721    1    1
  NOTICE   main:            5    61370346          80  0.0000000E+00
  NOTICE   main:            5    61370352         160  0.0000000E+00
  NOTICE   main:            5    61370358         240  0.0000000E+00
  NOTICE   main:            5    61370364         320  0.0000000E+00
  NOTICE   main:            5    61370370         400  0.0000000E+00
  NOTICE   main:            5    61370376         480  0.0000000E+00
  NOTICE   main:            5    61370382         560  0.0000000E+00
  NOTICE   main:            5    61370388         640  0.0000000E+00
  NOTICE   main:            5    61370394         720  0.0000000E+00
  NOTICE   main:            5    61370400         800  0.0000000E+00
  NOTICE advrng: (kg ,xyr,xy1) -  1 1440  721    1    1
  NOTICE   main:            6    61370406          80  0.0000000E+00
  NOTICE   main:            6    61370412         160  0.0000000E+00
  NOTICE   main:            6    61370418         240  0.0000000E+00
  NOTICE   main:            6    61370424         320  0.0000000E+00
  NOTICE   main:            6    61370430         400  0.0000000E+00
  NOTICE   main:            6    61370436         480  0.0000000E+00
  NOTICE   main:            6    61370442         560  0.0000000E+00
  NOTICE   main:            6    61370448         640  0.0000000E+00
  NOTICE   main:            6    61370454         720  0.0000000E+00
  NOTICE condsk:  1 output at    9   7   9   0
  NOTICE   main:            6    61370460         800  0.0000000E+00
  NOTICE advrng: (kg ,xyr,xy1) -  1 1440  721    1    1
  NOTICE metpos: (mtime,ftime) -     61370640    61370280
  NOTICE metinp: kwbc 1  481 1440  721    1    1 61370640 16  9  7 12
  NOTICE   main:            7    61370472           4  2.0864317E+14
  NOTICE   main:            7    61370484           9  3.6034071E+14
  NOTICE   main:            7    61370496          15  5.1250660E+14
  NOTICE   main:            7    61370508          22  6.6510260E+14
  NOTICE   main:            7    61370520          30  8.1816657E+14
  NOTICE advrng: (kg ,xyr,xy1) -  1 1440  721    1    1
  NOTICE   main:            8    61370532          32  9.6810241E+14
  NOTICE   main:            8    61370544          43  1.1164780E+15
  NOTICE   main:            8    61370556          57  1.2601646E+15
  NOTICE   main:            8    61370568          72  1.4026577E+15
  NOTICE   main:            8    61370580          88  1.5433169E+15
  NOTICE advrng: (kg ,xyr,xy1) -  1 1440  721    1    1
  NOTICE   main:            9    61370590          94  1.6485794E+15
  NOTICE   main:            9    61370600         114  1.7500071E+15
  NOTICE   main:            9    61370610         136  1.8520786E+15
  NOTICE   main:            9    61370620         161  1.9550027E+15
  NOTICE   main:            9    61370630         186  2.0630674E+15
  NOTICE condsk:  1 output at    9   7  12   0
  NOTICE   main:            9    61370640         214  2.1686021E+15
  NOTICE advrng: (kg ,xyr,xy1) -  1 1440  721    1    1
  NOTICE metpos: (mtime,ftime) -     61370820    61370460
  NOTICE metinp: kwbc 1  641 1440  721    1    1 61370820 16  9  7 15
  NOTICE   main:           10    61370650         205  2.3322266E+15
  NOTICE   main:           10    61370660         234  2.4383152E+15
  NOTICE   main:           10    61370670         268  2.5423227E+15
  NOTICE   main:           10    61370680         303  2.6429049E+15
  NOTICE   main:           10    61370690         341  2.7483480E+15
  NOTICE   main:           10    61370700         382  2.8555976E+15
  NOTICE advrng: (kg ,xyr,xy1) -  1 1440  721    1    1
  NOTICE   main:           11    61370710         414  2.9541719E+15
  NOTICE   main:           11    61370720         460  3.0499985E+15
  NOTICE   main:           11    61370730         510  3.1449189E+15
  NOTICE   main:           11    61370740         564  3.2416180E+15
  NOTICE   main:           11    61370750         625  3.3379031E+15
  NOTICE   main:           11    61370760         691  3.4358101E+15
  NOTICE advrng: (kg ,xyr,xy1) -  1 1440  721    1    1
  NOTICE   main:           12    61370766         749  3.4764437E+15
  NOTICE   main:           12    61370772         821  3.5181505E+15
  NOTICE   main:           12    61370778         901  3.5582934E+15
  NOTICE   main:           12    61370784         985  3.5963570E+15
  NOTICE   main:           12    61370790        1073  3.6326586E+15
  NOTICE   main:           12    61370796        1168  3.6694837E+15
  NOTICE   main:           12    61370802        1263  3.7072123E+15
  NOTICE   main:           12    61370808        1362  3.7448518E+15
  NOTICE   main:           12    61370814        1469  3.7833009E+15
  NOTICE condsk:  1 output at    9   7  15   0
  NOTICE   main:           12    61370820        1578  3.8211830E+15
                                                        Index  Height %Mass
                                                            5   430.0  5.30
                                                            4   265.0  9.49
                                                            3   140.0 39.18
                                                            2    55.0 34.39
                                                            1    10.0 11.63
  NOTICE advrng: (kg ,xyr,xy1) -  1 1440  721    1    1
  NOTICE metpos: (mtime,ftime) -     61371000    61370640
  NOTICE metinp: kwbc 1  801 1440  721    1    1 61371000 16  9  7 18
  NOTICE   main:           13    61370826        1560  3.9102418E+15
  NOTICE   main:           13    61370832        1677  3.9463035E+15
  NOTICE   main:           13    61370838        1797  3.9838613E+15
  NOTICE   main:           13    61370844        1918  4.0209487E+15
  NOTICE   main:           13    61370850        2048  4.0583194E+15
  NOTICE   main:           13    61370856        2182  4.0947805E+15
  NOTICE   main:           13    61370862        2318  4.1305484E+15
  NOTICE   main:           13    61370868        2460  4.1647584E+15
  NOTICE   main:           13    61370874        2606  4.2009816E+15
  NOTICE   main:           13    61370880        2756  4.2372880E+15
  NOTICE advrng: (kg ,xyr,xy1) -  1 1440  721    1    1
  NOTICE   main:           14    61370886        2895  4.2744451E+15
  NOTICE   main:           14    61370892        3059  4.3114621E+15
  NOTICE   main:           14    61370898        3227  4.3485814E+15
  NOTICE   main:           14    61370904        3401  4.3841233E+15
  NOTICE   main:           14    61370910        3578  4.4196813E+15
  NOTICE   main:           14    61370916        3757  4.4565850E+15
  NOTICE   main:           14    61370922        3941  4.4927530E+15
  NOTICE   main:           14    61370928        4130  4.5293506E+15
  NOTICE   main:           14    61370934        4324  4.5663303E+15
  NOTICE   main:           14    61370940        4521  4.6028209E+15
  NOTICE advrng: (kg ,xyr,xy1) -  1 1440  721    1    1
  NOTICE   main:           15    61370945        4701  4.6320396E+15
  NOTICE   main:           15    61370950        4908  4.6605791E+15
  NOTICE   main:           15    61370955        5120  4.6909692E+15
  NOTICE   main:           15    61370960        5334  4.7224314E+15
  NOTICE   main:           15    61370965        5550  4.7528419E+15
  NOTICE   main:           15    61370970        5769  4.7833657E+15
  NOTICE   main:           15    61370975        5996  4.8133661E+15
  NOTICE   main:           15    61370980        6229  4.8437310E+15
  NOTICE   main:           15    61370985        6467  4.8741066E+15
  NOTICE   main:           15    61370990        6708  4.9046035E+15
  NOTICE   main:           15    61370995        6959  4.9348702E+15
  NOTICE condsk:  1 output at    9   7  18   0
  NOTICE   main:           15    61371000        7217  4.9657907E+15
  NOTICE advrng: (kg ,xyr,xy1) -  1 1440  721    1    1
  NOTICE metpos: (mtime,ftime) -     61371180    61370820
  NOTICE metinp: kwbc 1  961 1440  721    1    1 61371180 16  9  7 21
  NOTICE   main:           16    61371004        7195  5.0252282E+15
  NOTICE   main:           16    61371008        7436  5.0489161E+15
  NOTICE   main:           16    61371012        7690  5.0718807E+15
  NOTICE   main:           16    61371016        7953  5.0953924E+15
  NOTICE   main:           16    61371020        8223  5.1186916E+15
  NOTICE   main:           16    61371024        8503  5.1411150E+15
  NOTICE   main:           16    61371028        8791  5.1626634E+15
  NOTICE   main:           16    61371032        9087  5.1835960E+15
  NOTICE   main:           16    61371036        9389  5.2044175E+15
  NOTICE   main:           16    61371040        9704  5.2248997E+15
 WARNING depsus: exceeding puff limit
  NOTICE   main:           16    61371044       10000  5.2462956E+15
 WARNING emspnt: exceeding puff limit
 WARNING depsus: exceeding puff limit
  NOTICE   main:           16    61371048       10000  5.2310474E+15
 WARNING emspnt: exceeding puff limit
 WARNING depsus: exceeding puff limit
  NOTICE   main:           16    61371052       10000  5.2177002E+15
 WARNING emspnt: exceeding puff limit
 WARNING depsus: exceeding puff limit
 ...
  WARNING emspnt: exceeding puff limit
 WARNING depsus: exceeding puff limit
  NOTICE   main:          165    61380000       10000  4.4689068E+14
 WARNING particle array limit reached; expand MAXPAR
and here is default_conc file

Code: Select all

16 09 07 03
1
35.41 126.42 50.0
165
0
10000.0
1
D:/과제/01_방사능/04.소프트웨어/방사능과제_건기원/api2arl_window/api2arl_Win_OS_051719/api2arl_Win_OS_051719/
160907_7days_incl_rain
1
C137
1.05E+15
165
16 09 07 03 00
1
35.55 127.08
0.02 0.02
3.3 3.04
./
Hanbit160907
1
0
16 09 07 03 00
16 09 13 00 00
00 03 00
1
0.4 1.9 1.0
0.002 136.96 0.0 0.0 0.0
0.0 3.2E-05 5.0E+05
10960.0
1.0E-06
Please let me know why i could not get concentration value which is not 0.
Also please let me know why i only got contour image like attachment.
And please let me know what should i change to get a colorful dispersion image and meaningful dispersion model results.


Thank you all advanced.
Attachments
hysplit contour image.gif

sonny.zinn
Posts: 20
Joined: May 8th, 2019, 1:31 pm
Registered HYSPLIT User: Yes

Re: Concentration matter and Contour display matter.

Post by sonny.zinn » February 6th, 2020, 7:30 am

Could you attach the SETUP.CFG and CONTROL files as well? Together with MESSAGE.TXT that you already posted they provide useful information for a diagnosis.

sol
Posts: 6
Joined: February 5th, 2020, 9:24 pm
Registered HYSPLIT User: No

Re: Concentration matter and Contour display matter.

Post by sol » February 7th, 2020, 9:21 am

Thank you for your reply.
I attached whole message file, setup file and control file.

Thank you advanced.
Attachments
MESSAGE.txt
(269.53 KiB) Downloaded 2 times
SETUP.txt
(561 Bytes) Downloaded 4 times
CONTROL.txt
(415 Bytes) Downloaded 3 times

Post Reply