Changeset 9314


Ignore:
Timestamp:
2018-02-09T12:03:12+01:00 (2 years ago)
Author:
dancopsey
Message:

Added logical to turn on extra diagnostic printing.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/UKMO/dev_r5518_new_runoff_coupling/NEMOGCM/NEMO/OPA_SRC/SBC/cpl_rnf_1d.F90

    r9313 r9314  
    3838   INTEGER, PUBLIC            :: nn_cpl_river   ! Maximum number of rivers being passed through the coupler 
    3939   INTEGER, PUBLIC            :: runoff_id      ! OASIS coupling id used in oasis_get command 
     40   LOGICAL                    :: ln_print_river_info  ! Diagnostic prints of river coupling information 
    4041    
    4142CONTAINS 
     
    5960      REAL(wp), POINTER, DIMENSION(:,:)         ::   river_number        ! 2D array containing the river outflow numbers 
    6061       
    61       NAMELIST/nam_cpl_rnf_1d/file_riv_number, nn_cpl_river 
     62      NAMELIST/nam_cpl_rnf_1d/file_riv_number, nn_cpl_river, ln_print_river_info 
    6263      !!---------------------------------------------------------------------- 
    6364 
     
    8182      IF(lwp) WRITE(numout,*) '          Namelist nam_cpl_rnf_1d : Coupled runoff using 1D array' 
    8283      IF(lwp) WRITE(numout,*) '             Input file that contains river numbers = ',file_riv_number 
     84      IF(lwp) WRITE(numout,*) '             Maximum number of rivers to couple = ',nn_cpl_river 
     85      IF(lwp) WRITE(numout,*) '             Print river information = ',ln_print_river_info 
    8386      IF(lwp) WRITE(numout,*) ' ' 
    8487       
     
    97100        DO jj = 1, jpj 
    98101          rivers%river_number(ii,jj) = INT(river_number(ii,jj)) 
    99            
    100102          IF ( rivers%river_number(ii,jj) > max_river ) THEN 
    101103            max_river = rivers%river_number(ii,jj) 
    102           END IF 
    103            
    104           IF ( ii == 59 .AND. jj == 29 ) THEN 
    105             WRITE(numout,*) 'Amazon grid point river number (float) = ',river_number(ii,jj) 
    106             WRITE(numout,*) 'Amazon grid point river number (int) = ',INT(river_number(ii,jj)) 
    107             WRITE(numout,*) 'Amazon grid point river number (store) = ',rivers%river_number(ii,jj) 
    108             WRITE(numout,*) 'max_river at this stage = ',max_river 
    109             WRITE(numout,*) 'Amazon grid box area = ',e1e2t(ii,jj) 
    110104          END IF 
    111105        END DO 
     
    113107       
    114108      ! Print out the largest river number 
    115       WRITE(numout,*) 'Maximum river number = ',max_river 
     109      WRITE(numout,*) 'Maximum river number in input file = ',max_river 
    116110       
    117111      ! Get the area of each river outflow 
     
    128122      ! Use mpp_sum to add together river areas on other processors 
    129123      CALL mpp_sum( rivers%river_area, nn_cpl_river ) 
    130       WRITE(numout,*) 'Area of river number 1 is ',rivers%river_area(1) 
    131       WRITE(numout,*) 'Area of river number 10 is ',rivers%river_area(10) 
    132        
     124      IF ( ln_print_river_info ) THEN 
     125        WRITE(numout,*) 'Area of rivers 1 to 10 are ',rivers%river_area(1:10) 
     126      END IF 
    133127       
    134128   END SUBROUTINE cpl_rnf_1d_init 
     
    153147      LOGICAL  ::   llaction               ! Has the get worked? 
    154148       
    155       IF ( ln_ctl ) THEN 
     149      IF ( ln_print_river_info ) THEN 
    156150         WRITE(numout,*)' Getting data from 1D river runoff coupling ' 
    157151      ENDIF 
     
    163157                   
    164158      ! Output coupling info 
    165       IF ( ln_ctl ) THEN 
     159      IF ( ln_print_river_info ) THEN 
    166160         WRITE(numout,*)' narea = ', narea 
    167161         WRITE(numout,*)' kstep = ', kstep 
     
    193187      END IF 
    194188       
    195       IF ( ln_ctl ) WRITE(numout,*)' River runoff flux of AMAZON (pe 351) is ', rnf(59,29)          
     189      IF ( ln_print_river_info ) WRITE(numout,*)' River runoff flux of AMAZON (pe 351) is ', rnf(59,29)          
    196190    
    197191   END SUBROUTINE cpl_rnf_1d_rcv 
Note: See TracChangeset for help on using the changeset viewer.