close
프로필 배경
프로필 로고

해솔 | 해처럼 밝고 소나무처럼 바르게

GrADS · 2019. 9. 3. fullscreen 넓게보기

[GrADS, ShellScript] 그라즈 아스키 (ASCII) 형식인 복사전달모델 (RTM) 결과를 이용한 가시화

목차

 정보

  • 업무명    : 그라즈 아스키 (ASCII) 형식인 복사전달모델 (RTM) 결과를 이용한 가시화

  • 작성자    : 이상호

  • 작성일    : 2019-09-03

  • 설   명     :

  • 수정이력 :

 

 내용

[특징]

  • 복사전달모델 결과를 이해하기 위해 가시화 도구가 필요하며 이 프로그램은 이러한 목적을 달성하기 위해 고안된 소프트웨어

 

[기능]

  • 아스키 형식인 입력 자료를 이용하여 매핑

  • 컬러바 생성

 

[활용 자료]

  • 자료 : 복사전달모델 결과

  • 기간 : 2015년 04월 01일 - 2015년 04월 01일

  • 해상도 : 매 6시간

 

[자료 처리 방안 및 활용 분석 기법]

  • 없음

 

[사용법]

  • GrADS 가시화에 필요한 라이브러리를 동일 디렉터리 위치

  • 소스 코드를 실행 (csh Visualization_Using_ASCII_Format_of_Radiation_Transfer_Model_Result.csh)

  • 가시화 결과를 확인

 

[사용 OS]

  • Linux

 

[사용 언어]

  • GrADS

  • ShellScript (csh)

 

 소스 코드

[명세]

  • 매핑을 위한 라이브러리

function main(args)

fname1=subwrd(args,1)
latco=subwrd(args,2)
lonco=subwrd(args,3)
valco=subwrd(args,4)
colco=subwrd(args,5)
* xi=subwrd(args,5)
* xd=subwrd(args,6)

color='14 4 13 3 10 7 12 8 2 6'
palgrey()

'q dims'
rec2=sublin(result,2)
rec3 =sublin(result,3)
xlon0=subwrd(rec2,6)
xlon1=subwrd(rec2,8)
ylat0=subwrd(rec3,6)
ylat1=subwrd(rec3,8)

ret=read(fname1)

'q ll2xy 'xlon0' 'ylat0
x0=subwrd(result,1)
y0=subwrd(result,2)
'q ll2xy  'xlon1' 'ylat1
x1=subwrd(result,1)
y1=subwrd(result,2)
xs=x1-x0
ys=y1-y0

dx2=xs*0.05/(xlon1-xlon0) 
dy2=ys*0.05/(ylat1-ylat0)

while (1)
  ret = read(fname1)
  rc = sublin(ret,1)
  if (rc>0) 
    if (rc!=2) 
      say 'File I/O Error'
      return
    endif
    break
  endif
  rec = sublin(ret,2)
  time= subwrd(rec,valco)
  ylat= subwrd(rec,latco)
  xlon= subwrd(rec,lonco)
  xco=subwrd(rec,colco)
  if (xlon < 0);xlon=xlon+360;endif;
  
  if (xlon >= xlon0 & xlon <= xlon1 & ylat >= ylat0 & ylat <= ylat1)
     x=x0+xs*(xlon-xlon0)/(xlon1-xlon0)
     y=y0+ys*(ylat-ylat0)/(ylat1-ylat0)

if (time < 0.000 ); xcol= 16 ;endif;

if (time >=      0.000 & time <     10.000 ); xcol= 17 ;endif;
if (time >=     10.000 & time <     20.000 ); xcol= 18 ;endif;
if (time >=     20.000 & time <     30.000 ); xcol= 19 ;endif;
if (time >=     30.000 & time <     40.000 ); xcol= 20 ;endif;
if (time >=     40.000 & time <     50.000 ); xcol= 21 ;endif;
if (time >=     50.000 & time <     60.000 ); xcol= 22 ;endif;
if (time >=     60.000 & time <     70.000 ); xcol= 23 ;endif;
if (time >=     70.000 & time <     80.000 ); xcol= 24 ;endif;
if (time >=     80.000 & time <     90.000 ); xcol= 25 ;endif;
if (time >=     90.000 & time <    100.000 ); xcol= 26 ;endif;
if (time >=    100.000 & time <    110.000 ); xcol= 27 ;endif;
if (time >=    110.000 & time <    120.000 ); xcol= 28 ;endif;
if (time >=    120.000 & time <    130.000 ); xcol= 29 ;endif;
if (time >=    130.000 & time <    140.000 ); xcol= 30 ;endif;
if (time >=    140.000 & time <    150.000 ); xcol= 31 ;endif;
if (time >=    150.000 & time <    160.000 ); xcol= 32 ;endif;
if (time >=    160.000 & time <    170.000 ); xcol= 33 ;endif;
if (time >=    170.000 & time <    180.000 ); xcol= 34 ;endif;
if (time >=    180.000 & time <    190.000 ); xcol= 35 ;endif;
if (time >=    190.000 & time <    200.000 ); xcol= 36 ;endif;
if (time >=    200.000 & time <    210.000 ); xcol= 37 ;endif;
if (time >=    210.000 & time <    220.000 ); xcol= 38 ;endif;
if (time >=    220.000 & time <    230.000 ); xcol= 39 ;endif;
if (time >=    230.000 & time <    240.000 ); xcol= 40 ;endif;
if (time >=    240.000 & time <    250.000 ); xcol= 41 ;endif;
if (time >=    250.000 & time <    260.000 ); xcol= 42 ;endif;
if (time >=    260.000 & time <    270.000 ); xcol= 43 ;endif;
if (time >=    270.000 & time <    280.000 ); xcol= 44 ;endif;
if (time >=    280.000 & time <    290.000 ); xcol= 45 ;endif;
if (time >=    290.000 & time <    300.000 ); xcol= 46 ;endif;
if (time >=    300.000 & time <    310.000 ); xcol= 47 ;endif;
if (time >=    310.000 & time <    320.000 ); xcol= 48 ;endif;
if (time >=    320.000 & time <    330.000 ); xcol= 49 ;endif;
if (time >=    330.000 & time <    340.000 ); xcol= 50 ;endif;
if (time >=    340.000 & time <    350.000 ); xcol= 51 ;endif;
if (time >=    350.000 & time <    360.000 ); xcol= 52 ;endif;
if (time >=    360.000 & time <    370.000 ); xcol= 53 ;endif;
if (time >=    370.000 & time <    380.000 ); xcol= 54 ;endif;
if (time >=    380.000 & time <    390.000 ); xcol= 55 ;endif;
if (time >=    390.000 & time <    400.000 ); xcol= 56 ;endif;
if (time >=    400.000 & time <    410.000 ); xcol= 57 ;endif;
if (time >=    410.000 & time <    420.000 ); xcol= 58 ;endif;
if (time >=    420.000 & time <    430.000 ); xcol= 59 ;endif;
if (time >=    430.000 & time <    440.000 ); xcol= 60 ;endif;
if (time >=    440.000 & time <    450.000 ); xcol= 61 ;endif;
if (time >=    450.000 & time <    460.000 ); xcol= 62 ;endif;
if (time >=    460.000 & time <    470.000 ); xcol= 63 ;endif;
if (time >=    470.000 & time <    480.000 ); xcol= 64 ;endif;
if (time >=    480.000 & time <    490.000 ); xcol= 65 ;endif;
if (time >=    490.000 & time <    500.000 ); xcol= 66 ;endif;
if (time >=    500.000 & time <    510.000 ); xcol= 67 ;endif;
if (time >=    510.000 & time <    520.000 ); xcol= 68 ;endif;
if (time >=    520.000 & time <    530.000 ); xcol= 69 ;endif;
if (time >=    530.000 & time <    540.000 ); xcol= 70 ;endif;
if (time >=    540.000 & time <    550.000 ); xcol= 71 ;endif;
if (time >=    550.000 & time <    560.000 ); xcol= 72 ;endif;
if (time >=    560.000 & time <    570.000 ); xcol= 73 ;endif;
if (time >=    570.000 & time <    580.000 ); xcol= 74 ;endif;
if (time >=    580.000 & time <    590.000 ); xcol= 75 ;endif;
if (time >=    590.000 & time <    600.000 ); xcol= 76 ;endif;
if (time >=    600.000 & time <    610.000 ); xcol= 77 ;endif;
if (time >=    610.000 & time <    620.000 ); xcol= 78 ;endif;
if (time >=    620.000 & time <    630.000 ); xcol= 79 ;endif;
if (time >=    630.000 & time <    640.000 ); xcol= 80 ;endif;
if (time >=    640.000 & time <    650.000 ); xcol= 81 ;endif;
if (time >=    650.000 & time <    660.000 ); xcol= 82 ;endif;
if (time >=    660.000 & time <    670.000 ); xcol= 83 ;endif;
if (time >=    670.000 & time <    680.000 ); xcol= 84 ;endif;
if (time >=    680.000 & time <    690.000 ); xcol= 85 ;endif;
if (time >=    690.000 & time <    700.000 ); xcol= 86 ;endif;
if (time >=    700.000 & time <    710.000 ); xcol= 87 ;endif;
if (time >=    710.000 & time <    720.000 ); xcol= 88 ;endif;
if (time >=    720.000 & time <    730.000 ); xcol= 89 ;endif;
if (time >=    730.000 & time <    740.000 ); xcol= 90 ;endif;
if (time >=    740.000 & time <    750.000 ); xcol= 91 ;endif;
if (time >=    750.000 & time <    760.000 ); xcol= 92 ;endif;
if (time >=    760.000 & time <    770.000 ); xcol= 93 ;endif;
if (time >=    770.000 & time <    780.000 ); xcol= 94 ;endif;
if (time >=    780.000 & time <    790.000 ); xcol= 95 ;endif;
if (time >=    790.000 & time <    800.000 ); xcol= 96 ;endif;
if (time >=    800.000 & time <    810.000 ); xcol= 97 ;endif;
if (time >=    810.000 & time <    820.000 ); xcol= 98 ;endif;
if (time >=    820.000 & time <    830.000 ); xcol= 99 ;endif;
if (time >=    830.000 & time <    840.000 ); xcol=100 ;endif;
if (time >=    840.000 & time <    850.000 ); xcol=101 ;endif;
if (time >=    850.000 & time <    860.000 ); xcol=102 ;endif;
if (time >=    860.000 & time <    870.000 ); xcol=103 ;endif;
if (time >=    870.000 & time <    880.000 ); xcol=104 ;endif;
if (time >=    880.000 & time <    890.000 ); xcol=105 ;endif;
if (time >=    890.000 & time <    900.000 ); xcol=106 ;endif;

