[gradsusr] Question about displaying minutes properly on Grads image
Nimrod Micael
nmicael at aggies.ncat.edu
Fri May 24 09:06:34 EDT 2013
I am trying to display on where it says Forecast Valid:
Date: Month: Yr Time: with the Z at the end of it. But there's
issues with it displaying the minutes.
Can you find out what the problem is? Also I would like the output file
where the printim image to be in the format of:
Date_Month_Yr_Time: with the Z at the end of it.
My Code:
''reinit'
'clear'
'open 20110416_d03.ctl'
'set display color white'
'set map 1 1 5'
'set annot 1 6'
'set xlopts 1 6 .17'
'set ylopts 1 6 .16'
'set mpdset hires'
*'set lon -105 -60'
*'set lat 24 54'
'set clab on'
'set grid off'
'set timelab off'
'set gxout grfill'
***************************************************
* Begin loop for levels and time
**************************************************
tt=170
while (tt<=190)
'clear'
'set grads off'
'set t 'tt
'q time'
time1=subwrd(result,3)
say time1
hour=substr(time1,1,2)
say hour
day=substr(time1,5,2)
say day
mon=substr(time1,6,3)
if (JAN = mon) ; mon = 01 ; endif
if (FEB = mon) ; mon = 02 ; endif
if (MAR = mon) ; mon = 03 ; endif
if (APR = mon) ; mon = 04 ; endif
if (MAY = mon) ; mon = 05 ; endif
if (JUN = mon) ; mon = 06 ; endif
if (JUL = mon) ; mon = 07 ; endif
if (AUG = mon) ; mon = 08 ; endif
if (SEP = mon) ; mon = 09 ; endif
if (OCT = mon) ; mon = 10 ; endif
if (NOV = mon) ; mon = 11 ; endif
if (DEC = mon) ; mon = 12 ; endif
say mon
year=substr(time1,9,5)
say year
date = year''mon''day'_'hour''
say date
time1=subwrd(result,3)
say time1
timeA=substr(time1,3,11)
timeB=substr(time1,7,16)
timeC=substr(time1,1,6)
'q dims'
lev1=subwrd(result,37)
********************************************************
*
* Color Shades
*******************************************************
*Black To white
'set rgb 59 0 0 0'
'set rgb 60 15 15 15'
'set rgb 61 30 30 30'
'set rgb 62 45 45 45'
'set rgb 63 60 60 60'
'set rgb 64 75 75 75'
'set rgb 65 90 90 90'
'set rgb 66 105 105 105'
'set rgb 67 120 120 120'
'set rgb 68 135 135 135'
'set rgb 69 150 150 150'
'set rgb 70 165 165 165'
'set rgb 71 180 180 180'
'set rgb 72 195 195 195'
'set rgb 73 210 210 210'
'set rgb 74 225 225 225'
'set rgb 75 245 245 245'
'set rgb 76 255 255 255'
*light yellow to dark red
'set rgb 21 255 250 170'
'set rgb 22 255 232 120'
'set rgb 23 255 192 60'
'set rgb 24 255 160 0'
'set rgb 25 255 96 0'
'set rgb 26 255 50 0'
'set rgb 27 225 20 0'
'set rgb 28 192 0 0'
'set rgb 29 165 0 0'
*light blue to dark blue
'set rgb 41 225 255 255'
'set rgb 42 180 240 250'
'set rgb 43 150 210 250'
'set rgb 44 120 185 250'
'set rgb 45 80 165 245'
'set rgb 46 60 150 245'
'set rgb 47 40 130 240'
'set rgb 48 30 110 235'
'set rgb 49 20 100 210'
*light green to dark green
'set rgb 31 230 255 225'
'set rgb 32 200 255 190'
'set rgb 33 180 250 170'
'set rgb 34 150 245 140'
'set rgb 35 120 245 115'
'set rgb 36 80 240 80'
'set rgb 37 55 210 60'
'set rgb 38 30 180 30'
'set rgb 39 15 160 15'
*set rgb 39 5 150 5
*light pink to dark rose
'set rgb 101 255 230 230'
'set rgb 102 255 200 200'
'set rgb 103 248 160 160'
'set rgb 104 230 140 140'
'set rgb 105 230 112 112'
'set rgb 106 230 80 80'
'set rgb 107 200 60 60'
'set rgb 108 180 40 40'
'set rgb 109 164 32 32'
'set xlint 3'
********************************
* Winds Speed Plot (shaded)
********************************
'set gxout shaded'
'set clevs 0 8 12 16 20 24 28 32 36 40 44 48 52 56 60 64 68 72 76 80'
'set ccols 76 41 43 44 35 36 37 38 21 22 23 24 25 26 27 28 106 104 102 42
44 46'
'd max_dbz'
**************************************
*
* Settings for contour plot
***************************************
'set gxout contour'
'set clab masked'
'set cthick 12'
'set clevs 32'
'set ccolor 1'
'set gxout vector'
'set cthick 5'
'set ccolor 1'
'set arrscl 0.15'
'd skip(u10m,4,4);v10m'
*********************************************
* Draw Title, Run Cbar, Draw Map
**********************************************
'draw map'
'run cbarn'
'draw shp c_02ap13'
'set t 1'
'q time'
time2=subwrd(result,3)
'set string 1 tl 5'
'set strsiz 0.17'
'draw string 1.9 8.2 `14km, GFS Initialization: 'time2''
'draw string 1.9 8 `1Reflectivity and Surface Winds'
'set string 1 tr 5'
'set strsiz 0.12'
'draw string 10.9 0.2 `1WRFV3.3.1'
'set string 1 tl 6'
'set strsiz 0.17'
'draw string 1.9 8.4 `1Forecast Valid:'timeB' 'timeC''
***********************************************
* Print picture file
***********************************************
'printim 'timeB '_reflectivity_sfc_wind.png x1000 y800'
"
Thanks!
--
Nimrod Micael
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20130524/848f5434/attachment-0003.html
More information about the gradsusr
mailing list