[MITgcm-support] how to specify the vertical and horizontal diffusivity at every gridding point

[Martin Losch]

Guang,

for the vertical diffusivity, you set
#define ALLOW_3D_DIFFKR
in CPP_OPTIONS.h and recompile
then you use the parameter diffKrFile (in data PARM05, I guess) to  
specify the filename with the 3D diffusivities (just as initial  
conditions for, say, theta)
Unfortunately, there is no code for variable horizonal diffusivities  
(as far as I know). You'll have to add it yourself. See model/src/ 
ini_mixing.F and all other occurences of the above CPP flags to see  
how you can to define such a new field and then apply it to pkg/ 
generic_advdiff/gad_diff_x.F and gad_diff_y.F instead of diffKh

Martin

On Jan 9, 2010, at 3:50 AM, 张正光 wrote:

> Dear MITgcm support list,
>
>     I want to run a experiment which I can specify the vertical and  
> horizontal diffusivity at every gridding point,so I can investigate  
> the influence of the mixing hotpoint on the circulation.
>     I've got some clue on the internet, for example, the 3D  
> diffusivity file  "diffkrfile" and the imformation in the  
> CPP_OPTION.h:
>        C o Allow full 3D specification of vertical diffusivity
>        #undef ALLOW_3D_DIFFKR
>     but I can't figure out how to exactly use them, can you give me  
> a detailed instruction about how to specify the vertical and  
> horizontal diffusivity at every gridding point?
>     Thank you very much!
>     call me  Guang is ok.
>
>  Guang
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support