[MITgcm-support] Read theta(salt)&eta from outer files when the model restarted

[Martin Losch]

I don't think, that I understand what you want to do. Do you really  
want to re-initialize the model in on-going integrations (without  
stopping and restarting from a restart/pickup file)? That's not going  
to work to well, because with every re-initialization there's going to  
be shock that the model has to deal with -> I would expect  
instabiliities and many waves. But if you really want to try then  
you'll have to take the code from S/R ini_fields and subroutines  
called from within and place them somewhere within S/R forward_step  
(maybe load_fields_driver/external_fields_load) with the appropriate  
instruction, as to when the fields are to be read (see  
external_fields_load.F for an example).

Instead you could restore your model to some theta/salinity fields  
(but not eta), e.g. with the help of the rbcs package.

Good luck,

Martin

On Aug 14, 2009, at 5:22 AM, 名 徐 wrote:

> Dear Martin:
>     Thank you for your advice.Now another question,if I want to   
> replace the fields for theta(salt)&eta in model computing process  
> instead of doing it via reading&writing pickup files，so I needn't  
> stop&restart the model.How can I do OR which routine  can been added  
> the
> postprocessing subroutine.
>
>
> --- 09年8月13日，周四, Martin Losch <Martin.Losch at awi.de> 写 
> 道：
>
> 发件人: Martin Losch <Martin.Losch at awi.de>
> 主题: Re: [MITgcm-support] Read theta(salt)&eta from outer files  
> when the model restarted
> 收件人: mitgcm-support at mitgcm.org
> 日期: 2009年8月13日,周四,下午2:23
>
> the simplest way is to modify the pickup files: replace the fields  
> for theta(salt)&eta with with your versions. The tendency terms will  
> still be the old ones, but may that's not much of a problem.
>
> Otherwise you could try to change the code: reverse the order of  
> subroutine calls in ini_fields.F but I don't guarantee that this  
> will not give unexpected results (notice the double negative).
>
> Martin
>
> On Aug 13, 2009, at 4:38 AM, 名 徐 wrote:
>
> > Dear all:
> >      I want to read the data of three-dimensional theta(salt)&eta  
> from outer files when the model restarted.For example, the model has  
> the output files every three days,then the output T(S,Eta)have a  
> postprocessing to generate the new  three-dimensional  
> theta(salt)&eta data with other programes，the model read the new  
> theta(salt)&eta data as input data instead of the theta(salt)&eta  
> data store in pickup file when it restart. But other parameters  
> still read from pickup file.Who can tell me how can I do.
> >      Thanks!
> >
> >
> >
> >
> > 好玩贺卡等你发，邮箱贺卡全新上线！ 
> _______________________________________________
> > MITgcm-support mailing list
> > MITgcm-support at mitgcm.org
> > http://mitgcm.org/mailman/listinfo/mitgcm-support
>
>
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support
>
> 好玩贺卡等你发，邮箱贺卡全新上线！ 
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support