[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
[GEMPAK #RJR-299468]: GDPVSF Status
- Subject: [GEMPAK #RJR-299468]: GDPVSF Status
- Date: Fri, 22 Feb 2013 15:52:45 -0700
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