[GEMPAK #RJR-299468]: GDPVSF Status

Chris and Tim,

Sorry for the delay on this, it's a difficult problem and I'm working on it 
alone.  GDPVSF was not an original GEMPAK program, it was developed by a user 
and added to the baseline before I took over the GEMPAK project at Unidata, and 
the code was out of date before the last update.  The level interpolation is 
the problem here and I will put some time towards it this weekend.

Michael


> 
> Michael:
> 
> Thanks very much for all your assistance to my student on this.  I think
> he's figured out how to use the dynamic tropopause routines in GEMPAK just
> fine, save resolving this problem.  I've told him to be patient as you
> work through this issue.
> 
> Best,
> 
> --Chris
> 
> 
> _________________________________________________________________
> 
> Christopher L. Castro, Ph.D.
> Associate Professor
> Department of Atmospheric Sciences
> University of Arizona
> Physics and Atmospheric Sciences Bldg., Rm 520
> 1118 East Fourth Street
> Tucson, AZ  85721-0081  USA
> 
> Office: (520) 626-5617
> Fax:  (520) 621-6833
> E-mail: address@hidden
> Web page: www.atmo.arizona.edu/personalpages/castro/castro.htm
> 
> On Wed, 13 Mar 2013, Unidata GEMPAK Support wrote:
> 
> > Hi Tim,
> >
> > I'm not there yet but I think I know what to try.  I need to block off some 
> > time for this.
> >
> > Michael
> >
> >
> >> Hi Michael,
> >>
> >> Thanks for responding to my questions about the errors I was having with 
> >> GDPVSF. Have you been able to determine the cause of these issues or 
> >> devise a solution for them? Let me know if you have made any progress on 
> >> fixing the code or if you need any more data samples.
> >>
> >> Thanks,
> >> Tim Lahmers
> >>
> >> On Feb 22, 2013, at 3:52 PM, Unidata GEMPAK Support wrote:
> >>
> >>> in gdpvlv.f, funny enough it says:
> >>>
> >>> C*      Note: the following code is not always quite correct. -JN
> >>>        DO WHILE ( i .lt. nlev .and.
> >>>     +           (istrt .eq. 0 .or. istop .eq. 0 ) )
> >>>          IF ( (ystrt .ge. rlvl ( i-1 ) .and. ystrt .lt. rlvl ( i ) )
> >>>     +                               .or.
> >>>     +         (ystrt .le. rlvl ( i-1 ) .and. ystrt .gt. rlvl ( i ) ) )
> >>>     +      istrt = i - 1
> >>>          IF ( (ystop .gt. rlvl ( i-1 ) .and. ystop .le. rlvl ( i ) )
> >>>     +                               .or.
> >>>     +         (ystop .lt. rlvl ( i-1 ) .and. ystop .ge. rlvl ( i ) ) )
> >>>     +      istop = i
> >>>          i = i + 1
> >>>        END DO
> >>>
> >>> this block is not setting istart and istop so later on the program 
> >>> assumes all levels are to be used.
> >>>
> >>> more later.
> >>>
> >>> Michael
> >>>
> >>>
> >>>>
> >>>> Hi Tim,
> >>>>
> >>>> Thanks for providing as much detail as you did.  I have run three tests 
> >>>> of different levels with GDVINT and GDPVSF on the NAM file to confirm 
> >>>> that the grids look different and I suspect there may be something in 
> >>>> the code which matches levels in the input to those specified by STARTL 
> >>>> and STOPL.  I'll update you when I know more.
> >>>>
> >>>> Michael
> >>>>
> >>>>
> >>>>
> >>>>
> >>>>> Hi Michael,
> >>>>>
> >>>>> Thank you for update on the GDPVSF tool. I am having some more problems 
> >>>>> with the GDPVSF tool, and I think there may be a coding issue. The data 
> >>>>> that I am using are in pressure coordinates, so I have to use GDVINT to 
> >>>>> interpolate the dataset to theta coordinates. When I make this 
> >>>>> conversion and run GDPVSF, the output that is plotted tends to have 
> >>>>> gaps in the dataset. The amount of data that comes through seems to be 
> >>>>> dependent upon the range of theta values that I interpolate the 
> >>>>> pressure data to. I have tried to do this with both archived 
> >>>>> operational NAM data and NARR data that I converted from GRIB files. In 
> >>>>> this particular example, I downloaded the file 2011081500_nam212.gem 
> >>>>> from the following link. 
> >>>>> http://mtarchive.geol.iastate.edu/2011/08/15/gempak/model/ The file was 
> >>>>> re-saved with the .grd extension. The settings that I am using for 
> >>>>> GDVINT are as follows:
> >>>>>
> >>>>> GDFILE = 2011081500_nam212.grd
> >>>>> GDOUTF = dtnam_test.grd
> >>>>> GDATTIM = f00
> >>>>> GLEVEL = 280-380-5
> >>>>> GVCORD = pres/thta
> >>>>> MAXGRD = 4000
> >>>>> GAREA = 15;-130;50;-70
> >>>>> VCOORD  = pres
> >>>>>
> >>>>> Upon interpolating the pressure coordinate data to theta coordinates, I 
> >>>>> then run GDPVSF with the following settings:
> >>>>>
> >>>>> GDFILE   =  dtnam_test.grd
> >>>>> GDOUTF   =  dtnam_test.grd
> >>>>> GFUNC    =  abs(pvor(pres,wnd))
> >>>>> GDATTIM  =  f00
> >>>>> GVCORD   =  thta
> >>>>> STARTL   =  380
> >>>>> STOPL    =  280
> >>>>> DESIRE   =  0.00000015
> >>>>> GDOUTL   =  15
> >>>>> GVOUTC   =  pvab
> >>>>> GPACK    =
> >>>>> GLIST    =  thta;pres
> >>>>> PMAX     =  800
> >>>>>
> >>>>> This process produced a DT surface with holes at higher values of theta 
> >>>>> (see attached file). I also tried to rerun all of these tools with an 
> >>>>> increased range of potential temperature (280-420k). The end result 
> >>>>> yields a DT surface with  most of the values over the CONUS missing. Do 
> >>>>> you know what is causing this inconsistent output with these dynamic 
> >>>>> tropopause surfaces?
> >>>>>
> >>>>> Thanks,
> >>>>> Tim Lahmers
> >>>>>
> >>>>>
> >>>>>
> >>>>>
> >>>>> On Feb 5, 2013, at 9:55 AM, Unidata GEMPAK Support wrote:
> >>>>>
> >>>>>> Hi Tim,
> >>>>>>
> >>>>>> I've added the fix for GDPVSf in GEMPAK 6.8.0 which can be downloaded 
> >>>>>> at http://www.unidata.ucar.edu/downloads/gempak/ or from github at 
> >>>>>> https://github.com/Unidata/gempak/tree/master/gempak/source/programs/gd/gdpvsf
> >>>>>>
> >>>>>> if you don't want to reinstall 6.8.0, you can download the contents of 
> >>>>>> that gdpvsf direvctory and in $GEMPAK/source/programs/gd/gdpvsf/ run 
> >>>>>> "rm $OS_LIB/gdpvsf.a; make everything"
> >>>>>>
> >>>>>> -Michael
> >>>>>>
> >>>>>>> Quick update.
> >>>>>>>
> >>>>>>> I've got the program working on the hrcbob data set once again.  Yet 
> >>>>>>> to test it on the nam and your other GD input, but if that seems okay 
> >>>>>>> I'll send you the update tomorrow.
> >>>>>>>
> >>>>>>> Michael
> >>>>>>>
> >>>>>>>
> >>>>>>>> Hi Tim,
> >>>>>>>>
> >>>>>>>> GDPVSF is a mess right now, yes.  I'm trying to fix it this morning 
> >>>>>>>> and seeing some problems with how the grid files are opened and 
> >>>>>>>> read, likely due to not keeping up with changes in the grid 
> >>>>>>>> diagnostic library.  I'll let you know what I find.
> >>>>>>>>
> >>>>>>>> Michael James
> >>>>>>>> Unidata
> >>>>>>>>
> >>>>>>>>
> >>>>>>>>> Hi GEMPAK Support,
> >>>>>>>>>
> >>>>>>>>> I am currently working on a project that requires me to use the 
> >>>>>>>>> GDPVSF
> >>>>>>>>> function in GEMPAK 6.7.0  to compute Dynamic Tropopause surfaces. 
> >>>>>>>>> When I
> >>>>>>>>> run GDPVSF, there is no output file created; however, the program 
> >>>>>>>>> does
> >>>>>>>>> not list any errors. I read from the below link from your support 
> >>>>>>>>> site
> >>>>>>>>> that GDPVSF has not been working properly due to changes from past
> >>>>>>>>> updates.
> >>>>>>>>>
> >>>>>>>>> http://www.unidata.ucar.edu/support/help/MailArchives/gempak/msg06131.html
> >>>>>>>>>
> >>>>>>>>> Do you know if the problems I am having are related to the same 
> >>>>>>>>> issues
> >>>>>>>>> described above with GDPVSF? What is the status of the GDPVSF 
> >>>>>>>>> function,
> >>>>>>>>> and do you know when a working version of it will be available? The
> >>>>>>>>> input parameters that I am using are listed at the bottom of this
> >>>>>>>>> message. These parameters originally came from a tutorial 
> >>>>>>>>> concerning the
> >>>>>>>>> function that I found from a university website (listed below); 
> >>>>>>>>> however,
> >>>>>>>>> I used a different input file.
> >>>>>>>>>
> >>>>>>>>> http://www.atmos.albany.edu/ftp/gempak/gdpvsf.readme
> >>>>>>>>>
> >>>>>>>>> Thanks,
> >>>>>>>>> Tim Lahmers
> >>>>>>>>>
> >>>>>>>>> GDFILE   =  2012081700_nam212.gem
> >>>>>>>>> GDOUTF   =  DT_test1.grd
> >>>>>>>>> GFUNC    =  mul(avor(obs),stap)
> >>>>>>>>> GDATTIM  =  last
> >>>>>>>>> GVCORD   =  thta
> >>>>>>>>> STARTL   =  260
> >>>>>>>>> STOPL    =  350
> >>>>>>>>> DESIRE   =  0.00000015
> >>>>>>>>> GDOUTL   =  15
> >>>>>>>>> GVOUTC   =  pvbl
> >>>>>>>>> GPACK    =
> >>>>>>>>> GLIST    =  uwnd;vwnd;pres
> >>>>>>>>> PMAX     =  700
> >>>>>>>>>
> >>>>>>>>>
> >>>>>>>>>
> >>>>>>>>
> >>>>>>>
> >>>>>>
> >>>>>> Ticket Details
> >>>>>> ===================
> >>>>>> Ticket ID: RJR-299468
> >>>>>> Department: Support GEMPAK
> >>>>>> Priority: Normal
> >>>>>> Status: Open
> >>>>>>
> >>>>>
> >>>>>
> >>>>>
> >>>>> Hi Michael,
> >>>>>
> >>>>> Thank you for update on the GDPVSF tool. I am having some more problems 
> >>>>> with the GDPVSF tool, and I think there may be a coding issue. The data 
> >>>>> that I am using are in pressure coordinates, so I have to use GDVINT to 
> >>>>> interpolate the dataset to theta coordinates. When I make this 
> >>>>> conversion and run GDPVSF, the output that is plotted tends to have 
> >>>>> gaps in the dataset. The amount of data that comes through seems to be 
> >>>>> dependent upon the range of theta values that I interpolate the 
> >>>>> pressure data to. I have tried to do this with both archived 
> >>>>> operational NAM data and NARR data that I converted from GRIB files. In 
> >>>>> this particular example, I downloaded the file 2011081500_nam212.gem 
> >>>>> from the following link. 
> >>>>> http://mtarchive.geol.iastate.edu/2011/08/15/gempak/model/ The file was 
> >>>>> re-saved with the .grd extension. The settings that I am using for 
> >>>>> GDVINT are as follows:
> >>>>>
> >>>>> GDFILE = 2011081500_nam212.grd
> >>>>> GDOUTF = dtnam_test.grd
> >>>>> GDATTIM = f00
> >>>>> GLEVEL = 280-380-5
> >>>>> GVCORD = pres/thta
> >>>>> MAXGRD = 4000
> >>>>> GAREA = 15;-130;50;-70
> >>>>> VCOORD  = pres
> >>>>>
> >>>>> Upon interpolating the pressure coordinate data to theta coordinates, I 
> >>>>> then run GDPVSF with the following settings:
> >>>>>
> >>>>> GDFILE   =  dtnam_test.grd
> >>>>> GDOUTF   =  dtnam_test.grd
> >>>>> GFUNC    =  abs(pvor(pres,wnd))
> >>>>> GDATTIM  =  f00
> >>>>> GVCORD   =  thta
> >>>>> STARTL   =  380
> >>>>> STOPL    =  280
> >>>>> DESIRE   =  0.00000015
> >>>>> GDOUTL   =  15
> >>>>> GVOUTC   =  pvab
> >>>>> GPACK    =
> >>>>> GLIST    =  thta;pres
> >>>>> PMAX     =  800
> >>>>>
> >>>>> This process produced a DT surface with holes at higher values of theta 
> >>>>> (see attached file). I also tried to rerun all of these tools with an 
> >>>>> increased range of potential temperature (280-420k). The end result 
> >>>>> yields a DT surface with  most of the values over the CONUS missing. Do 
> >>>>> you know what is causing this inconsistent output with these dynamic 
> >>>>> tropopause surfaces?
> >>>>>
> >>>>> Thanks,
> >>>>> Tim Lahmers
> >>>>>
> >>>>>
> >>>>>
> >>>>>
> >>>>> On Feb 5, 2013, at 9:55 AM, Unidata GEMPAK Support wrote:
> >>>>>
> >>>>>> Hi Tim,
> >>>>>>
> >>>>>> I've added the fix for GDPVSf in GEMPAK 6.8.0 which can be downloaded 
> >>>>>> at http://www.unidata.ucar.edu/downloads/gempak/ or from github at 
> >>>>>> https://github.com/Unidata/gempak/tree/master/gempak/source/programs/gd/gdpvsf
> >>>>>>
> >>>>>> if you don't want to reinstall 6.8.0, you can download the contents of 
> >>>>>> that gdpvsf direvctory and in $GEMPAK/source/programs/gd/gdpvsf/ run 
> >>>>>> "rm $OS_LIB/gdpvsf.a; make everything"
> >>>>>>
> >>>>>> -Michael
> >>>>>>
> >>>>>>> Quick update.
> >>>>>>>
> >>>>>>> I've got the program working on the hrcbob data set once again.  Yet 
> >>>>>>> to test it on the nam and your other GD input, but if that seems okay 
> >>>>>>> I'll send you the update tomorrow.
> >>>>>>>
> >>>>>>> Michael
> >>>>>>>
> >>>>>>>
> >>>>>>>> Hi Tim,
> >>>>>>>>
> >>>>>>>> GDPVSF is a mess right now, yes.  I'm trying to fix it this morning 
> >>>>>>>> and seeing some problems with how the grid files are opened and 
> >>>>>>>> read, likely due to not keeping up with changes in the grid 
> >>>>>>>> diagnostic library.  I'll let you know what I find.
> >>>>>>>>
> >>>>>>>> Michael James
> >>>>>>>> Unidata
> >>>>>>>>
> >>>>>>>>
> >>>>>>>>> Hi GEMPAK Support,
> >>>>>>>>>
> >>>>>>>>> I am currently working on a project that requires me to use the 
> >>>>>>>>> GDPVSF
> >>>>>>>>> function in GEMPAK 6.7.0  to compute Dynamic Tropopause surfaces. 
> >>>>>>>>> When I
> >>>>>>>>> run GDPVSF, there is no output file created; however, the program 
> >>>>>>>>> does
> >>>>>>>>> not list any errors. I read from the below link from your support 
> >>>>>>>>> site
> >>>>>>>>> that GDPVSF has not been working properly due to changes from past
> >>>>>>>>> updates.
> >>>>>>>>>
> >>>>>>>>> http://www.unidata.ucar.edu/support/help/MailArchives/gempak/msg06131.html
> >>>>>>>>>
> >>>>>>>>> Do you know if the problems I am having are related to the same 
> >>>>>>>>> issues
> >>>>>>>>> described above with GDPVSF? What is the status of the GDPVSF 
> >>>>>>>>> function,
> >>>>>>>>> and do you know when a working version of it will be available? The
> >>>>>>>>> input parameters that I am using are listed at the bottom of this
> >>>>>>>>> message. These parameters originally came from a tutorial 
> >>>>>>>>> concerning the
> >>>>>>>>> function that I found from a university website (listed below); 
> >>>>>>>>> however,
> >>>>>>>>> I used a different input file.
> >>>>>>>>>
> >>>>>>>>> http://www.atmos.albany.edu/ftp/gempak/gdpvsf.readme
> >>>>>>>>>
> >>>>>>>>> Thanks,
> >>>>>>>>> Tim Lahmers
> >>>>>>>>>
> >>>>>>>>> GDFILE   =  2012081700_nam212.gem
> >>>>>>>>> GDOUTF   =  DT_test1.grd
> >>>>>>>>> GFUNC    =  mul(avor(obs),stap)
> >>>>>>>>> GDATTIM  =  last
> >>>>>>>>> GVCORD   =  thta
> >>>>>>>>> STARTL   =  260
> >>>>>>>>> STOPL    =  350
> >>>>>>>>> DESIRE   =  0.00000015
> >>>>>>>>> GDOUTL   =  15
> >>>>>>>>> GVOUTC   =  pvbl
> >>>>>>>>> GPACK    =
> >>>>>>>>> GLIST    =  uwnd;vwnd;pres
> >>>>>>>>> PMAX     =  700
> >>>>>>>>>
> >>>>>>>>>
> >>>>>>>>>
> >>>>>>>>
> >>>>>>>
> >>>>>>
> >>>>>> Ticket Details
> >>>>>> ===================
> >>>>>> Ticket ID: RJR-299468
> >>>>>> Department: Support GEMPAK
> >>>>>> Priority: Normal
> >>>>>> Status: Open
> >>>>>>
> >>>>>
> >>>>>
> >>>>>
> >>>>> Hi Michael,
> >>>>>
> >>>>> Thank you for update on the GDPVSF tool. I am having some more problems 
> >>>>> with the GDPVSF tool, and I think there may be a coding issue. The data 
> >>>>> that I am using are in pressure coordinates, so I have to use GDVINT to 
> >>>>> interpolate the dataset to theta coordinates. When I make this 
> >>>>> conversion and run GDPVSF, the output that is plotted tends to have 
> >>>>> gaps in the dataset. The amount of data that comes through seems to be 
> >>>>> dependent upon the range of theta values that I interpolate the 
> >>>>> pressure data to. I have tried to do this with both archived 
> >>>>> operational NAM data and NARR data that I converted from GRIB files. In 
> >>>>> this particular example, I downloaded the file 2011081500_nam212.gem 
> >>>>> from the following link. 
> >>>>> http://mtarchive.geol.iastate.edu/2011/08/15/gempak/model/ The file was 
> >>>>> re-saved with the .grd extension. The settings that I am using for 
> >>>>> GDVINT are as follows:
> >>>>>
> >>>>> GDFILE = 2011081500_nam212.grd
> >>>>> GDOUTF = dtnam_test.grd
> >>>>> GDATTIM = f00
> >>>>> GLEVEL = 280-380-5
> >>>>> GVCORD = pres/thta
> >>>>> MAXGRD = 4000
> >>>>> GAREA = 15;-130;50;-70
> >>>>> VCOORD  = pres
> >>>>>
> >>>>> Upon interpolating the pressure coordinate data to theta coordinates, I 
> >>>>> then run GDPVSF with the following settings:
> >>>>>
> >>>>> GDFILE   =  dtnam_test.grd
> >>>>> GDOUTF   =  dtnam_test.grd
> >>>>> GFUNC    =  abs(pvor(pres,wnd))
> >>>>> GDATTIM  =  f00
> >>>>> GVCORD   =  thta
> >>>>> STARTL   =  380
> >>>>> STOPL    =  280
> >>>>> DESIRE   =  0.00000015
> >>>>> GDOUTL   =  15
> >>>>> GVOUTC   =  pvab
> >>>>> GPACK    =
> >>>>> GLIST    =  thta;pres
> >>>>> PMAX     =  800
> >>>>>
> >>>>> This process produced a DT surface with holes at higher values of theta 
> >>>>> (see attached file). I also tried to rerun all of these tools with an 
> >>>>> increased range of potential temperature (280-420k). The end result 
> >>>>> yields a DT surface with  most of the values over the CONUS missing. Do 
> >>>>> you know what is causing this inconsistent output with these dynamic 
> >>>>> tropopause surfaces?
> >>>>>
> >>>>> Thanks,
> >>>>> Tim Lahmers
> >>>>>
> >>>>>
> >>>>>
> >>>>>
> >>>>> On Feb 5, 2013, at 9:55 AM, Unidata GEMPAK Support wrote:
> >>>>>
> >>>>>> Hi Tim,
> >>>>>>
> >>>>>> I've added the fix for GDPVSf in GEMPAK 6.8.0 which can be downloaded 
> >>>>>> at http://www.unidata.ucar.edu/downloads/gempak/ or from github at 
> >>>>>> https://github.com/Unidata/gempak/tree/master/gempak/source/programs/gd/gdpvsf
> >>>>>>
> >>>>>> if you don't want to reinstall 6.8.0, you can download the contents of 
> >>>>>> that gdpvsf direvctory and in $GEMPAK/source/programs/gd/gdpvsf/ run 
> >>>>>> "rm $OS_LIB/gdpvsf.a; make everything"
> >>>>>>
> >>>>>> -Michael
> >>>>>>
> >>>>>>> Quick update.
> >>>>>>>
> >>>>>>> I've got the program working on the hrcbob data set once again.  Yet 
> >>>>>>> to test it on the nam and your other GD input, but if that seems okay 
> >>>>>>> I'll send you the update tomorrow.
> >>>>>>>
> >>>>>>> Michael
> >>>>>>>
> >>>>>>>
> >>>>>>>> Hi Tim,
> >>>>>>>>
> >>>>>>>> GDPVSF is a mess right now, yes.  I'm trying to fix it this morning 
> >>>>>>>> and seeing some problems with how the grid files are opened and 
> >>>>>>>> read, likely due to not keeping up with changes in the grid 
> >>>>>>>> diagnostic library.  I'll let you know what I find.
> >>>>>>>>
> >>>>>>>> Michael James
> >>>>>>>> Unidata
> >>>>>>>>
> >>>>>>>>
> >>>>>>>>> Hi GEMPAK Support,
> >>>>>>>>>
> >>>>>>>>> I am currently working on a project that requires me to use the 
> >>>>>>>>> GDPVSF
> >>>>>>>>> function in GEMPAK 6.7.0  to compute Dynamic Tropopause surfaces. 
> >>>>>>>>> When I
> >>>>>>>>> run GDPVSF, there is no output file created; however, the program 
> >>>>>>>>> does
> >>>>>>>>> not list any errors. I read from the below link from your support 
> >>>>>>>>> site
> >>>>>>>>> that GDPVSF has not been working properly due to changes from past
> >>>>>>>>> updates.
> >>>>>>>>>
> >>>>>>>>> http://www.unidata.ucar.edu/support/help/MailArchives/gempak/msg06131.html
> >>>>>>>>>
> >>>>>>>>> Do you know if the problems I am having are related to the same 
> >>>>>>>>> issues
> >>>>>>>>> described above with GDPVSF? What is the status of the GDPVSF 
> >>>>>>>>> function,
> >>>>>>>>> and do you know when a working version of it will be available? The
> >>>>>>>>> input parameters that I am using are listed at the bottom of this
> >>>>>>>>> message. These parameters originally came from a tutorial 
> >>>>>>>>> concerning the
> >>>>>>>>> function that I found from a university website (listed below); 
> >>>>>>>>> however,
> >>>>>>>>> I used a different input file.
> >>>>>>>>>
> >>>>>>>>> http://www.atmos.albany.edu/ftp/gempak/gdpvsf.readme
> >>>>>>>>>
> >>>>>>>>> Thanks,
> >>>>>>>>> Tim Lahmers
> >>>>>>>>>
> >>>>>>>>> GDFILE   =  2012081700_nam212.gem
> >>>>>>>>> GDOUTF   =  DT_test1.grd
> >>>>>>>>> GFUNC    =  mul(avor(obs),stap)
> >>>>>>>>> GDATTIM  =  last
> >>>>>>>>> GVCORD   =  thta
> >>>>>>>>> STARTL   =  260
> >>>>>>>>> STOPL    =  350
> >>>>>>>>> DESIRE   =  0.00000015
> >>>>>>>>> GDOUTL   =  15
> >>>>>>>>> GVOUTC   =  pvbl
> >>>>>>>>> GPACK    =
> >>>>>>>>> GLIST    =  uwnd;vwnd;pres
> >>>>>>>>> PMAX     =  700
> >>>>>>>>>
> >>>>>>>>>
> >>>>>>>>>
> >>>>>>>>
> >>>>>>>
> >>>>>>
> >>>>>> Ticket Details
> >>>>>> ===================
> >>>>>> Ticket ID: RJR-299468
> >>>>>> Department: Support GEMPAK
> >>>>>> Priority: Normal
> >>>>>> Status: Open
> >>>>>>
> >>>>>
> >>>>>
> >>>>>
> >>>>> Hi Michael,
> >>>>>
> >>>>> Thank you for update on the GDPVSF tool. I am having some more problems 
> >>>>> with the GDPVSF tool, and I think there may be a coding issue. The data 
> >>>>> that I am using are in pressure coordinates, so I have to use GDVINT to 
> >>>>> interpolate the dataset to theta coordinates. When I make this 
> >>>>> conversion and run GDPVSF, the output that is plotted tends to have 
> >>>>> gaps in the dataset. The amount of data that comes through seems to be 
> >>>>> dependent upon the range of theta values that I interpolate the 
> >>>>> pressure data to. I have tried to do this with both archived 
> >>>>> operational NAM data and NARR data that I converted from GRIB files. In 
> >>>>> this particular example, I downloaded the file 2011081500_nam212.gem 
> >>>>> from the following link. 
> >>>>> http://mtarchive.geol.iastate.edu/2011/08/15/gempak/model/ The file was 
> >>>>> re-saved with the .grd extension. The settings that I am using for 
> >>>>> GDVINT are as follows:
> >>>>>
> >>>>> GDFILE = 2011081500_nam212.grd
> >>>>> GDOUTF = dtnam_test.grd
> >>>>> GDATTIM = f00
> >>>>> GLEVEL = 280-380-5
> >>>>> GVCORD = pres/thta
> >>>>> MAXGRD = 4000
> >>>>> GAREA = 15;-130;50;-70
> >>>>> VCOORD  = pres
> >>>>>
> >>>>> Upon interpolating the pressure coordinate data to theta coordinates, I 
> >>>>> then run GDPVSF with the following settings:
> >>>>>
> >>>>> GDFILE   =  dtnam_test.grd
> >>>>> GDOUTF   =  dtnam_test.grd
> >>>>> GFUNC    =  abs(pvor(pres,wnd))
> >>>>> GDATTIM  =  f00
> >>>>> GVCORD   =  thta
> >>>>> STARTL   =  380
> >>>>> STOPL    =  280
> >>>>> DESIRE   =  0.00000015
> >>>>> GDOUTL   =  15
> >>>>> GVOUTC   =  pvab
> >>>>> GPACK    =
> >>>>> GLIST    =  thta;pres
> >>>>> PMAX     =  800
> >>>>>
> >>>>> This process produced a DT surface with holes at higher values of theta 
> >>>>> (see attached file). I also tried to rerun all of these tools with an 
> >>>>> increased range of potential temperature (280-420k). The end result 
> >>>>> yields a DT surface with  most of the values over the CONUS missing. Do 
> >>>>> you know what is causing this inconsistent output with these dynamic 
> >>>>> tropopause surfaces?
> >>>>>
> >>>>> Thanks,
> >>>>> Tim Lahmers
> >>>>>
> >>>>>
> >>>>>
> >>>>>
> >>>>> On Feb 5, 2013, at 9:55 AM, Unidata GEMPAK Support wrote:
> >>>>>
> >>>>>> Hi Tim,
> >>>>>>
> >>>>>> I've added the fix for GDPVSf in GEMPAK 6.8.0 which can be downloaded 
> >>>>>> at http://www.unidata.ucar.edu/downloads/gempak/ or from github at 
> >>>>>> https://github.com/Unidata/gempak/tree/master/gempak/source/programs/gd/gdpvsf
> >>>>>>
> >>>>>> if you don't want to reinstall 6.8.0, you can download the contents of 
> >>>>>> that gdpvsf direvctory and in $GEMPAK/source/programs/gd/gdpvsf/ run 
> >>>>>> "rm $OS_LIB/gdpvsf.a; make everything"
> >>>>>>
> >>>>>> -Michael
> >>>>>>
> >>>>>>> Quick update.
> >>>>>>>
> >>>>>>> I've got the program working on the hrcbob data set once again.  Yet 
> >>>>>>> to test it on the nam and your other GD input, but if that seems okay 
> >>>>>>> I'll send you the update tomorrow.
> >>>>>>>
> >>>>>>> Michael
> >>>>>>>
> >>>>>>>
> >>>>>>>> Hi Tim,
> >>>>>>>>
> >>>>>>>> GDPVSF is a mess right now, yes.  I'm trying to fix it this morning 
> >>>>>>>> and seeing some problems with how the grid files are opened and 
> >>>>>>>> read, likely due to not keeping up with changes in the grid 
> >>>>>>>> diagnostic library.  I'll let you know what I find.
> >>>>>>>>
> >>>>>>>> Michael James
> >>>>>>>> Unidata
> >>>>>>>>
> >>>>>>>>
> >>>>>>>>> Hi GEMPAK Support,
> >>>>>>>>>
> >>>>>>>>> I am currently working on a project that requires me to use the 
> >>>>>>>>> GDPVSF
> >>>>>>>>> function in GEMPAK 6.7.0  to compute Dynamic Tropopause surfaces. 
> >>>>>>>>> When I
> >>>>>>>>> run GDPVSF, there is no output file created; however, the program 
> >>>>>>>>> does
> >>>>>>>>> not list any errors. I read from the below link from your support 
> >>>>>>>>> site
> >>>>>>>>> that GDPVSF has not been working properly due to changes from past
> >>>>>>>>> updates.
> >>>>>>>>>
> >>>>>>>>> http://www.unidata.ucar.edu/support/help/MailArchives/gempak/msg06131.html
> >>>>>>>>>
> >>>>>>>>> Do you know if the problems I am having are related to the same 
> >>>>>>>>> issues
> >>>>>>>>> described above with GDPVSF? What is the status of the GDPVSF 
> >>>>>>>>> function,
> >>>>>>>>> and do you know when a working version of it will be available? The
> >>>>>>>>> input parameters that I am using are listed at the bottom of this
> >>>>>>>>> message. These parameters originally came from a tutorial 
> >>>>>>>>> concerning the
> >>>>>>>>> function that I found from a university website (listed below); 
> >>>>>>>>> however,
> >>>>>>>>> I used a different input file.
> >>>>>>>>>
> >>>>>>>>> http://www.atmos.albany.edu/ftp/gempak/gdpvsf.readme
> >>>>>>>>>
> >>>>>>>>> Thanks,
> >>>>>>>>> Tim Lahmers
> >>>>>>>>>
> >>>>>>>>> GDFILE   =  2012081700_nam212.gem
> >>>>>>>>> GDOUTF   =  DT_test1.grd
> >>>>>>>>> GFUNC    =  mul(avor(obs),stap)
> >>>>>>>>> GDATTIM  =  last
> >>>>>>>>> GVCORD   =  thta
> >>>>>>>>> STARTL   =  260
> >>>>>>>>> STOPL    =  350
> >>>>>>>>> DESIRE   =  0.00000015
> >>>>>>>>> GDOUTL   =  15
> >>>>>>>>> GVOUTC   =  pvbl
> >>>>>>>>> GPACK    =
> >>>>>>>>> GLIST    =  uwnd;vwnd;pres
> >>>>>>>>> PMAX     =  700
> >>>>>>>>>
> >>>>>>>>>
> >>>>>>>>>
> >>>>>>>>
> >>>>>>>
> >>>>>>
> >>>>>> Ticket Details
> >>>>>> ===================
> >>>>>> Ticket ID: RJR-299468
> >>>>>> Department: Support GEMPAK
> >>>>>> Priority: Normal
> >>>>>> Status: Open
> >>>>>>
> >>>>>
> >>>>>
> >>>>>
> >>>>
> >>>
> >>>
> >>> Ticket Details
> >>> ===================
> >>> Ticket ID: RJR-299468
> >>> Department: Support GEMPAK
> >>> Priority: Normal
> >>> Status: Open
> >>>
> >>
> >>
> >
> >
> > Ticket Details
> > ===================
> > Ticket ID: RJR-299468
> > Department: Support GEMPAK
> > Priority: Normal
> > Status: Open
> >
> 
> 


Ticket Details
===================
Ticket ID: RJR-299468
Department: Support GEMPAK
Priority: Normal
Status: Open