***************************************************** load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl" load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl" load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl" ;************************************************************ begin a = addfile("/home/bittner/Uni/Bachelorarbeit/Daten/Dezember15/ruc2_252_20101215_1600_000.grb" , "r") b = addfile("/home/bittner/Uni/Bachelorarbeit/Daten/Januar07/ruc2_252_20110107_1600_000.grb" , "r") c = addfile("/home/bittner/Uni/Bachelorarbeit/Daten/Juli21/ruc2_252_20110721_1500_000.grb" , "r") d = addfile("/home/bittner/Uni/Bachelorarbeit/Daten/August24/ruc2_252_20110824_1500_000.grb" , "r") Feuchte1 = a->R_H_252_ISBL ;Temperatur auf Isobarenlevel, 3D-Variable (Longnitude, Latitude, Höhe) Feuchte2 = b->R_H_252_ISBL Feuchte3 = c->R_H_252_ISBL Feuchte4 = d->R_H_252_ISBL hohe = a->lv_ISBL2 ;***********************************Plot************************************************************************* wks_type = "pdf" wks = gsn_open_wks(wks_type,"Feuchteprofil_24.08.2011 ") res = True res@tiMainString ="profil des 24.08.2010" ;Titel res@tiYAxisString ="Höhe in hPa" res@tiXAxisString ="Temperatur in K" res@xyLineColors =(/"red"/) ;Linienfarbe res@trYReverse =True ;reverser Y-axis plot= gsn_csm_xy(wks, Feuchte4(:,126,213), hohe,res) res = True res@tiYAxisString ="Höhe in hPa" res@tiXAxisString ="Temperatur in K" res@xyLineColors =(/"red"/) ;Linienfarbe res@trYReverse =True ;reverser Y-axis plot= gsn_csm_xy(wks, Feuchte3(:,126,213), hohe,res)