Diabatic Heating

Sat Aug 4 01:20:22 EDT 2007

   YAHOO.Shortcuts.hasSensitiveText = false; YAHOO.Shortcuts.doUlt = false; YAHOO.Shortcuts.location = "us"; YAHOO.Shortcuts.lang = "us"; YAHOO.Shortcuts.document_id = 0; YAHOO.Shortcuts.document_type = ""; YAHOO.Shortcuts.document_title = ""; YAHOO.Shortcuts.document_publish_date = ""; YAHOO.Shortcuts.document_author = ""; YAHOO.Shortcuts.annotationSet = { "lw_1186204683_0": { "text": "airtemp.nc", "extended": 0, "startchar": 856, "endchar": 865, "start": 856, "end": 865, "extendedFrom": "", "weight": 1, "type": ["shortcuts:/us/place/virtual/web_site"], "category": ["IDENTIFIER"] ,"context": " MEAN) of air temprature,u,v and omega.  \x27sdfopen airtemp.nc\x27 \x27sdfopen uwind.nc\x27 \x27sdfopen vwind.nc\x27 \x27sdfopen"}, "lw_1186204683_1": { "text": "uwind.nc", "extended": 0, "startchar": 888, "endchar": 895, "start": 888, "end": 895, "extendedFrom": "", "weight": 1, "type": ["shortcuts:/us/place/virtual/web_site"], "category": ["IDENTIFIER"] ,"context": " and omega.
 \x27sdfopen airtemp.nc\x27 \x27sdfopen uwind.nc\x27 \x27sdfopen vwind.nc\x27 \x27sdfopen omega.nc\x27   \x27PP ="}, "lw_1186204683_2": { "text": "vwind.nc", "extended": 0, "startchar": 918, "endchar": 925, "start": 918, "end": 925, "extendedFrom": "", "weight": 1, "type": ["shortcuts:/us/place/virtual/web_site"], "category": ["IDENTIFIER"] ,"context": " airtemp.nc\x27 \x27sdfopen uwind.nc\x27 \x27sdfopen vwind.nc\x27 \x27sdfopen omega.nc\x27   \x27PP = lev\x27 \x27set lat 0"}, "lw_1186204683_3": { "text": "omega.nc", "extended": 0, "startchar": 948, "endchar": 955, "start": 948, "end": 955, "extendedFrom": "", "weight": 1, "type": ["shortcuts:/us/place/virtual/web_site"], "category": ["IDENTIFIER"] ,"context": " uwind.nc\x27 \x27sdfopen vwind.nc\x27 \x27sdfopen omega.nc\x27   \x27PP = lev\x27 \x27set lat 0 50\x27 \x27set lon 400"}, "lw_1186204683_4": { "text": "Tehran, Iran", "extended": 0, "startchar": 2599, "endchar": 2610, "start": 2599, "end": 2610, "extendedFrom": "",
 "weight": 0.999122, "type": ["shortcuts:/us/instance/place/ir/town"], "category": ["PLACE"] ,"context": " Dept. Tarbiat Modares University Tehran, Iran        Need a vacation? Get great deals to" , "metaData": { "geoArea": "85.9871", "geoCountry": "Iran", "geoIsoCountryCode": "IR", "geoLocation": "(51.408829, 35.702591)", "geoName": "Tehran", "geoPlaceType": "Town", "geoState": "Tehran", "geoTown": "Tehran", "type": "shortcuts:/us/instance/place/ir/town" } }, "lw_1186204683_5": { "text": "Get great deals to amazing places ", "extended": 0, "startchar": 2863, "endchar": 2897, "start": 2863, "end": 2897, "extendedFrom": "", "weight": 1, "type": ["shortcuts:/us/instance/identifier/hyperlink/http"], "category": ["IDENTIFIER"] ,"context": " University Tehran, Iran        Need a vacation? Get great deals to amazing places on Yahoo! Travel." , "metaData": { "linkHref":
 "http://us.rd.yahoo.com/evt=48256/*http://travel.yahoo.com/;_ylc=X3oDMTFhN2hucjlpBF9TAzk3NDA3NTg5BHBvcwM1BHNlYwNncm91cHMEc2xrA2VtYWlsLW5jbQ--", "linkProtocol": "http", "linkRel": "nofollow", "linkTarget": "_blank" } }, "lw_1186204683_6": { "text": "Yahoo! Travel", "extended": 0, "startchar": 2905, "endchar": 2917, "start": 2905, "end": 2917, "extendedFrom": "", "weight": 0.462173, "type": ["shortcuts:/us/instance/organization/company/yahoo_property"], "category": ["ORGANIZATION"] ,"context": " a vacation? Get great deals to amazing places on Yahoo! Travel." , "metaData": { "yprop_name": "Yahoo! Travel", "yprop_url": "http://travel.yahoo.com/" } } };  YAHOO.Shortcuts.overlaySpaceId = "97546169";  YAHOO.Shortcuts.hostSpaceId = "97546168";       Dear GrADS-users,
I am trying to write a script for calculating diabatic heating.
Actually, it can be computed from the thermodynamic equation.
It is consist of three following terms:
''(dtheta/dt)+(u*dtheta/dx+v*dtheta/dy)-(omega*dtheta/dp)''
I have written a script and it seems that I have succeeded to compute terms 1 and 2 but somthing is wrong with term3.
  As I am not sure about their accuracy, would somebody please help me and check if the script is correct for terms 1 and 2?
and also,
I have really appreciate any help to make a correction of term3.
  Sorry for long note and thank you very much for your time.
Azar Zarrin

'reinit'
*****opening four files (NCEP/NCAR DAILY MEAN) of air temprature,u,v and omega.
'sdfopen airtemp.nc'
'sdfopen uwind.nc'
'sdfopen vwind.nc'
'sdfopen omega.nc'
  'PP = lev'
'set lat 0 50'
'set lon 400 480'
'set lev 1000 100'
'set t 2'
  "define Term1 = (air(t-1)*pow(1000/PP(t-1),0.286)-air*pow(1000/PP,0.286)))/24"
'set lat 32.5'
'set zlog on'
'd Term1'
  say 'hit return to continue and plot'
pull cmd
'c'
  'set lat 0 50'
'set lon 400 480'
'set lev 1000 100'
'set t 2'
'define pt=air+273.15*pow(1000/PP,0.286)'
"dptx =(cdiff(pt,x))"
"dpty =(cdiff(pt,y))"
'set dfile 2'
'u = uwnd'
"dux =(cdiff(u,x))"
'set dfile 3'
'v=vwnd'
"dvy =(cdiff(v,y))"
'dx = cdiff(lon,x)*cos(lat*3.14/180)*111325'
"dy =cdiff(lat,y)*55445"
'define Term2 = (pt*(dux/dx))+(pt*(dvy/dy))'
'set lat 32.5'
'set zlog on'
'd Term2'

say 'hit return to continue and plot'
pull cmd
'c'
  'set lat 0 50'
'set lon 400 480'
'set lev 1000 100'
'set t 2'
'set dfile 1'
"define Term3=(air(z-1)*pow(1000/PP(z-1),0.286)-air(z+1)*pow(1000/PP(z+1),0.286)))/(PP(z-1)-PP(z+1))"
'set lat 32.5'
'set zlog on'
'd Term3'
say error.....
Cannot plot data - all undefined values.
  ----------------------------------------------------------------
Azar Zarrin
Geography Dept.
Tarbiat Modares University
Tehran, Iran

---------------------------------
Ready for the edge of your seat? Check out tonight's top picks on Yahoo! TV.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20070803/3303d72b/attachment.html 