Changeset 9788


Ignore:
Timestamp:
2018-06-13T12:12:50+02:00 (2 years ago)
Author:
jchanut
Message:

Fixes top and agrif compilation issues, #1972

Location:
NEMO/trunk/src
Files:
13 edited

Legend:

Unmodified
Added
Removed
  • NEMO/trunk/src/NST/agrif_top_interp.F90

    r9598 r9788  
    4848   END SUBROUTINE Agrif_trc 
    4949 
    50  
    5150   SUBROUTINE interptrn( ptab, i1, i2, j1, j2, k1, k2, n1, n2, before, nb, ndir ) 
    5251      !!---------------------------------------------------------------------- 
    53       !!                   *** ROUTINE interptrn *** 
     52      !!                  *** ROUTINE interptrn *** 
    5453      !!---------------------------------------------------------------------- 
    5554      REAL(wp), DIMENSION(i1:i2,j1:j2,k1:k2,n1:n2), INTENT(inout) ::   ptab 
     
    5756      LOGICAL                                     , INTENT(in   ) ::   before 
    5857      INTEGER                                     , INTENT(in   ) ::   nb , ndir 
    59       !! 
    60       INTEGER ::   ji, jj, jk, jn   ! dummy loop indices 
    61       INTEGER ::   imin, imax, jmin, jmax 
    62       LOGICAL ::   ll_west, ll_east, ll_north, ll_south 
     58      ! 
     59      INTEGER  ::   ji, jj, jk, jn, iref, jref, ibdy, jbdy   ! dummy loop indices 
     60      INTEGER  ::   imin, imax, jmin, jmax, N_in, N_out 
    6361      REAL(wp) ::   zrhox, z1, z2, z3, z4, z5, z6, z7 
    64       !!---------------------------------------------------------------------- 
    65       ! 
    66       INTEGER  ::   ji, jj, jk, jn   ! dummy loop indices 
    67       INTEGER  ::   imin, imax, jmin, jmax, N_in, N_out 
    68       REAL(wp) ::   zrhox , zalpha1, zalpha2, zalpha3 
    69       REAL(wp) ::   zalpha4, zalpha5, zalpha6, zalpha7 
    7062      LOGICAL :: western_side, eastern_side,northern_side,southern_side 
    7163      ! vertical interpolation: 
     
    7365      REAL(wp), DIMENSION(k1:k2,n1:n2-1) :: tabin 
    7466      REAL(wp), DIMENSION(k1:k2) :: h_in 
    75       REAL(wp), DIMENSION(1:jpk) :: h_out(1:jpk) 
    76       REAL(wp) :: h_diff, zrhoxy 
    77  
    78       zrhoxy = Agrif_rhox()*Agrif_rhoy() 
    79       IF (before) THEN          
    80          DO jn = 1,jpts 
     67      REAL(wp), DIMENSION(1:jpk) :: h_out 
     68      REAL(wp) :: h_diff 
     69 
     70      IF( before ) THEN          
     71         DO jn = 1,jptra 
    8172            DO jk=k1,k2 
    8273               DO jj=j1,j2 
     
    8475                       ptab(ji,jj,jk,jn) = trn(ji,jj,jk,jn) 
    8576                 END DO 
    86                END DO 
    87             END DO 
    88          END DO 
     77              END DO 
     78           END DO 
     79        END DO 
     80 
    8981# if defined key_vertical 
    9082        DO jk=k1,k2 
    9183           DO jj=j1,j2 
    9284              DO ji=i1,i2 
    93                  ptab(ji,jj,jk,jpts+1) = tmask(ji,jj,jk) * e3t_n(ji,jj,jk)  
     85                 ptab(ji,jj,jk,jptra+1) = tmask(ji,jj,jk) * e3t_n(ji,jj,jk)  
    9486              END DO 
    9587           END DO 
    9688        END DO 
    9789# endif 
    98  
    99       ELSE 
    100          ! 
    101          western_side  = (nb == 1).AND.(ndir == 1) 
    102          eastern_side  = (nb == 1).AND.(ndir == 2) 
    103          southern_side = (nb == 2).AND.(ndir == 1) 
    104          northern_side = (nb == 2).AND.(ndir == 2) 
     90      ELSE  
     91 
     92         western_side  = (nb == 1).AND.(ndir == 1)   ;   eastern_side  = (nb == 1).AND.(ndir == 2) 
     93         southern_side = (nb == 2).AND.(ndir == 1)   ;   northern_side = (nb == 2).AND.(ndir == 2) 
    10594 
    10695# if defined key_vertical               
     
    127116               ENDDO 
    128117               IF (N_in > 0) THEN 
    129                   h_diff = sum(h_out(1:N_out))-sum(h_in(1:N_in)) 
    130118                  DO jn=1,jptra 
    131119                     call reconstructandremap(tabin(1:N_in,jn),h_in,ptab_child(ji,jj,1:N_out,jn),h_out,N_in,N_out) 
     
    137125         ptab_child(i1:i2,j1:j2,1:jpk,1:jptra) = ptab(i1:i2,j1:j2,1:jpk,1:jptra) 
    138126# endif 
    139  
    140127         ! 
    141          zrhox = Agrif_Rhox() 
    142          !  
    143          zalpha1 = ( zrhox - 1. ) * 0.5 
    144          zalpha2 = 1. - zalpha1 
    145          !  
    146          zalpha3 = ( zrhox - 1. ) / ( zrhox + 1. ) 
    147          zalpha4 = 1. - zalpha3 
    148          !  
    149          zalpha6 = 2. * ( zrhox - 1. ) / ( zrhox + 1. ) 
    150          zalpha7 =    - ( zrhox - 1. ) / ( zrhox + 3. ) 
    151          zalpha5 = 1. - zalpha6 - zalpha7 
    152          ! 
    153          imin = i1 
    154          imax = i2 
    155          jmin = j1 
    156          jmax = j2 
    157          !  
    158          ! Remove CORNERS 
    159          IF((nbondj == -1).OR.(nbondj == 2)) jmin = 3 
    160          IF((nbondj == +1).OR.(nbondj == 2)) jmax = nlcj-2 
    161          IF((nbondi == -1).OR.(nbondi == 2)) imin = 3 
    162          IF((nbondi == +1).OR.(nbondi == 2)) imax = nlci-2         
    163          ! 
    164          IF( eastern_side) THEN 
    165             DO jn = 1, jptra 
    166                tra(nlci,j1:j2,1:jpk,jn) = zalpha1 * ptab_child(nlci,j1:j2,1:jpk,jn) + zalpha2 * ptab_child(nlci-1,j1:j2,1:jpk,jn) 
    167                DO jk = 1, jpkm1 
    168                   DO jj = jmin,jmax 
    169                      IF( umask(nlci-2,jj,jk) == 0.e0 ) THEN 
    170                         tra(nlci-1,jj,jk,jn) = tra(nlci,jj,jk,jn) * tmask(nlci-1,jj,jk) 
    171                      ELSE 
    172                         tra(nlci-1,jj,jk,jn)=(zalpha4*tra(nlci,jj,jk,jn)+zalpha3*tra(nlci-2,jj,jk,jn))*tmask(nlci-1,jj,jk) 
    173                         IF( un(nlci-2,jj,jk) > 0.e0 ) THEN 
    174                            tra(nlci-1,jj,jk,jn)=( zalpha6*tra(nlci-2,jj,jk,jn)+zalpha5*tra(nlci,jj,jk,jn) &  
    175                                  + zalpha7*tra(nlci-3,jj,jk,jn) ) * tmask(nlci-1,jj,jk) 
    176                         ENDIF 
    177                      END DO 
    178                   END DO 
    179                END DO 
    180             ENDDO 
     128         DO jn=1, jptra 
     129            tra(i1:i2,j1:j2,1:jpk,jn)=ptab_child(i1:i2,j1:j2,1:jpk,jn)*tmask(i1:i2,j1:j2,1:jpk)  
     130         END DO 
     131 
     132         IF ( .NOT.lk_agrif_clp ) THEN  
     133            ! 
     134            zrhox = Agrif_Rhox() 
     135            z1 = ( zrhox - 1. ) * 0.5 
     136            z3 = ( zrhox - 1. ) / ( zrhox + 1. ) 
     137            z6 = 2. * ( zrhox - 1. ) / ( zrhox + 1. ) 
     138            z7 =    - ( zrhox - 1. ) / ( zrhox + 3. ) 
     139            ! 
     140            z2 = 1. - z1 
     141            z4 = 1. - z3 
     142            z5 = 1. - z6 - z7 
     143            ! 
     144            imin = i1 ; imax = i2 
     145            jmin = j1 ; jmax = j2 
     146            !  
     147            ! Remove CORNERS 
     148            IF((nbondj == -1).OR.(nbondj == 2)) jmin = 2 + nbghostcells 
     149            IF((nbondj == +1).OR.(nbondj == 2)) jmax = nlcj - nbghostcells - 1 
     150            IF((nbondi == -1).OR.(nbondi == 2)) imin = 1 + nbghostcells 
     151            IF((nbondi == +1).OR.(nbondi == 2)) imax = nlci - nbghostcells - 1       
     152            ! 
     153            IF( eastern_side ) THEN 
     154               ibdy = nlci-nbghostcells 
     155               DO jn = 1, jptra 
     156                  tra(ibdy+1,jmin:jmax,k1:k2,jn) = z1 * ptab_child(ibdy+1,jmin:jmax,k1:k2,jn) + z2 * ptab_child(ibdy,jmin:jmax,k1:k2,jn) 
     157                  DO jk = 1, jpkm1 
     158                     DO jj = jmin,jmax 
     159                        IF( umask(ibdy-1,jj,jk) == 0._wp ) THEN 
     160                           tra(ibdy,jj,jk,jn) = tra(ibdy+1,jj,jk,jn) * tmask(ibdy,jj,jk) 
     161                        ELSE 
     162                           tra(ibdy,jj,jk,jn)=(z4*tra(ibdy+1,jj,jk,jn)+z3*tra(ibdy-1,jj,jk,jn))*tmask(ibdy,jj,jk) 
     163                           IF( un(ibdy-1,jj,jk) > 0._wp ) THEN 
     164                              tra(ibdy,jj,jk,jn)=( z6*tra(ibdy-1,jj,jk,jn)+z5*tra(ibdy+1,jj,jk,jn) &  
     165                                                 + z7*tra(ibdy-2,jj,jk,jn) ) * tmask(ibdy,jj,jk) 
     166                           ENDIF 
     167                        ENDIF 
     168                     END DO 
     169                  END DO 
     170                  ! Restore ghost points: 
     171                  tra(ibdy+1,jmin:jmax,k1:k2,jn) = ptab_child(ibdy+1,jmin:jmax,k1:k2,jn) * tmask(ibdy+1,jmin:jmax,k1:k2) 
     172               END DO 
     173            ENDIF 
     174            !  
     175            IF( northern_side ) THEN 
     176               jbdy = nlcj-nbghostcells          
     177               DO jn = 1, jptra 
     178                  tra(imin:imax,jbdy+1,k1:k2,jn) = z1 * ptab_child(imin:imax,jbdy+1,k1:k2,jn) + z2 * ptab_child(imin:imax,jbdy,k1:k2,jn) 
     179                  DO jk = 1, jpkm1 
     180                     DO ji = imin,imax 
     181                        IF( vmask(ji,jbdy-1,jk) == 0._wp ) THEN 
     182                           tra(ji,jbdy,jk,jn) = tra(ji,jbdy+1,jk,jn) * tmask(ji,jbdy,jk) 
     183                        ELSE 
     184                           tra(ji,jbdy,jk,jn)=(z4*tra(ji,jbdy+1,jk,jn)+z3*tra(ji,jbdy-1,jk,jn))*tmask(ji,jbdy,jk)         
     185                           IF (vn(ji,jbdy-1,jk) > 0._wp ) THEN 
     186                              tra(ji,jbdy,jk,jn)=( z6*tra(ji,jbdy-1,jk,jn)+z5*tra(ji,jbdy+1,jk,jn)  & 
     187                                                 + z7*tra(ji,jbdy-2,jk,jn) ) * tmask(ji,jbdy,jk) 
     188                           ENDIF 
     189                        ENDIF 
     190                     END DO 
     191                  END DO 
     192                  ! Restore ghost points: 
     193                  tra(imin:imax,jbdy+1,k1:k2,jn) = ptab_child(imin:imax,jbdy+1,k1:k2,jn) * tmask(imin:imax,jbdy+1,k1:k2) 
     194               END DO 
     195            ENDIF 
     196            ! 
     197            IF( western_side ) THEN     
     198               ibdy = 1+nbghostcells        
     199               DO jn = 1, jptra 
     200                  tra(ibdy-1,jmin:jmax,k1:k2,jn) = z1 * ptab_child(ibdy-1,jmin:jmax,k1:k2,jn) + z2 * ptab_child(ibdy,jmin:jmax,k1:k2,jn) 
     201                  DO jk = 1, jpkm1 
     202                     DO jj = jmin,jmax 
     203                        IF( umask(ibdy,jj,jk) == 0._wp ) THEN 
     204                           tra(ibdy,jj,jk,jn) = tra(ibdy-1,jj,jk,jn) * tmask(ibdy,jj,jk) 
     205                        ELSE 
     206                           tra(ibdy,jj,jk,jn)=(z4*tra(ibdy-1,jj,jk,jn)+z3*tra(ibdy+1,jj,jk,jn))*tmask(ibdy,jj,jk)         
     207                           IF( un(ibdy,jj,jk) < 0._wp ) THEN 
     208                              tra(ibdy,jj,jk,jn)=(z6*tra(ibdy+1,jj,jk,jn)+z5*tra(ibdy-1,jj,jk,jn)+z7*tra(ibdy+2,jj,jk,jn))*tmask(ibdy,jj,jk) 
     209                           ENDIF 
     210                        ENDIF 
     211                     END DO 
     212                  END DO 
     213                  ! Restore ghost points: 
     214                  tra(ibdy-1,jmin:jmax,k1:k2,jn) = ptab_child(ibdy-1,jmin:jmax,k1:k2,jn) * tmask(ibdy-1,jmin:jmax,k1:k2) 
     215               END DO 
     216            ENDIF 
     217            ! 
     218            IF( southern_side ) THEN   
     219               jbdy=1+nbghostcells         
     220               DO jn = 1, jptra 
     221                  tra(imin:imax,jbdy-1,k1:k2,jn) = z1 * ptab_child(imin:imax,jbdy-1,k1:k2,jn) + z2 * ptab_child(imin:imax,jbdy,k1:k2,jn) 
     222                  DO jk = 1, jpk       
     223                     DO ji=imin,imax 
     224                        IF( vmask(ji,jbdy,jk) == 0._wp ) THEN 
     225                           tra(ji,jbdy,jk,jn)=tra(ji,jbdy-1,jk,jn) * tmask(ji,jbdy,jk) 
     226                        ELSE 
     227                           tra(ji,jbdy,jk,jn)=(z4*tra(ji,jbdy-1,jk,jn)+z3*tra(ji,jbdy+1,jk,jn))*tmask(ji,jbdy,jk) 
     228                           IF( vn(ji,jbdy,jk) < 0._wp ) THEN 
     229                              tra(ji,jbdy,jk,jn)=(z6*tra(ji,jbdy+1,jk,jn)+z5*tra(ji,jbdy-1,jk,jn)+z7*tra(ji,jbdy+2,jk,jn))*tmask(ji,jbdy,jk) 
     230                           ENDIF 
     231                        ENDIF 
     232                     END DO 
     233                  END DO 
     234                  ! Restore ghost points: 
     235                  tra(imin:imax,jbdy-1,k1:k2,jn) = tra(imin:imax,jbdy-1,k1:k2,jn) * tmask(imin:imax,jbdy-1,k1:k2) 
     236               END DO 
     237            ENDIF 
     238            ! 
    181239         ENDIF 
    182          !  
    183          IF( northern_side ) THEN             
    184             DO jn = 1, jptra 
    185                tra(i1:i2,nlcj,1:jpk,jn) = zalpha1 * ptab_child(i1:i2,nlcj,1:jpk,jn) + zalpha2 * ptab_child(i1:i2,nlcj-1,1:jpk,jn) 
    186                DO jk = 1, jpkm1 
    187                   DO ji = imin,imax 
    188                      IF( vmask(ji,nlcj-2,jk) == 0.e0 ) THEN 
    189                         tra(ji,nlcj-1,jk,jn) = tra(ji,nlcj,jk,jn) * tmask(ji,nlcj-1,jk) 
    190                      ELSE 
    191                         tra(ji,nlcj-1,jk,jn)=(zalpha4*tra(ji,nlcj,jk,jn)+zalpha3*tra(ji,nlcj-2,jk,jn))*tmask(ji,nlcj-1,jk)         
    192                         IF (vn(ji,nlcj-2,jk) > 0.e0 ) THEN 
    193                            tra(ji,nlcj-1,jk,jn)=( zalpha6*tra(ji,nlcj-2,jk,jn)+zalpha5*tra(ji,nlcj,jk,jn)  & 
    194                                  + zalpha7*tra(ji,nlcj-3,jk,jn) ) * tmask(ji,nlcj-1,jk) 
    195                         ENDIF 
    196                      END DO 
    197                   END DO 
    198                END DO 
    199             ENDDO 
    200          ENDIF 
    201          ! 
    202          IF( western_side) THEN             
    203             DO jn = 1, jptra 
    204                tra(1,j1:j2,1:jpk,jn) = zalpha1 * ptab_child(1,j1:j2,1:jpk,jn) + zalpha2 * ptab_child(2,j1:j2,1:jpk,jn) 
    205                DO jk = 1, jpkm1 
    206                   DO jj = jmin,jmax 
    207                      IF( umask(2,jj,jk) == 0.e0 ) THEN 
    208                         tra(2,jj,jk,jn) = tra(1,jj,jk,jn) * tmask(2,jj,jk) 
    209                      ELSE 
    210                         tra(2,jj,jk,jn)=(zalpha4*tra(1,jj,jk,jn)+zalpha3*tra(3,jj,jk,jn))*tmask(2,jj,jk)         
    211                         IF( un(2,jj,jk) < 0.e0 ) THEN 
    212                            tra(2,jj,jk,jn)=(zalpha6*tra(3,jj,jk,jn)+zalpha5*tra(1,jj,jk,jn)+zalpha7*tra(4,jj,jk,jn))*tmask(2,jj,jk) 
    213                         ENDIF 
    214                      END DO 
    215                   END DO 
    216                END DO 
    217             END DO 
    218          ENDIF 
    219          ! 
    220          IF( southern_side ) THEN            
    221             DO jn = 1, jptra 
    222                tra(i1:i2,1,1:jpk,jn) = zalpha1 * ptab_child(i1:i2,1,1:jpk,jn) + zalpha2 * ptab_child(i1:i2,2,1:jpk,jn) 
    223                DO jk=1,jpkm1 
    224                   DO ji=imin,imax 
    225                      IF( vmask(ji,2,jk) == 0.e0 ) THEN 
    226                         tra(ji,2,jk,jn)=tra(ji,1,jk,jn) * tmask(ji,2,jk) 
    227                      ELSE 
    228                         tra(ji,2,jk,jn)=(zalpha4*tra(ji,1,jk,jn)+zalpha3*tra(ji,3,jk,jn))*tmask(ji,2,jk) 
    229                         IF( vn(ji,2,jk) < 0.e0 ) THEN 
    230                            tra(ji,2,jk,jn)=(zalpha6*tra(ji,3,jk,jn)+zalpha5*tra(ji,1,jk,jn)+zalpha7*tra(ji,4,jk,jn))*tmask(ji,2,jk) 
    231                         ENDIF 
    232                      END DO 
    233                   END DO 
    234                END DO 
    235             ENDIF 
    236             ! 
    237             ! Treatment of corners 
    238             IF( ll_east .AND.((nbondj == -1).OR.(nbondj == 2)) )   tra(nlci-1,   2  ,:,:) = ptab(nlci-1,   2  ,:,:)   ! East south 
    239             IF( ll_east .AND.((nbondj ==  1).OR.(nbondj == 2)) )   tra(nlci-1,nlcj-1,:,:) = ptab(nlci-1,nlcj-1,:,:)   ! East north 
    240             IF( ll_west .AND.((nbondj == -1).OR.(nbondj == 2)) )   tra(   2  ,   2  ,:,:) = ptab(   2  ,   2  ,:,:)   ! West south 
    241             IF( ll_west .AND.((nbondj ==  1).OR.(nbondj == 2)) )   tra(   2  ,nlcj-1,:,:) = ptab(   2  ,nlcj-1,:,:)   ! West north 
    242             ! 
    243          ENDIF 
    244          ! 
    245          ! Treatment of corners 
    246          !  
    247          ! East south 
    248          IF ((eastern_side).AND.((nbondj == -1).OR.(nbondj == 2))) THEN 
    249             tra(nlci-1,2,:,:) = ptab_child(nlci-1,2,:,:) 
    250          ENDIF 
    251          ! East north 
    252          IF ((eastern_side).AND.((nbondj == 1).OR.(nbondj == 2))) THEN 
    253             tra(nlci-1,nlcj-1,:,:) = ptab_child(nlci-1,nlcj-1,:,:) 
    254          ENDIF 
    255          ! West south 
    256          IF ((western_side).AND.((nbondj == -1).OR.(nbondj == 2))) THEN 
    257             tra(2,2,:,:) = ptab_child(2,2,:,:) 
    258          ENDIF 
    259          ! West north 
    260          IF ((western_side).AND.((nbondj == 1).OR.(nbondj == 2))) THEN 
    261             tra(2,nlcj-1,:,:) = ptab_child(2,nlcj-1,:,:) 
    262          ENDIF 
    263          ! 
    264240      ENDIF 
    265241      ! 
  • NEMO/trunk/src/NST/agrif_user.F90

    r9780 r9788  
    624624      ENDIF 
    625625 
    626       ENDIF 
    627626      ! Check passive tracer cell 
    628627      IF( nn_dttrc .NE. 1 ) THEN 
  • NEMO/trunk/src/TOP/PISCES/P2Z/p2zexp.F90

    r9598 r9788  
    253253 
    254254   !!====================================================================== 
    255 END MODULE  p2zexp 
     255END MODULE p2zexp 
  • NEMO/trunk/src/TOP/PISCES/P2Z/p2zopt.F90

    r9598 r9788  
    201201 
    202202   !!====================================================================== 
    203 END MODULE  p2zopt 
     203END MODULE p2zopt 
  • NEMO/trunk/src/TOP/PISCES/P2Z/p2zsed.F90

    r9598 r9788  
    152152 
    153153   !!====================================================================== 
    154 END MODULE  p2zsed 
     154END MODULE p2zsed 
  • NEMO/trunk/src/TOP/PISCES/P4Z/p4zche.F90

    r9598 r9788  
    831831 
    832832   !!====================================================================== 
    833 END MODULE  p4zche 
     833END MODULE p4zche 
  • NEMO/trunk/src/TOP/PISCES/P4Z/p5zice.F90

    r9598 r9788  
    600600   END FUNCTION p5z_lim_alloc 
    601601   !!====================================================================== 
    602 END MODULE  p5zice 
     602END MODULE p5zice 
  • NEMO/trunk/src/TOP/PISCES/P4Z/p5zmeso.F90

    r9598 r9788  
    427427 
    428428   !!====================================================================== 
    429 END MODULE  p5zmeso 
     429END MODULE p5zmeso 
  • NEMO/trunk/src/TOP/PISCES/P4Z/p5zmicro.F90

    r9598 r9788  
    369369 
    370370   !!====================================================================== 
    371 END MODULE  p5zmicro 
     371END MODULE p5zmicro 
  • NEMO/trunk/src/TOP/PISCES/P4Z/p5zmort.F90

    r9598 r9788  
    310310 
    311311   !!====================================================================== 
    312 END MODULE  p5zmort 
     312END MODULE p5zmort 
  • NEMO/trunk/src/TOP/PISCES/P4Z/p5zprod.F90

    r9598 r9788  
    618618   END FUNCTION p5z_prod_alloc 
    619619   !!====================================================================== 
    620 END MODULE  p5zprod 
     620END MODULE p5zprod 
  • NEMO/trunk/src/TOP/TRP/trcrad.F90

    r9598 r9788  
    1313   !!   trc_rad    : correction of negative concentrations 
    1414   !!---------------------------------------------------------------------- 
     15   USE par_trc             ! need jptra, number of passive tracers  
    1516   USE oce_trc             ! ocean dynamics and tracers variables 
    1617   USE trc                 ! ocean passive tracers variables 
  • NEMO/trunk/src/TOP/trcbc.F90

    r9598 r9788  
    3939   TYPE(FLD), SAVE, PUBLIC, ALLOCATABLE, DIMENSION(:)  :: sf_trccbc    ! structure of data input CBC (file informations, fields read) 
    4040   REAL(wp) , SAVE, PUBLIC, ALLOCATABLE, DIMENSION(:)  :: rf_trofac    ! multiplicative factor for OBCtracer values 
    41    TYPE(FLD), SAVE, PUBLIC, ALLOCATABLE, DIMENSION(:), TARGET  :: sf_trcobc    ! structure of data input OBC (file informations, fields read) 
     41   TYPE(FLD), SAVE, PUBLIC, ALLOCATABLE, DIMENSION(:)  :: sf_trcobc    ! structure of data input OBC (file informations, fields read) 
    4242   TYPE(MAP_POINTER), ALLOCATABLE, DIMENSION(:) :: nbmap_ptr   ! array of pointers to nbmap 
    4343 
Note: See TracChangeset for help on using the changeset viewer.