Changeset 5236 for branches/2015/dev_r5204_CNRS_PISCES_dcy/NEMOGCM/NEMO/OPA_SRC/SBC/sbccpl.F90

Timestamp:

2015-04-24T14:08:11+02:00 (9 years ago)

Author:

cetlod

Message:

NEMOGCM_dev_r5204_CNRS_PISCES_dcy : update routines according to the new strategy, see ticket #1484

File:

• branches/2015/dev_r5204_CNRS_PISCES_dcy/NEMOGCM/NEMO/OPA_SRC/SBC/sbccpl.F90 (modified) (4 diffs)

branches/2015/dev_r5204_CNRS_PISCES_dcy/NEMOGCM/NEMO/OPA_SRC/SBC/sbccpl.F90

@@ -135 +135 @@
 
    REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:)   ::   albedo_oce_mix     ! ocean albedo sent to atmosphere (mix clear/overcast sky)
-
-   REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:)   :: qsr_tot_tmp     ! arrays containing consecutives qsr in a day
-   REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:,:) :: qsr_ice_tmp     !  ===                        === 
 
    INTEGER , ALLOCATABLE, SAVE, DIMENSION(    :) ::   nrcvinfo           ! OASIS info argument
@@ -586 +583 @@
       ENDIF
       !
-      n_cpl_qsr = INT( 86400 / ( cpl_freq( 'O_QsrOce' ) + cpl_freq( 'O_QsrMix' ) + cpl_freq( 'S_QsrOce' ) + cpl_freq( 'S_QsrMix' ) ) )
-      !
-      IF( ln_dm2dc .AND. n_cpl_qsr /= 1 )   &
+      ncpl_qsr_freq = INT( 86400 / ( cpl_freq( 'O_QsrOce' ) + cpl_freq( 'O_QsrMix' ) + cpl_freq( 'S_QsrOce' ) + cpl_freq( 'S_QsrMix' ) ) )
+      !
+      IF( ln_dm2dc .AND. ncpl_qsr_freq /= 1 )   &
          &   CALL ctl_stop( 'sbc_cpl_init: diurnal cycle reconstruction (ln_dm2dc) needs daily couping for solar radiation' )
 
@@ -1290 +1287 @@
       ENDIF
       !
-      IF( l_trcdm2dc )  CALL sbc_cpl_qsr_mean( it )    ! computation of daily mean qsr for biogeochemical model if needed
-
-      !                                                      ! ========================= !
       SELECT CASE( TRIM( sn_rcv_dqnsdt%cldes ) )             !          d(qns)/dt        !
       !                                                      ! ========================= !
@@ -1327 +1321 @@
    END SUBROUTINE sbc_cpl_ice_flx
 
-   SUBROUTINE sbc_cpl_qsr_mean( kt )
-      !!----------------------------------------------------------------------
-      !!             ***  ROUTINE sbc_cpl_mean  ***
-      !!
-      !! ** Purpose :   Compute daily mean qsr for biogeochmeical model in case of diurnal cycle
-      !!
-      !!----------------------------------------------------------------------
-      INTEGER, INTENT(in) ::   kt
-      INTEGER  :: jn
-
-      IF( kt == nit000 ) THEN
-         ALLOCATE( qsr_tot_tmp(jpi,jpj,n_cpl_qsr), qsr_ice_tmp(jpi,jpj,jpl,n_cpl_qsr) ) 
-         DO jn = 1, n_cpl_qsr
-            qsr_tot_tmp(:,:  ,jn) = qsr_tot(:,:  )
-            qsr_ice_tmp(:,:,:,jn) = qsr_ice(:,:,:)
-         ENDDO
-         qsr_mean    (:,:  ) = qsr_tot(:,:  ) 
-         qsr_ice_mean(:,:,:) = qsr_ice(:,:,:)
-      ENDIF
-      !
-      IF( kt /= nit000 .AND.  nrcvinfo(jpr_qsroce) == OASIS_Rcv ) THEN !  => need to be done only when we receive the field
-         DO jn = 1, n_cpl_qsr - 1
-            qsr_tot_tmp(:,:  ,jn) = qsr_tot_tmp(:,:  ,jn+1)
-            qsr_ice_tmp(:,:,:,jn) = qsr_ice_tmp(:,:,:,jn+1)
-         ENDDO
-         qsr_tot_tmp(:,:  ,n_cpl_qsr ) = qsr_tot(:,:  )
-         qsr_ice_tmp(:,:,:,n_cpl_qsr ) = qsr_ice(:,:,:)
-         !
-         qsr_mean    (:,:  ) = SUM( qsr_tot_tmp(:,:,:  ), 3 ) / n_cpl_qsr 
-         qsr_ice_mean(:,:,:) = SUM( qsr_ice_tmp(:,:,:,:), 4 ) / n_cpl_qsr 
-         !
-      ENDIF
-      !
-   END SUBROUTINE sbc_cpl_qsr_mean
-   
-   
+
    SUBROUTINE sbc_cpl_snd( kt )
       !!----------------------------------------------------------------------