if (time >= 900.0 ); xcol= 199 ;endif;


**    xcl=color2(time,xi,xd)
**    say xcl
**    xcl=subwrd(color,xcol)
*
*
     'set line 'xcol
*    'set line 'xco

    'draw recf 'x-dx2' 'y-dy2' 'x+dx2' 'y+dy2
*    'draw recf 'x'-'dx2' 'y'-'dy2' 'x'+'dx2' 'y'+'dy2
*    'draw mark 3 'x' 'y' 0.1'
 endif
endwhile

return


function palgrey()

***************** Rainbow ********************
* 13 color
*
*'set rgb 16 64 64 64'
*
*'set rgb 17 240   0 130'
*'set rgb 18 250  60  60'
*'set rgb 19 240 130  40'
*'set rgb 20 230 175  45'
*'set rgb 21 230 220  50'
*'set rgb 22 160 230  50'
*'set rgb 23   0 220   0'
*'set rgb 24   0 210 140'
*'set rgb 25   0 200 200'
*'set rgb 26   0 160 255'
*'set rgb 27  30  60 255'
*'set rgb 28 130   0 220'
*'set rgb 29 160   0 200'
*
*'set rgb 99 255 255 255'

***************** KMA Radar ********************
* color 30
*
*'set rgb  16  64  64  64'
*
*'set rgb  17 135 217 255'
*'set rgb  18  62 193 255'
*'set rgb  19   7 171 255'
*'set rgb  20   0 141 255'
*'set rgb  21   0 119 255'
*'set rgb  22 105 252 105'
*'set rgb  23  30 242 105'
*'set rgb  24   0 213   0'
*'set rgb  25   0 164   0'
*'set rgb  26   0 128   0'
*'set rgb  27 255 241 111'
*'set rgb  28 255 226  86'
*'set rgb  29 255 208  57'
*'set rgb  30 255 188  30'
*'set rgb  31 255 170   9'
*'set rgb  32 255 156   0'
*'set rgb  33 255 139  27'
*'set rgb  34 255 128  81'
*'set rgb  35 255 110 110'
*'set rgb  36 246  94 103'
*'set rgb  37 232  74  86'
*'set rgb  38 217  52  62'
*'set rgb  39 201  31  37'
*'set rgb  40 188  13  15'
*'set rgb  41 179   0   0'
*'set rgb  42 197 125 255'
*'set rgb  43 194  62 255'
*'set rgb  44 173   7 255'
*'set rgb  45 146   0 255'
*'set rgb  46 127   0 191'
*
*'set rgb 99 255 255 255'


************* Rainbow *************************
* 39 color

