[RegCNET] New NCEP files
Steve Hostetler
steve at coas.oregonstate.edu
Sun Feb 13 20:35:04 CET 2011
Can someone provide the 194712...200209 binary GISST file or the code
used to convert the ascii BADC files? I have the raw files from BADC
that I am using but for some reason my gridded binary files don't look
right even though it should be straight forward to convert them.
Ultimately I want to convert the binary to netCDF, which I will make
available if you are interested in it
Thanks
Steve
__________________________
Steve Hostetler
US Geological Survey
Department of Geosciences
Wilkinson Hall Room 104
Oregon State University
Corvallis, OR 97331
tel +1 541 737-8928
fac +1 541 737-1200
On 1/24/2011 8:17 AM, bixq wrote:
>
> Dear all:
>
> Thanks to Steve Hostetler of OSU for him find out the problem of RegCM
> related to NNRP1 data of year 2010 and his suggestion of solution to
> the problem.
>
> The following two versions of RegCM3 has been updated to fix the problem
> related to 2010 NNRP1 data and/or monthly OISST probelms.
>
> http://www.ictp.it/~pubregcm/RegCM3/regcm.tar.gz
>
> http://www.ictp.it/~pubregcm/RegCM3/regcm_beta_20090916.tar.gz
>
> Dr. Ha of Vietnam also mentioned the monthly OISST problem to me before,
> the reason is that CDC/NOAA occasionally change the order of fields in
> their netCDF monthly OISST files. This problem has been solved in the
> above two version.
>
> In RegCM4.0, there is also problem related to 2010 NNRP1 data, so version
> 4.0 should also be updated soon.
>
> Best regards !
> Xunqiang Bi
>
>
> On Thu, 20 Jan 2011, Steve Hostetler wrote:
>
>> In an attempt to save you from dealing with user panic/work:
>>
>> Beginning with 2010, NCEP (RI for now, RII is not updated yet) has
>> started writing some of their files (eg air.2010.nc and the weekly
>> Oi_sst) in a different order than previous files (eg air.2009.nc).
>> This causes errors in reading the offset, scale and values in the V3
>> SST_1DEG.f and V3/V4 ICBC.f codes because they are hardwired to read
>> a certain varid (5 in ICBC) 5. ( V4 mod_sst_1deg.f90 has been fixed
>> for this issue.) In both codes, calls to nf_inq_varid() need to be
>> made to get the correct varid in order to free them from being
>> hardwired. The ICBC code needs nf_inq_varid() in several places. In
>> addition, V3/V4 SST_1DEG.f/mod_1deg.f90 need to have a trap to stop
>> them in case they attempt to read past the number of records in the
>> OI_SST file. I just use:
>>
>> status=nf_inq_unlimdim(inet,recid)
>> status=nf_inq_dim(inet,recid,recname,nrecs)
>> if (it .gt. nrecs) stop '>>>> NUMBER OF SST TIME SLICES EXCEEDED'
>>
>> The netCDF error trap would work just as well. As it is, if it does
>> attempt to make SSTs for dates that exceed the SST data, all
>> temperatures are set to the last good values.
>>
>> Steve
>>
>> __________________________
>> Steve Hostetler
>> US Geological Survey
>> Department of Geosciences
>> Wilkinson Hall Room 104
>> Oregon State University
>> Corvallis, OR 97331
>> tel +1 541 737-8928
>> fac +1 541 737-1200
>>
>>
More information about the RegCNET
mailing list