[Wrf-users] Meaning of fields in record8 section of namelist.obsproc ?

Sky Patrick patricksky852 at gmail.com
Fri May 15 22:36:58 MDT 2009


Hi Oscar,

For converting the observation data to LITTLE_R, we may need to know the
format of LITTLE_R first. And the format specification can found in the
following page:

http://www.mmm.ucar.edu/mm5/documents/MM5_tut_Web_notes/OA/OA.htm

Following is the thing is my observation that may need to take care for
using LITTLE_R with obsproc :
1. I am also new to Fortran, and see the there have Fortran I/O format
column in this article. This field is analogy to the C printf formatting
descriptor. e.g.

F20.5  = %20.5f
I10 = %10d
A40 = %40s

2. After realized the LITTLE_R format, another thing worth to mention is
that, the 'platform' field in the report header seem important for
obsproc.exe, as the following post said:

http://mailman.ucar.edu/pipermail/wrf-users/2008/001009.html

3. In the report header, the fields '
- id', 'name', 'platform', 'source' need to have left alignment.
- another fields need right alignment

For converting my data to LITTLE_R, i just write a program to scan my
observation file and convert it to LITTLE_R.

*************************************************************************************************
Back to my issue,
1. what is the meaning of NESTIX, NESTJX, NEXTI and NEXT fields graphically
in the namelist.obsproc ?
2. For specify the domain of time in namelist.obsproc, what is the rule to
fill the value of time_window_min, time_analysis and time_window_max ?
Becasse i have an observation in LITTLE_R format, the range of date of the
observation in this file is, let say 2005-02-28_15:00:00 to
2005-02-29_17:00:00, so i set

&record2
 time_window_min  = '2005-02-28_15:00:00',
 time_analysis    = '2005-02-29_16:00:00',
 time_window_max  = '2005-02-29_17:00:00',
/

but there have a line in command prompt shown the following:


Screwy NDATE: 2005-02-29_16:00:00
Screwy NDATE: 2005-02-29_17:00:00

 Error: The time window [2005-02-28_15:00:00, 2005-02-29_17:00:00] does not
menc 2005-02-29_16:00:00


Finally, no 3DVAR was generated. What is the problem of that, and is it have
another field need to set for the time domain ?

Thank you very much,

Regards,
Patrick

On Sat, May 16, 2009 at 2:12 AM, Oscar Chimborazo <ochimborazo at gmail.com>wrote:

> Hi Patrick...
>
> Sorry I don't know the answer to your question.. and I apologize for
> writting to you in this way ...I am Oscar Chimborazo and I am a new user of
> WRFDA... is it possible that you can tell me how can I convert my own
> observation dat to LITTLE_R format? I have been looking for that information
> but without luck. If you can help me with this I would appreciate it very
> much.
>
> Best Regards
>
> Oscar Chimborazo
>
> On Mon, May 11, 2009 at 10:21 AM, Sky Patrick <patricksky852 at gmail.com>wrote:
>
>> Hi,
>>
>> I am using the WRFDA package, and converted my own observation data to
>> LITTLE_R format and try to process it using the obsproc.exe. So i edit the
>> file namelist.obsproc to make some configuration. obsproc.exe can read my
>> LITTLE_R file, but the result at the bottom of this email was printed and no
>> file like 3dvar file was generated. From the runtime log message of
>> obsproc.exe, it shown all of the data are out of domain. I suspect that the
>> record8 setting of namelist.obsproc was set incorrectly. What is the meaning
>> of all of this field graphically ? Especially for the field, NESTIX, NESTJX,
>> NEXTI and NEXTJ, because i don't catch the meaning of this from the ARW user
>> guide. Thank very much ~
>>
>> Regards,
>> Patrick
>>
>>  READ NAMELIST FILE: namelist.obsproc
>>  ------------------
>> &RECORD1
>>  OBS_GTS_FILENAME="/home/user/data/test/output.ob
>> ",
>>  OBS_ERR_FILENAME="obserr.txt
>> ",
>>  FG_FORMAT="MM5",  /
>>
>> &RECORD2
>>  TIME_WINDOW_MIN="2005-03-27_23:00:00",
>>  TIME_ANALYSIS="2005-03-29_12:00:00",
>>  TIME_WINDOW_MAX="2005-03-30_23:59:59",  /
>>
>> &RECORD3
>>  MAX_NUMBER_OF_OBS=     400000,
>>  FATAL_IF_EXCEED_MAX_OBS=T,  /
>>
>> &RECORD4
>>  QC_TEST_VERT_CONSISTENCY=T,
>>  QC_TEST_CONVECTIVE_ADJ=T,
>>  QC_TEST_ABOVE_LID=T,
>>  REMOVE_ABOVE_LID=T,
>>  DOMAIN_CHECK_H=T,
>>  THINING_SATOB=F,
>>  THINING_SSMI=F,
>>  THINING_QSCAT=F,  /
>>
>> &RECORD5
>>  PRINT_GTS_READ=T,
>>  PRINT_GPSPW_READ=T,
>>  PRINT_RECOVERP=T,
>>  PRINT_DUPLICATE_LOC=T,
>>  PRINT_DUPLICATE_TIME=T,
>>  PRINT_RECOVERH=T,
>>  PRINT_QC_VERT=T,
>>  PRINT_QC_CONV=T,
>>  PRINT_QC_LID=T,
>>  PRINT_UNCOMPLETE=T,  /
>>
>> &RECORD6
>>  PTOP= 1000.00000    ,
>>  PS0= 100000.000    ,
>>  TS0=  290.00000    ,
>>  TLP=  50.000000    ,
>>  PIS0=  20000.000    ,
>>  TIS0=  215.00000    ,
>>  BASE_PRES= 100000.000    ,
>>  BASE_TEMP=  290.00000    ,
>>  BASE_LAPSE=  50.000000    ,
>>  BASE_STRAT_TEMP=  215.00000    ,
>>  BASE_TROPO_PRES=  20000.000    ,  /
>>
>> &RECORD7
>>  IPROJ=          1,
>>  PHIC=  22.263700    ,
>>  XLONC=  113.91190    ,
>>  TRUELAT1=  20.000000    ,
>>  TRUELAT2=  60.000000    ,
>>  MOAD_CEN_LAT=  22.263700    ,
>>  STANDARD_LON=  113.91190    ,  /
>>
>> &RECORD8
>>  IDD=          1,
>>  MAXNES=          1,
>>  NESTIX=         60,        200,        136,        181,        211,
>>   5*0          ,
>>  NESTJX=         90,        200,        181,        196,        211,
>>   5*0          ,
>>  DIS=  60.000000    , 10.0000000    ,  3.3000000    , 2*1.1000000      ,
>> 5*0.0000000      ,
>>
>>  NUMC= 2*1          ,          2,          3,          4, 5*0          ,
>>
>>  NESTI=          1,         40,         28,         35,         45,
>>   5*0          ,
>>  NESTJ=          1,         60,         25,         65,         55,
>>   5*0          ,  /
>>
>> &RECORD9
>>  PREPBUFR_OUTPUT_FILENAME="prepbufr_output_filename
>> ",
>>  PREPBUFR_TABLE_FILENAME="prepbufr_table_filename
>> ",
>>  OUTPUT_OB_FORMAT=          2,
>>  USE_FOR="3DVAR",
>>  NUM_SLOTS_PAST=          3,
>>  NUM_SLOTS_AHEAD=          3,
>>  WRITE_SYNOP=T,
>>  WRITE_SHIP=T,
>>  WRITE_METAR=T,
>>  WRITE_BUOY=T,
>>  WRITE_PILOT=T,
>>  WRITE_SOUND=T,
>>  WRITE_AMDAR=T,
>>  WRITE_SATEM=T,
>>  WRITE_SATOB=T,
>>  WRITE_AIREP=T,
>>  WRITE_GPSPW=T,
>>  WRITE_GPSZTD=T,
>>  WRITE_GPSREF=T,
>>  WRITE_GPSEPH=T,
>>  WRITE_SSMT1=T,
>>  WRITE_SSMT2=T,
>>  WRITE_SSMI=T,
>>  WRITE_TOVS=T,
>>  WRITE_QSCAT=T,
>>  WRITE_PROFL=T,
>>  WRITE_BOGUS=T,
>>  WRITE_AIRS=T,
>>  WRITE_TAMDAR=T,  /
>>
>>  COARSE GRID CENTER IS AT X =      0.0 KM AND Y =  -8863.1 KM.
>> DIS(1)=    60.000
>>
>> ------------------------------------------------------------------------------
>> READ GTS OBSERVATIONS IN FILE /home/user/data/test/output.ob
>>
>> Diagnostics in file obs_gts_read.diag
>>
>> IO_ERROR=-1     1 FM-35 TEMP                                 22.264N
>> 113.912E  19000101000000       0.000
>> Have reached end of observations file.
>>
>>
>> ------------------------------------------------------------------------------
>> GTS OBSERVATIONS READ:
>>
>>  SYNOP reports:     0
>>  SHIPS reports:     0
>>  BUOYS reports:     0
>>  BOGUS reports:     0
>>  METAR reports:     0
>>  PILOT reports:     0
>>  SOUND reports:  8640
>>  AMDAR reports:     0
>>  SATEM reports:     0
>>  SATOB reports:     0
>>  GPSPW reports:     0
>>  GPSZD reports:     0
>>  GPSRF reports:     0
>>  GPSEP reports:     0
>>  AIREP reports:     0
>> TAMDAR reports:     0
>>  SSMT1 reports:     0
>>  SSMT2 reports:     0
>>  SSMI  reports:     0
>>  TOVS  reports:     0
>>  QSCAT reports:     0
>>  PROFL reports:     0
>>  AIRST reports:     0
>>  OTHER reports:     0
>>  Total reports:  8640
>>
>> Number of observations read:              8640
>> Number of empty observations:                0
>> Number of observations out of domain:     8640
>> Number of observations for ingestion:        0
>>
>>
>> AIRCRAFT DATA: Total=     0  Above cut_height=     0
>> STOP 99999
>>
>> _______________________________________________
>> Wrf-users mailing list
>> Wrf-users at ucar.edu
>> http://mailman.ucar.edu/mailman/listinfo/wrf-users
>>
>>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/wrf-users/attachments/20090516/e38ac72c/attachment.html 


More information about the Wrf-users mailing list