*'set rgb 16  64  64  64'
*
*'set rgb 17 153 102 255'
*'set rgb 18 131  87 255'
*'set rgb 19 109  73 255'
*'set rgb 20  87  58 255'
*'set rgb 21  66  44 255'
*'set rgb 22  44  29 255'
*'set rgb 23   0   0 255'
*'set rgb 24   0  72 255'
*'set rgb 25   0 117 255'
*'set rgb 26   0 156 255'
*'set rgb 27   0 192 255'
*'set rgb 28   0 213 255'
*'set rgb 29   0 232 255'
*'set rgb 30   0 246 255'
*'set rgb 31   0 255 249'
*'set rgb 32   0 255 213'
*'set rgb 33   0 255 193'
*'set rgb 34   0 255 159'
*'set rgb 35   0 255 126'
*'set rgb 36   0 255  77'
*'set rgb 37   0 255  15'
*'set rgb 38  50 255   0'
*'set rgb 39 114 255   0'
*'set rgb 40 152 255   0'
*'set rgb 41 181 255   0'
*'set rgb 42 227 255   0'
*'set rgb 43 243 255   0'
*'set rgb 44 255 253   0'
*'set rgb 45 255 232   0'
*'set rgb 46 255 211   0'
*'set rgb 47 255 189   0'
*'set rgb 48 255 168   0'
*'set rgb 49 255 146   0'
*'set rgb 50 255 124   0'
*'set rgb 51 255 102   0'
*'set rgb 52 255  87   0'
*'set rgb 53 255  73   0'
*'set rgb 54 255  58   0'
*'set rgb 55 255  29   0'
*'set rgb 56 255   0   0'
*
*'set rgb 99 255 255 255'
*
************* Diff ****************************
* 37 color
*
*'set rgb 16   0   0   0'
*
*'set rgb 17   0   0 255'
*'set rgb 18  38  38 255'
*'set rgb 19  51  51 255'
*'set rgb 20  63  63 255'
*'set rgb 21  76  76 255'
*'set rgb 22  89  89 255'
*'set rgb 23 102 102 255'
*'set rgb 24 114 114 255'
*'set rgb 25 127 127 255'
*'set rgb 26 140 140 255'
*'set rgb 27 153 153 255'
*'set rgb 28 165 165 255'
*'set rgb 29 178 178 255'
*'set rgb 30 191 191 255'
*'set rgb 31 204 204 255'
*'set rgb 32 216 216 255'
*'set rgb 33 229 229 255'
*'set rgb 34 242 242 255'
*'set rgb 35 255 255 255'
*'set rgb 36 255 243 243'
*'set rgb 37 255 230 230'
*'set rgb 38 255 217 217'
*'set rgb 39 255 205 205'
*'set rgb 40 255 192 192'
*'set rgb 41 255 179 179'
*'set rgb 42 255 166 166'
*'set rgb 43 255 154 154'
*'set rgb 44 255 141 141'
*'set rgb 45 255 128 128'
*'set rgb 46 255 115 115'
*'set rgb 47 255 102 102'
*'set rgb 48 255  90  90'
*'set rgb 49 255  77  77'
*'set rgb 50 255  64  64'
*'set rgb 51 255  51  51'
*'set rgb 52 255  38  38'
*'set rgb 53 255   0   0'
*
*'set rgb 99 255 255 255'

**************  Forest  ****************
* 65 color
*
*'set rgb 16  64  64  64'
*
*'set rgb 17   0   0  30'
*'set rgb 18   5  26   2'
*'set rgb 19  10  32   4'
*'set rgb 20  15  38   6'
*'set rgb 21  20  44   8'
*'set rgb 22  25  50  10'
*'set rgb 23  30  56  12'
*'set rgb 24  35  62  14'
*'set rgb 25  40  68  16'
*'set rgb 26  45  74  18'
*'set rgb 27  50  80  20'
*'set rgb 28  55  86  22'
*'set rgb 29  60  92  24'
*'set rgb 30  65  98  26'
*'set rgb 31  70 104  28'
*'set rgb 32  75 110  30'
*'set rgb 33  80 116  32'
*'set rgb 34  85 122  34'
*'set rgb 35  90 128  36'
*'set rgb 36  95 134  38'
*'set rgb 37 100 140  40'
*'set rgb 38 105 146  42'
*'set rgb 39 110 152  44'
*'set rgb 40 115 158  46'
*'set rgb 41 120 164  48'
*'set rgb 42 125 170  50'
*'set rgb 43 130 176  52'
*'set rgb 44 135 182  54'
*'set rgb 45 140 188  56'
*'set rgb 46 145 194  58'
*'set rgb 47 150 200  60'
*'set rgb 48 155 206  62'
*'set rgb 49 255 250 110'
*'set rgb 50 249 238  98'
*'set rgb 51 243 226  86'
*'set rgb 52 237 214  74'
*'set rgb 53 231 202  62'
*'set rgb 54 225 190  50'
*'set rgb 55 219 178  38'
*'set rgb 56 213 166  26'
*'set rgb 57 207 154  14'
*'set rgb 58 201 142   2'
*'set rgb 59 195 130   0'
*'set rgb 60 189 118   0'
*'set rgb 61 183 106   0'
*'set rgb 62 177  94   0'
*'set rgb 63 171  82   0'
*'set rgb 64 165  70   0'
*'set rgb 65 159  58   0'
*'set rgb 66 153  46   0'
*'set rgb 67 147  34   0'
*'set rgb 68 141  22   0'
*'set rgb 69 135  10   0'
*'set rgb 70 129   0   0'
*'set rgb 71 123   0   0'
*'set rgb 72 117   0   0'
*'set rgb 73 111   0   0'
*'set rgb 74 105   0   0'
*'set rgb 75  99   0   0'
*'set rgb 76  93   0   0'
*'set rgb 77  87   0   0'
*'set rgb 78  81   0   0'
*'set rgb 79  75   0   0'
*'set rgb 80  69   0   0'
*'set rgb 81  20   0   0'
*
*'set rgb 99 255 255 255'
*
**************  Rainbow [Red ~ Blue] ****************
* color 100 

'set rgb  16  64  64  64'

'set rgb  17   0   0 130'
'set rgb  18   0   0 134'
'set rgb  19   0   0 146'
'set rgb  20   0   0 154'
'set rgb  21   0   0 166'
'set rgb  22   0   0 174'
'set rgb  23   0   0 186'
'set rgb  24   0   0 198'
'set rgb  25   0   0 206'
'set rgb  26   0   0 218'
'set rgb  27   0   0 226'
'set rgb  28   0   0 238'
'set rgb  29   0   0 250'
'set rgb  30   0   0 254'
'set rgb  31   0  11 254'
'set rgb  32   0  19 254'
'set rgb  33   0  31 254'
'set rgb  34   0  39 254'
'set rgb  35   0  51 254'
'set rgb  36   0  63 254'
'set rgb  37   0  71 254'
'set rgb  38   0  83 254'
'set rgb  39   0  91 254'
'set rgb  40   0 103 254'
'set rgb  41   0 115 254'
'set rgb  42   0 123 254'
'set rgb  43   0 134 254'
'set rgb  44   0 142 254'
'set rgb  45   0 154 254'
'set rgb  46   0 162 254'
'set rgb  47   0 174 254'
'set rgb  48   0 186 254'
'set rgb  49   0 194 254'
'set rgb  50   0 206 254'
'set rgb  51   0 214 254'
'set rgb  52   0 226 254'
'set rgb  53   0 238 254'
'set rgb  54   0 246 254'
'set rgb  55   0 254 254'
'set rgb  56   7 254 246'
'set rgb  57  19 254 234'
'set rgb  58  31 254 222'
'set rgb  59  39 254 214'
'set rgb  60  51 254 202'
'set rgb  61  59 254 194'
'set rgb  62  71 254 181'
'set rgb  63  79 254 174'
'set rgb  64  91 254 162'
'set rgb  65 103 254 150'
'set rgb  66 111 254 142'
'set rgb  67 123 254 130'
'set rgb  68 130 254 123'
'set rgb  69 142 254 111'
'set rgb  70 154 254  99'
'set rgb  71 162 254  91'
'set rgb  72 174 254  79'
'set rgb  73 182 254  71'
'set rgb  74 194 254  59'
'set rgb  75 202 254  51'
'set rgb  76 214 254  39'
'set rgb  77 226 254  27'
'set rgb  78 234 254  19'
'set rgb  79 246 254   7'
'set rgb  80 254 254   0'
'set rgb  81 254 242   0'
'set rgb  82 254 230   0'
'set rgb  83 254 222   0'
'set rgb  84 254 210   0'
'set rgb  85 254 202   0'
'set rgb  86 254 190   0'
'set rgb  87 254 178   0'
'set rgb  88 254 170   0'
'set rgb  89 254 158   0'
'set rgb  90 254 150   0'
'set rgb  91 254 138   0'
'set rgb  92 254 130   0'
'set rgb  93 254 119   0'
'set rgb  94 254 107   0'
'set rgb  95 254  99   0'
'set rgb  96 254  87   0'
'set rgb  97 254  79   0'
'set rgb  98 254  67   0'
'set rgb  99 254  55   0'
'set rgb 100 254  47   0'
'set rgb 101 254  35   0'
'set rgb 102 254  27   0'
'set rgb 103 254  15   0'
'set rgb 104 254   7   0'
'set rgb 105 249   0   0'
'set rgb 106 236   0   0'
'set rgb 107 227   0   0'
'set rgb 108 214   0   0'
'set rgb 109 205   0   0'
'set rgb 110 192   0   0'
'set rgb 111 179   0   0'
'set rgb 112 170   0   0'
'set rgb 113 157   0   0'
'set rgb 114 148   0   0'
'set rgb 115 135   0   0'
'set rgb 116 130   0   0'

