Problem using ana_m2obc

[fragoso]

Dear Rommers,

I´m trying to model a large and flat river in the Amazon basin. I have information of flow and velocity in the western boundary, wich I fit to an analytical expression and I´m trying to use ana_m2obc to force the model. When compiling the model I receive this error messages:

analytical.f90:304.39:

& IminS, ImaxS, JminS, JmaxS, &
1
Error: Symbol 'imaxs' at (1) has no IMPLICIT type
analytical.f90:304.32:

& IminS, ImaxS, JminS, JmaxS, &
1
Error: Symbol 'imins' at (1) has no IMPLICIT type
analytical.f90:304.53:

& IminS, ImaxS, JminS, JmaxS, &
1
Error: Symbol 'jmaxs' at (1) has no IMPLICIT type
analytical.f90:304.46:

& IminS, ImaxS, JminS, JmaxS, &
1
Error: Symbol 'jmins' at (1) has no IMPLICIT type

My .h file is like below:

#undef UV_ADV
#define UV_COR
#undef UV_C4ADVECTION
#define UV_LDRAG
#undef TS_C4HADVECTION
#undef TS_C4VADVECTION
#undef SALINITY
#undef NONLIN_EOS

#define UV_VIS2
#undef MIX_S_UV
#undef VISC_GRID

#undef TS_DIF2
#undef MIX_S_TS
#undef DIFF_GRID


/* model configuration options */
#undef SOLVE3D
#define MASKING
#undef DIAGNOSTICS_UV
#undef DIAGNOSTICS_TS
#define CURVGRID

/* vertical turbulent mixing scheme */


#undef MY25_MIXING

#ifdef MY25_MIXING
# define N2S2_HORAVG /* use if Large et al. (1994) interior closure */
# define KANTHA_CLAYSON /* use if Kantha and Clayson stability function */
#endif

#define ANA_INITIAL

#define ANA_SMFLUX
#define ANA_SRFLUX


/* lateral boundary conditions */

#define SOUTHERN WALL


#define EAST_FSCHAPMAN
#define EAST_M2RADIATION
#define EAST_TRADIATION


#define WEST_FSCHAPMAN
#define WEST_TRADIATION
#define WEST_M2CLAMPED

#define NORTH_FSCHAPMAN
#define NORTH_M2FLATHER
#define NORTH_TRADIATION

#define ANA_FSOBC
#define ANA_M2OBC
#define ANA_TOBC

Any hints?

Thank you.

Mauricio
PROOCEANO - BRASIL

[kate]

Did you start from one of the ana_m2obc.h example files? Did you mix-and-match between ROMS versions? Is it complaining about lines in ana_m2obc or in ana_m2obc_tile? Did you include "tile.h" as in the example?

See https://www.myroms.org/wiki/index.php/P ... ode_syntax for what IminS and friends are.

[fragoso]

Hi Kate. Thanks for the reply.

I´m running serial, so I guess it´s complaining about ana_m2obc.h.

I used ana_m2obc I downloaded along with ROMS code and just modified after the "else" part:

Code: Select all

#else
...

# ifdef WEST_M2OBC
      IF (WESTERN_EDGE) THEN
        DO j=JstrR,JendR
          BOUNDARY(ng)%ubar_west(j)=-0.65*(((j-1)*0.0133)-1)**2.+0.65
        END DO
        DO j=Jstr,JendR
          BOUNDARY(ng)%vbar_west(j)=0.0_r8
        END DO
      END IF
# endif  

...

Thanks again,
Mauricio.