Ignore:
Timestamp:
2018-01-18T14:47:37+01:00 (3 years ago)
Author:
jpalmier
Message:

JPALM — GMED 364 — apply Richard review suggestion

Location:
branches/NERC/dev_r5518_GO6_conserv_check_up/NEMOGCM/NEMO/TOP_SRC
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • branches/NERC/dev_r5518_GO6_conserv_check_up/NEMOGCM/NEMO/TOP_SRC/trcini.F90

    r9208 r9260  
    218218      IF ( lwp ) WRITE(numout,*) 
    219219      IF ( lwp ) WRITE(numout,*)    ' Elemental cycle totals: ' 
    220       IF ( lwp ) CALL flush(numout) 
    221220      ! nitrogen 
    222221      ztot3d(:,:,:) = trn(:,:,:,jpphn) + trn(:,:,:,jpphd) + trn(:,:,:,jpzmi) + & 
     
    227226      cycletot(1)   = zsum3d + zsum2d 
    228227      IF ( lwp ) WRITE(numout,9010) 'nitrogen', zsum3d, zsum2d, cycletot(1) 
    229       IF ( lwp ) CALL flush(numout) 
    230228      ! silicon 
    231229      ztot3d(:,:,:) = trn(:,:,:,jppds) + trn(:,:,:,jpsil) 
     
    235233      cycletot(2)   = zsum3d + zsum2d 
    236234      IF ( lwp ) WRITE(numout,9010) 'silicon', zsum3d, zsum2d, cycletot(2) 
    237       IF ( lwp ) CALL flush(numout) 
    238235      ! iron 
    239236      ztot3d(:,:,:) = ((trn(:,:,:,jpphn) + trn(:,:,:,jpphd) + trn(:,:,:,jpzmi) + & 
     
    244241      cycletot(3)   = zsum3d + zsum2d 
    245242      IF ( lwp ) WRITE(numout,9010) 'iron', zsum3d, zsum2d, cycletot(3) 
    246       IF ( lwp ) CALL flush(numout) 
    247243      ! carbon (uses fixed C:N ratios on plankton tracers) 
    248244      ztot3d(:,:,:) = (trn(:,:,:,jpphn) * xthetapn)  + (trn(:,:,:,jpphd) * xthetapd)  +  & 
     
    254250      cycletot(4)   = zsum3d + zsum2d 
    255251      IF ( lwp ) WRITE(numout,9010) 'carbon', zsum3d, zsum2d, cycletot(4) 
    256       IF ( lwp ) CALL flush(numout) 
    257252      ! alkalinity (note benthic correction) 
    258253      ztot3d(:,:,:) = trn(:,:,:,jpalk) 
     
    262257      cycletot(5)   = zsum3d + zsum2d 
    263258      IF ( lwp ) WRITE(numout,9010) 'alkalinity', zsum3d, zsum2d, cycletot(5) 
    264       IF ( lwp ) CALL flush(numout) 
    265259      ! oxygen (note no benthic) 
    266260      ztot3d(:,:,:) = trn(:,:,:,jpoxy) 
     
    270264      cycletot(6)   = zsum3d + zsum2d 
    271265      IF ( lwp ) WRITE(numout,9010) 'oxygen', zsum3d, zsum2d, cycletot(6) 
    272       IF ( lwp ) CALL flush(numout) 
    273266      ! Check 
    274267      zsum3d        = glob_sum( cvol(:,:,:) ) 
    275268      zsum2d        = glob_sum( carea(:,:) )       
    276       IF ( lwp ) WRITE(numout,*) 
    277       IF ( lwp ) WRITE(numout,*) ' check : cvol    : ', zsum3d 
    278       IF ( lwp ) WRITE(numout,*) ' check : carea   : ', zsum2d 
    279       IF ( lwp ) WRITE(numout,*) 
    280       IF ( lwp ) CALL flush(numout) 
    281       ! 
    282 # endif 
    283  
    284       IF(lwp) WRITE(numout,*) 
    285       IF(lwp) WRITE(numout,*) 'trc_init : passive tracer set up completed' 
    286       IF(lwp) WRITE(numout,*) '~~~~~~~' 
    287       IF(lwp) CALL flush(numout) 
     269      IF ( lwp ) THEN 
     270         WRITE(numout,*) 
     271         WRITE(numout,*) ' check : cvol    : ', zsum3d 
     272         WRITE(numout,*) ' check : carea   : ', zsum2d 
     273         WRITE(numout,*) 
     274      ENDIF 
     275      ! 
     276# endif 
     277 
     278      IF(lwp) THEN  
     279          WRITE(numout,*) 
     280          WRITE(numout,*) 'trc_init : passive tracer set up completed' 
     281          WRITE(numout,*) '~~~~~~~' 
     282      ENDIF  
    288283# if defined key_debug_medusa 
    289284         CALL trc_rst_stat 
  • branches/NERC/dev_r5518_GO6_conserv_check_up/NEMOGCM/NEMO/TOP_SRC/trcrst.F90

    r9208 r9260  
    844844      zsum3d        = glob_sum( zvol(:,:,:) ) 
    845845      zsum2d        = glob_sum( zarea(:,:) ) 
    846       IF ( lwp ) WRITE(numout,*) 
    847       IF ( lwp ) WRITE(numout,*) ' check : cvol    : ', zsum3d 
    848       IF ( lwp ) WRITE(numout,*) ' check : carea   : ', zsum2d 
    849       IF ( lwp ) WRITE(numout,*) 
    850       IF ( lwp ) CALL flush(numout) 
     846      IF ( lwp ) THEN  
     847         WRITE(numout,*) 
     848         WRITE(numout,*) ' check : cvol    : ', zsum3d 
     849         WRITE(numout,*) ' check : carea   : ', zsum2d 
     850         WRITE(numout,*) 
     851      ENDIF 
    851852      ! 
    8528539010  FORMAT(' element:',a10,                     & 
Note: See TracChangeset for help on using the changeset viewer.