Ignore:
Timestamp:
2014-01-17T17:56:32+01:00 (7 years ago)
Author:
jchanut
Message:

Restore AGRIF and BDY compatibility, see ticket #1133

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2013/dev_MERGE_2013/NEMOGCM/NEMO/OPA_SRC/BDY/bdydta.F90

    r4333 r4354  
    2020   !!    bdy_dta_init   : initialise arrays etc for reading of external data 
    2121   !!---------------------------------------------------------------------- 
    22    USE wrk_nemo        ! Memory Allocation 
    2322   USE timing          ! Timing 
    2423   USE oce             ! ocean dynamics and tracers 
     
    5150   LOGICAL,           DIMENSION(jp_bdy) ::   ln_full_vel_array ! =T => full velocities in 3D boundary conditions 
    5251                                                               ! =F => baroclinic velocities in 3D boundary conditions 
    53  
     52!$AGRIF_DO_NOT_TREAT 
    5453   TYPE(FLD), PUBLIC, ALLOCATABLE, DIMENSION(:), TARGET ::   bf        ! structure of input fields (file informations, fields read) 
    55  
     54!$AGRIF_END_DO_NOT_TREAT 
    5655   TYPE(MAP_POINTER), ALLOCATABLE, DIMENSION(:) :: nbmap_ptr   ! array of pointers to nbmap 
    5756 
     
    103102         ! Calculate depth-mean currents 
    104103         !----------------------------- 
    105          CALL wrk_alloc(jpi,jpj,pun2d,pvn2d)  
    106  
    107          pun2d(:,:) = 0.e0 
    108          pvn2d(:,:) = 0.e0 
    109          DO ik = 1, jpkm1   !! Vertically integrated momentum trends 
    110              pun2d(:,:) = pun2d(:,:) + fse3u(:,:,ik) * umask(:,:,ik) * un(:,:,ik) 
    111              pvn2d(:,:) = pvn2d(:,:) + fse3v(:,:,ik) * vmask(:,:,ik) * vn(:,:,ik) 
    112          END DO 
    113          pun2d(:,:) = pun2d(:,:) * hur(:,:) 
    114          pvn2d(:,:) = pvn2d(:,:) * hvr(:,:) 
    115104          
    116105         DO ib_bdy = 1, nb_bdy 
     
    135124                     ii = idx_bdy(ib_bdy)%nbi(ib,igrd) 
    136125                     ij = idx_bdy(ib_bdy)%nbj(ib,igrd) 
    137                      dta_bdy(ib_bdy)%u2d(ib) = pun2d(ii,ij) * umask(ii,ij,1)          
     126                     dta_bdy(ib_bdy)%u2d(ib) = un_b(ii,ij) * umask(ii,ij,1)          
    138127                  END DO  
    139128               END IF 
     
    143132                     ii = idx_bdy(ib_bdy)%nbi(ib,igrd) 
    144133                     ij = idx_bdy(ib_bdy)%nbj(ib,igrd) 
    145                      dta_bdy(ib_bdy)%v2d(ib) = pvn2d(ii,ij) * vmask(ii,ij,1)          
     134                     dta_bdy(ib_bdy)%v2d(ib) = vn_b(ii,ij) * vmask(ii,ij,1)          
    146135                  END DO  
    147136               END IF 
     
    156145                        ii = idx_bdy(ib_bdy)%nbi(ib,igrd) 
    157146                        ij = idx_bdy(ib_bdy)%nbj(ib,igrd) 
    158                         dta_bdy(ib_bdy)%u3d(ib,ik) =  ( un(ii,ij,ik) - pun2d(ii,ij) ) * umask(ii,ij,ik)          
     147                        dta_bdy(ib_bdy)%u3d(ib,ik) =  ( un(ii,ij,ik) - un_b(ii,ij) ) * umask(ii,ij,ik)          
    159148                     END DO 
    160149                  END DO  
     
    166155                        ii = idx_bdy(ib_bdy)%nbi(ib,igrd) 
    167156                        ij = idx_bdy(ib_bdy)%nbj(ib,igrd) 
    168                         dta_bdy(ib_bdy)%v3d(ib,ik) =  ( vn(ii,ij,ik) - pvn2d(ii,ij) ) * vmask(ii,ij,ik)          
     157                        dta_bdy(ib_bdy)%v3d(ib,ik) =  ( vn(ii,ij,ik) - vn_b(ii,ij) ) * vmask(ii,ij,ik)          
    169158                        END DO 
    170159                  END DO  
     
    262251         ENDDO ! ib_bdy 
    263252 
    264          CALL wrk_dealloc(jpi,jpj,pun2d,pvn2d)  
    265253 
    266254      ENDIF ! kt .eq. nit000 
     
    919907   !!============================================================================== 
    920908END MODULE bdydta 
     909 
     910 
Note: See TracChangeset for help on using the changeset viewer.