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: 27
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 27 times
SETUP.txt
(561 Bytes) Downloaded 28 times
CONTROL.txt
(415 Bytes) Downloaded 27 times

Post Reply