[MITgcm-support] Tidal flow over topography

[Martin Losch]

I don’t know what you are doing with pkg/mypackage. I was under the impression (but might be wrong) that this is a template for buildung your own package. All you need to do is modify obcs_calc.F either in place (pkg/obcs/obcs_calc.F) or make a copy in your “code” directory and edit that.

Your fortran compile error appears, because you are defining uVel etc. twice, once as a DUMMY argument in obcs_calc.F and once DYNVARS.h (have a look at obcs_calc.f to see this). 

Do you really want to use the dynamic variable etaN (which is updated each timestep) to compute u on the boundary? In that case you need to either pass etaN, just like uVel, etc as a DUMMY variable to obcs_calc.F, or rename all dummy variables (uVel -> uVelLoc or similar, etc.) within obcs_calc.F and then include DYNVARS.h.

Or do you want to prescribe the surface elevation along the boundary so that it is consistent with u0*sin(omega*t)? Then you need to asign values just like u0, omega etc. also to OBXeta or similar.

Martin

> On 8. May 2018, at 15:22, subhajit kar <subhajitkar19 at gmail.com> wrote:
> 
> Thank you Martin. The problem has now been resolved.
> 
> The problem I have been trying to model is tidal flow over a 3-D Gaussian topography.
> 
> I have been modeling tidal flow as -  u = u0*sin(omega*t) - c/H*etaN in the western boundary. Here, u0, omega, c and H are constant and etaN is the free-surface elevation.
> 
> To model this I have modified the obcs_calc.F in pkg/mypackage module.
> 
> However, etaN variable is defined in DYNVARS.h and while including this library I am getting errors 
> 
> 
> 
> Real*8  uVel (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
>                   1
> Error: COMMON attribute conflicts with DUMMY attribute in ‘uvel’ at (1)
> obcs_calc.f:3154:18:
> 
>        Real*8  vVel (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
>                   1
> Error: COMMON attribute conflicts with DUMMY attribute in ‘vvel’ at (1)
> obcs_calc.f:3155:18:
> 
>        Real*8  wVel (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
>                   1
> Error: COMMON attribute conflicts with DUMMY attribute in ‘wvel’ at (1)
> 
> 
> While not included, I am getting -  etaN has no implicit type.
> 
> 
> Is there anything wrong I am doing and also how to get rid of it.
> 
> Thanks,
> Subhajit Kar  
> 
> 
> On Tue, May 8, 2018 at 1:06 PM, Martin Losch <Martin.Losch at awi.de> wrote:
> Hi,
> you need to provide more details about what you actually did to the code in order to get specific help.
> 
> This error you could have googled. I found this in stackoverflow,
> https://stackoverflow.com/questions/26887689/error-non-numeric-character-in-statement-at-label-1
> 
> you compiler assumes fixed source form with file suffixes .f, but you start your code in column 4 (
> 
> C234567
>    LOGICAL OBCSisON
> 
> In fixed source form, the first 6 columns are served for special labels, comment characters, etc.
> 
> Martin
> 
> 
>> On 8. May 2018, at 08:59, subhajit kar <subhajitkar19 at gmail.com> wrote:
>> 
>> Hello MITgcm users,
>> 
>>   
>>         I am trying model a toy-problem of Internal Gravity Waves generation through a tidal flow.
>> 
>>         I have modified the obcs_calc.F in the pkg/mypackage. I am getting several errors - 
>> 
>> 
>> 
>>    obcs_calc.f:3008:3:
>> 
>>    LOGICAL OBCSisON
>>    1
>> Error: Non-numeric character in statement label at (1)
>> obcs_calc.f:3008:3:
>> 
>>    LOGICAL OBCSisON
>>    1
>> Error: Unclassifiable statement at (1)
>> 
>> Error: Non-numeric character in statement label at (1)
>> obcs_calc.f:3009:3:
>> 
>>    COMMON /OBCS_PACKAGE/ OBCSisON
>>    1
>> Error: Unclassifiable statement at (1)
>> obcs_calc.f:3033:1:
>> 
>>  & OB_Jnorth,OB_Jsouth,OB_Ieast,OB_Iwest,
>>  1
>> Error: Non-numeric character in statement label at (1)
>> 
>> ......
>> ......
>> ......
>> 
>> 
>> Does anyone know how to remove this errors?
>> 
>> Thanks for the help.
>> 
>> 
>> Regrads,
>> Subhajit Kar
>> 
>> _______________________________________________
>> MITgcm-support mailing list
>> MITgcm-support at mitgcm.org
>> http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support
> 
> 
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support
> 
> 
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support