hirlam

Hi,

I am running harmonie in the souhern hemnisphere (42W to 9W and 72S to 78S) but the boundary files extracted from the archive are 180deg out from the specified domain, ie 142E to 166E.

Has anyone come across this problem before or have any hints for how I can resolve this issue?

Cheers,

Klara

I would start checking the domain you have defined by drawing some of the fields in your climatic file.

Anyway I think there is a variable LARGE_EC_BD that if it is set to yes' it will request a global mars field.

Javier

The climate files are created for the wrong domain (180 deg out from my specified domain) so the mars retrivals are using the climate file domain. So now the question becomes what is creating the wrong domain for the climate files. I think the problem is somewhere in gl but not sure where.

Klara,

Where can I find you experiment on ecgate?

Ulf

Klara,

Please try the change below in gl/grb/reglambert2.f and gl_grib_api/grb/reglambert2.f. It produces the right results, but I have to check it more carefully before I make it default. Remember to remove your old ECMWF boundaries from the disk before you restart your experiment.

Ulf

svn diff grb/reglambert2.f
Index: grb/reglambert2.f
===================================================================
--- grb/reglambert2.f (revision 10733)
+++ grb/reglambert2.f (working copy)
@@ -123,9 +123,9 @@
x=alon(i,j)
y=alat(i,j)
rho=sign(1._jprb,n)*sqrt(x*x+(rho0-y)*(rho0-y))
- theta=atan2d(x,(rho0-y))
+ theta=atan2d(x,sign(1._jprb,n)*(rho0-y))

- lon(i,j)=theta/n+lov
+ lon(i,j)=theta/abs(n)+lov
lat(i,j)=(2*atand((rearth*F/rho)**(1/n)))-90._jprb

enddo