[MITgcm-support] Obcs problems
Dimitris Menemenlis
menemenlis at jpl.nasa.gov
Sat May 9 10:41:53 EDT 2009
Chuncheng, I have no experience with the useOBCSbalance flag.
I always balance offline.
There was some discussion in the support list back in December:
http://forge.csail.mit.edu/pipermail/mitgcm-support/2008-December/005841.html
Not clear if useOBCSbalance works or if it is broken.
Dimitris
On May 8, 2009, at 7:50 PM, ouc.edu.cn wrote:
> Hi,Dimitris,
> Thank you very much for your reply.I have solved the problem
> already,but there are still some problems that confuse me.MITgcm has
> the code for obcs balance(in obcs_calc.F),and the method is the same
> as the one you offered me.I just set useOBCSbalance=.TRUE. in
> data.obcs,however,it doesn't work.The code for obcs balance doesn't
> appear in obcs_balance.F,and the sea level keeps rising or
> decreasing with forward time step.So,What's wrong with this problem?
> should I activate some 'switches' in other files?
> Thank you very much.
> Best regards,
> Chuncheng
>
>
> 在2009-05-09,"Dimitris Menemenlis" <dmenemenlis at gmail.com> 写
> 道:
> >Chuncheng, I am so sorry, your message got buried and went
> >unanswered. If this is something that you are still having problems
> >with, let me know and I will try to answer your questions. Regards,
> >Dimitris
> >
> >On Dec 25, 2008, at 6:36 AM, ouc.edu.cn wrote:
> >
> >> Hi,Dimitris,
> >> Thank you for helping me in solving the problem of obcs
> >> balance.You offered me a program which is http://ecco2.jpl.nasa.gov/data1/arctic/run_template_cube81/mk_run_template.m
> >> However,I read the program,and understand the algorithm.But I'm
> >> confused with follows:
> >>
> >> for t=1:(nt+1), mydisp(t)
> >> for g=1:length(genBC)
> >> switch genBC{g}
> >> case 'E'
> >> tmp=readbin([pout 'OBEu_' nme '_' dim '.bin'],[ny nz],
> >> 1,'real*4',t-1);
> >> OBE(t)=sum(sum(tmp.*OBEmask));.......
> >>
> >>
> >> (1)So,what's "nt" and "mydisp"?And what's "pout" and "t-1" stand
> for?
> >>
> >> ...
> >> tmp(find(OBWmask))=tmp(find(OBWmask))+ ...
> >> mean(OBE(t)-OBS(t)+OBN(t)-OBW(t))/sum(sum(OBWmask));
> >> writebin([pout 'OBWu_' nme '_' dim '.balance'],tmp,
> >> 1,'real*4',t-1);
> >> ...
> >> (2) Here tmp is a 1D arrray,but we want to output file with the
> size
> >> of ny*nz.So what's the problem?
> >>
> >> Thank you very much.
> >> Best regards!
> >> Chuncheng
> >>
> >>
> >> [广告] 重奖 悬赏kfc3v3 球衣
> >
> >Dimitris Menemenlis <DMenemenlis at gmail.com>
> >5056 Oakwood Ave, La Canada, CA 91011-2450, USA
> >tel/fax: 818-790-6735; cell: 818-625-6498
> >
>
>
> 穿越地震带 纪念汶川地震一周年
More information about the MITgcm-support
mailing list