Changeset 12547
- Timestamp:
- 2020-03-13T11:12:08+01:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/trunk/src/OCE/BDY/bdydta.F90
r12396 r12547 92 92 INTEGER :: ii, ij, ik, igrd, ipl ! local integers 93 93 INTEGER, DIMENSION(jpbgrd) :: ilen1 94 INTEGER, DIMENSION(:), POINTER :: nblen, nblenrim ! short cuts95 94 TYPE(OBC_DATA) , POINTER :: dta_alias ! short cut 96 95 TYPE(FLD), DIMENSION(:), POINTER :: bf_alias … … 108 107 DO jbdy = 1, nb_bdy 109 108 ! 110 nblen => idx_bdy(jbdy)%nblen111 nblenrim => idx_bdy(jbdy)%nblenrim112 !113 109 IF( nn_dyn2d_dta(jbdy) == 0 ) THEN 114 ilen1(:) = nblen(:)115 110 IF( dta_bdy(jbdy)%lneed_ssh ) THEN 116 111 igrd = 1 117 DO ib = 1, i len1(igrd)112 DO ib = 1, idx_bdy(jbdy)%nblenrim(igrd) ! ssh is allocated and used only on the rim 118 113 ii = idx_bdy(jbdy)%nbi(ib,igrd) 119 114 ij = idx_bdy(jbdy)%nbj(ib,igrd) … … 123 118 IF( dta_bdy(jbdy)%lneed_dyn2d) THEN 124 119 igrd = 2 125 DO ib = 1, ilen1(igrd)120 DO ib = 1, SIZE(dta_bdy(jbdy)%u2d) ! u2d is used only on the rim except if ln_full_vel = T, see bdy_dta_init 126 121 ii = idx_bdy(jbdy)%nbi(ib,igrd) 127 122 ij = idx_bdy(jbdy)%nbj(ib,igrd) … … 129 124 END DO 130 125 igrd = 3 131 DO ib = 1, ilen1(igrd)126 DO ib = 1, SIZE(dta_bdy(jbdy)%v2d) ! v2d is used only on the rim except if ln_full_vel = T, see bdy_dta_init 132 127 ii = idx_bdy(jbdy)%nbi(ib,igrd) 133 128 ij = idx_bdy(jbdy)%nbj(ib,igrd) … … 138 133 ! 139 134 IF( nn_dyn3d_dta(jbdy) == 0 ) THEN 140 ilen1(:) = nblen(:)141 135 IF( dta_bdy(jbdy)%lneed_dyn3d ) THEN 142 136 igrd = 2 143 DO ib = 1, i len1(igrd)137 DO ib = 1, idx_bdy(jbdy)%nblen(igrd) 144 138 DO ik = 1, jpkm1 145 139 ii = idx_bdy(jbdy)%nbi(ib,igrd) … … 149 143 END DO 150 144 igrd = 3 151 DO ib = 1, i len1(igrd)145 DO ib = 1, idx_bdy(jbdy)%nblen(igrd) 152 146 DO ik = 1, jpkm1 153 147 ii = idx_bdy(jbdy)%nbi(ib,igrd) … … 160 154 161 155 IF( nn_tra_dta(jbdy) == 0 ) THEN 162 ilen1(:) = nblen(:)163 156 IF( dta_bdy(jbdy)%lneed_tra ) THEN 164 157 igrd = 1 165 DO ib = 1, i len1(igrd)158 DO ib = 1, idx_bdy(jbdy)%nblen(igrd) 166 159 DO ik = 1, jpkm1 167 160 ii = idx_bdy(jbdy)%nbi(ib,igrd) … … 176 169 #if defined key_si3 177 170 IF( nn_ice_dta(jbdy) == 0 ) THEN ! set ice to initial values 178 ilen1(:) = nblen(:)179 171 IF( dta_bdy(jbdy)%lneed_ice ) THEN 180 172 igrd = 1 181 173 DO jl = 1, jpl 182 DO ib = 1, i len1(igrd)174 DO ib = 1, idx_bdy(jbdy)%nblen(igrd) 183 175 ii = idx_bdy(jbdy)%nbi(ib,igrd) 184 176 ij = idx_bdy(jbdy)%nbj(ib,igrd) … … 338 330 DO jbdy = 1, nb_bdy ! Tidal component added in ts loop 339 331 IF ( nn_dyn2d_dta(jbdy) .GE. 2 ) THEN 340 nblen => idx_bdy(jbdy)%nblen 341 nblenrim => idx_bdy(jbdy)%nblenrim 342 IF( cn_dyn2d(jbdy) == 'frs' ) THEN ; ilen1(:)=nblen(:) 343 ELSE ; ilen1(:)=nblenrim(:) 332 IF( cn_dyn2d(jbdy) == 'frs' ) THEN ; ilen1(:)=idx_bdy(jbdy)%nblen(:) 333 ELSE ; ilen1(:)=idx_bdy(jbdy)%nblenrim(:) 344 334 ENDIF 345 335 IF ( dta_bdy(jbdy)%lneed_ssh ) dta_bdy_s(jbdy)%ssh(1:ilen1(1)) = dta_bdy(jbdy)%ssh(1:ilen1(1))
Note: See TracChangeset
for help on using the changeset viewer.