امروز:پنج شنبه, ۹ فروردین , ۱۴۰۳
زمان انتشار : دوشنبه, آبان 26ام, 1393 | پرینت مطلب |رضا ابراهیمی| بازديد: 743 بار

چگونه جبهه را درگردس شناسایی کنیم

اسکریپت زیر نحوه ترسیم جبهه را در نرم افزار گردس نشان می دهد

تعداد روزها

n=156

مشخصات جغرافیایی

lonrange='0 80'
latrange='0 60'
level=1000

reinit
مسیر دهی

inpath1='H:\Data\air\air.'
inpath2='H:\Data\uwnd\uwnd.'
inpath3='H:\Data\vwnd\vwnd.'

inext='.nc'

oupath='C:\Era\Dust\'

dates='C:\Era\DUst\DustDates.txt'

i=1
while i<=n
say i
rec=read(dates)
rec=sublin(rec,2)
y.i=subwrd(rec,1)
m.i=subwrd(rec,2)
d.i=subwrd(rec,3)

infile1=inpath1%y.i%inext
infile2=inpath2%y.i%inext
infile3=inpath3%y.i%inext

'sdfopen ' infile1
'sdfopen ' infile2
'sdfopen ' infile3

نام گذاری ماهها

if m.i=1 ; mon=jan ; endif
if m.i=2 ; mon=feb ; endif
if m.i=3 ; mon=mar ; endif
if m.i=4 ; mon=apr ; endif
if m.i=5 ; mon=may ; endif
if m.i=6 ; mon=jnn ; endif
if m.i=7 ; mon=jul ; endif
if m.i=8 ; mon=aug ; endif
if m.i=9 ; mon=sep ; endif
if m.i=10 ; mon=oct ; endif
if m.i=11 ; mon=nov ; endif
if m.i=12 ; mon=dec ; endif

tar=y.i%'_'%m.i%'_'%d.i
now='00z'%d.i%mon%y.i

'set lon 'lonrange
'set lat ' latrange
'set lev 'level
'set time 'now

'th=air*pow((1000/lev),(287/1005))'

معادلات

'define dx = cdiff(lon,x)*(3.1416/180)*cos(lat*3.1416/180)* 6.37e6*2'
'define dy = cdiff(lat,y)*(3.1416/180)*6.37e6*2'
'define dthx = cdiff(th,x)'
'define dthy = cdiff(th,y) '
'define dux = cdiff(uwnd.2,x)'
'define duy = cdiff(uwnd.2,y) '
'define dvx = cdiff(vwnd.3,x)'
'define dvy = cdiff(vwnd.3,y) '

'dth=mag(dthx/dx,dthy/dy)'

معادلات وزش


'A=(-1*dthx/dx) * ( ((dux/dx)*(dthx/dx)) + ((dvx/dx)*(dthy/dy)) )'
'B=(dthy/dy) * ( ((duy/dy)*(dthx/dx)) + ((dvy/dy)*(dthy/dy)) )'

'f=(1/dth)*(A-B)*3600*24*100000'
* f is in K/100Km/day

نمایه های ترسیم

'set gxout contour'
'c'
'set csmooth on'
'set map 10 10 10'
'set ccolor 1'
'set cmin 0.1'

'd f'



مطالب مرتبط

avatar

نویسنده: رضا ابراهیمی

reza ebrahimi Ph.D Climatology yazd University