gouravsaha1995
New member
Hi all,
This is Gourav here. I am trying to implement WRF observational nudging following the details given in this guide: https://www2.mmm.ucar.edu/wrf/users/docs/ObsNudgingGuide.pdf. I am trying to use local meteorological data (MSLP, wind speed, wind direction, temperature and humidity - which are the minimum data fields required to implement FDDA in WRF) by creating the desired input file OBS_DOMAIN101 directly using a python script and following the necessary formatting options and QC flags (as instructed in the PDF guide). After making the necessary changes to the namelist.input file, I am able to activate the observational nudging option and the desired lines are being printed in the log file showing that the file (OBS_DOMAIN101) is being read:
OBS NUDGING: Reading new obs for time window TBACK = -3.000 TFORWD = 3.000 for grid = 1
opening first fdda obs file, fonc=01 inest= 1
ifon= 1
Saving index of first ob after end of current time window for nest = 1:
ob index = 99, time of ob = 6.0000 hrs, end of time window = 3.0000 hrs
****** CALL IN4DOB AT KTAU = 0 AND XTIME = 0.00: NSTA = 98 ******
*****NUDGING INDIVIDUAL OBS ON MESH # 1 WITH RINXY = 0.100E+03 KM, RINSIG = 0.100E+00 AND TWINDO (HALF-PERIOD) = 0.180E+03 MIN
NUDGING IND. OBS WINDS WITH GIV = 0.500E-01
NUDGING IND. OBS TEMPERATURE WITH GIT = 0.500E-01
NUDGING IND. OBS MOISTURE WITH GIQ = 0.200E-01
ABOVE THE SURFACE LAYER, OBS NUDGING IS PERFORMED ON PRESSURE LEVELS,
WHERE RINXY VARIES LINEARLY FROM 0.100E+03 KM AT THE SURFACE TO 0.200E+03 KM AT 500.00 MB AND ABOVE
IN THE SURFACE LAYER, WXY IS A FUNCTION OF DPSMX = 75.00 MB WITH DCON = 0.133E+00 - SEE SUBROUTINE NUDOB
FOR EFFICIENCY, THE OBS NUDGING FREQUENCY FOR MESH # 1 IS 2 CGM TIMESTEPS
d01 2025-02-03_00:00:00 Input data is acceptable to use: wrfbdy_d01
Timing for processing lateral boundary for domain 1: 0.76873 elapsed seconds
Tile Strategy is not specified. Assuming 1D-Y
WRF TILE 1 IS 1 IE 253 JS 1 JE 313
WRF NUMBER OF TILES = 1
OBS NUDGING is requested on a total of 1 domain(s).
++++++CALL ERROB AT KTAU = 0 AND INEST = 1: NSTA = 98 ++++++
REPORTING OBS MASS-PT LOCS FOR NEST 1 AT XTIME= 0.0 MINUTES
FREQ= 1, MAX= 150 LOCS, NEWLY READ OBS ONLY, -999 => OBS OFF PROC
OBS# I J K OBS LAT OBS LON XLAT(I,J) XLONG(I,J) TIME(hrs) OBS STATION ID
1 189.560 85.519 1.000 8.033 93.533 8.033 93.533 3.00 43382
2 104.981 87.214 1.000 8.733 77.816 8.733 77.816 3.00 43376
3 100.289 88.798 1.000 9.000 76.933 9.000 76.933 3.00 43354
4 185.273 91.287 1.000 9.150 92.816 9.150 92.816 3.00 43367
5 98.429 91.931 1.000 9.567 76.566 9.567 76.566 3.00 43355
OBS NUDGING FOR IN,J,KTAU,XTIME,IVAR,IPL: 1 10 0 0.00 3 3 rindx= 4.8
OBS NUDGING FOR IN,J,KTAU,XTIME,IVAR,IPL: 1 10 0 0.00 4 4 rindx= 4.8
OBS NUDGING FOR IN,J,KTAU,XTIME,IVAR,IPL: 1 10 0 0.00 1 1 rindx= 4.8
OBS NUDGING FOR IN,J,KTAU,XTIME,IVAR,IPL: 1 10 0 0.00 2 2 rindx= 4.8
d01 2025-02-03_00:00:00 ----------------------------------------
d01 2025-02-03_00:00:00 W-DAMPING BEGINS AT W-COURANT NUMBER = 1.000000
d01 2025-02-03_00:00:00 ----------------------------------------
Timing for main: time 2025-02-03_00:01:30 on domain 1: 113.03083 elapsed seconds
However, when I am trying the generate the anomaly plots (output with data assimilation - output without data assimilation), the plots show zero values everywhere, proving the fact that even though the file is being opened and read, the data is not being assimilated properly by the model. I am sharing the OBS_DOMAIN101 file, my current namelist.input file and my rsl.out.0000 log file for better understanding. Can anyone please reply to this thread as early as possible and help me diagnose this issue?
With regards,
Gourav Saha
Email: gs27@iitbbs.ac.in
P.S I am trying to upload the files and it shows: "The upload failed because the temporary directory was missing. The site administrator will need to resolve this before any files can be uploaded"
This is Gourav here. I am trying to implement WRF observational nudging following the details given in this guide: https://www2.mmm.ucar.edu/wrf/users/docs/ObsNudgingGuide.pdf. I am trying to use local meteorological data (MSLP, wind speed, wind direction, temperature and humidity - which are the minimum data fields required to implement FDDA in WRF) by creating the desired input file OBS_DOMAIN101 directly using a python script and following the necessary formatting options and QC flags (as instructed in the PDF guide). After making the necessary changes to the namelist.input file, I am able to activate the observational nudging option and the desired lines are being printed in the log file showing that the file (OBS_DOMAIN101) is being read:
OBS NUDGING: Reading new obs for time window TBACK = -3.000 TFORWD = 3.000 for grid = 1
opening first fdda obs file, fonc=01 inest= 1
ifon= 1
Saving index of first ob after end of current time window for nest = 1:
ob index = 99, time of ob = 6.0000 hrs, end of time window = 3.0000 hrs
****** CALL IN4DOB AT KTAU = 0 AND XTIME = 0.00: NSTA = 98 ******
*****NUDGING INDIVIDUAL OBS ON MESH # 1 WITH RINXY = 0.100E+03 KM, RINSIG = 0.100E+00 AND TWINDO (HALF-PERIOD) = 0.180E+03 MIN
NUDGING IND. OBS WINDS WITH GIV = 0.500E-01
NUDGING IND. OBS TEMPERATURE WITH GIT = 0.500E-01
NUDGING IND. OBS MOISTURE WITH GIQ = 0.200E-01
ABOVE THE SURFACE LAYER, OBS NUDGING IS PERFORMED ON PRESSURE LEVELS,
WHERE RINXY VARIES LINEARLY FROM 0.100E+03 KM AT THE SURFACE TO 0.200E+03 KM AT 500.00 MB AND ABOVE
IN THE SURFACE LAYER, WXY IS A FUNCTION OF DPSMX = 75.00 MB WITH DCON = 0.133E+00 - SEE SUBROUTINE NUDOB
FOR EFFICIENCY, THE OBS NUDGING FREQUENCY FOR MESH # 1 IS 2 CGM TIMESTEPS
d01 2025-02-03_00:00:00 Input data is acceptable to use: wrfbdy_d01
Timing for processing lateral boundary for domain 1: 0.76873 elapsed seconds
Tile Strategy is not specified. Assuming 1D-Y
WRF TILE 1 IS 1 IE 253 JS 1 JE 313
WRF NUMBER OF TILES = 1
OBS NUDGING is requested on a total of 1 domain(s).
++++++CALL ERROB AT KTAU = 0 AND INEST = 1: NSTA = 98 ++++++
REPORTING OBS MASS-PT LOCS FOR NEST 1 AT XTIME= 0.0 MINUTES
FREQ= 1, MAX= 150 LOCS, NEWLY READ OBS ONLY, -999 => OBS OFF PROC
OBS# I J K OBS LAT OBS LON XLAT(I,J) XLONG(I,J) TIME(hrs) OBS STATION ID
1 189.560 85.519 1.000 8.033 93.533 8.033 93.533 3.00 43382
2 104.981 87.214 1.000 8.733 77.816 8.733 77.816 3.00 43376
3 100.289 88.798 1.000 9.000 76.933 9.000 76.933 3.00 43354
4 185.273 91.287 1.000 9.150 92.816 9.150 92.816 3.00 43367
5 98.429 91.931 1.000 9.567 76.566 9.567 76.566 3.00 43355
OBS NUDGING FOR IN,J,KTAU,XTIME,IVAR,IPL: 1 10 0 0.00 3 3 rindx= 4.8
OBS NUDGING FOR IN,J,KTAU,XTIME,IVAR,IPL: 1 10 0 0.00 4 4 rindx= 4.8
OBS NUDGING FOR IN,J,KTAU,XTIME,IVAR,IPL: 1 10 0 0.00 1 1 rindx= 4.8
OBS NUDGING FOR IN,J,KTAU,XTIME,IVAR,IPL: 1 10 0 0.00 2 2 rindx= 4.8
d01 2025-02-03_00:00:00 ----------------------------------------
d01 2025-02-03_00:00:00 W-DAMPING BEGINS AT W-COURANT NUMBER = 1.000000
d01 2025-02-03_00:00:00 ----------------------------------------
Timing for main: time 2025-02-03_00:01:30 on domain 1: 113.03083 elapsed seconds
However, when I am trying the generate the anomaly plots (output with data assimilation - output without data assimilation), the plots show zero values everywhere, proving the fact that even though the file is being opened and read, the data is not being assimilated properly by the model. I am sharing the OBS_DOMAIN101 file, my current namelist.input file and my rsl.out.0000 log file for better understanding. Can anyone please reply to this thread as early as possible and help me diagnose this issue?
With regards,
Gourav Saha
Email: gs27@iitbbs.ac.in
P.S I am trying to upload the files and it shows: "The upload failed because the temporary directory was missing. The site administrator will need to resolve this before any files can be uploaded"
Last edited: