problem in MM5 pregrid:date not found

Dear all,
The problem in regrid is as follow. Does anyone solve it? Do I have to change compiler PGI from 6.1 to 6? but 6.0 can not be used in Fedora Core 5.

REGRID/pregrid> ./pregrid.csh

Processing for SourceType = 3D



Starting ON84 processing for type 3D

cd /home/mm5/MM5V3/REGRID/pregrid/on84

ln -s /home/mm5/MM5V3/REGRID/pregrid/NCEP_ON84.9303 ON84FILEA

\



** pregrid version number: 0.28 **





HSTART, HEND = 1993-03-13_00:00:00 1993-03-14_00:00:00
Opening: ON84FILEA


Subroutine DATINT:
Interpolating 3-d files to fill in any missing data…


Looking for data at time 1993-03-13_00

Date not found: 1993-03-13_00:00:00
FORTRAN STOP
foreach: No match.

Hi zhangbl,

I’m not an expert in running pregrid, but it looks like you haven’t configured pregrid.csh to use the correct dates for your data file.

  • Mat

Thank you,mkcolg,
Such problem also occurs when reading other kind of data.
Some MM5 user also complained about this. Someone said it is a bug of PGI 6.1, because they said there was no such problen in PGI 6.
I am sure I have configured the date correctly.
I think PGI has problem when read in such dataset.
If you are expert in PGI, I think you can find the reseaon.
Thank you again.

Hi zhangbl,

I’m wondering if this is related to a problem caused by a program error in pregrid. See: HERE

If not, can you document the steps you’ve taken to get this error so we can recreate it here? I’d need all steps including other preprocessing steps (i.e. Terrain), flags you used to build pregridder, what changes you’ve made to pregid.csh, what data set you used (and where I can obtain it), and what OS your running on.

Thanks,
Mat

dear mkcolg,

This is an example of using Regrid of MM5, because this is not using the result of Terrain, so I think the problem lies in reading the data.
From the result of ./pregrid.csh, it has found the data, but could not read.

The Os is Fedora Core 5
the data (NCEP_ON84.9303) can be found here:

http://www.mmm.ucar.edu/mm5/On-Line-Tutorial/regrid/regrid_data.html

Many people suggest me to use 6.0.2, so I think there are pronlems in 6.1.6

the same error can also be found here:
http://mailman.ucar.edu/pipermail/mm5-users/2006/001007.html
http://mailman.ucar.edu/pipermail/mm5-users/2006/001008.html

Hi zhangbl,

I was able to successfully compile and run pregrid.csh with 6.1-6 once I fixed the source code error in the on84 version of pregridder. It’s only by luck that it works with the 6.0 release.

From my earlier post:

The Fortran character string ‘on84flnm’ has length of 9 and does not end with a blank or NULL. This is fine except that the copen routine in cio.c does expect a blank or NULL to end the string. I think the best solution is to append a ‘\0’ to the end of ‘on84flnm’ when COPEN is called. Also, if you look at other pregridder directories, such as grib.misc’s rd_grib.F at line 134, you’ll see that this is exactly what is done.

To fix, change line 31 to ron84_record.F

call COPEN(iunit, ifd, trim(on84flnm)//char(0), 1, ierr, -1)

Please give this a try.

Thanks,
Mat

Dear mkcolg,

You are great.
I finished compiling REGRID of MM5.

Thank you very much!

Maybe this will happen when reading other data.
Best Regards.

Hi zhangbl,

I’m glad this solution worked. Would your mind posting this to the UCAR mailing list? I’d rather not have the MM5 user community think the PGI 6.1 release is “buggy” when the real cause is a source code problem.

Thanks!
Mat

Hi Mat,

I have send e-mail to MM5 user.
Thank you again for your help

zhangbl

Hi Mat,

I’m having a similar problem with grib.misc/pregrid_grib.exe.

The data I’m using is from http://dss.ucar.edu/datasets/ds090.0/data/pgbf00-grb2d/

Parameters I used in pregrid.csh:

   set DataDir = /home/mark/Download/MM5V3/REGRID
   set SRC3D = GRIB  # Many GRIB-format datasets
   set InFiles = ( ${DataDir}/pgb.f* )
   set SRCSST = $SRC3D
   set InSST = ( )
   set SRCSNOW = $SRC3D
   set InSnow = ()
   set SRCSOIL = $SRC3D
   set InSoil = ()

   set VT3D = ( grib.misc/Vtable.NNRP3D )
   set VTSST = ( grib.misc/Vtable.NNRPSST )
   set VTSNOW = ( grib.misc/Vtable.NNRPSNOW )
   set VTSOIL = ( grib.misc/Vtable.NNRPSOIL )

Following is the error message:

HSTART, HEND = 2006-02-15_00:00:00    2006-03-15_00:00:00

###############################################################################

Inventory for date = ^@^@^@^@^@^@^@^@^@^@ ^@^@^@^@^@^@^@^@

PRES   SKINTEMP
-------------------------------------------------------------------------------
-------------------------------------------------------------------------------

###############################################################################

Inventory for date = ^@^@^@^@^@^@^@^@^@^@ ^@^@^@^@^@^@^@^@

PRES   SKINTEMP
-------------------------------------------------------------------------------
-------------------------------------------------------------------------------
 First pass done, doing a reprocess
 nfiles =             1
 ntimes =           113
 filedates = ^@^@^@^@^@^@^@^@^@^@^@^@^@:00:00.0000
 hstart, hend = 2006-02-15_00:00:00  2006-03-15_06:00:00
 filedates(n) = ^@^@^@^@^@^@^@^@^@^@^@^@^@:00:00.0000

**********
Subroutine DATINT:
  Interpolating 3-d files to fill in any missing data...
**********


  Looking for data at time 2006-02-15_00
 Date not found: 2006-02-15_00:00:00.0000
foreach: No match.

BTW, I’m running the model on Linux x86_64 (Fedora 11) with PGI 9.0.

Any help is appreciated.

Thanking you in advance.

Surya

Hi,

Could an error come from not having SST input files? If you run ON84 data set, you will see there are SST input files. I would check with MM5 folks on how to get SST data or configure pregrid.csh to get it right.

Hongyon