Changeset 10753
- Timestamp:
- 2019-03-14T16:34:25+01:00 (5 years ago)
- Location:
- branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/ICB
- Files:
-
- 9 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/ICB/icbclv.F90
r8046 r10753 56 56 INTEGER :: imx ! temporary integer for max berg class 57 57 LOGICAL, SAVE :: ll_first_call = .TRUE. 58 LOGICAL :: ll_write 58 59 !!---------------------------------------------------------------------- 59 60 ! … … 68 69 ! nn_coupled_iceshelf_fluxes uninitialised unless lk_oasis=true 69 70 IF( nn_coupled_iceshelf_fluxes .gt. 0 ) THEN 70 71 ll_write = ( MOD( kt, sn_cfctl%ptimincr ) == 0 ) .OR. ( kt == nitend ) .AND. lwp 71 72 ! Adjust total calving rates so that sum of iceberg calving and iceshelf melting in the northern 72 73 ! and southern hemispheres equals rate of increase of mass of greenland and antarctic ice sheets … … 80 81 81 82 ! check 82 IF(l wp) WRITE(numout, *) 'Greenland iceberg calving climatology (kg/s) : ',zgreenland_calving_sum83 IF(ll_write) WRITE(numout, *) 'Greenland iceberg calving climatology (kg/s) : ',zgreenland_calving_sum 83 84 zgreenland_calving_sum = SUM( berg_grid%calving(:,:) * greenland_icesheet_mask(:,:) ) 84 85 IF( lk_mpp ) CALL mpp_sum( zgreenland_calving_sum ) 85 IF(l wp) WRITE(numout, *) 'Greenland iceberg calving adjusted value (kg/s) : ',zgreenland_calving_sum86 IF(ll_write) WRITE(numout, *) 'Greenland iceberg calving adjusted value (kg/s) : ',zgreenland_calving_sum 86 87 87 88 zantarctica_calving_sum = SUM( berg_grid%calving(:,:) * antarctica_icesheet_mask(:,:) ) … … 92 93 93 94 ! check 94 IF(l wp) WRITE(numout, *) 'Antarctica iceberg calving climatology (kg/s) : ',zantarctica_calving_sum95 IF(ll_write) WRITE(numout, *) 'Antarctica iceberg calving climatology (kg/s) : ',zantarctica_calving_sum 95 96 zantarctica_calving_sum = SUM( berg_grid%calving(:,:) * antarctica_icesheet_mask(:,:) ) 96 97 IF( lk_mpp ) CALL mpp_sum( zantarctica_calving_sum ) 97 IF(l wp) WRITE(numout, *) 'Antarctica iceberg calving adjusted value (kg/s) : ',zantarctica_calving_sum98 IF(ll_write) WRITE(numout, *) 'Antarctica iceberg calving adjusted value (kg/s) : ',zantarctica_calving_sum 98 99 99 100 ENDIF … … 218 219 CALL lbc_lnk( berg_grid%stored_heat, 'T', 1._wp ) 219 220 ! 220 IF( nn_verbose_level > 0 .AND. icntmax > 1 ) WRITE(numicb,*) 'icb_clv: icnt=', icnt,' on', narea221 IF( nn_verbose_level > 0 .AND. icntmax > 1 .AND. numicb.NE.-1) WRITE(numicb,*) 'icb_clv: icnt=', icnt,' on', narea 221 222 ! 222 223 END SUBROUTINE icb_clv -
branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/ICB/icbdia.F90
r6491 r10753 311 311 'returned',calving_ret_net) 312 312 ENDIF 313 WRITE( numicb, '("calved by class = ",i6,20(",",i6))') (nbergs_calved_by_class(ik),ik=1,nclasses) 314 IF ( nspeeding_tickets > 0 ) WRITE( numicb, '("speeding tickets issued = ",i6)') nspeeding_tickets 313 IF(numicb.NE.-1) THEN 314 WRITE( numicb, '("calved by class = ",i6,20(",",i6))') (nbergs_calved_by_class(ik),ik=1,nclasses) 315 IF ( nspeeding_tickets > 0 ) WRITE( numicb, '("speeding tickets issued = ",i6)') nspeeding_tickets 316 ENDIF 315 317 316 318 nbergs_start = nbergs_end … … 418 420 stored_start = SUM( berg_grid%stored_ice(:,:,:) ) 419 421 IF( lk_mpp ) CALL mpp_sum( stored_start ) 420 WRITE(numicb,'(a,es13.6,a)') 'icb_dia_income: initial stored mass=',stored_start,' kg'422 IF(numicb.NE.-1) WRITE(numicb,'(a,es13.6,a)') 'icb_dia_income: initial stored mass=',stored_start,' kg' 421 423 ! 422 424 stored_heat_start = SUM( berg_grid%stored_heat(:,:) ) 423 425 IF( lk_mpp ) CALL mpp_sum( stored_heat_start ) 424 WRITE(numicb,'(a,es13.6,a)') 'icb_dia_income: initial stored heat=',stored_heat_start,' J'426 IF(numicb.NE.-1) WRITE(numicb,'(a,es13.6,a)') 'icb_dia_income: initial stored heat=',stored_heat_start,' J' 425 427 ENDIF 426 428 ! … … 495 497 ! 496 498 IF ( PRESENT(kbergs) ) THEN 497 WRITE(numicb,100) cd_budgetstr // ' state:', & 499 IF(numicb.NE.-1) & 500 & WRITE(numicb,100) cd_budgetstr // ' state:', & 498 501 cd_startstr // ' start', pstartval, cd_budgetunits, & 499 502 cd_endstr // ' end', pendval, cd_budgetunits, & … … 501 504 '# of bergs', kbergs 502 505 ELSE 503 WRITE(numicb,100) cd_budgetstr // ' state:', & 506 IF(numicb.NE.-1) & 507 & WRITE(numicb,100) cd_budgetstr // ' state:', & 504 508 cd_startstr // ' start', pstartval, cd_budgetunits, & 505 509 cd_endstr // ' end', pendval, cd_budgetunits, & … … 517 521 !!---------------------------------------------------------------------- 518 522 ! 519 WRITE(numicb,200) cd_budgetstr // ' check:', & 523 IF(numicb.NE.-1) & 524 & WRITE(numicb,200) cd_budgetstr // ' check:', & 520 525 cd_startstr, pstartval, cd_budgetunits, & 521 526 cd_endstr, pendval, cd_budgetunits, & … … 538 543 & MAX( 1.e-30, MAX( abs( pendval - pstartval ) , ABS( pinval - poutval ) ) ) 539 544 540 WRITE(numicb,200) cd_budgetstr // ' budget:', & 545 IF(numicb.NE.-1) & 546 & WRITE(numicb,200) cd_budgetstr // ' budget:', & 541 547 & cd_instr // ' in', pinval, cd_budgetunits, & 542 548 & cd_outstr // ' out', poutval, cd_budgetunits, & … … 554 560 INTEGER, INTENT(in) :: pstartval, pendval 555 561 ! 556 WRITE(numicb,100) cd_budgetstr // ' state:', & 562 IF(numicb.NE.-1) & 563 & WRITE(numicb,100) cd_budgetstr // ' state:', & 557 564 & cd_startstr // ' start', pstartval, & 558 565 & cd_endstr // ' end', pendval, & … … 571 578 !!---------------------------------------------------------------------- 572 579 ! 573 WRITE(numicb,200) cd_budgetstr // ' budget:', & 580 IF(numicb.NE.-1) & 581 & WRITE(numicb,200) cd_budgetstr // ' budget:', & 574 582 cd_instr // ' in', pinval, & 575 583 cd_outstr // ' out', poutval, & -
branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/ICB/icbdyn.F90
r6486 r10753 372 372 ! ! check the speed and acceleration limits 373 373 IF( ABS( zuveln ) > pp_vel_lim .OR. ABS( zvveln ) > pp_vel_lim ) & 374 WRITE( numicb,'("pe=",i3,x,a)') narea,'Dump triggered by excessive velocity'374 WRITE(*,'("pe=",i3,x,a)') narea,'Dump triggered by excessive velocity' 375 375 IF( ABS( pax ) > pp_accel_lim .OR. ABS( pay ) > pp_accel_lim ) & 376 WRITE( numicb,'("pe=",i3,x,a)') narea,'Dump triggered by excessive acceleration'376 WRITE(*,'("pe=",i3,x,a)') narea,'Dump triggered by excessive acceleration' 377 377 ! 378 378 END SUBROUTINE icb_accel -
branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/ICB/icbini.F90
r6823 r10753 78 78 ! ! note that we choose to do this on all processors since we cannot 79 79 ! ! predict where icebergs will be ahead of time 80 CALL ctl_opn( numicb, 'icebergs.stat', 'REPLACE', 'FORMATTED', 'SEQUENTIAL', -1, numout, lwp, narea ) 80 numicb=-1 81 IF(nn_verbose_level>0 .AND. nprint>0) THEN 82 CALL ctl_opn( numicb, 'icebergs.stat', 'REPLACE', 'FORMATTED', 'SEQUENTIAL', -1, numout, lwp, narea ) 83 ELSE 84 IF(lwp) CALL ctl_opn( numicb, 'icebergs.stat', 'REPLACE', 'FORMATTED', 'SEQUENTIAL', -1, numout, lwp, narea ) 85 ENDIF 81 86 82 87 ! set parameters (mostly from namelist) … … 187 192 ENDIF 188 193 ! 189 IF( nn_verbose_level > 0 ) THEN194 IF( nn_verbose_level > 0 .AND. numicb.NE.-1) THEN 190 195 WRITE(numicb,*) 'processor ', narea 191 196 WRITE(numicb,*) 'jpi, jpj ', jpi, jpj … … 241 246 CALL iom_close( inum ) ! close file 242 247 ! 243 WRITE(numicb,*) 244 WRITE(numicb,*) ' calving read in a file' 248 IF(numicb.NE.-1) THEN 249 WRITE(numicb,*) 250 WRITE(numicb,*) ' calving read in a file' 251 ENDIF 245 252 ALLOCATE( sf_icb(1), STAT=istat1 ) ! Create sf_icb structure (calving) 246 253 ALLOCATE( sf_icb(1)%fnow(jpi,jpj,1), STAT=istat2 ) … … 336 343 ibergs = icb_utl_count() 337 344 IF( lk_mpp ) CALL mpp_sum(ibergs) 338 WRITE(numicb,'(a,i6,a)') 'diamonds, icb_ini_gen: ',ibergs,' were generated'345 IF(numicb.NE.-1) WRITE(numicb,'(a,i6,a)') 'diamonds, icb_ini_gen: ',ibergs,' were generated' 339 346 ! 340 347 END SUBROUTINE icb_ini_gen … … 368 375 READ ( numnam_cfg, namberg, IOSTAT = ios, ERR = 902 ) 369 376 902 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namberg in configuration namelist', lwp ) 370 IF(lwm ) WRITE ( numond, namberg )377 IF(lwm .AND. nprint>3) WRITE ( numond, namberg ) 371 378 #else 372 379 IF(lwp) THEN -
branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/ICB/icblbc.F90
r6486 r10753 197 197 ENDIF 198 198 199 IF( nn_verbose_level >= 2 ) THEN199 IF( nn_verbose_level >= 2 .AND. numicb .NE. -1) THEN 200 200 WRITE(numicb,*) 'processor west : ', ipe_W 201 201 WRITE(numicb,*) 'processor east : ', ipe_E … … 240 240 this => this%next 241 241 ibergs_to_send_e = ibergs_to_send_e + 1 242 IF( nn_verbose_level >= 4 ) THEN242 IF( nn_verbose_level >= 4 .AND. numicb .NE. -1) THEN 243 243 WRITE(numicb,*) 'bergstep ',nktberg,' packing berg ',tmpberg%number(:),' for transfer to east' 244 244 CALL flush( numicb ) … … 253 253 this => this%next 254 254 ibergs_to_send_w = ibergs_to_send_w + 1 255 IF( nn_verbose_level >= 4 ) THEN255 IF( nn_verbose_level >= 4 .AND. numicb.NE.-1) THEN 256 256 WRITE(numicb,*) 'bergstep ',nktberg,' packing berg ',tmpberg%number(:),' for transfer to west' 257 257 CALL flush( numicb ) … … 267 267 END DO 268 268 ENDIF 269 IF( nn_verbose_level >= 3 ) THEN269 IF( nn_verbose_level >= 3 .AND. numicb .NE. -1) THEN 270 270 WRITE(numicb,*) 'bergstep ',nktberg,' send ew: ', ibergs_to_send_e, ibergs_to_send_w 271 271 CALL flush(numicb) … … 301 301 ibergs_rcvd_from_w = INT( zwebergs(2) ) 302 302 END SELECT 303 IF( nn_verbose_level >= 3 ) THEN303 IF( nn_verbose_level >= 3 .AND. numicb .NE. -1) THEN 304 304 WRITE(numicb,*) 'bergstep ',nktberg,' recv ew: ', ibergs_rcvd_from_w, ibergs_rcvd_from_e 305 305 CALL flush(numicb) … … 315 315 IF( ibergs_to_send_e > 0 .AND. l_isend ) CALL mpi_wait( iml_req1, iml_stat, iml_err ) 316 316 DO i = 1, ibergs_rcvd_from_e 317 IF( nn_verbose_level >= 4 ) THEN317 IF( nn_verbose_level >= 4 .AND. numicb .NE. -1) THEN 318 318 WRITE(numicb,*) 'bergstep ',nktberg,' unpacking berg ',INT(ibuffer_e%data(16,i)),' from east' 319 319 CALL flush( numicb ) … … 335 335 IF( ibergs_to_send_e > 0 .AND. l_isend ) CALL mpi_wait( iml_req3, iml_stat, iml_err ) 336 336 DO i = 1, ibergs_rcvd_from_e 337 IF( nn_verbose_level >= 4 ) THEN337 IF( nn_verbose_level >= 4 .AND. numicb .NE. -1) THEN 338 338 WRITE(numicb,*) 'bergstep ',nktberg,' unpacking berg ',INT(ibuffer_e%data(16,i)),' from east' 339 339 CALL flush( numicb ) … … 342 342 END DO 343 343 DO i = 1, ibergs_rcvd_from_w 344 IF( nn_verbose_level >= 4 ) THEN344 IF( nn_verbose_level >= 4 .AND. numicb .NE. -1) THEN 345 345 WRITE(numicb,*) 'bergstep ',nktberg,' unpacking berg ',INT(ibuffer_w%data(16,i)),' from west' 346 346 CALL flush( numicb ) … … 356 356 IF( ibergs_to_send_w > 0 .AND. l_isend ) CALL mpi_wait( iml_req4, iml_stat, iml_err ) 357 357 DO i = 1, ibergs_rcvd_from_w 358 IF( nn_verbose_level >= 4 ) THEN358 IF( nn_verbose_level >= 4 .AND. numicb .NE. -1) THEN 359 359 WRITE(numicb,*) 'bergstep ',nktberg,' unpacking berg ',INT(ibuffer_w%data(16,i)),' from west' 360 360 CALL flush( numicb ) … … 378 378 this => this%next 379 379 ibergs_to_send_n = ibergs_to_send_n + 1 380 IF( nn_verbose_level >= 4 ) THEN380 IF( nn_verbose_level >= 4 .AND. numicb .NE. -1) THEN 381 381 WRITE(numicb,*) 'bergstep ',nktberg,' packing berg ',tmpberg%number(:),' for transfer to north' 382 382 CALL flush( numicb ) … … 388 388 this => this%next 389 389 ibergs_to_send_s = ibergs_to_send_s + 1 390 IF( nn_verbose_level >= 4 ) THEN390 IF( nn_verbose_level >= 4 .AND. numicb .NE. -1) THEN 391 391 WRITE(numicb,*) 'bergstep ',nktberg,' packing berg ',tmpberg%number(:),' for transfer to south' 392 392 CALL flush( numicb ) … … 399 399 END DO 400 400 ENDIF 401 if( nn_verbose_level >= 3 ) then401 if( nn_verbose_level >= 3 .AND. numicb .NE. -1) then 402 402 write(numicb,*) 'bergstep ',nktberg,' send ns: ', ibergs_to_send_n, ibergs_to_send_s 403 403 call flush(numicb) … … 432 432 ibergs_rcvd_from_s = INT( zsnbergs(2) ) 433 433 END SELECT 434 if( nn_verbose_level >= 3 ) then434 if( nn_verbose_level >= 3 .AND. numicb .NE. -1) then 435 435 write(numicb,*) 'bergstep ',nktberg,' recv ns: ', ibergs_rcvd_from_s, ibergs_rcvd_from_n 436 436 call flush(numicb) … … 446 446 IF( ibergs_to_send_n > 0 .AND. l_isend ) CALL mpi_wait( iml_req1, iml_stat, iml_err ) 447 447 DO i = 1, ibergs_rcvd_from_n 448 IF( nn_verbose_level >= 4 ) THEN448 IF( nn_verbose_level >= 4 .AND. numicb .NE. -1) THEN 449 449 WRITE(numicb,*) 'bergstep ',nktberg,' unpacking berg ',INT(ibuffer_n%data(16,i)),' from north' 450 450 CALL flush( numicb ) … … 466 466 IF( ibergs_to_send_n > 0 .AND. l_isend ) CALL mpi_wait( iml_req3, iml_stat, iml_err ) 467 467 DO i = 1, ibergs_rcvd_from_n 468 IF( nn_verbose_level >= 4 ) THEN468 IF( nn_verbose_level >= 4 .AND. numicb .NE. -1) THEN 469 469 WRITE(numicb,*) 'bergstep ',nktberg,' unpacking berg ',INT(ibuffer_n%data(16,i)),' from north' 470 470 CALL flush( numicb ) … … 473 473 END DO 474 474 DO i = 1, ibergs_rcvd_from_s 475 IF( nn_verbose_level >= 4 ) THEN475 IF( nn_verbose_level >= 4 .AND. numicb .NE. -1) THEN 476 476 WRITE(numicb,*) 'bergstep ',nktberg,' unpacking berg ',INT(ibuffer_s%data(16,i)),' from south' 477 477 CALL flush( numicb ) … … 487 487 IF( ibergs_to_send_s > 0 .AND. l_isend ) CALL mpi_wait( iml_req4, iml_stat, iml_err ) 488 488 DO i = 1, ibergs_rcvd_from_s 489 IF( nn_verbose_level >= 4 ) THEN489 IF( nn_verbose_level >= 4 .AND. numicb .NE. -1) THEN 490 490 WRITE(numicb,*) 'bergstep ',nktberg,' unpacking berg ',INT(ibuffer_s%data(16,i)),' from south' 491 491 CALL flush( numicb ) … … 495 495 END SELECT 496 496 497 IF( nn_verbose_level > 0 ) THEN497 IF( nn_verbose_level > 0 .AND. numicb .NE. -1) THEN 498 498 ! compare the number of icebergs on this processor from the start to the end 499 499 ibergs_end = icb_utl_count() … … 535 535 IF( npolj /= 0 .AND. jpni > 1 ) CALL icb_lbc_mpp_nfld( ) 536 536 537 IF( nn_verbose_level > 0 537 IF( nn_verbose_level > 0) THEN 538 538 i = 0 539 539 this => first_berg … … 547 547 ijne .GT. mjg(nicbej)) THEN 548 548 i = i + 1 549 WRITE(numicb,*) 'berg lost in halo: ', this%number(:),iine,ijne 550 WRITE(numicb,*) ' ', nimpp, njmpp 551 WRITE(numicb,*) ' ', nicbdi, nicbei, nicbdj, nicbej 552 CALL flush( numicb ) 549 IF(numicb .NE. -1) THEN 550 WRITE(numicb,*) 'berg lost in halo: ', this%number(:),iine,ijne 551 WRITE(numicb,*) ' ', nimpp, njmpp 552 WRITE(numicb,*) ' ', nicbdi, nicbei, nicbdj, nicbej 553 CALL flush( numicb ) 554 ENDIF 553 555 ENDIF 554 556 this => this%next … … 659 661 IF( nicbfldproc(jjn) .eq. INT(znbergs(1)) ) EXIT 660 662 END DO 661 IF( jjn .GT. jpni ) write( numicb,*) 'ICB ERROR'663 IF( jjn .GT. jpni ) write(*,*) 'ICB ERROR' 662 664 nicbfldexpect(jjn) = INT( znbergs(2) ) 663 665 !IF ( nicbfldexpect(jjn) .GT. 0) write(numicb,*) 'ICB expecting ',nicbfldexpect(jjn),' from ', nicbfldproc(jjn) … … 714 716 tmpberg => this 715 717 ibergs_to_send = ibergs_to_send + 1 716 IF( nn_verbose_level >= 4 ) THEN718 IF( nn_verbose_level >= 4 .AND. numicb.NE.-1) THEN 717 719 WRITE(numicb,*) 'bergstep ',nktberg,' packing berg ',tmpberg%number(:),' for north fold' 718 720 CALL flush( numicb ) … … 727 729 END DO 728 730 ENDIF 729 if( nn_verbose_level >= 3 ) then731 if( nn_verbose_level >= 3 .AND. numicb.NE.-1) then 730 732 write(numicb,*) 'bergstep ',nktberg,' send nfld: ', ibergs_to_send 731 733 call flush(numicb) … … 758 760 ! 759 761 DO jk = 1, ibergs_to_rcv 760 IF( nn_verbose_level >= 4 ) THEN762 IF( nn_verbose_level >= 4 .AND. numicb.NE.-1) THEN 761 763 WRITE(numicb,*) 'bergstep ',nktberg,' unpacking berg ',INT(ibuffer_f%data(16,jk)),' from north fold' 762 764 CALL flush( numicb ) -
branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/ICB/icbrst.F90
r9321 r10753 134 134 ! Sanity check 135 135 jn = icb_utl_count() 136 IF (n n_verbose_level >= 0) &136 IF (nprint > 0 .AND. lwp) & 137 137 WRITE(numout,'(2(a,i5))') 'icebergs, read_restart_bergs: # bergs =',jn,' on PE',narea-1 138 138 IF( lk_mpp ) THEN … … 141 141 CALL mpp_sum(jn) 142 142 ENDIF 143 IF(lwp) WRITE(numout,'(a,i5,a,i5,a)') 'icebergs, read_restart_bergs: there were',ibergs_in_file, & 143 IF(lwp .AND. nprint > 1) & 144 & WRITE(numout,'(a,i5,a,i5,a)') 'icebergs, read_restart_bergs: there were',ibergs_in_file, & 144 145 & ' bergs in the restart file and', jn,' bergs have been read' 145 146 ! … … 147 148 CALL iom_close( ncid ) 148 149 ! 149 IF(nn_timing == 2) CALL timing_stop('iom_rstget')150 IF( lwp .and. n n_verbose_level>= 0) WRITE(numout,'(a)') 'icebergs, read_restart_bergs: completed'150 IF(nn_timing == 1) CALL timing_stop('iom_rstget') 151 IF( lwp .and. nprint >= 0) WRITE(numout,'(a)') 'icebergs, read_restart_bergs: completed' 151 152 ! 152 153 END SUBROUTINE icb_rst_read … … 190 191 WRITE(cl_filename,'(A,"_icebergs_",A,"_restart.nc")') TRIM(cexper), TRIM(ADJUSTL(clkt)) 191 192 ENDIF 192 IF (nn_verbose_level >= 0) WRITE(numout,'(2a)') 'icebergs, write_restart: creating ',TRIM(cl_path)//TRIM(cl_filename) 193 IF (nprint > 0 .AND. lwp) & 194 WRITE(numout,'(2a)') 'icebergs, write_restart: creating ',TRIM(cl_path)//TRIM(cl_filename) 193 195 194 196 nret = NF90_CREATE(TRIM(cl_path)//TRIM(cl_filename), NF90_CLOBBER, ncid) … … 324 326 ENDIF 325 327 ENDDO 326 IF( lwp ) WRITE(numout,*) 'file: ',TRIM(cl_path)//TRIM(cl_filename),' var: stored_ice written'328 IF( lwp .AND. nprint > 1) WRITE(numout,*) 'file: ',TRIM(cl_path)//TRIM(cl_filename),' var: stored_ice written' 327 329 328 330 nret = NF90_PUT_VAR( ncid, nkountid, num_bergs(:) ) … … 331 333 nret = NF90_PUT_VAR( ncid, nsheatid, berg_grid%stored_heat(:,:) ) 332 334 IF (nret .ne. NF90_NOERR) CALL ctl_stop('icebergs, write_restart: nf_put_var stored_heat failed') 333 IF( lwp ) WRITE(numout,*) 'file: ',TRIM(cl_path)//TRIM(cl_filename),' var: stored_heat written'335 IF( lwp .AND. nprint > 1) WRITE(numout,*) 'file: ',TRIM(cl_path)//TRIM(cl_filename),' var: stored_heat written' 334 336 335 337 nret = NF90_PUT_VAR( ncid, ncalvid , src_calving(:,:) ) … … 337 339 nret = NF90_PUT_VAR( ncid, ncalvhid, src_calving_hflx(:,:) ) 338 340 IF (nret .ne. NF90_NOERR) CALL ctl_stop('icebergs, write_restart: nf_put_var calving_hflx failed') 339 IF( lwp ) WRITE(numout,*) 'file: ',TRIM(cl_path)//TRIM(cl_filename),' var: calving written'341 IF( lwp .AND. nprint > 1) WRITE(numout,*) 'file: ',TRIM(cl_path)//TRIM(cl_filename),' var: calving written' 340 342 341 343 IF ( ASSOCIATED(first_berg) ) THEN -
branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/ICB/icbstp.F90
r6486 r10753 88 88 89 89 ! write out time 90 IF( ll_verbose ) WRITE(numicb,9100) nktberg, ndastp, nsec_day90 IF( ll_verbose .AND. numicb.NE.-1) WRITE(numicb,9100) nktberg, ndastp, nsec_day 91 91 9100 FORMAT('kt= ',i8, ' day= ',i8,' secs=',i8) 92 92 … … 166 166 IF( nn_sample_rate .GT. 0 ) CALL icb_trj_end() 167 167 168 IF(lwp) WRITE(numout,'(a,i6)') 'icebergs: icb_end complete', narea 169 CALL flush( numicb ) 170 CLOSE( numicb ) 168 IF(lwp .AND. nprint>3) WRITE(numout,'(a,i6)') 'icebergs: icb_end complete', narea 169 IF(numicb.NE.-1) THEN 170 CALL flush( numicb ) 171 CLOSE( numicb ) 172 ENDIF 171 173 ! 172 174 END SUBROUTINE icb_end -
branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/ICB/icbtrj.F90
r7733 r10753 83 83 ELSE ; WRITE(cl_filename,'("trajectory_icebergs_",A,"-",A ,".nc")') TRIM(ADJUSTL(cldate_ini)), TRIM(ADJUSTL(cldate_end)) 84 84 ENDIF 85 IF ( lwp .AND. n n_verbose_level >= 0) WRITE(numout,'(2a)') 'icebergs, icb_trj_init: creating ',TRIM(cl_filename)85 IF ( lwp .AND. nprint > 2) WRITE(numout,'(2a)') 'icebergs, icb_trj_init: creating ',TRIM(cl_filename) 86 86 87 87 iret = NF90_CREATE(TRIM(cl_filename), NF90_CLOBBER, ntrajid) … … 246 246 this=>this%next 247 247 END DO 248 IF( lwp . and. nn_verbose_level> 0 ) WRITE(numout,*) 'trajectory write to frame ', jn248 IF( lwp .AND. nprint > 0 ) WRITE(numout,*) 'trajectory write to frame ', jn 249 249 num_traj = jn 250 250 ! -
branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/ICB/icbutl.F90
r6486 r10753 661 661 ! 662 662 pt => berg%current_point 663 IF(numicb.NE.-1) THEN 663 664 WRITE(numicb, 9200) kt, berg%number(1), & 664 665 pt%xi, pt%yj, pt%lon, pt%lat, pt%uvel, pt%vvel, & 665 666 pt%uo, pt%vo, pt%ua, pt%va, pt%ui, pt%vi 666 667 CALL flush( numicb ) 668 ENDIF 667 669 9200 FORMAT(5x,i5,2x,i10,6(2x,2f10.4)) 668 670 ! … … 685 687 ! 686 688 this => first_berg 687 IF( ASSOCIATED(this) ) THEN 688 WRITE(numicb,'(a," pe=(",i3,")")' ) cd_label, narea 689 WRITE(numicb,'(a8,4x,a6,12x,a5,15x,a7,19x,a3,17x,a5,17x,a5,17x,a5)' ) & 689 IF(numicb.NE.-1) THEN 690 IF( ASSOCIATED(this) ) THEN 691 WRITE(numicb,'(a," pe=(",i3,")")' ) cd_label, narea 692 WRITE(numicb,'(a8,4x,a6,12x,a5,15x,a7,19x,a3,17x,a5,17x,a5,17x,a5)' ) & 690 693 & 'timestep', 'number', 'xi,yj','lon,lat','u,v','uo,vo','ua,va','ui,vi' 694 ENDIF 691 695 ENDIF 692 696 DO WHILE( ASSOCIATED(this) ) … … 697 701 inbergs = ibergs 698 702 IF( lk_mpp ) CALL mpp_sum(inbergs) 699 IF( ibergs > 0 ) WRITE(numicb,'(a," there are",i5," bergs out of",i6," on PE ",i4)') &703 IF( ibergs > 0 .AND. numicb.NE.-1) WRITE(numicb,'(a," there are",i5," bergs out of",i6," on PE ",i4)') & 700 704 & cd_label, ibergs, inbergs, narea 701 705 !
Note: See TracChangeset
for help on using the changeset viewer.