Ignore:
Timestamp:
2015-12-07T14:32:39+01:00 (5 years ago)
Author:
timgraham
Message:

Tidying of DIU code

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2015/dev_MetOffice_merge_2015/NEMOGCM/NEMO/OPA_SRC/DIU/diurnal_bulk.F90

    r5676 r6010  
    8181   END SUBROUTINE diurnal_sst_bulk_init 
    8282    
    83    SUBROUTINE diurnal_sst_takaya_step(psolflux, pqflux, ptauflux, prho, p_rdt,& 
    84             &                  pla, pthick, pcoolthick, pmu, ld_calcfrac, & 
     83   SUBROUTINE diurnal_sst_takaya_step(kt, psolflux, pqflux, ptauflux, prho, p_rdt,& 
     84            &                  pla, pthick, pcoolthick, pmu, & 
    8585            &                  p_fvel_bkginc, p_hflux_bkginc) 
    8686      !!---------------------------------------------------------------------- 
     
    123123      REAL(wp), INTENT(IN) :: p_rdt                         ! time-step 
    124124       
    125       LOGICAL, OPTIONAL, & 
    126       &     INTENT(IN)  :: ld_calcfrac   ! Set TRUE to recalculate the  
    127                                          ! solar fraction 
    128        
    129125      ! Local variables  
    130126      REAL(wp), DIMENSION(jpi,jpj) :: z_fvel              ! friction velocity      
     
    136132       
    137133      INTEGER :: ji,jj 
     134      INTEGER, INTENT(IN) :: kt                           ! time step 
    138135 
    139136      ! Set optional arguments to their defaults 
     
    158155         zla(:,:) = pla(:,:) 
    159156      ENDIF 
    160       IF ( .NOT. PRESENT(ld_calcfrac) ) THEN 
    161          ll_calcfrac = .FALSE. 
    162       ELSE 
    163          ll_calcfrac = ld_calcfrac 
    164       ENDIF      
    165157       
    166158      ! If not done already, calculate the solar fraction 
    167       IF (ll_calcfrac ) THEN 
     159      IF ( kt==nit000 ) THEN 
    168160         DO jj = 1,jpj 
    169161            DO ji = 1, jpi 
Note: See TracChangeset for help on using the changeset viewer.