MODULE cpl_rnf_1d !!====================================================================== !! *** MODULE cpl_rnf_1d *** !! Ocean forcing: River runoff passed from the atmosphere using !! 1D array. One value per river. !!===================================================================== !! History : ?.? ! 2018-01 (D. Copsey) Initial setup !!---------------------------------------------------------------------- !!---------------------------------------------------------------------- !! cpl_rnf_1d_init : runoffs initialisation !!---------------------------------------------------------------------- #if defined key_oasis3 USE mod_oasis ! OASIS3-MCT module #endif USE timing ! Timing USE in_out_manager ! I/O units USE lib_mpp ! MPP library USE iom USE wrk_nemo ! Memory allocation USE dom_oce ! Domain sizes (for grid box area e1e2t) USE sbc_oce ! Surface boundary condition: ocean fields USE cpl_oasis3 ! Coupling information (for n_rivers and runoff_id) IMPLICIT NONE PRIVATE PUBLIC cpl_rnf_1d_init ! routine called in nemo_init PUBLIC cpl_rnf_1d_rcv ! routine called in sbccpl.F90 TYPE, PUBLIC :: RIVERS_DATA !: Storage for river outflow data INTEGER, POINTER, DIMENSION(:,:) :: river_number !: River outflow number REAL(wp), POINTER, DIMENSION(:) :: river_area ! 1D array listing areas of each river outflow (m2) END TYPE RIVERS_DATA TYPE(RIVERS_DATA), PUBLIC, TARGET :: rivers !: River data CONTAINS SUBROUTINE cpl_rnf_1d_init !!---------------------------------------------------------------------- !! *** SUBROUTINE cpl_rnf_1d_init *** !! !! ** Purpose : - Read in file for river outflow numbers. !! Calculate 2D area of river outflow points. !! Called from nemo_init (nemogcm.F90). !! !!---------------------------------------------------------------------- !! namelist variables !!------------------- CHARACTER(len=80) :: file_riv_number !: Filename for river numbers INTEGER :: ios ! Local integer output status for namelist read INTEGER :: inum INTEGER :: ii, jj !: Loop indices INTEGER :: max_river REAL(wp), POINTER, DIMENSION(:,:) :: river_number ! 2D array containing the river outflow numbers NAMELIST/nam_cpl_rnf_1d/file_riv_number !!---------------------------------------------------------------------- IF( nn_timing == 1 ) CALL timing_start('cpl_rnf_1d_init') IF(lwp) WRITE(numout,*) IF(lwp) WRITE(numout,*) 'cpl_rnf_1d_init : initialization of river runoff coupling' IF(lwp) WRITE(numout,*) '~~~~~~~~~~~~' REWIND(numnam_cfg) ! Read the namelist READ ( numnam_ref, nam_cpl_rnf_1d, IOSTAT = ios, ERR = 901) 901 IF( ios /= 0 ) CALL ctl_nam ( ios , 'nam_cpl_rnf_1d in reference namelist', lwp ) READ ( numnam_cfg, nam_cpl_rnf_1d, IOSTAT = ios, ERR = 902 ) 902 IF( ios /= 0 ) CALL ctl_nam ( ios , 'nam_cpl_rnf_1d in configuration namelist', lwp ) IF(lwm) WRITE ( numond, nam_cpl_rnf_1d ) ! ! Parameter control and print IF(lwp) WRITE(numout,*) ' ' IF(lwp) WRITE(numout,*) ' Namelist nam_cpl_rnf_1d : Coupled runoff using 1D array' IF(lwp) WRITE(numout,*) ' Input file that contains river numbers = ',file_riv_number IF(lwp) WRITE(numout,*) ' ' ! Assign space for river numbers ALLOCATE( rivers%river_number( jpi, jpj ) ) CALL wrk_alloc( jpi, jpj, river_number ) ! Read the river numbers from netcdf file CALL iom_open (file_riv_number , inum ) CALL iom_get ( inum, jpdom_data, 'river_number', river_number ) CALL iom_close( inum ) ! Convert from a real array to an integer array max_river=0 DO ii = 1, jpi DO jj = 1, jpj rivers%river_number(ii,jj) = INT(river_number(ii,jj)) IF ( rivers%river_number(ii,jj) > max_river ) THEN max_river = rivers%river_number(ii,jj) END IF IF ( ii == 59 .AND. jj == 29 ) THEN WRITE(numout,*) 'Amazon grid point river number (float) = ',river_number(ii,jj) WRITE(numout,*) 'Amazon grid point river number (int) = ',INT(river_number(ii,jj)) WRITE(numout,*) 'Amazon grid point river number (store) = ',rivers%river_number(ii,jj) WRITE(numout,*) 'max_river at this stage = ',max_river WRITE(numout,*) 'Amazon grid box area = ',e1e2t(ii,jj) END IF END DO END DO ! Print out the largest river number WRITE(numout,*) 'Maximum river number = ',max_river ! Get the area of each river outflow ALLOCATE( rivers%river_area( n_rivers ) ) rivers%river_area(:) = 0.0 DO ii = 1, jpi DO jj = 1, jpj IF ( rivers%river_number(ii,jj) > 0 .AND. rivers%river_number(ii,jj) <= n_rivers ) THEN rivers%river_area(rivers%river_number(ii,jj)) = rivers%river_area(rivers%river_number(ii,jj)) + e1e2t(ii,jj) END IF END DO END DO ! Use mpp_sum to add together river areas on other processors CALL mpp_sum( rivers%river_area, n_rivers ) WRITE(numout,*) 'Area of river number 1 is ',rivers%river_area(1) WRITE(numout,*) 'Area of river number 10 is ',rivers%river_area(10) END SUBROUTINE cpl_rnf_1d_init SUBROUTINE cpl_rnf_1d_rcv( kstep) !!---------------------------------------------------------------------- !! *** SUBROUTINE cpl_rnf_1d_rcv *** !! !! ** Purpose : - Get river outflow from 1D array (passed from the !! atmosphere) and transfer it to the 2D NEMO runoff !! field. !! Called from sbc_cpl_rcv (sbccpl.F90). !! !!---------------------------------------------------------------------- INTEGER , INTENT(in ) :: kstep ! ocean time-step in seconds INTEGER :: kinfo ! OASIS3 info argument REAL(wp) :: runoff_1d(n_rivers) ! River runoff. One value per river. INTEGER :: ii, jj ! Loop indices LOGICAL :: llaction ! Has the get worked? IF ( ln_ctl ) THEN WRITE(numout,*)' Getting data from 1D river runoff coupling ' ENDIF ! Get the river runoff sent by the atmosphere CALL oasis_get ( runoff_id, kstep, runoff_1d, kinfo ) llaction = kinfo == OASIS_Recvd .OR. kinfo == OASIS_FromRest .OR. & & kinfo == OASIS_RecvOut .OR. kinfo == OASIS_FromRestOut ! Output coupling info IF ( ln_ctl ) THEN WRITE(numout,*)' narea = ', narea WRITE(numout,*)' kstep = ', kstep WRITE(numout,*)' River runoff of river number 1 is ', runoff_1d(1) WRITE(numout,*)' River runoff of river number 10 is ', runoff_1d(10) WRITE(numout,*)' kinfo = ', kinfo WRITE(numout,*)' llaction = ', llaction WRITE(numout,*)' OASIS_Recvd = ',OASIS_Recvd WRITE(numout,*)' OASIS_FromRest = ',OASIS_FromRest WRITE(numout,*)' OASIS_RecvOut = ',OASIS_RecvOut WRITE(numout,*)' OASIS_FromRestOut = ',OASIS_FromRestOut WRITE(numout,*)'-------' ENDIF IF ( llaction ) THEN ! Convert the 1D total runoff per river to 2D runoff flux by ! dividing by the area of each runoff zone. DO ii = 1, jpi DO jj = 1, jpj IF ( rivers%river_number(ii,jj) > 0 .AND. rivers%river_number(ii,jj) <= n_rivers ) THEN rnf(ii,jj) = runoff_1d(rivers%river_number(ii,jj)) / rivers%river_area(rivers%river_number(ii,jj)) ELSE rnf(ii,jj) = 0.0 END IF END DO END DO END IF IF ( ln_ctl ) WRITE(numout,*)' River runoff flux of AMAZON (pe 351) is ', rnf(59,29) END SUBROUTINE cpl_rnf_1d_rcv END MODULE cpl_rnf_1d