기상 기후 강의 노트

https://www.daleseo.com/python-enumerate/

https://truman.tistory.com/88

http://melonicedlatte.com/2021/03/14/181700.html

scp -r ./data/GDAS/201601 {username}@XXX.XXX.XXX.103:/data/GDAS

출처: https://yobbicorgi.tistory.com/3 [상엽’s Python 블로그]

Wrfchemi_ 가 해당날짜에 없을 수 있다.

해결

  http://mailman.ucar.edu/pipermail/wrf-users/2011/002486.html

아래와 같이 만들어 주면,

Ob too early 라는 메세지는 사라진다.  In rsl.out.0000

그러나, n =     21  unknown ob of type FM-88 SATOB 는 여전히 남는다. 여튼 결과는 ws10m 가 별 변화가 없다.

Sounds like you have 2 issues:

1. OBSGRID.EXE is not writing output at the times you want

2. WRF-Chem is bombing when nudging

I posted a bit on problem 1 - I was misinterpreting the DOCs (or the

DOCs are misleading, depending on your viewpoint).  See

http://mailman.ucar.edu/pipermail/wrf-users/2011/002424.html. 

Basically, set &share::interval_seconds to match your  input GRIB data

interval (must match the interval you ran METGRID.EXE with).  This is

21600 if your GRIB files come every 6 hours.

Then make sure your little_r format files contain data every

&record7::intf4d (from your email, it looks sorta like you want to nudge

every 3 hours).  The main point is that OBSGRID.EXE will look for a new

file every intf4d seconds, and will not look in a file that's already

been opened (from a previous interval).

I'm nearly certain that you must combine both your surface and upper-air

data into the same file, but that's pretty easy - you can just "cat"

them together. 

I suggest you fix problem 1 first, so you have non-empty OBS_DOMAIN1??

files for the in-between hours, before addressing the bombing of WRF.

Bart Brashers

From: wrf-users-bounces at ucar.edu [mailto:wrf-users-bounces at ucar.edu] On

Behalf Of Chris Klich

Sent: Sunday, October 30, 2011 12:31 PM

To: wrf-users at ucar.edu

Subject: [Wrf-users] Help with Obs Nudging Crash

Hi all, I have recently been trying to use obs-nudging with my WRF run.

I am currently running WRF-Chem version 3.2.1 due to a long-term project

and unable to upgrade to 3.3.

I currently have both upper air and surface data in little-r format

(separate however).  I convert these files into the proper naming format

for OBSGRID, run obsgrid and move the files to the /run directory,

concatenate all the OBS_DOMAIN files into OBS_DOMAIN101, and then run

real.exe and wrf.exe.  However, when OBSGRID runs, the output files are

every 3 hours, while my met_em files and observation files are every 6.

These 3 hour intermediate times (3, 9, 15, 21) are empty when the

OBS_DOMAIN and plotobs and qc files are outputted, while the normal

intervals of 0, 6, 12, 18 are several megabytes.  Even in this case, I

had tried concatenating all files to OBS_DOMAIN101 and running WRF.

However, it seems the entire run finishes only when I use just upper air

data specifically.  When I run using just surface obs or concatenating

surface to upper air, the run crashes, and it seems at a very random

time.  When using just surface it crashed after about 12:49 into the

run, while running with the concatenated upper air/surface, it crashed

after about 3 days, 12:49.  I can't seem to figure out what is causing

this.  The end of my rsl.error.0000 file looks like:

Timing for main (dt=135.00): time 2008-05-20_12:45:13 on domain   1:

1.08750 elapsed seconds.

Timing for main (dt=135.00): time 2008-05-20_12:47:28 on domain   1:

1.05240 elapsed seconds.

 OBS NUDGING: Reading new obs for time window TBACK =   12.125 TFORWD =

13.458

  for grid =  1

 OBS NUDGING:     1 previously read obs are now too old for the current

window a

 nd have been removed.

 ****** CALL IN4DOB AT KTAU =   326 AND XTIME =     767.48:  NSTA =

43 ****

 **

 ++++++CALL ERROB AT KTAU =   326 AND INEST =  1:  NSTA =    43 ++++++

  OBS NUDGING FOR IN,J,KTAU,XTIME,IVAR,IPL:  1 10   326   767.48  3  3

rindx= 5.

 3

  OBS NUDGING FOR IN,J,KTAU,XTIME,IVAR,IPL:  1 10   326   767.48  4  4

rindx= 5.

 3

  OBS NUDGING FOR IN,J,KTAU,XTIME,IVAR,IPL:  1 10   326   767.48  1  1

rindx= 5.

 3

  OBS NUDGING FOR IN,J,KTAU,XTIME,IVAR,IPL:  1 10   326   767.48  2  2

rindx= 5.

 3

Timing for main (dt=135.00): time 2008-05-20_12:49:43 on domain   1:

1.07960 elapsed seconds.

forrtl: error (78): process killed (SIGTERM)

Image              PC                Routine            Line

Source            

libc.so.6          00000033C34DD1E3  Unknown               Unknown

Unknown

libopen-pal.so.0   00002B47160653DD  Unknown               Unknown

Unknown

libopen-pal.so.0   00002B471606213D  Unknown               Unknown

Unknown

libopen-pal.so.0   00002B471605555C  Unknown               Unknown

Unknown

libmpi.so.0        00002B4715B43AA8  Unknown               Unknown

Unknown

libmpi.so.0        00002B4715B721AC  Unknown               Unknown

Unknown

