- Timestamp:
- 2017-12-19T09:26:25+01:00 (6 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2017/dev_merge_2017/NEMOGCM/NEMO/OPA_SRC/SBC/sbcdcy.F90
r6140 r9124 16 16 USE dom_oce ! ocean space and time domain 17 17 USE sbc_oce ! Surface boundary condition: ocean fields 18 ! 18 19 USE in_out_manager ! I/O manager 19 20 USE lib_mpp ! MPP library 20 USE timing ! Timing21 21 22 22 IMPLICIT NONE … … 63 63 !! Part 1: a diurnally forced OGCM. Climate Dynamics 29:6, 575-590. 64 64 !!---------------------------------------------------------------------- 65 LOGICAL , OPTIONAL, INTENT(in) :: l_mask! use the routine for night mask computation65 LOGICAL , OPTIONAL , INTENT(in) :: l_mask ! use the routine for night mask computation 66 66 REAL(wp), DIMENSION(jpi,jpj), INTENT(in) :: pqsrin ! input daily QSR flux 67 REAL(wp), DIMENSION(jpi,jpj) :: zqsrout ! output QSR flux with diurnal cycle 67 68 !! 68 69 INTEGER :: ji, jj ! dummy loop indices … … 73 74 REAL(wp) :: ztmp, ztmp1, ztmp2, ztest 74 75 REAL(wp) :: ztmpm, ztmpm1, ztmpm2 75 REAL(wp), DIMENSION(jpi,jpj) :: zqsrout ! output QSR flux with diurnal cycle76 76 !---------------------------statement functions------------------------ 77 77 REAL(wp) :: fintegral, pt1, pt2, paaa, pbbb, pccc ! dummy statement function arguments … … 80 80 & - paaa * pt1 - zinvtwopi * pbbb * SIN(pccc + ztwopi * pt1) 81 81 !!--------------------------------------------------------------------- 82 !83 IF( nn_timing == 1 ) CALL timing_start('sbc_dcy')84 82 ! 85 83 ! Initialization … … 199 197 DO jj = 1, jpj 200 198 DO ji = 1, jpi 201 ztmpm = 0. 0199 ztmpm = 0._wp 202 200 IF( ABS(rab(ji,jj)) < 1. ) THEN ! day duration is less than 24h 203 201 ! … … 241 239 END DO 242 240 ! 243 IF 241 IF( PRESENT(l_mask) .AND. l_mask ) THEN 244 242 zqsrout(:,:) = float(imask_night(:,:)) 245 243 ENDIF 246 244 ! 247 IF( nn_timing == 1 ) CALL timing_stop('sbc_dcy')248 !249 245 END FUNCTION sbc_dcy 250 246
Note: See TracChangeset
for help on using the changeset viewer.