Hi Klemens, Uh, long time ago since I used Qpack-1. And even longer time since I used the CONTABS option of Qpack-1.
I rememebr that you easily ended up with numerical problems with CONTABS. But you write that you get a dimension error. By the way, what is the error message? Right now I have nothing more to add. If you can send a (complete) test case that I can run, that is probably the simplest way to sort this out. Bye, Patrick On 2012-05-25 7:21 PM, Klemens Hocke wrote: > Dear all, > > I tried to use the Qpack1 option 3 for retrieval of continuum absorption: > > Q.CONTABS_DO = 3; > Q.CONTABS_ORDER = 1; > Q.CONTABS_LIMITS = [-1 -1]; > > Q.CONTABS_COVMAT = 'sx.continua.aa'; > Q.CONTABS_KGRID = 'grid_strato.aa' ; > > Q.CONTABS_REF_SPECIES = [ '"O2-MPM93", "H2O-MPM93", "N2-SelfContMPM93"']; > > > I selected Q.CONTABS_KGRID same as Q.SPECIES_KGRIDS. Is this okay? > > While the Sx matrix looks fine for contabs_do=2 , Sx is totally wrong > for option 3 (leading later to a dimension-error in qpcls). > > I guess, the problem is that I don't know how Q.CONTABS_COVMAT is defined. > Is it pressure in the first column, error of O2 in the second column , > error of H2O-MPM3 in the third column and so on? > > > sx.continua.aa from above looks at the moment : > > 2 > 2 1 > 0 > 0 > 2 3 > 6 1 0 > -4 1 0 > > Thank you and kind regards, > Klemens > > > > > Dr. Klemens Hocke > Institute of Applied Physics and Oeschger Centre for Climate Change Research > University of Bern > Sidlerstrasse 5 > CH-3012 Bern > E-mail: [email protected] <mailto:[email protected]> > Phone: +41 316318958 > Fax: +41 316313765 > Web: http://www.iap.unibe.ch > http://www.oeschger.unibe.ch > > > > _______________________________________________ > qpack mailing list > [email protected] > https://www.sat.ltu.se/mailman/listinfo/qpack _______________________________________________ qpack mailing list [email protected] https://www.sat.ltu.se/mailman/listinfo/qpack