wrf.exe            0000000001C64282  Unknown               Unknown

Unknown

wrf.exe            0000000000F373F4  Unknown               Unknown

Unknown

wrf.exe            00000000014F73B5  Unknown               Unknown

Unknown

wrf.exe            0000000000E8EC4F  Unknown               Unknown

Unknown

wrf.exe            0000000000D8DB6D  Unknown               Unknown

Unknown

wrf.exe            000000000052ECDF  Unknown               Unknown

Unknown

wrf.exe            00000000004BC733  Unknown               Unknown

Unknown

wrf.exe            00000000004BC6E7  Unknown               Unknown

Unknown

wrf.exe            00000000004BC67C  Unknown               Unknown

Unknown

libc.so.6          00000033C341EC9D  Unknown               Unknown

Unknown

wrf.exe            00000000004BC579  Unknown               Unknown

Unknown

In addition, when I use upper air data, I get thousands and thousands of

lines in my rsl.error.0000 file that read something like:

n=        9061 unknown ob of type FM-88 SATOB

Is this normal or is there an option I am missing where it will read

these correctly?

Any help for either of these problems would be greatly appreciated as I

need to figure this out ASAP.

I've also attached both namelist.input and namelist.oa for any help that

may provide.

Thank you.

--

Christopher Klich

Graduate Student

Florida State University

(908) 208-9743

Pasted from <http://mailman.ucar.edu/pipermail/wrf-users/2011/002486.html>

2021년 1월 18일 월요일

오후 4:08

d01 2010-05-16_00:00:00 med_initialdata_input: calling input_input

 i_parent_start from namelist.input file =             4

 i_parent_start from gridded input file  =             1

 j_parent_start from namelist.input file =             7

 j_parent_start from gridded input file  =             1

-------------- FATAL CALLED ---------------

FATAL CALLED FROM FILE:  <stdin>  LINE:     540

Nest start locations do not match: namelist.input vs gridded input file

-------------------------------------------

해결

Ncdump 를 이용해서, global attribute 를 수정해 줘야 한다.

http://nco.sourceforge.net/nco.html#ncatted-netCDF-Attribute-Editor

Ncatted 가 wrfinput_d0X 내의  global attributes  편집 가능하게 한다.  아래처럼

 ncatted -a GRID_ID,global,m,i,2 wrfinput_d02

 ncatted -a I_PARENT_START,global,m,i,4 wrfinput_d02

 ncatted -a J_PARENT_START,global,m,i,7 wrfinput_d02

 ncatted -a PARENT_GRID_RATIO,global,m,i,3 wrfinput_d02

 ncatted -a DT,global,m,f,40.0 wrfinput_d02

 ncatted -a PARENT_ID,global,m,i,1 wrfinput_d02

ncatted -a GRID_ID,global,m,i,3 wrfinput_d03

 ncatted -a I_PARENT_START,global,m,i,32 wrfinput_d03

 ncatted -a J_PARENT_START,global,m,i,20 wrfinput_d03

 ncatted -a PARENT_GRID_RATIO,global,m,i,3 wrfinput_d03

 ncatted -a DT,global,m,f,13.3333 wrfinput_d03

 ncatted -a PARENT_ID,global,m,i,2 wrfinput_d03

d01 2010-05-15_00:00:00 open_aux_u : opening ./VPRM_input/archive/VPRM_input_d01 for reading. DATASET DATASET=AUXINPUT15

d01 2010-05-15_00:00:00 calling wrf_open_for_read_begin in open_u_dataset

d01 2010-05-15_00:00:00 calling insub in open_u_dataset

            4  input_wrf: wrf_get_next_time current_date: 2010-05-15_00:00:00 Status =          -102

-------------- FATAL CALLED ---------------

FATAL CALLED FROM FILE:  <stdin>  LINE:     930

 ... Could not find matching time in input file

-------------------------------------------

application called MPI_Abort(MPI_COMM_WORLD, 1) - process 0

• Mediation_integrate.f90

 CASE ( AUXINPUT15_ALARM )

   CALL open_aux_u( grid, config_flags, stream, AUXINPUT15_ALARM,       &

                    config_flags%auxinput15_inname, grid%auxinput15_oid, &

                    input_auxinput15, ierr )

   CALL input_auxinput15 ( grid%auxinput15_oid, grid , config_flags , ierr )

 CASE ( AUXINPUT15_ALARM )

     IF ( grid%nframes(stream) >= config_flags%frames_per_auxinput15 ) THEN

       CALL close_dataset ( grid%auxinput15_oid , config_flags , "DATASET=AUXINPUT15" )

       grid%auxinput15_oid = 0

       grid%nframes(stream) = 0

     ENDIF

Pasted from <http://meteo.edu.vn/~trungnq/Download/Save_WRFCHEM_for_Sri/WRFV3/share/mediation_integrate.f90>

-------------- FATAL CALLED ---------------

FATAL CALLED FROM FILE:  <stdin>  LINE:     930

 ... Could not find matching time in input file wrfchemi_d01_2016-05-02_00:00:00

-------------------------------------------

application called MPI_Abort(MPI_COMM_WORLD, 1) - process 0

Sol>

namelist.input 내 auxinput5_interval_m = 360 라인을 comment out 시킨다. 이 라인이 존재하면, 동일한 wrfchemi_d01_2016-05-02_00:00:00 파일에서 6시간 간격의  emission 시간을 찾는다. 없애면 default 로 매 파일 마다 매시간 데이터를 불러온다.

• VPRM input 읽을 때 문제