WRF precipitation data assimilation

J

janac4

Member
Hello everyone,
I'm familiar with WRF model, but now I would like to start with WRFDA. I have rain gauge hourly data and want to know the possible options and workflows to assimilate this information to adjust WRF model forecast.
In the User Guide there is a section called "Precipitation Data Assimilation in WRFDA 4DVAR", but is based on NCEP Stage IV radar and gauge precipitation data. I want to use measurements from other sources and parts of the world, is it possible?
I don't see a clear path to follow since it says: Users wishing to use their own observations instead of NCEP Stage IV are responsible for writing a Fortran main program and call subroutine writerainobs (in write_rainobs.f90) to generate their own precipitation data.
So any additional information would be very helpful!
 
Hello janac4,

Did you succeed?I am currently working on the exact same challenge. I also want to assimilate my own local rain gauge data using WRFDA 4DVAR instead of relying on the NCEP Stage IV data, and I'm looking at the same write_rainobs.f90 subroutine constraint.

Since it has been a while since you posted this, I was wondering if you were successful in getting this to work? If so, I would love to ask you a few questions:

How exactly did you process and format your rain gauge data to generate the custom precipitation data?

Did you assimilate any other types of meteorological data alongside your rainfall data?

I would be incredibly grateful if you could share your workflow, any scripts, or insights you figured out along the way.

Thanks in advance!
 
lygedasego said:
Hello janac4,

Did you succeed?I am currently working on the exact same challenge. I also want to assimilate my own local rain gauge data using WRFDA 4DVAR instead of relying on the NCEP Stage IV data, and I'm looking at the same write_rainobs.f90 subroutine constraint.

Since it has been a while since you posted this, I was wondering if you were successful in getting this to work? If so, I would love to ask you a few questions:

How exactly did you process and format your rain gauge data to generate the custom precipitation data?

Did you assimilate any other types of meteorological data alongside your rainfall data?

I would be incredibly grateful if you could share your workflow, any scripts, or insights you figured out along the way.

Thanks in advance!
Click to expand...
Hi,

unfortunately I haven’t managed to make any progress on this yet.
I basically got stuck at the same point you mention—the documentation is quite focused on NCEP Stage IV and it’s not very clear how to properly adapt it to external observations without diving into custom Fortran development.

So at the moment I don’t have a working solution or scripts to share.
 
You must log in or register to reply here.
Back
Top