[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