- Timestamp:
- 2020-06-07T18:26:09+02:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
utils/tools_dev_r12970_AGRIF_CMEMS/DOMAINcfg/src/agrif_domzgr.F90
r13024 r13056 1 MODULE agrif_domzgr 2 3 USE agrif_profiles 4 USE dom_oce 5 6 IMPLICIT NONE 7 PRIVATE 8 9 PUBLIC :: agrif_create_bathy_meter 10 11 12 CONTAINS 13 1 14 #if defined key_agrif 2 subroutine agrif_domzgr3 end subroutine agrif_domzgr4 15 5 subroutine agrif_create_bathy_meter 6 use agrif_profiles 7 external :: init_bathy 16 SUBROUTINE agrif_create_bathy_meter 8 17 9 callAgrif_Init_variable(bathy_id, procname = init_bathy)18 CALL Agrif_Init_variable(bathy_id, procname = init_bathy) 10 19 11 end subroutineagrif_create_bathy_meter20 END SUBROUTINE agrif_create_bathy_meter 12 21 13 SUBROUTINE init_bathy( ptab, i1, i2, j1, j2, before, nb,ndir) 14 use dom_oce 22 SUBROUTINE init_bathy( ptab, i1, i2, j1, j2, before, nb,ndir) 15 23 !!---------------------------------------------------------------------- 16 24 !! *** ROUTINE interpsshn *** … … 20 28 LOGICAL , INTENT(in ) :: before 21 29 INTEGER , INTENT(in ) :: nb , ndir 22 LOGICAL :: western_side, eastern_side,northern_side,southern_side23 30 ! 24 31 !!---------------------------------------------------------------------- 25 32 INTEGER :: ji,jj 26 ! 27 western_side = (nb == 1).AND.(ndir == 1) 28 eastern_side = (nb == 1).AND.(ndir == 2) 29 southern_side = (nb == 2).AND.(ndir == 1) 30 northern_side = (nb == 2).AND.(ndir == 2) 33 31 34 IF( before) THEN 32 35 ptab(i1:i2,j1:j2) = bathy(i1:i2,j1:j2) 33 34 do jj=j1,j2 35 do ji=i1,i2 36 ptab(ji,jj) = SUM ( e3t_0(ji,jj, 1:mbkt(ji,jj) ) ) * ssmask(ji,jj) 37 enddo 38 enddo 36 DO jj=j1,j2 37 DO ji=i1,i2 38 ptab(ji,jj) = SUM( e3t_0(ji,jj, 1:mbkt(ji,jj) ) ) * ssmask(ji,jj) 39 END DO 40 END DO 39 41 ELSE 40 42 bathy(i1:i2,j1:j2)=ptab … … 42 44 ! 43 45 END SUBROUTINE init_bathy 46 44 47 #else 45 subroutine agrif_domzgr_empty 46 end subroutine agrif_domzgr_empty 48 SUBROUTINE agrif_create_bathy_meter 49 END SUBROUTINE agrif_create_bathy_meter 47 50 #endif 51 END MODULE agrif_domzgr
Note: See TracChangeset
for help on using the changeset viewer.