- Timestamp:
- 2016-12-01T11:30:29+01:00 (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2016/dev_merge_2016/NEMOGCM/NEMO/OPA_SRC/DOM/dommsk.F90
r6140 r7412 24 24 USE oce ! ocean dynamics and tracers 25 25 USE dom_oce ! ocean space and time domain 26 !26 USE bdy_oce 27 27 USE in_out_manager ! I/O manager 28 28 USE lbclnk ! ocean lateral boundary conditions (or mpp link) 29 29 USE lib_mpp ! 30 USE iom 30 31 USE wrk_nemo ! Memory allocation 31 32 USE timing ! Timing … … 88 89 !! are defined with the proper value at lateral domain boundaries. 89 90 !! 90 !! In case of open boundaries (l k_bdy=T):91 !! In case of open boundaries (ln_bdy=T): 91 92 !! - tmask is set to 1 on the points to be computed bay the open 92 93 !! boundaries routines. … … 102 103 INTEGER :: iif, iil, ii0, ii1, ii ! local integers 103 104 INTEGER :: ijf, ijl, ij0, ij1 ! - - 104 INTEGER :: ios 105 INTEGER :: ios, inum 105 106 INTEGER :: isrow ! index for ORCA1 starting row 106 107 INTEGER , POINTER, DIMENSION(:,:) :: imsk … … 108 109 !! 109 110 NAMELIST/namlbc/ rn_shlat, ln_vorlat 111 NAMELIST/nambdy/ ln_bdy ,nb_bdy, ln_coords_file, cn_coords_file, & 112 & ln_mask_file, cn_mask_file, cn_dyn2d, nn_dyn2d_dta, & 113 & cn_dyn3d, nn_dyn3d_dta, cn_tra, nn_tra_dta, & 114 & ln_tra_dmp, ln_dyn3d_dmp, rn_time_dmp, rn_time_dmp_out, & 115 & cn_ice_lim, nn_ice_lim_dta, & 116 & rn_ice_tem, rn_ice_sal, rn_ice_age, & 117 & ln_vol, nn_volctl, nn_rimwidth, nb_jpk_bdy 110 118 !!--------------------------------------------------------------------- 111 119 ! … … 155 163 END DO 156 164 165 REWIND( numnam_ref ) ! Namelist nambdy in reference namelist :Unstructured open boundaries 166 READ ( numnam_ref, nambdy, IOSTAT = ios, ERR = 903) 167 903 IF( ios /= 0 ) CALL ctl_nam ( ios , 'nambdy in reference namelist', lwp ) 168 169 REWIND( numnam_cfg ) ! Namelist nambdy in configuration namelist :Unstructured open boundaries 170 READ ( numnam_cfg, nambdy, IOSTAT = ios, ERR = 904 ) 171 904 IF( ios /= 0 ) CALL ctl_nam ( ios , 'nambdy in configuration namelist', lwp ) 172 IF(lwm) WRITE ( numond, nambdy ) 173 174 IF( ln_bdy .AND. ln_mask_file ) THEN ! correct for bdy mask 175 CALL iom_open( cn_mask_file, inum ) 176 CALL iom_get ( inum, jpdom_data, 'bdy_msk', bdytmask(:,:) ) 177 CALL iom_close( inum ) 178 179 ! Mask corrections 180 ! ---------------- 181 DO jk = 1, jpkm1 182 DO jj = 1, jpj 183 DO ji = 1, jpi 184 tmask(ji,jj,jk) = tmask(ji,jj,jk) * bdytmask(ji,jj) 185 END DO 186 END DO 187 END DO 188 ENDIF 189 157 190 ! (ISF) define barotropic mask and mask the ice shelf point 158 191 ssmask(:,:)=tmask(:,:,1) ! at this stage ice shelf is not masked
Note: See TracChangeset
for help on using the changeset viewer.