Changeset 2 for Roms_tools


Ignore:
Timestamp:
05/04/11 15:27:05 (13 years ago)
Author:
cholod
Message:

patch_bry_20_09_2010.tar : Pb of array dim for bry forcing in parallel mode

Location:
Roms_tools/Roms_Agrif
Files:
6 edited

Legend:

Unmodified
Added
Removed

  • Roms_tools/Roms_Agrif/boundary.h

    r1 r2  
    11! 
    2 !====================================================================== 
    3 ! ROMS_AGRIF is a branch of ROMS developped at IRD and INRIA, in France 
    4 ! The two other branches from UCLA (Shchepetkin et al)  
    5 ! and Rutgers University (Arango et al) are under MIT/X style license. 
    6 ! ROMS_AGRIF specific routines (nesting) are under CeCILL-C license. 
    7 !  
    8 ! ROMS_AGRIF website : http://roms.mpl.ird.fr 
    9 !====================================================================== 
     2! $Id: boundary.h,v 1.4 2005/10/26 15:18:52 pmarches Exp $ 
    103! 
    114#ifndef ANA_BRY 
    … …  
    119112# if defined OBC_WEST || defined AGRIF_OBC_WEST 
    120113#  ifdef Z_FRC_BRY 
    121       real zetabry_west(0:Mm+1), zetabry_west_dt(0:Mm+1,2) 
     114      real zetabry_west(GLOBAL_1D_ARRAYETA), 
     115     &    zetabry_west_dt(GLOBAL_1D_ARRAYETA,2) 
    122116      common /bry_zeta_west/ zetabry_west, zetabry_west_dt 
    123117#  endif 
    124118#  ifdef M2_FRC_BRY 
    125       real ubarbry_west(0:Mm+1), ubarbry_west_dt(0:Mm+1,2), 
    126      &     vbarbry_west(0:Mm+1), vbarbry_west_dt(0:Mm+1,2) 
     119      real ubarbry_west(GLOBAL_1D_ARRAYETA), 
     120     &    ubarbry_west_dt(GLOBAL_1D_ARRAYETA,2) 
     121     &    ,vbarbry_west(GLOBAL_1D_ARRAYETA), 
     122     &    vbarbry_west_dt(GLOBAL_1D_ARRAYETA,2) 
    127123      common /bry_ubar_west/ ubarbry_west, ubarbry_west_dt, 
    128124     &                       vbarbry_west, vbarbry_west_dt 
    … …  
    130126#  ifdef SOLVE3D 
    131127#   ifdef M3_FRC_BRY 
    132       real ubry_west(0:Mm+1,N), ubry_west_dt(0:Mm+1,N,2), 
    133      &     vbry_west(0:Mm+1,N), vbry_west_dt(0:Mm+1,N,2) 
     128      real ubry_west(GLOBAL_1D_ARRAYETA,N), 
     129     &    ubry_west_dt(GLOBAL_1D_ARRAYETA,N,2) 
     130     &    ,vbry_west(GLOBAL_1D_ARRAYETA,N), 
     131     &    vbry_west_dt(GLOBAL_1D_ARRAYETA,N,2) 
    134132      common /bry_u_west/ ubry_west, ubry_west_dt, 
    135133     &                    vbry_west, vbry_west_dt 
    136134#   endif 
    137135#   ifdef T_FRC_BRY 
    138       real tbry_west(0:Mm+1,N,NT), tbry_west_dt(0:Mm+1,N,2,NT) 
     136      real tbry_west(GLOBAL_1D_ARRAYETA,N,NT), 
     137     &    tbry_west_dt(GLOBAL_1D_ARRAYETA,N,2,NT) 
    139138      common /bry_t_west/ tbry_west, tbry_west_dt 
    140139#   endif 
    … …  
    144143# if defined OBC_EAST || defined AGRIF_OBC_EAST 
    145144#  ifdef Z_FRC_BRY 
    146       real zetabry_east(0:Mm+1), zetabry_east_dt(0:Mm+1,2) 
     145      real zetabry_east(GLOBAL_1D_ARRAYETA), 
     146     &    zetabry_east_dt(GLOBAL_1D_ARRAYETA,2) 
    147147      common /bry_zeta_east/ zetabry_east, zetabry_east_dt 
    148148#  endif 
    149149#  ifdef M2_FRC_BRY 
    150       real ubarbry_east(0:Mm+1), ubarbry_east_dt(0:Mm+1,2), 
    151      &     vbarbry_east(0:Mm+1), vbarbry_east_dt(0:Mm+1,2) 
     150      real ubarbry_east(GLOBAL_1D_ARRAYETA), 
     151     &    ubarbry_east_dt(GLOBAL_1D_ARRAYETA,2) 
     152     &    ,vbarbry_east(GLOBAL_1D_ARRAYETA), 
     153     &    vbarbry_east_dt(GLOBAL_1D_ARRAYETA,2) 
    152154      common /bry_ubar_east/ ubarbry_east, ubarbry_east_dt, 
    153155     &                       vbarbry_east, vbarbry_east_dt 
    … …  
    155157#  ifdef SOLVE3D  
    156158#   ifdef M3_FRC_BRY 
    157       real ubry_east(0:Mm+1,N), ubry_east_dt(0:Mm+1,N,2), 
    158      &     vbry_east(0:Mm+1,N), vbry_east_dt(0:Mm+1,N,2) 
     159      real ubry_east(GLOBAL_1D_ARRAYETA,N), 
     160     &    ubry_east_dt(GLOBAL_1D_ARRAYETA,N,2) 
     161     &    ,vbry_east(GLOBAL_1D_ARRAYETA,N), 
     162     &    vbry_east_dt(GLOBAL_1D_ARRAYETA,N,2) 
    159163      common /bry_u_east/ ubry_east, ubry_east_dt, 
    160164     &                    vbry_east, vbry_east_dt 
    161165#   endif 
    162166#   ifdef T_FRC_BRY 
    163       real tbry_east(0:Mm+1,N,NT), tbry_east_dt(0:Mm+1,N,2,NT) 
     167      real tbry_east(GLOBAL_1D_ARRAYETA,N,NT), 
     168     &    tbry_east_dt(GLOBAL_1D_ARRAYETA,N,2,NT) 
    164169      common /bry_t_east/ tbry_east, tbry_east_dt 
    165170#   endif 
    … …  
    169174# if defined OBC_SOUTH || defined AGRIF_OBC_SOUTH 
    170175#  ifdef Z_FRC_BRY  
    171       real zetabry_south(0:Lm+1), zetabry_south_dt(0:Lm+1,2) 
     176      real zetabry_south(GLOBAL_1D_ARRAYXI), 
     177     &    zetabry_south_dt(GLOBAL_1D_ARRAYXI,2) 
    172178      common /bry_zeta_south/ zetabry_south, zetabry_south_dt 
    173179#  endif 
    174180#  ifdef M2_FRC_BRY 
    175       real ubarbry_south(0:Lm+1), ubarbry_south_dt(0:Lm+1,2), 
    176      &     vbarbry_south(0:Lm+1), vbarbry_south_dt(0:Lm+1,2) 
     181      real ubarbry_south(GLOBAL_1D_ARRAYXI), 
     182     &    ubarbry_south_dt(GLOBAL_1D_ARRAYXI,2) 
     183     &    ,vbarbry_south(GLOBAL_1D_ARRAYXI), 
     184     &    vbarbry_south_dt(GLOBAL_1D_ARRAYXI,2) 
    177185      common /bry_ubar_south/ ubarbry_south, ubarbry_south_dt, 
    178186     &                        vbarbry_south, vbarbry_south_dt 
    … …  
    180188#  ifdef SOLVE3D 
    181189#   ifdef M3_FRC_BRY 
    182       real ubry_south(0:Lm+1,N), ubry_south_dt(0:Lm+1,N,2), 
    183      &     vbry_south(0:Lm+1,N), vbry_south_dt(0:Lm+1,N,2) 
     190      real ubry_south(GLOBAL_1D_ARRAYXI,N), 
     191     &    ubry_south_dt(GLOBAL_1D_ARRAYXI,N,2) 
     192     &    ,vbry_south(GLOBAL_1D_ARRAYXI,N), 
     193     &    vbry_south_dt(GLOBAL_1D_ARRAYXI,N,2) 
    184194      common /bry_u_south/ ubry_south, ubry_south_dt, 
    185195     &                     vbry_south, vbry_south_dt 
    186196#   endif 
    187197#   ifdef T_FRC_BRY 
    188       real tbry_south(0:Lm+1,N,NT), tbry_south_dt(0:Lm+1,N,2,NT) 
     198      real tbry_south(GLOBAL_1D_ARRAYXI,N,NT), 
     199     &    tbry_south_dt(GLOBAL_1D_ARRAYXI,N,2,NT) 
    189200      common /bry_t_south/ tbry_south, tbry_south_dt 
    190201#   endif 
    … …  
    194205# if defined OBC_NORTH || defined AGRIF_OBC_NORTH 
    195206#  ifdef Z_FRC_BRY 
    196       real zetabry_north(0:Lm+1), zetabry_north_dt(0:Lm+1,2) 
     207      real zetabry_north(GLOBAL_1D_ARRAYXI), 
     208     &    zetabry_north_dt(GLOBAL_1D_ARRAYXI,2) 
    197209      common /bry_zeta_north/ zetabry_north, zetabry_north_dt 
    198210#  endif 
    199211#  ifdef M2_FRC_BRY 
    200       real ubarbry_north(0:Lm+1), ubarbry_north_dt(0:Lm+1,2), 
    201      &     vbarbry_north(0:Lm+1), vbarbry_north_dt(0:Lm+1,2) 
     212      real ubarbry_north(GLOBAL_1D_ARRAYXI), 
     213     &    ubarbry_north_dt(GLOBAL_1D_ARRAYXI,2) 
     214     &    ,vbarbry_north(GLOBAL_1D_ARRAYXI), 
     215     &    vbarbry_north_dt(GLOBAL_1D_ARRAYXI,2) 
    202216      common /bry_ubar_north/ ubarbry_north, ubarbry_north_dt, 
    203217     &                        vbarbry_north, vbarbry_north_dt 
    … …  
    205219#  ifdef SOLVE3D 
    206220#   ifdef M3_FRC_BRY 
    207       real ubry_north(0:Lm+1,N), ubry_north_dt(0:Lm+1,N,2), 
    208      &     vbry_north(0:Lm+1,N), vbry_north_dt(0:Lm+1,N,2) 
     221      real ubry_north(GLOBAL_1D_ARRAYXI,N), 
     222     &    ubry_north_dt(GLOBAL_1D_ARRAYXI,N,2) 
     223     &    ,vbry_north(GLOBAL_1D_ARRAYXI,N), 
     224     &    vbry_north_dt(GLOBAL_1D_ARRAYXI,N,2) 
    209225      common /bry_u_north/ ubry_north, ubry_north_dt, 
    210226     &                     vbry_north, vbry_north_dt 
    211227#   endif 
    212228#   ifdef T_FRC_BRY 
    213       real tbry_north(0:Lm+1,N,NT), tbry_north_dt(0:Lm+1,N,2,NT) 
     229      real tbry_north(GLOBAL_1D_ARRAYXI,N,NT), 
     230     &    tbry_north_dt(GLOBAL_1D_ARRAYXI,N,2,NT) 
    214231      common /bry_t_north/ tbry_north, tbry_north_dt 
    215232#   endif 

  • Roms_tools/Roms_Agrif/get_bry.F

    r1 r2  
    11! 
    2 !====================================================================== 
    3 ! ROMS_AGRIF is a branch of ROMS developped at IRD and INRIA, in France 
    4 ! The two other branches from UCLA (Shchepetkin et al)  
    5 ! and Rutgers University (Arango et al) are under MIT/X style license. 
    6 ! ROMS_AGRIF specific routines (nesting) are under CeCILL-C license. 
    7 !  
    8 ! ROMS_AGRIF website : http://roms.mpl.ird.fr 
    9 !====================================================================== 
     2! $Id: get_bry.F,v 1.3 2005/10/07 15:32:39 pmarches Exp $ 
    103! 
    114#include "cppdefs.h" 
    … …  
    229222#ifdef OBC_WEST 
    230223#  ifdef Z_FRC_BRY 
    231             ierr=nf_read_bry_EW (zetabry_west_dt(0,itbry), bry_id, 
    232      &                           zetabry_west_id, bry_rec, r2dvar) 
     224            ierr=nf_read_bry_EW (zetabry_west_dt(START_1D_ARRAYETA,itbry 
     225     &          ), bry_id,zetabry_west_id, bry_rec, r2dvar) 
    233226            if (ierr.ne.nf_noerr) write(stdout,2) 'zeta_west' 
    234227            ierr_all=ierr_all+ierr 
    235228#  endif 
    236229#  ifdef M2_FRC_BRY 
    237             ierr=nf_read_bry_EW (ubarbry_west_dt(0,itbry), bry_id, 
    238      &                           ubarbry_west_id, bry_rec, u2dvar) 
     230            ierr=nf_read_bry_EW (ubarbry_west_dt(START_1D_ARRAYETA,itbry 
     231     &          ), bry_id,ubarbry_west_id, bry_rec, u2dvar) 
    239232            if (ierr.ne.nf_noerr) write(stdout,2) 'ubar_west' 
    240233            ierr_all=ierr_all+ierr 
    241             ierr=nf_read_bry_EW (vbarbry_west_dt(0,itbry), bry_id, 
    242      &                           vbarbry_west_id, bry_rec, v2dvar) 
     234            ierr=nf_read_bry_EW (vbarbry_west_dt(START_1D_ARRAYETA,itbry 
     235     &          ), bry_id,vbarbry_west_id, bry_rec, v2dvar) 
    243236            if (ierr.ne.nf_noerr) write(stdout,2) 'vbar_west' 
    244237            ierr_all=ierr_all+ierr 
    … …  
    246239#  ifdef SOLVE3D 
    247240#   ifdef M3_FRC_BRY 
    248             ierr=nf_read_bry_EW (ubry_west_dt(0,1,itbry), bry_id, 
    249      &                           ubry_west_id,  bry_rec,  u3dvar) 
     241            ierr=nf_read_bry_EW (ubry_west_dt(START_1D_ARRAYETA,1,itbry) 
     242     &          , bry_id,ubry_west_id,  bry_rec,  u3dvar) 
    250243            if (ierr.ne.nf_noerr)  write(stdout,2)  'u_west' 
    251244            ierr_all=ierr_all+ierr 
    252             ierr=nf_read_bry_EW (vbry_west_dt(0,1,itbry), bry_id, 
    253      &                           vbry_west_id,  bry_rec,  v3dvar) 
     245            ierr=nf_read_bry_EW (vbry_west_dt(START_1D_ARRAYETA,1,itbry) 
     246     &          , bry_id,vbry_west_id,  bry_rec,  v3dvar) 
    254247            if (ierr.ne.nf_noerr)  write(stdout,2)  'v_west' 
    255248            ierr_all=ierr_all+ierr 
    … …  
    257250#   ifdef T_FRC_BRY 
    258251            do itrc=1,isalt 
    259              ierr=nf_read_bry_EW (tbry_west_dt(0,1,itbry,itrc),  
    260      &               bry_id,tbry_west_id(itrc),  bry_rec, r3dvar) 
     252              ierr=nf_read_bry_EW (tbry_west_dt(START_1D_ARRAYETA,1 
     253     &            ,itbry,itrc),bry_id,tbry_west_id(itrc),  bry_rec, 
     254     &            r3dvar) 
    261255              if (ierr.ne.nf_noerr)  write(stdout,2)  'trc_west' 
    262256              ierr_all=ierr_all+ierr 
    … …  
    268262# ifdef OBC_EAST 
    269263#  ifdef Z_FRC_BRY 
    270             ierr=nf_read_bry_EW (zetabry_east_dt(0,itbry), bry_id, 
    271      &                           zetabry_east_id, bry_rec, r2dvar) 
     264            ierr=nf_read_bry_EW (zetabry_east_dt(START_1D_ARRAYETA,itbry 
     265     &          ), bry_id,zetabry_east_id, bry_rec, r2dvar) 
    272266            if (ierr.ne.nf_noerr) write(stdout,2) 'zeta_east' 
    273267            ierr_all=ierr_all+ierr 
    274268#  endif 
    275269#  ifdef M2_FRC_BRY 
    276             ierr=nf_read_bry_EW (ubarbry_east_dt(0,itbry), bry_id, 
    277      &                           ubarbry_east_id, bry_rec, u2dvar) 
     270            ierr=nf_read_bry_EW (ubarbry_east_dt(START_1D_ARRAYETA,itbry 
     271     &          ), bry_id,ubarbry_east_id, bry_rec, u2dvar) 
    278272            if (ierr.ne.nf_noerr) write(stdout,2) 'ubar_east' 
    279273            ierr_all=ierr_all+ierr 
    280             ierr=nf_read_bry_EW (vbarbry_east_dt(0,itbry), bry_id, 
    281      &                           vbarbry_east_id, bry_rec, v2dvar) 
     274            ierr=nf_read_bry_EW (vbarbry_east_dt(START_1D_ARRAYETA,itbry 
     275     &          ), bry_id,vbarbry_east_id, bry_rec, v2dvar) 
    282276            if (ierr.ne.nf_noerr) write(stdout,2) 'vbar_east' 
    283277            ierr_all=ierr_all+ierr 
    … …  
    285279#  ifdef SOLVE3D 
    286280#   ifdef M3_FRC_BRY 
    287             ierr=nf_read_bry_EW (ubry_east_dt(0,1,itbry), bry_id, 
    288      &                           ubry_east_id,  bry_rec,  u3dvar) 
     281            ierr=nf_read_bry_EW (ubry_east_dt(START_1D_ARRAYETA,1,itbry) 
     282     &          , bry_id,ubry_east_id,  bry_rec,  u3dvar) 
    289283            if (ierr.ne.nf_noerr)  write(stdout,2)  'u_east' 
    290284            ierr_all=ierr_all+ierr 
    291             ierr=nf_read_bry_EW (vbry_east_dt(0,1,itbry), bry_id,  
    292      &                           vbry_east_id,  bry_rec,  v3dvar) 
     285            ierr=nf_read_bry_EW (vbry_east_dt(START_1D_ARRAYETA,1,itbry) 
     286     &          , bry_id,vbry_east_id,  bry_rec,  v3dvar) 
    293287            if (ierr.ne.nf_noerr)  write(stdout,2)  'v_east' 
    294288            ierr_all=ierr_all+ierr 
    … …  
    296290#   ifdef T_FRC_BRY 
    297291            do itrc=1,isalt 
    298               ierr=nf_read_bry_EW (tbry_east_dt(0,1,itbry,itrc),  
    299      &               bry_id,tbry_east_id(itrc),  bry_rec, r3dvar) 
     292              ierr=nf_read_bry_EW (tbry_east_dt(START_1D_ARRAYETA,1 
     293     &            ,itbry,itrc),bry_id,tbry_east_id(itrc),  bry_rec, 
     294     &            r3dvar) 
    300295              if (ierr.ne.nf_noerr)  write(stdout,2)  'trc_east' 
    301296              ierr_all=ierr_all+ierr 
    … …  
    307302# ifdef OBC_SOUTH 
    308303#  ifdef Z_FRC_BRY 
    309             ierr=nf_read_bry_NS (zetabry_south_dt(0,itbry), bry_id, 
    310      &                          zetabry_south_id, bry_rec,  r2dvar) 
     304            ierr=nf_read_bry_NS (zetabry_south_dt(START_1D_ARRAYXI,itbry 
     305     &          ), bry_id,zetabry_south_id, bry_rec,  r2dvar) 
    311306            if (ierr.ne.nf_noerr) write(stdout,2) 'zeta_south' 
    312307            ierr_all=ierr_all+ierr 
    313308#  endif 
    314309#  ifdef M2_FRC_BRY 
    315             ierr=nf_read_bry_NS (ubarbry_south_dt(0,itbry), bry_id, 
    316      &                           ubarbry_south_id, bry_rec, u2dvar) 
     310            ierr=nf_read_bry_NS (ubarbry_south_dt(START_1D_ARRAYXI,itbry 
     311     &          ), bry_id,ubarbry_south_id, bry_rec, u2dvar) 
    317312            if (ierr.ne.nf_noerr) write(stdout,2) 'ubar_south' 
    318313            ierr_all=ierr_all+ierr 
    319             ierr=nf_read_bry_NS (vbarbry_south_dt(0,itbry), bry_id, 
    320      &                           vbarbry_south_id, bry_rec, v2dvar) 
     314            ierr=nf_read_bry_NS (vbarbry_south_dt(START_1D_ARRAYXI,itbry 
     315     &          ), bry_id,vbarbry_south_id, bry_rec, v2dvar) 
    321316            if (ierr.ne.nf_noerr) write(stdout,2) 'vbar_south' 
    322317            ierr_all=ierr_all+ierr 
    … …  
    324319#  ifdef SOLVE3D 
    325320#   ifdef M3_FRC_BRY 
    326             ierr=nf_read_bry_NS (ubry_south_dt(0,1,itbry), bry_id, 
    327      &                           ubry_south_id,  bry_rec,  u3dvar) 
     321            ierr=nf_read_bry_NS (ubry_south_dt(START_1D_ARRAYXI,1,itbry) 
     322     &          , bry_id,ubry_south_id,  bry_rec,  u3dvar) 
    328323            if (ierr.ne.nf_noerr)  write(stdout,2)  'u_south' 
    329324            ierr_all=ierr_all+ierr 
    330325 
    331             ierr=nf_read_bry_NS (vbry_south_dt(0,1,itbry), bry_id, 
    332      &                           vbry_south_id,  bry_rec,  v3dvar) 
     326            ierr=nf_read_bry_NS (vbry_south_dt(START_1D_ARRAYXI,1,itbry) 
     327     &          , bry_id,vbry_south_id,  bry_rec,  v3dvar) 
    333328            if (ierr.ne.nf_noerr)  write(stdout,2)  'v_south' 
    334329            ierr_all=ierr_all+ierr 
    … …  
    336331#   ifdef T_FRC_BRY 
    337332            do itrc=1,isalt 
    338               ierr=nf_read_bry_NS (tbry_south_dt(0,1,itbry,itrc),  
    339      &               bry_id,tbry_south_id(itrc),  bry_rec, r3dvar) 
     333              ierr=nf_read_bry_NS (tbry_south_dt(START_1D_ARRAYXI,1 
     334     &            ,itbry,itrc),bry_id,tbry_south_id(itrc),  bry_rec, 
     335     &            r3dvar) 
    340336              if (ierr.ne.nf_noerr) write(stdout,2) 'trc_south' 
    341337              ierr_all=ierr_all+ierr 
    … …  
    347343# ifdef OBC_NORTH 
    348344#  ifdef Z_FRC_BRY 
    349             ierr=nf_read_bry_NS (zetabry_north_dt(0,itbry), bry_id, 
    350      &                           zetabry_north_id, bry_rec, r2dvar) 
     345            ierr=nf_read_bry_NS (zetabry_north_dt(START_1D_ARRAYXI,itbry 
     346     &          ), bry_id,zetabry_north_id, bry_rec, r2dvar) 
    351347            if (ierr.ne.nf_noerr) write(stdout,2) 'zeta_north' 
    352348            ierr_all=ierr_all+ierr 
    353349#  endif 
    354350#  ifdef M2_FRC_BRY 
    355             ierr=nf_read_bry_NS (ubarbry_north_dt(0,itbry), bry_id, 
    356      &                           ubarbry_north_id, bry_rec, u2dvar) 
     351            ierr=nf_read_bry_NS (ubarbry_north_dt(START_1D_ARRAYXI,itbry 
     352     &          ), bry_id,ubarbry_north_id, bry_rec, u2dvar) 
    357353            if (ierr.ne.nf_noerr) write(stdout,2) 'ubar_north' 
    358354            ierr_all=ierr_all+ierr 
    359             ierr=nf_read_bry_NS (vbarbry_north_dt(0,itbry), bry_id,  
    360      &                           vbarbry_north_id, bry_rec, v2dvar) 
     355            ierr=nf_read_bry_NS (vbarbry_north_dt(START_1D_ARRAYXI,itbry 
     356     &          ), bry_id,vbarbry_north_id, bry_rec, v2dvar) 
    361357            if (ierr.ne.nf_noerr) write(stdout,2) 'vbar_north' 
    362358            ierr_all=ierr_all+ierr 
    … …  
    364360#  ifdef SOLVE3D 
    365361#   ifdef M3_FRC_BRY 
    366             ierr=nf_read_bry_NS (ubry_north_dt(0,1,itbry), bry_id, 
    367      &                           ubry_north_id,  bry_rec,  u3dvar) 
     362            ierr=nf_read_bry_NS (ubry_north_dt(START_1D_ARRAYXI,1,itbry) 
     363     &          , bry_id,ubry_north_id,  bry_rec,  u3dvar) 
    368364            if (ierr.ne.nf_noerr)  write(stdout,2)  'u_north' 
    369365            ierr_all=ierr_all+ierr 
    370             ierr=nf_read_bry_NS (vbry_north_dt(0,1,itbry), bry_id, 
    371      &                           vbry_north_id,  bry_rec,  v3dvar) 
     366            ierr=nf_read_bry_NS (vbry_north_dt(START_1D_ARRAYXI,1,itbry) 
     367     &          , bry_id,vbry_north_id,  bry_rec,  v3dvar) 
    372368            if (ierr.ne.nf_noerr)  write(stdout,2)  'v_north' 
    373369            ierr_all=ierr_all+ierr 
    … …  
    375371#   ifdef T_FRC_BRY 
    376372            do itrc=1,isalt 
    377               ierr=nf_read_bry_NS (tbry_north_dt(0,1,itbry,itrc),  
    378      &               bry_id,tbry_north_id(itrc),  bry_rec, r3dvar) 
     373              ierr=nf_read_bry_NS (tbry_north_dt(START_1D_ARRAYXI,1 
     374     &            ,itbry,itrc),bry_id,tbry_north_id(itrc),  bry_rec, 
     375     &            r3dvar) 
    379376              if (ierr.ne.nf_noerr)  write(stdout,2)  'temp_north' 
    380377              ierr_all=ierr_all+ierr 

  • Roms_tools/Roms_Agrif/get_bry_bio.F

    r1 r2  
    11! 
    2 !====================================================================== 
    3 ! ROMS_AGRIF is a branch of ROMS developped at IRD and INRIA, in France 
    4 ! The two other branches from UCLA (Shchepetkin et al)  
    5 ! and Rutgers University (Arango et al) are under MIT/X style license. 
    6 ! ROMS_AGRIF specific routines (nesting) are under CeCILL-C license. 
    7 !  
    8 ! ROMS_AGRIF website : http://roms.mpl.ird.fr 
    9 !====================================================================== 
     2! $Id$ 
    103! 
    114#include "cppdefs.h" 
    … …  
    533526                  call set_cycle (bry_id, bry_tid(itrc),  
    534527     &            ntbry1(itrc),bry_cycle1(itrc), bry_ncycle1(itrc),  
    535      $                    bry_rec1(itrc)) 
     528     &                    bry_rec1(itrc)) 
    536529                  bry_time1(1,itrc)=-1.E+20 
    537530                  bry_time1(2,itrc)=-1.E+20 
    … …  
    555548        if(got_tbry(itrc)) then 
    556549          do while (bry_time1(itbry1(itrc),itrc).lt.time+dt .and.  
    557      $             ierr.eq.nf_noerr) 
     550     &             ierr.eq.nf_noerr) 
    558551            call advance_cycle (bry_cycle1(itrc),ntbry1(itrc), 
    559      $           bry_ncycle1(itrc),bry_rec1(itrc)) 
     552     &           bry_ncycle1(itrc),bry_rec1(itrc)) 
    560553            if (ierr.eq.nf_noerr) then 
    561554            ierr=nf_get_var1_FTYPE  
    562      $           (bry_id, bry_tid(itrc), bry_rec1(itrc), cff) 
     555     &           (bry_id, bry_tid(itrc), bry_rec1(itrc), cff) 
    563556             if (ierr.eq.nf_noerr) then 
    564557               itbry1(itrc)=min(3-itbry1(itrc),ntbry1(itrc)) 
    565558               bry_time1(itbry1(itrc),itrc)=cff*day2sec +  
    566      $              bry_cycle1(itrc)*bry_ncycle1(itrc) 
     559     &              bry_cycle1(itrc)*bry_ncycle1(itrc) 
    567560               ierr_all=0 
    568561# ifdef OBC_WEST 
    569562 
    570                ierr=nf_read_bry_EW (tbry_west_dt(0,1,itbry1(itrc),itrc), 
    571      &              bry_id,tbry_west_id(itrc),  bry_rec1(itrc), r3dvar) 
     563               ierr=nf_read_bry_EW (tbry_west_dt(START_1D_ARRAYETA,1 
     564     &             ,itbry1(itrc),itrc),bry_id,tbry_west_id(itrc), 
     565     &             bry_rec1(itrc), r3dvar) 
    572566               if (ierr.ne.nf_noerr)  write(stdout,2)  'trc_west' 
    573567               ierr_all=ierr_all+ierr 
    … …  
    575569# ifdef OBC_EAST 
    576570 
    577               ierr=nf_read_bry_EW (tbry_east_dt(0,1,itbry1(itrc),itrc),  
    578      &              bry_id,tbry_east_id(itrc),  bry_rec1(itrc), r3dvar) 
     571               ierr=nf_read_bry_EW (tbry_east_dt(START_1D_ARRAYETA,1 
     572     &             ,itbry1(itrc),itrc),bry_id,tbry_east_id(itrc), 
     573     &             bry_rec1(itrc), r3dvar) 
    579574               if (ierr.ne.nf_noerr)  write(stdout,2)  'trc_east' 
    580575               ierr_all=ierr_all+ierr 
    581576# endif 
    582577# ifdef OBC_SOUTH 
    583               ierr=nf_read_bry_NS (tbry_south_dt(0,1,itbry1(itrc),itrc), 
    584      &              bry_id,tbry_south_id(itrc),  bry_rec1(itrc), r3dvar) 
     578               ierr=nf_read_bry_NS (tbry_south_dt(START_1D_ARRAYXI,1 
     579     &             ,itbry1(itrc),itrc),bry_id,tbry_south_id(itrc), 
     580     &             bry_rec1(itrc), r3dvar) 
    585581               if (ierr.ne.nf_noerr)  write(stdout,2)  'trc_south' 
    586582               ierr_all=ierr_all+ierr 
    587583# endif 
    588584# ifdef OBC_NORTH 
    589               ierr=nf_read_bry_NS (tbry_north_dt(0,1,itbry1(itrc),itrc), 
    590      &              bry_id,tbry_north_id(itrc),  bry_rec1(itrc), r3dvar) 
     585               ierr=nf_read_bry_NS (tbry_north_dt( START_1D_ARRAYXI,1 
     586     &             ,itbry1(itrc),itrc),bry_id,tbry_north_id(itrc), 
     587     &             bry_rec1(itrc), r3dvar) 
    591588               if (ierr.ne.nf_noerr)  write(stdout,2)  'trc_north' 
    592589               ierr_all=ierr_all+ierr 

  • Roms_tools/Roms_Agrif/init_arrays.F

    r1 r2  
    11! 
    2 !====================================================================== 
    3 ! ROMS_AGRIF is a branch of ROMS developped at IRD and INRIA, in France 
    4 ! The two other branches from UCLA (Shchepetkin et al)  
    5 ! and Rutgers University (Arango et al) are under MIT/X style license. 
    6 ! ROMS_AGRIF specific routines (nesting) are under CeCILL-C license. 
    7 !  
    8 ! ROMS_AGRIF website : http://roms.mpl.ird.fr 
    9 !====================================================================== 
     2! $Id: init_arrays.F,v 1.14 2005/10/10 13:40:17 pmarches Exp $ 
    103! 
    114#include "cppdefs.h" 
    … …  
    610603      if (SOUTHERN_EDGE) then 
    611604#  ifdef Z_FRC_BRY 
    612         do i=JstrR,IendR 
     605        do i=IstrR,IendR 
    613606          zetabry_south(i)=init 
    614607        enddo 
    615608#  endif 
    616609#  ifdef M2_FRC_BRY 
    617         do i=JstrR,IendR 
     610        do i=IstrR,IendR 
    618611          ubarbry_south(i)=init 
    619612          vbarbry_south(i)=init 
    … …  
    622615#  if defined SOLVE3D && (defined M3_FRC_BRY || defined T_FRC_BRY) 
    623616        do k=1,N 
    624           do i=JstrR,IendR 
     617          do i=IstrR,IendR 
    625618#   ifdef M3_FRC_BRY 
    626619            ubry_south(i,k)=init 
    … …  
    641634      if (NORTHERN_EDGE) then 
    642635#  ifdef Z_FRC_BRY 
    643         do i=JstrR,IendR 
     636        do i=IstrR,IendR 
    644637          zetabry_north(i)=init 
    645638        enddo 
    646639#  endif 
    647640#  ifdef M2_FRC_BRY 
    648         do i=JstrR,IendR 
     641        do i=IstrR,IendR 
    649642          ubarbry_north(i)=init 
    650643          vbarbry_north(i)=init 
    … …  
    653646#  if defined SOLVE3D && (defined M3_FRC_BRY || defined T_FRC_BRY) 
    654647        do k=1,N 
    655           do i=JstrR,IendR 
     648          do i=IstrR,IendR 
    656649#   ifdef M3_FRC_BRY 
    657650            ubry_north(i,k)=init 

  • Roms_tools/Roms_Agrif/nf_read_bry.F

    r1 r2  
    11! 
    2 !====================================================================== 
    3 ! ROMS_AGRIF is a branch of ROMS developped at IRD and INRIA, in France 
    4 ! The two other branches from UCLA (Shchepetkin et al)  
    5 ! and Rutgers University (Arango et al) are under MIT/X style license. 
    6 ! ROMS_AGRIF specific routines (nesting) are under CeCILL-C license. 
    7 !  
    8 ! ROMS_AGRIF website : http://roms.mpl.ird.fr 
    9 !====================================================================== 
     2! $Id: nf_read_bry.F,v 1.2 2005/10/06 12:56:21 pmarches Exp $ 
    103! 
    114#ifndef EAST_WEST  
    … …  
    3326#include "param.h" 
    3427#ifdef EAST_WEST  
    35       real A(0:Mm+1,N) 
     28      real A(GLOBAL_1D_ARRAYETA,N) 
    3629#else 
    37       real A(0:Lm+1,N)  
     30      real A(GLOBAL_1D_ARRAYXI,N)  
    3831#endif 
    3932      integer ncid, varid, record, type,  i,j,k, shft, ierr 

  • Roms_tools/Roms_Agrif/set_global_definitions.h

    r1 r2  
    11! 
    2 !====================================================================== 
    3 ! ROMS_AGRIF is a branch of ROMS developped at IRD and INRIA, in France 
    4 ! The two other branches from UCLA (Shchepetkin et al)  
    5 ! and Rutgers University (Arango et al) are under MIT/X style license. 
    6 ! ROMS_AGRIF specific routines (nesting) are under CeCILL-C license. 
    7 !  
    8 ! ROMS_AGRIF website : http://roms.mpl.ird.fr 
    9 !====================================================================== 
     2! $Id: set_global_definitions.h,v 1.8 2005/10/10 13:40:19 pmarches Exp $ 
    103! 
    114/* This is "global_definitions.h": It contains a set of predetermined 
    … …  
    182175# define GLOBAL_1D_ARRAYETA -1:Mm+2+padd_E 
    183176# define START_2D_ARRAY -1,-1 
     177# define START_1D_ARRAYXI -1 
     178# define START_1D_ARRAYETA -1 
    184179#else 
    185180# ifdef EW_PERIODIC 
    186181#   define GLOBAL_1D_ARRAYXI -1:Lm+2+padd_X 
     182#   define START_1D_ARRAYXI -1 
    187183#  ifdef NS_PERIODIC 
    188184#   define GLOBAL_2D_ARRAY -1:Lm+2+padd_X,-1:Mm+2+padd_E 
    189185#   define GLOBAL_1D_ARRAYETA -1:Mm+2+padd_E 
    190186#   define START_2D_ARRAY -1,-1 
     187#   define START_1D_ARRAYETA -1 
    191188#  else 
    192189#   define GLOBAL_2D_ARRAY -1:Lm+2+padd_X,0:Mm+1+padd_E 
    193190#   define GLOBAL_1D_ARRAYETA 0:Mm+1+padd_E 
    194191#   define START_2D_ARRAY -1,0 
     192#   define START_1D_ARRAYETA 0 
    195193#  endif 
    196194# else 
    197195#   define GLOBAL_1D_ARRAYXI 0:Lm+1+padd_X 
     196#   define START_1D_ARRAYXI 0 
    198197#  ifdef NS_PERIODIC 
    199198#   define GLOBAL_2D_ARRAY 0:Lm+1+padd_X,-1:Mm+2+padd_E 
    200199#   define GLOBAL_1D_ARRAYETA -1:Mm+2+padd_E 
    201200#   define START_2D_ARRAY 0,-1 
     201#   define START_1D_ARRAYETA -1 
    202202#  else 
    203203#   define GLOBAL_2D_ARRAY 0:Lm+1+padd_X,0:Mm+1+padd_E 
    204204#   define GLOBAL_1D_ARRAYETA 0:Mm+1+padd_E 
    205205#   define START_2D_ARRAY 0,0 
     206#   define START_1D_ARRAYETA 0 
    206207#  endif 
    207208# endif 
Note: See TracChangeset for help on using the changeset viewer.