Yes, it's me again. If you are reading that, obviously you would like to run REGRID using the Reanalysis GRIB files generated by LATS/GrADS couple. But I found a problem when doing that, regridder always give me back the message: "Extrapolating Point..." and my outputs looks crazy like a pool, just losting mass, and mass, and mass... I don't know if it is a bug, but I discovered that in pregrid phase the GRIB files were read fine, exception the DLAT (delta latitude) info. It was read with negative values... Here is the part of the program (REGRID/pregrid/grib.misc/rd_grib.F) where I found "my problem"... ^^^-----------^^^---------^^^---------^^^-------^^^-------^^^-------^^^------^^^ call build_hdate(hdate,year,month,day,hour,minute,second) call geth_newdate(hdate,hdate,3600*fcst) ! Store information about the grid on which the data is. ! This stuff gets stored in the MAP variable, as defined in module GRIDINFO if (ksec2(4).eq.0) then ! Lat/Lon grid map%igrid = 0 map%nx = infogrid(1) map%ny = infogrid(2) map%dx = ginfo(8) map%dy = ginfo(9) map%lat1 = ginfo(3) map%lon1 = ginfo(4) map%dy = ginfo(9) elseif (ksec2(4).eq.1) then ! Mercator grid map%igrid = 1 map%nx = infogrid(1) map%ny = infogrid(2) map%dx = ginfo(8) map%dy = ginfo(9) map%lat1 = ginfo(3) map%lon1 = ginfo(4) map%truelat1 = ginfo(11) elseif (ksec2(4).eq.3) then ! Lambert Conformal Grid ^^^-----------^^^---------^^^---------^^^-------^^^-------^^^-------^^^------^^^ And here is the fix of "my problem", just a "-1*" ahead ginfo(9) ^^^-----------^^^---------^^^---------^^^-------^^^-------^^^-------^^^------^^^ call build_hdate(hdate,year,month,day,hour,minute,second) call geth_newdate(hdate,hdate,3600*fcst) ! Store information about the grid on which the data is. ! This stuff gets stored in the MAP variable, as defined in module GRIDINFO if (ksec2(4).eq.0) then ! Lat/Lon grid map%igrid = 0 map%nx = infogrid(1) map%ny = infogrid(2) map%dx = ginfo(8) map%dy = ginfo(9) map%lat1 = ginfo(3) map%lon1 = ginfo(4) ! print *,'RMS::: ',map%lat1,map%dy map%dy = -1*ginfo(9) !RMS elseif (ksec2(4).eq.1) then ! Mercator grid map%igrid = 1 map%nx = infogrid(1) map%ny = infogrid(2) map%dx = ginfo(8) map%dy = ginfo(9) map%lat1 = ginfo(3) map%lon1 = ginfo(4) map%truelat1 = ginfo(11) elseif (ksec2(4).eq.3) then ! Lambert Conformal Grid ^^^-----------^^^---------^^^---------^^^-------^^^-------^^^-------^^^------^^^ This just happens in the level fields (3D and sea level pressure), so surface fields will do not pass by this one. I'm using grib.misc files as you noted, and here is my Vtables: Vtable.NCEP =========== GRIB | Level| Level| Level| REGRID | REGRID | REGRID | Code | Code | 1 | 2 | Name | Units | Description | -----+------+------+------+----------+----------+---------------------------------------------+ 11 | 100 | * | | T | K | Temperature | 33 | 100 | * | | U | m s{-1} | U | 34 | 100 | * | | V | m s{-1} | V | 52 | 100 | * | | RH | % | Relative Humidity | 7 | 100 | * | | HGT | m | Height | 11 | 105 | 2 | | T | K | Temperature | 51 | 105 | 2 | | SPECHUMD | kg kg{-1}| | | 105 | 2 | | RH | % | Relative Humidity | 33 | 105 | 10 | | U | m s{-1} | U | 34 | 105 | 10 | | V | m s{-1} | V | 2 | 102 | 0 | | PMSL | Pa | Sea-level Pressure | -----+------+------+------+----------+----------+---------------------------------------------+ Vtable.SST ========== GRIB | Level| Level| Level| REGRID | REGRID | REGRID | Code | Code | 1 | 2 | Name | Units | Description | -----+------+------+------+----------+----------+---------------------------------------------+ 11 | 1 | 0 | | SST | K | Sea-Surface Temperature | -----+------+------+------+----------+----------+---------------------------------------------+ I'm not sure how much people use this kind of data, but anyway, if one person cames to use that I'm going fell better because I spent one entire month until get this one.