﻿id	summary	reporter	owner	description	type	status	priority	milestone	component	version	resolution	keywords	cc
433	DMI ROPP9.0 wishlist (BAROCLIM etc)	Ian Culverwell	Ian Culverwell	"Stig Syndergaard (DMI) writes:

{{{
Hi Ian,

Sorry for this rather lengthy email. For some time now, I wanted to 
update you on the dmi_trunk ROPP work here at DMI, but didn't want 
to spam you with every little thing. So you get everything in a big chunk! 
Hope you don't mind.


We finally got BAROCLIM into the dmi_trunk:

https://trac.romsaf.org/ropp/changeset/4527/ropp_src/branches/dev/Share/dmi_trunk_8.0

Helge did most of the work to implement BAROCLIM. On top of this I 
added a new 2-parm fit approach for the stat.opt. (the one I 
developed in the struct.uncert. project some years back). I did 
this because the current 2-parm fit approach in ROPP has some minor 
flaws (sometimes not fitting correctly because of too few points to 
fit to) that I believe would overshadow any improvement we might 
see with BAROCLIM over MSIS when using the 2-parm fit.

Both are implemented as options, so one can choose either MSIS or 
BAROCLIM (in the configuration file), in combination with old or 
new 2-parm fit (or old 1-parm fit).

For the new 2-parm fit approach there are a couple of new
subroutines:

ropp_pp_fit_model_refraction_new.f90
ropp_pp_search_model_refraction_new.f90

with '_new' following the precedence of Chris B.'s implementation of
ropp_fm_refrac_1d_new.f90 etc...

In the .cf files the new method is called 'regular' and the old one 
'convoluted'. This because the old search and fit method is in fact 
convoluted in that it is fitting within the search algorithm. The 
new approach does first the search, then the fit. I'd be happy to 
explain more details if necessary.


The work on the 2-parm fit was done here:

https://trac.romsaf.org/ropp/changeset/4526/ropp_src/branches/dev/Share/ssy_8.0

The work on BAROCLIM was done here:

https://trac.romsaf.org/ropp/changeset/4334/ropp_src/branches/dev/Share/hjs_6.1
https://trac.romsaf.org/ropp/changeset/4482/ropp_src/branches/dev/Share/hjs_6.1
https://trac.romsaf.org/ropp/changeset/4518/ropp_src/branches/dev/Share/ssy_8.0
https://trac.romsaf.org/ropp/changeset/4519/ropp_src/branches/dev/Share/ssy_8.0
https://trac.romsaf.org/ropp/changeset/4520/ropp_src/branches/dev/Share/ssy_8.0

In the process of all this I found a small bug that in some 
compilations (when I didn't compile from scratch after a merge) 
seemed to make havoc:

https://trac.romsaf.org/ropp/changeset/4528/ropp_src/branches/dev/Share/dmi_trunk_8.0

There are also a few older changes by Hans, Joe, and myself that dates back to
summer last year (but after the 'deadline' of ROPP 8.0 wishes):

Changes to Tdry calculation by Hans:

https://trac.romsaf.org/ropp/changeset/4287/ropp_src/branches/dev/Share/dmi_trunk_6.1

Changes to background covariance calculation by Joe:

https://trac.romsaf.org/ropp/changeset/4331/ropp_src/branches/dev/Share/dmi_trunk_7.1

Changes related to eum2ropp by myself:

https://trac.romsaf.org/ropp/changeset/4278/ropp_src/branches/dev/Share/dmi_trunk_6.1
https://trac.romsaf.org/ropp/changeset/4412/ropp_src/branches/dev/Share/dmi_trunk_7.1

Hope you will consider taking all this on board for ROPP 9.0. We  should discuss
at next ROPP-GG. I understand that there will also be testing and documentation
that needs to be done. Of course, we  have tested BAROCLIM and the new 2-parm
fit here at DMI with a few profiles and some stats, but I suppose you would also
want to do some final testing when implemented in ROPP 9.0. I'd be happy to help
with these things.
}}}

Seems fairy nuff.
"	enhancement	closed	normal	9.1	ROPP (all)	8.0	fixed		
