- Timestamp:
- 2019-07-29T18:10:24+02:00 (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/UKMO/NEMO_4.0_GO8_coupled_iodef/src/OCE/step.F90
r11355 r11367 39 39 ! 40 40 USE iom ! xIOs server 41 USE ice 41 42 42 43 IMPLICIT NONE … … 94 95 IF( ln_timing ) CALL timing_start('stp') 95 96 96 IF(narea == 419) THEN 97 WRITE(9419,*) 'max sshn step before update IO = ',MAXVAL( ABS( sshn(:,:) ) ) 98 WRITE(9419,*) 'max ssha = ',MAXVAL( ABS( ssha(:,:) ) ) 99 WRITE(9419,*) 'max emp_b = ',MAXVAL( ABS( emp_b(:,:) ) ) 100 WRITE(9419,*) 'max emp = ',MAXVAL( ABS( emp(:,:) ) ) 97 IF(narea == 68) THEN 98 WRITE(968,*) 'in step before update IO: u_ice = ',SUM(u_ice) 99 WRITE(968,*) 'v_ice(54,0:2) = ',SUM(v_ice) 100 WRITE(968,*) 'h_i(54,0:2) = ',SUM(h_i) 101 WRITE(968,*) 'h_s(54,0:2) = ',SUM(h_s) 102 WRITE(968,*) 'v_i(54,0:2) = ',SUM(v_i) 103 WRITE(968,*) 'sv_i(54,0:2) = ',SUM(sv_i) 104 WRITE(968,*) 'e_s(54,0:2) = ',SUM(e_s) 105 WRITE(968,*) 'e_i(54,0:2) = ',SUM(e_i) 101 106 ENDIF 102 107 … … 115 120 IF( ln_crs ) CALL iom_setkt( kstp - nit000 + 1, TRIM(cxios_context)//"_crs" ) ! tell IOM we are at time step kstp 116 121 117 IF(narea == 419) THEN 118 WRITE(9419,*) 'max sshn step before external forcing = ',MAXVAL( ABS( sshn(:,:) ) ) 119 WRITE(9419,*) 'max ssha = ',MAXVAL( ABS( ssha(:,:) ) ) 120 WRITE(9419,*) 'max emp_b = ',MAXVAL( ABS( emp_b(:,:) ) ) 121 WRITE(9419,*) 'max emp = ',MAXVAL( ABS( emp(:,:) ) ) 122 IF(narea == 68) THEN 123 WRITE(968,*) 'in step before sbc: u_ice = ',SUM(u_ice) 124 WRITE(968,*) 'v_ice(54,0:2) = ',SUM(v_ice) 125 WRITE(968,*) 'h_i(54,0:2) = ',SUM(h_i) 126 WRITE(968,*) 'h_s(54,0:2) = ',SUM(h_s) 127 WRITE(968,*) 'v_i(54,0:2) = ',SUM(v_i) 128 WRITE(968,*) 'sv_i(54,0:2) = ',SUM(sv_i) 129 WRITE(968,*) 'e_s(54,0:2) = ',SUM(e_s) 130 WRITE(968,*) 'e_i(54,0:2) = ',SUM(e_i) 122 131 ENDIF 123 132 … … 139 148 IF( ln_sto_eos ) CALL sto_pts( tsn ) ! Random T/S fluctuations 140 149 141 IF(narea == 419) THEN 142 WRITE(9419,*) 'max sshn step before ocean physics = ',MAXVAL( ABS( sshn(:,:) ) ) 143 WRITE(9419,*) 'max ssha = ',MAXVAL( ABS( ssha(:,:) ) ) 144 WRITE(9419,*) 'max emp_b = ',MAXVAL( ABS( emp_b(:,:) ) ) 145 WRITE(9419,*) 'max emp = ',MAXVAL( ABS( emp(:,:) ) ) 150 IF(narea == 68) THEN 151 WRITE(968,*) 'in step before ocean physics: u_ice = ',SUM(u_ice) 152 WRITE(968,*) 'v_ice(54,0:2) = ',SUM(v_ice) 153 WRITE(968,*) 'h_i(54,0:2) = ',SUM(h_i) 154 WRITE(968,*) 'h_s(54,0:2) = ',SUM(h_s) 155 WRITE(968,*) 'v_i(54,0:2) = ',SUM(v_i) 156 WRITE(968,*) 'sv_i(54,0:2) = ',SUM(sv_i) 157 WRITE(968,*) 'e_s(54,0:2) = ',SUM(e_s) 158 WRITE(968,*) 'e_i(54,0:2) = ',SUM(e_i) 146 159 ENDIF 147 160 … … 180 193 IF( l_ldfdyn_time ) CALL ldf_dyn( kstp ) ! eddy viscosity coeff. 181 194 182 IF(narea == 419) THEN 183 WRITE(9419,*) 'max sshn step before ocean dynamics = ',MAXVAL( ABS( sshn(:,:) ) ) 184 WRITE(9419,*) 'max ssha = ',MAXVAL( ABS( ssha(:,:) ) ) 185 WRITE(9419,*) 'max emp_b = ',MAXVAL( ABS( emp_b(:,:) ) ) 186 WRITE(9419,*) 'max emp = ',MAXVAL( ABS( emp(:,:) ) ) 195 IF(narea == 68) THEN 196 WRITE(968,*) 'in step before dynamics: u_ice = ',SUM(u_ice) 197 WRITE(968,*) 'v_ice(54,0:2) = ',SUM(v_ice) 198 WRITE(968,*) 'h_i(54,0:2) = ',SUM(h_i) 199 WRITE(968,*) 'h_s(54,0:2) = ',SUM(h_s) 200 WRITE(968,*) 'v_i(54,0:2) = ',SUM(v_i) 201 WRITE(968,*) 'sv_i(54,0:2) = ',SUM(sv_i) 202 WRITE(968,*) 'e_s(54,0:2) = ',SUM(e_s) 203 WRITE(968,*) 'e_i(54,0:2) = ',SUM(e_i) 187 204 ENDIF 188 205 … … 241 258 & CALL dyn_asm_inc ( kstp ) ! apply dynamics assimilation increment 242 259 243 IF(narea == 419) THEN 244 WRITE(9419,*) 'max ssha step before bdy_dyn3d_dmp = ',MAXVAL( ABS( ssha(:,:) ) ) 260 IF(narea == 68) THEN 261 WRITE(968,*) 'in step before bdy_dyn3d: u_ice = ',SUM(u_ice) 262 WRITE(968,*) 'v_ice(54,0:2) = ',SUM(v_ice) 263 WRITE(968,*) 'h_i(54,0:2) = ',SUM(h_i) 264 WRITE(968,*) 'h_s(54,0:2) = ',SUM(h_s) 265 WRITE(968,*) 'v_i(54,0:2) = ',SUM(v_i) 266 WRITE(968,*) 'sv_i(54,0:2) = ',SUM(sv_i) 267 WRITE(968,*) 'e_s(54,0:2) = ',SUM(e_s) 268 WRITE(968,*) 'e_i(54,0:2) = ',SUM(e_i) 245 269 ENDIF 246 270 … … 287 311 CALL dyn_spg ( kstp ) ! surface pressure gradient 288 312 289 290 IF(narea == 419) THEN 291 WRITE(9419,*) 'max ssha step before div_hor = ',MAXVAL( ABS( ssha(:,:) ) ) 313 IF(narea == 68) THEN 314 WRITE(968,*) 'in step before div_hor: u_ice = ',SUM(u_ice) 315 WRITE(968,*) 'v_ice(54,0:2) = ',SUM(v_ice) 316 WRITE(968,*) 'h_i(54,0:2) = ',SUM(h_i) 317 WRITE(968,*) 'h_s(54,0:2) = ',SUM(h_s) 318 WRITE(968,*) 'v_i(54,0:2) = ',SUM(v_i) 319 WRITE(968,*) 'sv_i(54,0:2) = ',SUM(sv_i) 320 WRITE(968,*) 'e_s(54,0:2) = ',SUM(e_s) 321 WRITE(968,*) 'e_i(54,0:2) = ',SUM(e_i) 292 322 ENDIF 293 323 … … 307 337 CALL dyn_zdf ( kstp ) ! vertical diffusion 308 338 309 IF(narea == 419) THEN 310 WRITE(9419,*) 'max sshn step before cool skin = ',MAXVAL( ABS( sshn(:,:) ) ) 311 WRITE(9419,*) 'max ssha = ',MAXVAL( ABS( ssha(:,:) ) ) 339 IF(narea == 68) THEN 340 WRITE(968,*) 'in step before cool skin: u_ice = ',SUM(u_ice) 341 WRITE(968,*) 'v_ice(54,0:2) = ',SUM(v_ice) 342 WRITE(968,*) 'h_i(54,0:2) = ',SUM(h_i) 343 WRITE(968,*) 'h_s(54,0:2) = ',SUM(h_s) 344 WRITE(968,*) 'v_i(54,0:2) = ',SUM(v_i) 345 WRITE(968,*) 'sv_i(54,0:2) = ',SUM(sv_i) 346 WRITE(968,*) 'e_s(54,0:2) = ',SUM(e_s) 347 WRITE(968,*) 'e_i(54,0:2) = ',SUM(e_i) 312 348 ENDIF 313 349 … … 317 353 IF ( ln_diurnal ) CALL stp_diurnal( kstp ) 318 354 319 IF(narea == 419) THEN 320 WRITE(9419,*) 'max sshn step before diagnostics = ',MAXVAL( ABS( sshn(:,:) ) ) 321 WRITE(9419,*) 'max ssha = ',MAXVAL( ABS( ssha(:,:) ) ) 355 IF(narea == 68) THEN 356 WRITE(968,*) 'in step before diagnostics: u_ice = ',SUM(u_ice) 357 WRITE(968,*) 'v_ice(54,0:2) = ',SUM(v_ice) 358 WRITE(968,*) 'h_i(54,0:2) = ',SUM(h_i) 359 WRITE(968,*) 'h_s(54,0:2) = ',SUM(h_s) 360 WRITE(968,*) 'v_i(54,0:2) = ',SUM(v_i) 361 WRITE(968,*) 'sv_i(54,0:2) = ',SUM(sv_i) 362 WRITE(968,*) 'e_s(54,0:2) = ',SUM(e_s) 363 WRITE(968,*) 'e_i(54,0:2) = ',SUM(e_i) 322 364 ENDIF 323 365 … … 343 385 #endif 344 386 345 IF(narea == 419) THEN 346 WRITE(9419,*) 'max sshn step before active tracers = ',MAXVAL( ABS( sshn(:,:) ) ) 347 WRITE(9419,*) 'max ssha = ',MAXVAL( ABS( ssha(:,:) ) ) 387 IF(narea == 68) THEN 388 WRITE(968,*) 'in step before active tracers: u_ice = ',SUM(u_ice) 389 WRITE(968,*) 'v_ice(54,0:2) = ',SUM(v_ice) 390 WRITE(968,*) 'h_i(54,0:2) = ',SUM(h_i) 391 WRITE(968,*) 'h_s(54,0:2) = ',SUM(h_s) 392 WRITE(968,*) 'v_i(54,0:2) = ',SUM(v_i) 393 WRITE(968,*) 'sv_i(54,0:2) = ',SUM(sv_i) 394 WRITE(968,*) 'e_s(54,0:2) = ',SUM(e_s) 395 WRITE(968,*) 'e_i(54,0:2) = ',SUM(e_i) 348 396 ENDIF 349 397 … … 356 404 & ln_trainc ) CALL tra_asm_inc ( kstp ) ! apply tracer assimilation increment 357 405 358 IF(narea == 419) THEN 359 WRITE(9419,*) 'max sshn step before tra_sbc = ',MAXVAL( ABS( sshn(:,:) ) ) 360 WRITE(9419,*) 'max ssha = ',MAXVAL( ABS( ssha(:,:) ) ) 406 IF(narea == 68) THEN 407 WRITE(968,*) 'in step before tra_sbc: u_ice = ',SUM(u_ice) 408 WRITE(968,*) 'v_ice(54,0:2) = ',SUM(v_ice) 409 WRITE(968,*) 'h_i(54,0:2) = ',SUM(h_i) 410 WRITE(968,*) 'h_s(54,0:2) = ',SUM(h_s) 411 WRITE(968,*) 'v_i(54,0:2) = ',SUM(v_i) 412 WRITE(968,*) 'sv_i(54,0:2) = ',SUM(sv_i) 413 WRITE(968,*) 'e_s(54,0:2) = ',SUM(e_s) 414 WRITE(968,*) 'e_i(54,0:2) = ',SUM(e_i) 361 415 ENDIF 362 416 363 417 CALL tra_sbc ( kstp ) ! surface boundary condition 364 418 365 IF(narea == 419) THEN 366 WRITE(9419,*) 'max sshn step after tra_sbc = ',MAXVAL( ABS( sshn(:,:) ) ) 367 WRITE(9419,*) 'max ssha = ',MAXVAL( ABS( ssha(:,:) ) ) 419 IF(narea == 68) THEN 420 WRITE(968,*) 'in step before tra_adv: u_ice = ',SUM(u_ice) 421 WRITE(968,*) 'v_ice(54,0:2) = ',SUM(v_ice) 422 WRITE(968,*) 'h_i(54,0:2) = ',SUM(h_i) 423 WRITE(968,*) 'h_s(54,0:2) = ',SUM(h_s) 424 WRITE(968,*) 'v_i(54,0:2) = ',SUM(v_i) 425 WRITE(968,*) 'sv_i(54,0:2) = ',SUM(sv_i) 426 WRITE(968,*) 'e_s(54,0:2) = ',SUM(e_s) 427 WRITE(968,*) 'e_i(54,0:2) = ',SUM(e_i) 368 428 ENDIF 369 429 … … 389 449 IF( ln_zdfnpc ) CALL tra_npc ( kstp ) ! update after fields by non-penetrative convection 390 450 391 IF(narea == 419) THEN 392 WRITE(9419,*) 'max sshn step before boundary conditions = ',MAXVAL( ABS( sshn(:,:) ) ) 393 WRITE(9419,*) 'max ssha = ',MAXVAL( ABS( ssha(:,:) ) ) 451 IF(narea == 68) THEN 452 WRITE(968,*) 'in step before tra_nxt: u_ice = ',SUM(u_ice) 453 WRITE(968,*) 'v_ice(54,0:2) = ',SUM(v_ice) 454 WRITE(968,*) 'h_i(54,0:2) = ',SUM(h_i) 455 WRITE(968,*) 'h_s(54,0:2) = ',SUM(h_s) 456 WRITE(968,*) 'v_i(54,0:2) = ',SUM(v_i) 457 WRITE(968,*) 'sv_i(54,0:2) = ',SUM(sv_i) 458 WRITE(968,*) 'e_s(54,0:2) = ',SUM(e_s) 459 WRITE(968,*) 'e_i(54,0:2) = ',SUM(e_i) 394 460 ENDIF 395 461 … … 413 479 CALL tra_nxt ( kstp ) ! finalize (bcs) tracer fields at next time step and swap 414 480 415 IF(narea == 419) THEN 416 WRITE(9419,*) 'max sshn step before dyn_nxt = ',MAXVAL( ABS( sshn(:,:) ) ) 417 WRITE(9419,*) 'max ssha = ',MAXVAL( ABS( ssha(:,:) ) ) 481 IF(narea == 68) THEN 482 WRITE(968,*) 'in step before dyn_nxt: u_ice = ',SUM(u_ice) 483 WRITE(968,*) 'v_ice(54,0:2) = ',SUM(v_ice) 484 WRITE(968,*) 'h_i(54,0:2) = ',SUM(h_i) 485 WRITE(968,*) 'h_s(54,0:2) = ',SUM(h_s) 486 WRITE(968,*) 'v_i(54,0:2) = ',SUM(v_i) 487 WRITE(968,*) 'sv_i(54,0:2) = ',SUM(sv_i) 488 WRITE(968,*) 'e_s(54,0:2) = ',SUM(e_s) 489 WRITE(968,*) 'e_i(54,0:2) = ',SUM(e_i) 418 490 ENDIF 419 491 420 492 CALL dyn_nxt ( kstp ) ! finalize (bcs) velocities at next time step and swap (always called after tra_nxt) 421 493 422 IF(narea == 419) THEN 423 WRITE(9419,*) 'max sshn step before ssh_swp = ',MAXVAL( ABS( sshn(:,:) ) ) 424 WRITE(9419,*) 'max ssha = ',MAXVAL( ABS( ssha(:,:) ) ) 494 IF(narea == 68) THEN 495 WRITE(968,*) 'in step before ssh_swp: u_ice = ',SUM(u_ice) 496 WRITE(968,*) 'v_ice(54,0:2) = ',SUM(v_ice) 497 WRITE(968,*) 'h_i(54,0:2) = ',SUM(h_i) 498 WRITE(968,*) 'h_s(54,0:2) = ',SUM(h_s) 499 WRITE(968,*) 'v_i(54,0:2) = ',SUM(v_i) 500 WRITE(968,*) 'sv_i(54,0:2) = ',SUM(sv_i) 501 WRITE(968,*) 'e_s(54,0:2) = ',SUM(e_s) 502 WRITE(968,*) 'e_i(54,0:2) = ',SUM(e_i) 425 503 ENDIF 426 504 … … 435 513 ! 436 514 437 IF(narea == 419) THEN 438 WRITE(9419,*) 'max sshn step before dia_hsb = ',MAXVAL( ABS( sshn(:,:) ) ) 439 WRITE(9419,*) 'max ssha = ',MAXVAL( ABS( ssha(:,:) ) ) 515 IF(narea == 68) THEN 516 WRITE(968,*) 'in step before dia_hsb: u_ice = ',SUM(u_ice) 517 WRITE(968,*) 'v_ice(54,0:2) = ',SUM(v_ice) 518 WRITE(968,*) 'h_i(54,0:2) = ',SUM(h_i) 519 WRITE(968,*) 'h_s(54,0:2) = ',SUM(h_s) 520 WRITE(968,*) 'v_i(54,0:2) = ',SUM(v_i) 521 WRITE(968,*) 'sv_i(54,0:2) = ',SUM(sv_i) 522 WRITE(968,*) 'e_s(54,0:2) = ',SUM(e_s) 523 WRITE(968,*) 'e_i(54,0:2) = ',SUM(e_i) 440 524 ENDIF 441 525 … … 450 534 IF( ln_sto_eos ) CALL sto_rst_write( kstp ) ! write restart file for stochastic parameters 451 535 452 IF(narea == 419) THEN 453 WRITE(9419,*) 'max sshn step before agrif = ',MAXVAL( ABS( sshn(:,:) ) ) 536 IF(narea == 68) THEN 537 WRITE(968,*) 'in step before agrif: u_ice = ',SUM(u_ice) 538 WRITE(968,*) 'v_ice(54,0:2) = ',SUM(v_ice) 539 WRITE(968,*) 'h_i(54,0:2) = ',SUM(h_i) 540 WRITE(968,*) 'h_s(54,0:2) = ',SUM(h_s) 541 WRITE(968,*) 'v_i(54,0:2) = ',SUM(v_i) 542 WRITE(968,*) 'sv_i(54,0:2) = ',SUM(sv_i) 543 WRITE(968,*) 'e_s(54,0:2) = ',SUM(e_s) 544 WRITE(968,*) 'e_i(54,0:2) = ',SUM(e_i) 454 545 ENDIF 455 546 … … 468 559 !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< 469 560 470 IF(narea == 419) THEN 471 WRITE(9419,*) 'max sshn step before stp_ctl = ',MAXVAL( ABS( sshn(:,:) ) ) 561 IF(narea == 68) THEN 562 WRITE(968,*) 'in step before stp_ctl: u_ice = ',SUM(u_ice) 563 WRITE(968,*) 'v_ice(54,0:2) = ',SUM(v_ice) 564 WRITE(968,*) 'h_i(54,0:2) = ',SUM(h_i) 565 WRITE(968,*) 'h_s(54,0:2) = ',SUM(h_s) 566 WRITE(968,*) 'v_i(54,0:2) = ',SUM(v_i) 567 WRITE(968,*) 'sv_i(54,0:2) = ',SUM(sv_i) 568 WRITE(968,*) 'e_s(54,0:2) = ',SUM(e_s) 569 WRITE(968,*) 'e_i(54,0:2) = ',SUM(e_i) 472 570 ENDIF 473 571
Note: See TracChangeset
for help on using the changeset viewer.