- Timestamp:
- 2016-11-08T17:05:22+01:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2015/dev_r5003_MERCATOR6_CRS/NEMOGCM/NEMO/OPA_SRC/CRS/crsini.F90
r7207 r7210 223 223 CALL crs_dom_e3( e1v, e2v, zfse3v, p_sfc_2d_crs=e1v_crs , cd_type='V', p_mask=vmask, p_e3_crs=e3v_0_crs, p_e3_max_crs=e3v_max_0_crs) 224 224 225 WHERE(e3t_max_0_crs == 0._wp) e3t_max_0_crs=r_inf 226 WHERE(e3u_max_0_crs == 0._wp) e3u_max_0_crs=r_inf 227 WHERE(e3v_max_0_crs == 0._wp) e3v_max_0_crs=r_inf 228 WHERE(e3w_max_0_crs == 0._wp) e3w_max_0_crs=r_inf 225 DO jk = 1, jpk 226 DO ji = 1, jpi_crs 227 DO jj = 1, jpj_crs 228 IF( e3t_max_0_crs(ji,jj,jk) == 0._wp ) e3t_max_0_crs(ji,jj,jk) = e3t_1d(jk) 229 IF( e3w_max_0_crs(ji,jj,jk) == 0._wp ) e3w_max_0_crs(ji,jj,jk) = e3w_1d(jk) 230 IF( e3u_max_0_crs(ji,jj,jk) == 0._wp ) e3u_max_0_crs(ji,jj,jk) = e3t_1d(jk) 231 IF( e3v_max_0_crs(ji,jj,jk) == 0._wp ) e3v_max_0_crs(ji,jj,jk) = e3t_1d(jk) 232 ENDDO 233 ENDDO 234 ENDDO 229 235 230 236 #if defined key_vvl … … 281 287 CALL crs_dom_ope( gdept_0, 'MAX', 'T', tmask, gdept_0_crs, p_e3=zfse3t, psgn=1.0 ) 282 288 CALL crs_dom_ope( gdepw_0, 'MAX', 'W', tmask, gdepw_0_crs, p_e3=zfse3w, psgn=1.0 ) 289 290 DO jk = 1, jpk 291 DO ji = 1, jpi_crs 292 DO jj = 1, jpj_crs 293 IF( gdept_0_crs(ji,jj,jk) .LE. 0._wp ) gdept_0_crs(ji,jj,jk) = gdept_1d(jk) 294 IF( gdepw_0_crs(ji,jj,jk) .LE. 0._wp ) gdepw_0_crs(ji,jj,jk) = gdepw_1d(jk) 295 ENDDO 296 ENDDO 297 ENDDO 298 283 299 #if defined key_vvl 284 300 gdept_n_crs(:,:,:) = gdept_0_crs(:,:,:)
Note: See TracChangeset
for help on using the changeset viewer.