[MITgcm-support] Tidal flow over topography

Tue May 8 11:29:04 EDT 2018

Yes you are correct. I misspelled it.  It would be *pkg/obcs/obcs_calc.F.*

I was seeing some papers on Internal Gravity Waves generation through the
mentioned barotropic flow.

I will try with whatever you said.

Thanks,
Subhajit Kar

On Tue, May 8, 2018 at 8:49 PM, Martin Losch <Martin.Losch at awi.de> wrote:

> 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
>
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.mitgcm.org/pipermail/mitgcm-support/attachments/20180508/deafdd4a/attachment.html>