'set rgb 199   0   0   0'

return

 

  • 컬러바 생성을 위한 라이브러리

function colorbar (args)

vert=subwrd(args,1)
xmid=subwrd(args,2)
ymid=subwrd(args,3)
barw=subwrd(args,4)
sfc=subwrd(args,5)
xsiz=subwrd(args,6)
ysiz=subwrd(args,7)
col1=subwrd(args,8)
col2=subwrd(args,9)
val1=subwrd(args,10)
val2=subwrd(args,11)
cl1=subwrd(args,12)
cl2=subwrd(args,13)
cli=subwrd(args,14)

dum2=cli
if (cli<0) ;dum2=cli*-1;endif
dum1=cl2-cl1
if (cl1>cl2);dum1=cl1-cl2;endif
blnum=dum1/dum2+1

* col1 & col2 : should be determined when dot-alb.f
*               (min and max # of color scale) 
* val1 & val2 : actual maginitude corresponding colors, col1 & col2
*   if you set alb=0 to col#=40, and set alb=1.0 to col#=85, then
*      val1=0, val2=1, col1=40, col2=85
* cl1=val1, cl2=val2, I will change later, 
* cli=interval of # under colorbar

cnum=col2-col1+1

  stroff=0.05*sfc
  strxsiz=0.12*sfc
  strysiz=0.13*sfc

*
*	logic for setting the bar orientation with user overides
*
    if(vert = 0) ; vchk = 0 ; endif
    if(vert = 1) ; vchk = 1 ; endif
*
*	vertical bar
*

  if (vchk = 1 )

    ywid = ysiz / cnum
    xl = xmid - barw/2
    xr = xmid + barw/2
    yb = ymid-ysiz/2
    y1 = yb
    vert = 1

  else

*
*	horizontal bar
*

    xwid = xsiz / cnum
    yt = ymid + barw/2
    yb = ymid - barw/2
    xl = xmid-xsiz/2
    x1 = xl
    'set string 1 tc 5'
    vert = 0
  endif

*
*  Plot colorbar
*
  num = 1
  while (num<=cnum) 
    col = col1-1+num
    if (vert) 
      yt = yb + ywid
    else 
      xr = xl + xwid
    endif

    'set line 'col
    'draw recf 'xl' 'yb' 'xr' 'yt

    num = num + 1
    if (vert); yb = yt;
    else; xl = xr; endif;
  endwhile

 bli = 1

  'set strsiz 'strxsiz' 'strysiz
      if (vert)
        'set string 1 l 5'
         cl = cl1
         yp=y1 + ywid
         xp=xr+stroff
         yd = (ysiz - 2*ywid)/dum1*dum2
        while (bli <= blnum)
        'draw string 'xp' 'yp' 'cl
         cl = cl + cli
         yp = yp + yd
         bli=bli+1
        endwhile
      else
        'set string 1 tc 5'
         xp=x1 + xwid
         cl = cl1
         yp=yb-stroff
         xd = (xsiz - 2*xwid)/dum1*dum2
       while (bli <= blnum)
        'draw string 'xp' 'yp' 'cl
         cl = cl + cli
         xp = xp + xd
         bli=bli+1
        endwhile
      endif

return

 

  • 가시화

#!/bin/csh

#rm -rf FIG/*

set a = `ls DATA/2013??????`

foreach fn ($a)
set fn2 = `echo $fn | cut -c 6-15`

set yy = `echo $fn | cut -c 6-9`
set mm = `echo $fn | cut -c 10-11`
set dd = `echo $fn | cut -c 12-13`
set tt = `echo $fn | cut -c 14-15`
set fn3 = `echo ${yy}.${mm}.${dd} ${tt}UTC`

paste DATA/Geo.dat DATA/$fn2 > fort.1

cat >! 2.exec << EOF
reinit
open GrADS/Map.ctl
enable print 2.gx 
set font 5
set vpage 0  8.5 0 11
set parea 0 8.5 2.0 10.0
set mpdraw off
set xlopts 1 5.5 0.15
set ylopts 1 5.5 0.15
set lat 31. 45.
set lon 121. 135.5
set cmin 9999
set cint 2
set xlint 2.0
set ylint 2.0
set grid on
set grads off
draw title  $fn3 [Clear]
set strsiz 0.15 0.15
draw string 3.5 1.4 Surface Solar Radiation [Wm\`a-2\`n]
d t

run GrADS/Library_For_Mapping_ASCII_Format_File.gs fort.1  1 2 3 3
run GrADS/Library_For_Creating_Colorbar_From_ASCII_Format_File.gs   0   4.25 1.0   0.28 1.1   8.0 1.0  17 116   0 900 0 900 100
*run GrADS/Library_For_Creating_Colorbar_From_ASCII_Format_File.gs   1   8.0 6.0    0.25 1.0   1.0 8.0  16 116   0 900 0 900 50

set rgb 116 0 0 0
draw shp GrADS/Map/GSHHG/GSHHS_shp/f/GSHHS_f_L1.shp
*draw shp GrADS/Map/KOR_adm/KOR_adm0.shp
*draw shp GrADS/Map/CHN_adm/CHN_adm0.shp
*draw shp GrADS/Map/PRK_adm/PRK_adm0.shp
*draw shp GrADS/Map/JPN_adm/JPN_adm0.shp

print
disable print
!gxps -c -i 2.gx -o 2.ps
quit
EOF

grads -bpc 'exec 2.exec'
convert 2.ps 2.gif
#display  2.gif 

mv -f 2.gif FIG/${fn2}_SFC.gif
end

rm -f 2.ps  2.gx

 

[GitHub Gist]

  • 매핑을 위한 라이브러리

 

function main(args)
fname1=subwrd(args,1)
latco=subwrd(args,2)
lonco=subwrd(args,3)
valco=subwrd(args,4)
colco=subwrd(args,5)
* xi=subwrd(args,5)
* xd=subwrd(args,6)
color='14 4 13 3 10 7 12 8 2 6'
palgrey()
'q dims'
rec2=sublin(result,2)
rec3 =sublin(result,3)
xlon0=subwrd(rec2,6)
xlon1=subwrd(rec2,8)
ylat0=subwrd(rec3,6)
ylat1=subwrd(rec3,8)
ret=read(fname1)
'q ll2xy 'xlon0' 'ylat0
x0=subwrd(result,1)
y0=subwrd(result,2)
'q ll2xy 'xlon1' 'ylat1
x1=subwrd(result,1)
y1=subwrd(result,2)
xs=x1-x0
ys=y1-y0
dx2=xs*0.05/(xlon1-xlon0)
dy2=ys*0.05/(ylat1-ylat0)
while (1)
ret = read(fname1)
rc = sublin(ret,1)
if (rc>0)
if (rc!=2)
say 'File I/O Error'
return
endif
break
endif
rec = sublin(ret,2)
time= subwrd(rec,valco)
ylat= subwrd(rec,latco)
xlon= subwrd(rec,lonco)
xco=subwrd(rec,colco)
if (xlon < 0);xlon=xlon+360;endif;
if (xlon >= xlon0 & xlon <= xlon1 & ylat >= ylat0 & ylat <= ylat1)
x=x0+xs*(xlon-xlon0)/(xlon1-xlon0)
y=y0+ys*(ylat-ylat0)/(ylat1-ylat0)
if (time < 0.000 ); xcol= 16 ;endif;
if (time >= 0.000 & time < 10.000 ); xcol= 17 ;endif;
if (time >= 10.000 & time < 20.000 ); xcol= 18 ;endif;
if (time >= 20.000 & time < 30.000 ); xcol= 19 ;endif;
if (time >= 30.000 & time < 40.000 ); xcol= 20 ;endif;
if (time >= 40.000 & time < 50.000 ); xcol= 21 ;endif;
if (time >= 50.000 & time < 60.000 ); xcol= 22 ;endif;
if (time >= 60.000 & time < 70.000 ); xcol= 23 ;endif;
if (time >= 70.000 & time < 80.000 ); xcol= 24 ;endif;
if (time >= 80.000 & time < 90.000 ); xcol= 25 ;endif;
if (time >= 90.000 & time < 100.000 ); xcol= 26 ;endif;
if (time >= 100.000 & time < 110.000 ); xcol= 27 ;endif;
if (time >= 110.000 & time < 120.000 ); xcol= 28 ;endif;
if (time >= 120.000 & time < 130.000 ); xcol= 29 ;endif;
if (time >= 130.000 & time < 140.000 ); xcol= 30 ;endif;
if (time >= 140.000 & time < 150.000 ); xcol= 31 ;endif;
if (time >= 150.000 & time < 160.000 ); xcol= 32 ;endif;
if (time >= 160.000 & time < 170.000 ); xcol= 33 ;endif;
if (time >= 170.000 & time < 180.000 ); xcol= 34 ;endif;
if (time >= 180.000 & time < 190.000 ); xcol= 35 ;endif;
if (time >= 190.000 & time < 200.000 ); xcol= 36 ;endif;
if (time >= 200.000 & time < 210.000 ); xcol= 37 ;endif;
if (time >= 210.000 & time < 220.000 ); xcol= 38 ;endif;
if (time >= 220.000 & time < 230.000 ); xcol= 39 ;endif;
if (time >= 230.000 & time < 240.000 ); xcol= 40 ;endif;
if (time >= 240.000 & time < 250.000 ); xcol= 41 ;endif;
if (time >= 250.000 & time < 260.000 ); xcol= 42 ;endif;
if (time >= 260.000 & time < 270.000 ); xcol= 43 ;endif;
if (time >= 270.000 & time < 280.000 ); xcol= 44 ;endif;
if (time >= 280.000 & time < 290.000 ); xcol= 45 ;endif;
if (time >= 290.000 & time < 300.000 ); xcol= 46 ;endif;
if (time >= 300.000 & time < 310.000 ); xcol= 47 ;endif;
if (time >= 310.000 & time < 320.000 ); xcol= 48 ;endif;
if (time >= 320.000 & time < 330.000 ); xcol= 49 ;endif;
if (time >= 330.000 & time < 340.000 ); xcol= 50 ;endif;
if (time >= 340.000 & time < 350.000 ); xcol= 51 ;endif;
if (time >= 350.000 & time < 360.000 ); xcol= 52 ;endif;
if (time >= 360.000 & time < 370.000 ); xcol= 53 ;endif;
if (time >= 370.000 & time < 380.000 ); xcol= 54 ;endif;
if (time >= 380.000 & time < 390.000 ); xcol= 55 ;endif;
if (time >= 390.000 & time < 400.000 ); xcol= 56 ;endif;
if (time >= 400.000 & time < 410.000 ); xcol= 57 ;endif;
if (time >= 410.000 & time < 420.000 ); xcol= 58 ;endif;
if (time >= 420.000 & time < 430.000 ); xcol= 59 ;endif;
if (time >= 430.000 & time < 440.000 ); xcol= 60 ;endif;
if (time >= 440.000 & time < 450.000 ); xcol= 61 ;endif;
if (time >= 450.000 & time < 460.000 ); xcol= 62 ;endif;
if (time >= 460.000 & time < 470.000 ); xcol= 63 ;endif;
if (time >= 470.000 & time < 480.000 ); xcol= 64 ;endif;
if (time >= 480.000 & time < 490.000 ); xcol= 65 ;endif;
if (time >= 490.000 & time < 500.000 ); xcol= 66 ;endif;
if (time >= 500.000 & time < 510.000 ); xcol= 67 ;endif;
if (time >= 510.000 & time < 520.000 ); xcol= 68 ;endif;
if (time >= 520.000 & time < 530.000 ); xcol= 69 ;endif;
if (time >= 530.000 & time < 540.000 ); xcol= 70 ;endif;
if (time >= 540.000 & time < 550.000 ); xcol= 71 ;endif;
if (time >= 550.000 & time < 560.000 ); xcol= 72 ;endif;
if (time >= 560.000 & time < 570.000 ); xcol= 73 ;endif;
if (time >= 570.000 & time < 580.000 ); xcol= 74 ;endif;
if (time >= 580.000 & time < 590.000 ); xcol= 75 ;endif;
if (time >= 590.000 & time < 600.000 ); xcol= 76 ;endif;
if (time >= 600.000 & time < 610.000 ); xcol= 77 ;endif;
if (time >= 610.000 & time < 620.000 ); xcol= 78 ;endif;
if (time >= 620.000 & time < 630.000 ); xcol= 79 ;endif;
if (time >= 630.000 & time < 640.000 ); xcol= 80 ;endif;
if (time >= 640.000 & time < 650.000 ); xcol= 81 ;endif;
if (time >= 650.000 & time < 660.000 ); xcol= 82 ;endif;
if (time >= 660.000 & time < 670.000 ); xcol= 83 ;endif;
if (time >= 670.000 & time < 680.000 ); xcol= 84 ;endif;
if (time >= 680.000 & time < 690.000 ); xcol= 85 ;endif;
if (time >= 690.000 & time < 700.000 ); xcol= 86 ;endif;
if (time >= 700.000 & time < 710.000 ); xcol= 87 ;endif;
if (time >= 710.000 & time < 720.000 ); xcol= 88 ;endif;
if (time >= 720.000 & time < 730.000 ); xcol= 89 ;endif;
if (time >= 730.000 & time < 740.000 ); xcol= 90 ;endif;
if (time >= 740.000 & time < 750.000 ); xcol= 91 ;endif;
if (time >= 750.000 & time < 760.000 ); xcol= 92 ;endif;
if (time >= 760.000 & time < 770.000 ); xcol= 93 ;endif;
if (time >= 770.000 & time < 780.000 ); xcol= 94 ;endif;
if (time >= 780.000 & time < 790.000 ); xcol= 95 ;endif;
if (time >= 790.000 & time < 800.000 ); xcol= 96 ;endif;
if (time >= 800.000 & time < 810.000 ); xcol= 97 ;endif;
if (time >= 810.000 & time < 820.000 ); xcol= 98 ;endif;
if (time >= 820.000 & time < 830.000 ); xcol= 99 ;endif;
if (time >= 830.000 & time < 840.000 ); xcol=100 ;endif;
if (time >= 840.000 & time < 850.000 ); xcol=101 ;endif;
if (time >= 850.000 & time < 860.000 ); xcol=102 ;endif;
if (time >= 860.000 & time < 870.000 ); xcol=103 ;endif;
if (time >= 870.000 & time < 880.000 ); xcol=104 ;endif;
if (time >= 880.000 & time < 890.000 ); xcol=105 ;endif;
if (time >= 890.000 & time < 900.000 ); xcol=106 ;endif;
if (time >= 900.0 ); xcol= 199 ;endif;
** xcl=color2(time,xi,xd)
** say xcl
** xcl=subwrd(color,xcol)
*
*
'set line 'xcol
* 'set line 'xco
'draw recf 'x-dx2' 'y-dy2' 'x+dx2' 'y+dy2
* 'draw recf 'x'-'dx2' 'y'-'dy2' 'x'+'dx2' 'y'+'dy2
* 'draw mark 3 'x' 'y' 0.1'
endif
endwhile
return
function palgrey()
***************** Rainbow ********************
* 13 color
*
*'set rgb 16 64 64 64'
*
*'set rgb 17 240 0 130'
*'set rgb 18 250 60 60'
*'set rgb 19 240 130 40'
*'set rgb 20 230 175 45'
*'set rgb 21 230 220 50'
*'set rgb 22 160 230 50'
*'set rgb 23 0 220 0'
*'set rgb 24 0 210 140'
*'set rgb 25 0 200 200'
*'set rgb 26 0 160 255'
*'set rgb 27 30 60 255'
*'set rgb 28 130 0 220'
*'set rgb 29 160 0 200'
*
*'set rgb 99 255 255 255'
***************** KMA Radar ********************
* color 30
*
*'set rgb 16 64 64 64'
*
*'set rgb 17 135 217 255'
*'set rgb 18 62 193 255'
*'set rgb 19 7 171 255'
*'set rgb 20 0 141 255'
*'set rgb 21 0 119 255'
*'set rgb 22 105 252 105'
*'set rgb 23 30 242 105'
*'set rgb 24 0 213 0'
*'set rgb 25 0 164 0'
*'set rgb 26 0 128 0'
*'set rgb 27 255 241 111'
*'set rgb 28 255 226 86'
*'set rgb 29 255 208 57'
*'set rgb 30 255 188 30'
*'set rgb 31 255 170 9'
*'set rgb 32 255 156 0'
*'set rgb 33 255 139 27'
*'set rgb 34 255 128 81'
*'set rgb 35 255 110 110'
*'set rgb 36 246 94 103'
*'set rgb 37 232 74 86'
*'set rgb 38 217 52 62'
*'set rgb 39 201 31 37'
*'set rgb 40 188 13 15'
*'set rgb 41 179 0 0'
*'set rgb 42 197 125 255'
*'set rgb 43 194 62 255'
*'set rgb 44 173 7 255'
*'set rgb 45 146 0 255'
*'set rgb 46 127 0 191'
*
*'set rgb 99 255 255 255'
************* Rainbow *************************
* 39 color
*'set rgb 16 64 64 64'
*
*'set rgb 17 153 102 255'
*'set rgb 18 131 87 255'
*'set rgb 19 109 73 255'
*'set rgb 20 87 58 255'
*'set rgb 21 66 44 255'
*'set rgb 22 44 29 255'
*'set rgb 23 0 0 255'
*'set rgb 24 0 72 255'
*'set rgb 25 0 117 255'
*'set rgb 26 0 156 255'
*'set rgb 27 0 192 255'
*'set rgb 28 0 213 255'
*'set rgb 29 0 232 255'
*'set rgb 30 0 246 255'
*'set rgb 31 0 255 249'
*'set rgb 32 0 255 213'
*'set rgb 33 0 255 193'
*'set rgb 34 0 255 159'
*'set rgb 35 0 255 126'
*'set rgb 36 0 255 77'
*'set rgb 37 0 255 15'
*'set rgb 38 50 255 0'
*'set rgb 39 114 255 0'
*'set rgb 40 152 255 0'
*'set rgb 41 181 255 0'
*'set rgb 42 227 255 0'
*'set rgb 43 243 255 0'
*'set rgb 44 255 253 0'
*'set rgb 45 255 232 0'
*'set rgb 46 255 211 0'
*'set rgb 47 255 189 0'
*'set rgb 48 255 168 0'
*'set rgb 49 255 146 0'
*'set rgb 50 255 124 0'
*'set rgb 51 255 102 0'
*'set rgb 52 255 87 0'
*'set rgb 53 255 73 0'
*'set rgb 54 255 58 0'
*'set rgb 55 255 29 0'
*'set rgb 56 255 0 0'
*
*'set rgb 99 255 255 255'
*
************* Diff ****************************
* 37 color
*
*'set rgb 16 0 0 0'
*
*'set rgb 17 0 0 255'
*'set rgb 18 38 38 255'
*'set rgb 19 51 51 255'
*'set rgb 20 63 63 255'
*'set rgb 21 76 76 255'
*'set rgb 22 89 89 255'
*'set rgb 23 102 102 255'
*'set rgb 24 114 114 255'
*'set rgb 25 127 127 255'
*'set rgb 26 140 140 255'
*'set rgb 27 153 153 255'
*'set rgb 28 165 165 255'
*'set rgb 29 178 178 255'
*'set rgb 30 191 191 255'
*'set rgb 31 204 204 255'
*'set rgb 32 216 216 255'
*'set rgb 33 229 229 255'
*'set rgb 34 242 242 255'
*'set rgb 35 255 255 255'
*'set rgb 36 255 243 243'
*'set rgb 37 255 230 230'
*'set rgb 38 255 217 217'
*'set rgb 39 255 205 205'
*'set rgb 40 255 192 192'
*'set rgb 41 255 179 179'
*'set rgb 42 255 166 166'
*'set rgb 43 255 154 154'
*'set rgb 44 255 141 141'
*'set rgb 45 255 128 128'
*'set rgb 46 255 115 115'
*'set rgb 47 255 102 102'
*'set rgb 48 255 90 90'
*'set rgb 49 255 77 77'
*'set rgb 50 255 64 64'
*'set rgb 51 255 51 51'
*'set rgb 52 255 38 38'
*'set rgb 53 255 0 0'
*
*'set rgb 99 255 255 255'
************** Forest ****************
* 65 color
*
*'set rgb 16 64 64 64'
*
*'set rgb 17 0 0 30'
*'set rgb 18 5 26 2'
*'set rgb 19 10 32 4'
*'set rgb 20 15 38 6'
*'set rgb 21 20 44 8'
*'set rgb 22 25 50 10'
*'set rgb 23 30 56 12'
*'set rgb 24 35 62 14'
*'set rgb 25 40 68 16'
*'set rgb 26 45 74 18'
*'set rgb 27 50 80 20'
*'set rgb 28 55 86 22'
*'set rgb 29 60 92 24'
*'set rgb 30 65 98 26'
*'set rgb 31 70 104 28'
*'set rgb 32 75 110 30'
*'set rgb 33 80 116 32'
*'set rgb 34 85 122 34'
*'set rgb 35 90 128 36'
*'set rgb 36 95 134 38'
*'set rgb 37 100 140 40'
*'set rgb 38 105 146 42'
*'set rgb 39 110 152 44'
*'set rgb 40 115 158 46'
*'set rgb 41 120 164 48'
*'set rgb 42 125 170 50'
*'set rgb 43 130 176 52'
*'set rgb 44 135 182 54'
*'set rgb 45 140 188 56'
*'set rgb 46 145 194 58'
*'set rgb 47 150 200 60'
*'set rgb 48 155 206 62'
*'set rgb 49 255 250 110'
*'set rgb 50 249 238 98'
*'set rgb 51 243 226 86'
*'set rgb 52 237 214 74'
*'set rgb 53 231 202 62'
*'set rgb 54 225 190 50'
*'set rgb 55 219 178 38'
*'set rgb 56 213 166 26'
*'set rgb 57 207 154 14'
*'set rgb 58 201 142 2'
*'set rgb 59 195 130 0'
*'set rgb 60 189 118 0'
*'set rgb 61 183 106 0'
*'set rgb 62 177 94 0'
*'set rgb 63 171 82 0'
*'set rgb 64 165 70 0'
*'set rgb 65 159 58 0'
*'set rgb 66 153 46 0'
*'set rgb 67 147 34 0'
*'set rgb 68 141 22 0'
*'set rgb 69 135 10 0'
*'set rgb 70 129 0 0'
*'set rgb 71 123 0 0'
*'set rgb 72 117 0 0'
*'set rgb 73 111 0 0'
*'set rgb 74 105 0 0'
*'set rgb 75 99 0 0'
*'set rgb 76 93 0 0'
*'set rgb 77 87 0 0'
*'set rgb 78 81 0 0'
*'set rgb 79 75 0 0'
*'set rgb 80 69 0 0'
*'set rgb 81 20 0 0'
*
*'set rgb 99 255 255 255'
*
************** Rainbow [Red ~ Blue] ****************
* color 100
'set rgb 16 64 64 64'
'set rgb 17 0 0 130'
'set rgb 18 0 0 134'
'set rgb 19 0 0 146'
'set rgb 20 0 0 154'
'set rgb 21 0 0 166'
'set rgb 22 0 0 174'
'set rgb 23 0 0 186'
'set rgb 24 0 0 198'
'set rgb 25 0 0 206'
'set rgb 26 0 0 218'
'set rgb 27 0 0 226'
'set rgb 28 0 0 238'
'set rgb 29 0 0 250'
'set rgb 30 0 0 254'
'set rgb 31 0 11 254'
'set rgb 32 0 19 254'
'set rgb 33 0 31 254'
'set rgb 34 0 39 254'
'set rgb 35 0 51 254'
'set rgb 36 0 63 254'
'set rgb 37 0 71 254'
'set rgb 38 0 83 254'
'set rgb 39 0 91 254'
'set rgb 40 0 103 254'
'set rgb 41 0 115 254'
'set rgb 42 0 123 254'
'set rgb 43 0 134 254'
'set rgb 44 0 142 254'
'set rgb 45 0 154 254'
'set rgb 46 0 162 254'
'set rgb 47 0 174 254'
'set rgb 48 0 186 254'
'set rgb 49 0 194 254'
'set rgb 50 0 206 254'
'set rgb 51 0 214 254'
'set rgb 52 0 226 254'
'set rgb 53 0 238 254'
'set rgb 54 0 246 254'
'set rgb 55 0 254 254'
'set rgb 56 7 254 246'
'set rgb 57 19 254 234'
'set rgb 58 31 254 222'
'set rgb 59 39 254 214'
'set rgb 60 51 254 202'
'set rgb 61 59 254 194'
'set rgb 62 71 254 181'
'set rgb 63 79 254 174'
'set rgb 64 91 254 162'
'set rgb 65 103 254 150'
'set rgb 66 111 254 142'
'set rgb 67 123 254 130'
'set rgb 68 130 254 123'
'set rgb 69 142 254 111'
'set rgb 70 154 254 99'
'set rgb 71 162 254 91'
'set rgb 72 174 254 79'
'set rgb 73 182 254 71'
'set rgb 74 194 254 59'
'set rgb 75 202 254 51'
'set rgb 76 214 254 39'
'set rgb 77 226 254 27'
'set rgb 78 234 254 19'
'set rgb 79 246 254 7'
'set rgb 80 254 254 0'
'set rgb 81 254 242 0'
'set rgb 82 254 230 0'
'set rgb 83 254 222 0'
'set rgb 84 254 210 0'
'set rgb 85 254 202 0'
'set rgb 86 254 190 0'
'set rgb 87 254 178 0'
'set rgb 88 254 170 0'
'set rgb 89 254 158 0'
'set rgb 90 254 150 0'
'set rgb 91 254 138 0'
'set rgb 92 254 130 0'
'set rgb 93 254 119 0'
'set rgb 94 254 107 0'
'set rgb 95 254 99 0'
'set rgb 96 254 87 0'
'set rgb 97 254 79 0'
'set rgb 98 254 67 0'
'set rgb 99 254 55 0'
'set rgb 100 254 47 0'
'set rgb 101 254 35 0'
'set rgb 102 254 27 0'
'set rgb 103 254 15 0'
'set rgb 104 254 7 0'
'set rgb 105 249 0 0'
'set rgb 106 236 0 0'
'set rgb 107 227 0 0'
'set rgb 108 214 0 0'
'set rgb 109 205 0 0'
'set rgb 110 192 0 0'
'set rgb 111 179 0 0'
'set rgb 112 170 0 0'
'set rgb 113 157 0 0'
'set rgb 114 148 0 0'
'set rgb 115 135 0 0'
'set rgb 116 130 0 0'
'set rgb 199 0 0 0'
return
view raw Library_For_Mapping_ASCII_Format_File.gs hosted with ❤ by GitHub

 

 

  • 컬러바 생성을 위한 라이브러리

 

function colorbar (args)
vert=subwrd(args,1)
xmid=subwrd(args,2)
ymid=subwrd(args,3)
barw=subwrd(args,4)
sfc=subwrd(args,5)
xsiz=subwrd(args,6)
ysiz=subwrd(args,7)
col1=subwrd(args,8)
col2=subwrd(args,9)
val1=subwrd(args,10)
val2=subwrd(args,11)
cl1=subwrd(args,12)
cl2=subwrd(args,13)
cli=subwrd(args,14)
dum2=cli
if (cli<0) ;dum2=cli*-1;endif
dum1=cl2-cl1
if (cl1>cl2);dum1=cl1-cl2;endif
blnum=dum1/dum2+1
* col1 & col2 : should be determined when dot-alb.f
* (min and max # of color scale)
* val1 & val2 : actual maginitude corresponding colors, col1 & col2
* if you set alb=0 to col#=40, and set alb=1.0 to col#=85, then
* val1=0, val2=1, col1=40, col2=85
* cl1=val1, cl2=val2, I will change later,
* cli=interval of # under colorbar
cnum=col2-col1+1
stroff=0.05*sfc
strxsiz=0.12*sfc
strysiz=0.13*sfc
*
* logic for setting the bar orientation with user overides
*
if(vert = 0) ; vchk = 0 ; endif
if(vert = 1) ; vchk = 1 ; endif
*
* vertical bar
*
if (vchk = 1 )
ywid = ysiz / cnum
xl = xmid - barw/2
xr = xmid + barw/2
yb = ymid-ysiz/2
y1 = yb
vert = 1
else
*
* horizontal bar
*
xwid = xsiz / cnum
yt = ymid + barw/2
yb = ymid - barw/2
xl = xmid-xsiz/2
x1 = xl
'set string 1 tc 5'
vert = 0
endif
*
* Plot colorbar
*
num = 1
while (num<=cnum)
col = col1-1+num
if (vert)
yt = yb + ywid
else
xr = xl + xwid
endif
'set line 'col
'draw recf 'xl' 'yb' 'xr' 'yt
num = num + 1
if (vert); yb = yt;
else; xl = xr; endif;
endwhile
bli = 1
'set strsiz 'strxsiz' 'strysiz
if (vert)
'set string 1 l 5'
cl = cl1
yp=y1 + ywid
xp=xr+stroff
yd = (ysiz - 2*ywid)/dum1*dum2
while (bli <= blnum)
'draw string 'xp' 'yp' 'cl
cl = cl + cli
yp = yp + yd
bli=bli+1
endwhile
else
'set string 1 tc 5'
xp=x1 + xwid
cl = cl1
yp=yb-stroff
xd = (xsiz - 2*xwid)/dum1*dum2
while (bli <= blnum)
'draw string 'xp' 'yp' 'cl
cl = cl + cli
xp = xp + xd
bli=bli+1
endwhile
endif
return
view raw Library_For_Creating_Colorbar_From_ASCII_Format_File.gs hosted with ❤ by GitHub

 

 

  • 가시화

 

#!/bin/csh
#rm -rf FIG/*
set a = `ls DATA/2013??????`
foreach fn ($a)
set fn2 = `echo $fn | cut -c 6-15`
set yy = `echo $fn | cut -c 6-9`
set mm = `echo $fn | cut -c 10-11`
set dd = `echo $fn | cut -c 12-13`
set tt = `echo $fn | cut -c 14-15`
set fn3 = `echo ${yy}.${mm}.${dd} ${tt}UTC`
paste DATA/Geo.dat DATA/$fn2 > fort.1
cat >! 2.exec << EOF
reinit
open GrADS/Map.ctl
enable print 2.gx
set font 5
set vpage 0 8.5 0 11
set parea 0 8.5 2.0 10.0
set mpdraw off
set xlopts 1 5.5 0.15
set ylopts 1 5.5 0.15
set lat 31. 45.
set lon 121. 135.5
set cmin 9999
set cint 2
set xlint 2.0
set ylint 2.0
set grid on
set grads off
draw title $fn3 [Clear]
set strsiz 0.15 0.15
draw string 3.5 1.4 Surface Solar Radiation [Wm\`a-2\`n]
d t
run GrADS/Library_For_Mapping_ASCII_Format_File.gs fort.1 1 2 3 3
run GrADS/Library_For_Creating_Colorbar_From_ASCII_Format_File.gs 0 4.25 1.0 0.28 1.1 8.0 1.0 17 116 0 900 0 900 100
*run GrADS/Library_For_Creating_Colorbar_From_ASCII_Format_File.gs 1 8.0 6.0 0.25 1.0 1.0 8.0 16 116 0 900 0 900 50
set rgb 116 0 0 0
draw shp GrADS/Map/GSHHG/GSHHS_shp/f/GSHHS_f_L1.shp
*draw shp GrADS/Map/KOR_adm/KOR_adm0.shp
*draw shp GrADS/Map/CHN_adm/CHN_adm0.shp
*draw shp GrADS/Map/PRK_adm/PRK_adm0.shp
*draw shp GrADS/Map/JPN_adm/JPN_adm0.shp
print
disable print
!gxps -c -i 2.gx -o 2.ps
quit
EOF
grads -bpc 'exec 2.exec'
convert 2.ps 2.gif
#display 2.gif
mv -f 2.gif FIG/${fn2}_SFC.gif
end
rm -f 2.ps 2.gx
view raw Visualization_Using_ASCII_Format_Radiation_Transfer_Model_Result.csh hosted with ❤ by GitHub

 

 

 결과

  • 복사전달모델 결과를 이용한 가시화

etc-image-0
그림. 복사전달모델 결과를 이용한 가시화.

 

 참고 문헌

[논문]

  • 없음

[보고서]

  • 없음

[URL]

  • 없음

 

 문의사항

[기상학/프로그래밍 언어]

  • sangho.lee.1990@gmail.com

[해양학/천문학/빅데이터]

  • saimang0804@gmail.com
 
저작자표시 비영리 변경금지

'프로그래밍 언어 > GrADS' 카테고리의 다른 글

[GrADS] 그라즈 컬러 테이블 및 팔레트 소개  (0) 2020.04.07
[GrADS, ShellScript] 그라즈 HDF 형식인 Terra/MODIS 및 Aqua/MODIS 자료를 이용한 에어로졸 광학두께 가시화  (0) 2019.09.22
[GrADS, ShellScript] 그라즈 NetCDF 형식인 NOAA 월평균 해수면 온도 자료를 이용한 가시화  (0) 2019.09.07
[GrADS, ShellScript] 그라즈 Grib 형식인 ECMWF 모델 결과를 이용한 가시화  (0) 2019.09.07
[GrADS, ShellScript] 그라즈 중규모 기상모델 (WRF ARW)을 이용한 가시화  (0) 2019.09.02
프로그래밍 언어/GrADS 관련 글
더 보기

티스토리툴바