[MITgcm-support] Obcs problems

[Martin Losch]

Don't forget to define OBCS_ALLOW_BALANCE in OBCS_OPTIONS.h

And it looks like I (or someone else) fixed the issue with the loop  
ranges that I talk about in that thread that Dimitris has pointed out.

Martin
On May 9, 2009, at 4:41 PM, Dimitris Menemenlis wrote:

> 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
>> >
>>
>>
>> 穿越地震带 纪念汶川地震一周年
>
>
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support