چگونه جبهه را درگردس شناسایی کنیم
اسکریپت زیر نحوه ترسیم جبهه را در نرم افزار گردس نشان می دهد
تعداد روزها
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'