Changeset 15053 for NEMO/trunk/src/OCE
- Timestamp:
- 2021-06-24T17:39:38+02:00 (3 years ago)
- Location:
- NEMO/trunk/src/OCE
- Files:
-
- 7 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/trunk/src/OCE/ISF/isfcavmlt.F90
r15004 r15053 213 213 ! compute upward heat flux zhtflx and upward water flux zwflx 214 214 ! Resolution of a 3d equation from equation 24, 25 and 26 (note conduction through the ice has been added to Eq 24) 215 DO_2D( 1, 1, 1, 1)215 DO_2D( nn_hls, nn_hls, nn_hls, nn_hls ) 216 216 ! 217 217 ! compute coeficient to solve the 2nd order equation -
NEMO/trunk/src/OCE/ISF/isfcpl.F90
r15004 r15053 374 374 ! ----------------------------------------------------------------------------------------- 375 375 ! case we open a cell but no neigbour cells available to get an estimate of T and S 376 DO_3D( 1, 1, 1, 1, 1,jpk-1 )376 DO_3D( 0, 0, 0, 0, 1,jpk-1 ) 377 377 IF (tmask(ji,jj,jk) == 1._wp .AND. ts(ji,jj,jk,2,Kmm) == 0._wp) & 378 378 & CALL ctl_stop('STOP', 'failing to fill all new weet cell, & -
NEMO/trunk/src/OCE/ISF/isfdiags.F90
r15004 r15053 101 101 zvar3d(:,:,:) = 0._wp 102 102 ! 103 DO_2D( 1, 1, 1, 1)103 DO_2D( nn_hls, nn_hls, nn_hls, nn_hls ) 104 104 ikt = ktop(ji,jj) 105 105 ikb = kbot(ji,jj) -
NEMO/trunk/src/OCE/ISF/isfload.F90
r14064 r15053 94 94 ! !- Surface value + ice shelf gradient 95 95 pload(:,:) = 0._wp ! compute pressure due to ice shelf load 96 DO_2D( 1, 1, 1, 1)96 DO_2D( nn_hls, nn_hls, nn_hls, nn_hls ) 97 97 ikt = mikt(ji,jj) 98 98 ! -
NEMO/trunk/src/OCE/ISF/isftbl.F90
r14834 r15053 78 78 ! compute tbl property at T point 79 79 pvarout(1,:) = 0._wp 80 DO_2D( 0, 1, 1, 1)80 DO_2D( nn_hls-1, nn_hls, nn_hls, nn_hls ) 81 81 pvarout(ji,jj) = 0.5_wp * (zvarout(ji,jj) + zvarout(ji-1,jj)) 82 82 END_2D … … 99 99 ! pvarout is an averaging of wet point 100 100 pvarout(:,1) = 0._wp 101 DO_2D( 1, 1, 0, 1)101 DO_2D( nn_hls, nn_hls, nn_hls-1, nn_hls ) 102 102 pvarout(ji,jj) = 0.5_wp * (zvarout(ji,jj) + zvarout(ji,jj-1)) 103 103 END_2D … … 138 138 ! 139 139 ! compute tbl top.bottom level and thickness 140 DO_2D( 1, 1, 1, 1)140 DO_2D( nn_hls, nn_hls, nn_hls, nn_hls ) 141 141 ! 142 142 ! tbl top/bottom indices initialisation -
NEMO/trunk/src/OCE/TRA/eosbn2.F90
r14834 r15053 1078 1078 z1_T0 = 1._wp/40._wp 1079 1079 ! 1080 DO_2D( 1, 1, 1, 1)1080 DO_2D( nn_hls, nn_hls, nn_hls, nn_hls ) 1081 1081 ! 1082 1082 zt = ctmp (ji,jj) * z1_T0 … … 1142 1142 ! 1143 1143 z1_S0 = 1._wp / 35.16504_wp 1144 DO_2D( 1, 1, 1, 1)1144 DO_2D( nn_hls, nn_hls, nn_hls, nn_hls ) 1145 1145 zs= SQRT( ABS( psal(ji,jj) ) * z1_S0 ) ! square root salinity 1146 1146 ptf(ji,jj) = ((((1.46873e-03_wp*zs-9.64972e-03_wp)*zs+2.28348e-02_wp)*zs & … … 1249 1249 CASE( np_teos10, np_eos80 ) !== polynomial TEOS-10 / EOS-80 ==! 1250 1250 ! 1251 DO_3D( 1, 1, 1, 1, 1, jpkm1 )1251 DO_3D( nn_hls, nn_hls, nn_hls, nn_hls, 1, jpkm1 ) 1252 1252 ! 1253 1253 zh = gdept(ji,jj,jk,Kmm) * r1_Z0 ! depth … … 1308 1308 CASE( np_seos ) !== Vallis (2006) simplified EOS ==! 1309 1309 ! 1310 DO_3D( 1, 1, 1, 1, 1, jpkm1 )1310 DO_3D( nn_hls, nn_hls, nn_hls, nn_hls, 1, jpkm1 ) 1311 1311 zt = pts(ji,jj,jk,jp_tem) - 10._wp ! temperature anomaly (t-T0) 1312 1312 zs = pts (ji,jj,jk,jp_sal) - 35._wp ! abs. salinity anomaly (s-S0) -
NEMO/trunk/src/OCE/TRA/trabbl.F90
r14834 r15053 339 339 ENDIF 340 340 ! !* bottom variables (T, S, alpha, beta, depth, velocity) 341 DO_2D( 1, 1, 1, 1)341 DO_2D( nn_hls, nn_hls, nn_hls, nn_hls ) 342 342 ik = mbkt(ji,jj) ! bottom T-level index 343 343 zts (ji,jj,jp_tem) = ts(ji,jj,ik,jp_tem,Kbb) ! bottom before T and S
Note: See TracChangeset
for help on using the changeset viewer.