[MITgcm-support] obcs in CPP_OPTIONS.h?, more questions

[Martin Losch]

OB_Jsouth = 21*1,104*0,
the OB_Jsouth(1:21)=21*1 does not work with all compilers.
Martin

On 28 May 2008, at 12:04, Renske Gelderloos wrote:

> Hi Patrick,
>
> Thanks for the help, this indeed works a lot better :-). I have  
> another question about the open boundary conditions. I'm sure it's  
> a very trivial thing, but I can't find a quick way. The problem is  
> the following: in data.obcs you can specify the gridpoints where  
> you want an open boundary condition. I would like, for example, an  
> OB on the first 21 gridpoints (out of 125) at the southern edge of  
> the domain. I would expect that I should type something like  
> OB_Jsouth(1:21)=1 or maybe OB_Jsouth(1:21)=20*1, but this doesn't  
> work. The only way I get this working is by individually typing  
> OB_Jsouth(1)=1, OB_Jsouth(2)=1, etc. I'm sure there must be a  
> quicker way. Can anyone tell me?
>
> Renske
>
>
> Patrick Heimbach wrote:
>>
>> HI there,
>>
>> since "obcs" is a package, it is enabled in
>> code/packages.conf
>> rather than in
>> code/CPP_OPTIONS.h
>> (see e.g. manual, section 6.1.1 Package Inclusion/Exclusion)
>> http://mitgcm.org/r2_manual/latest/online_documents/node226.html
>>
>> In order to avoid confusion,
>> this practice is "enforced" by having genmake2 stop if it finds
>> such a #define in a CPP file.
>>
>> You can verify this by checking the contents of
>> verification/plume_on_slope/code/
>> (which, by the way, in newer versions of the code has changed into
>> verification/tutorial_plume_on_slope/
>> and is documented in section 3.16.3
>> http://mitgcm.org/r2_manual/latest/online_documents/node152.html
>> and as I am typing and double-checking this, I see
>> that the manual has the wrong instructions, since it wrongly says
>> "Open boundaries are enabled with #define ALLOW_OBCS in code/ 
>> CPP_OPTIONS.h").
>> Our apologies...
>>
>> -Patrick
>>
>>
>>
>> On May 27, 2008, at 10:10 AM, Renske Gelderloos wrote:
>>
>>> Hi,
>>>
>>> I've got a question about open boundary conditions. In the  
>>> tutorial (3.14.3) it says that you should add
>>>
>>> #define ALLOW_OBCS
>>>
>>> in code/CPP_OPTIONS.h. However, if I do so I cannot compile  
>>> anymore and get the following error (in genmake2):
>>>
>>> Searching for *OPTIONS.h files  in order to warn about the  
>>> presence of "#define "-type statements that are no longer allowed:
>>> found CPP_OPTIONS="../code//CPP_OPTIONS.h"
>>> Error: In ../code//CPP_OPTIONS.h there is an illegal line:  
>>> #define ALLOW_OBCS
>>>
>>> Furthermore, in the CPP_OPTIONS.h in verification/plume_on_slope/ 
>>> code I cannot find this switch. Has it changed in any way? Or am  
>>> I missing something else? Thanks for any advice!
>>>
>>> Renske
>>>
>>>
>>>
>>> _______________________________________________
>>> MITgcm-support mailing list
>>> MITgcm-support at mitgcm.org
>>> http://mitgcm.org/mailman/listinfo/mitgcm-support
>>
>> ---
>> Patrick Heimbach | heimbach at mit.edu | http://www.mit.edu/~heimbach
>> MIT | EAPS 54-1518 | 77 Massachusetts Ave | Cambridge MA 02139 USA
>> FON +1-617-253-5259 | FAX +1-617-253-4464 | SKYPE patrick.heimbach
>>
>>
>> _______________________________________________
>> 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