- Timestamp:
- 2019-07-26T11:32:13+02:00 (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/UKMO/NEMO_4.0_GO8_coupled_iodef/src/ICE/icethd_do.F90
r11082 r11355 207 207 END DO 208 208 209 IF(narea == 419) THEN 210 WRITE(9419,*) 'npti in icetd_do = ',npti 211 WRITE(9419,*) 'max qlead = ',MAXVAL( qlead(:,:) ) 212 WRITE(9419,*) 'min qlead = ',MINVAL( qlead(:,:) ) 213 WRITE(9419,*) 'qlead(48,31) = ', qlead(48,31) 214 WRITE(9419,*) 'tau_icebfr(48,31) = ',tau_icebfr(48,31) 215 ENDIF 216 209 217 ! Move from 2-D to 1-D vectors 210 218 IF ( npti > 0 ) THEN … … 230 238 CALL tab_2d_1d( npti, nptidx(1:npti), rn_amax_1d(1:npti) , rn_amax_2d ) 231 239 CALL tab_2d_1d( npti, nptidx(1:npti), sss_1d (1:npti) , sss_m ) 240 241 IF(narea == 419) THEN 242 WRITE(9419,*) 'wfx_opw(48,31) = ',wfx_opw(48,31) 243 WRITE(9419,*) 'max wfx_opw_1d icethd_do after to 1d = ',MAXVAL( ABS( wfx_opw_1d(:) ) ) 244 ENDIF 232 245 233 246 ! Convert units for ice internal energy … … 295 308 sfx_opw_1d(ji) = sfx_opw_1d(ji) - zv_newice(ji) * rhoi * zs_newice(ji) * r1_rdtice 296 309 END DO 310 311 IF(narea == 419) THEN 312 WRITE(9419,*) 'max wfx_opw_1d icethd_do after calculation = ',MAXVAL( ABS( wfx_opw_1d(:) ) ) 313 WRITE(9419,*) 'max zv_newice icethd_do after calculation = ',MAXVAL( ABS( zv_newice(:) ) ) 314 WRITE(9419,*) 'rhoi = ',rhoi 315 WRITE(9419,*) 'r1_rdtice = ',r1_rdtice 316 WRITE(9419,*) 'max qlead_1d icethd_do after calculation = ',MAXVAL( ABS( qlead_1d(:) ) ) 317 WRITE(9419,*) 'max ze_newice icethd_do after calculation = ',MAXVAL( ABS( ze_newice(:) ) ) 318 WRITE(9419,*) 'max t_bo_1d icethd_do after calculation = ',MAXVAL( ABS( t_bo_1d(:) ) ) 319 ENDIF 297 320 298 321 zv_frazb(1:npti) = 0._wp
Note: See TracChangeset
for help on using the changeset viewer.