Changeset 2228 for branches/DEV_r2106_LOCEAN2010/NEMO/OPA_SRC/SBC/sbcdcy.F90

Timestamp:

2010-10-12T16:33:06+02:00 (14 years ago)

Author:

smasson

Message:

update diurnal cycle from dev_r2174_DCY into DEV_r2106_LOCEAN2010

File:

• branches/DEV_r2106_LOCEAN2010/NEMO/OPA_SRC/SBC/sbcdcy.F90 (modified) (11 diffs)

branches/DEV_r2106_LOCEAN2010/NEMO/OPA_SRC/SBC/sbcdcy.F90

@@ -16 +16 @@
    USE phycst           ! ocean physics
    USE dom_oce          ! ocean space and time domain
+   USE sbc_oce          ! Surface boundary condition: ocean fields
    USE in_out_manager   ! I/O manager
 
    IMPLICIT NONE
    PRIVATE
-   INTEGER                      ::   nday_qsr                    ! day when parameters were computed
+   INTEGER, PUBLIC              ::   nday_qsr                    ! day when parameters were computed
    REAL(wp), DIMENSION(jpi,jpj) ::   raa , rbb  , rcc  , rab     ! parameters used to compute the diurnal cycle
    REAL(wp), DIMENSION(jpi,jpj) ::   rtmd, rdawn, rdusk, rscal   !     -       -         -           -      -
-   REAL(wp), DIMENSION(jpi,jpj) ::   qsr_daily                   ! to hold daily mean QSR
   
    PUBLIC   sbc_dcy     ! routine called by sbc
@@ -34 +34 @@
 CONTAINS
 
-      SUBROUTINE sbc_dcy( kt, pqsr )
+      FUNCTION sbc_dcy( pqsrin ) RESULT( zqsrout )
       !!----------------------------------------------------------------------
       !!                  ***  ROUTINE sbc_dcy  ***
@@ -48 +48 @@
       !!              Part 1: a diurnally forced OGCM. Climate Dynamics 29:6, 575-590.
       !!----------------------------------------------------------------------
-      INTEGER,                        INTENT(in   ) ::   kt     ! ocean time-step index
-      REAL(wp), DIMENSION(jpi,jpj,1), INTENT(inout) ::   pqsr   ! QSR flux with diurnal cycle
-      !!
-      INTEGER  ::   ji, jj                                      ! dummy loop indices
+      REAL(wp), DIMENSION(jpi,jpj), INTENT(in) ::   pqsrin    ! input daily QSR flux 
+      !!
+      INTEGER  ::   ji, jj                                       ! dummy loop indices
       REAL(wp) ::   ztwopi, zinvtwopi, zconvrad 
       REAL(wp) ::   zlo, zup, zlousd, zupusd
-      REAL(wp) ::   zdsws, zdecrad, ztx
+      REAL(wp) ::   zdsws, zdecrad, ztx, zsin, zcos
       REAL(wp) ::   ztmp, ztmp1, ztmp2, ztest
+      REAL(wp), DIMENSION(jpi,jpj) ::   zqsrout                  ! output QSR flux with diurnal cycle
       !---------------------------statement functions------------------------
-      REAL(wp) ::   fintegral, pt1, pt2, paaa, pbbb, pccc       ! dummy statement function arguments
+      REAL(wp) ::   fintegral, pt1, pt2, paaa, pbbb, pccc        ! dummy statement function arguments
       fintegral( pt1, pt2, paaa, pbbb, pccc ) =                         &
          &   paaa * pt2 + zinvtwopi * pbbb * SIN(pccc + ztwopi * pt2)   &
@@ -70 +70 @@
 
       ! When are we during the day (from 0 to 1)
-      zlo = MOD( rdt / rday * REAL( kt-nit000, wp ), 1.)
-      zup = zlo + rdt / rday
+      zlo = ( REAL(nsec_day, wp) - 0.5   * rdttra(1) ) / rday
+      zup = zlo + ( REAL(nn_fsbc, wp) * rdttra(1) ) / rday
 
       !                                          
-      IF( kt == nit000 ) THEN       ! first time step only               
+      IF( nday_qsr == -1 ) THEN       ! first time step only               
          IF(lwp) THEN
             WRITE(numout,*)
@@ -81 +81 @@
             WRITE(numout,*)
          ENDIF
-         nday_qsr = 0
          ! Compute rcc needed to compute the time integral of the diurnal cycle
          rcc(:,:) = zconvrad * glamt(:,:) - rpi
@@ -99 +98 @@
          nday_qsr = nday 
          ! number of days since the previous winter solstice (supposed to be always 21 December)         
-         zdsws = 11 + nday_year
+         zdsws = REAL(11 + nday_year, wp)
          ! declination of the earths orbit
          zdecrad = (-23.5 * zconvrad) * COS( zdsws * ztwopi / REAL(nyear_len(1),wp) )
-         ! save the daily QSR for nest hours of the day
-         qsr_daily(:,:) = pqsr(:,:,1)
          ! Compute A and B needed to compute the time integral of the diurnal cycle
         
+         zsin = SIN( zdecrad )   ;   zcos = COS( zdecrad )
          DO jj = 1, jpj
             DO ji = 1, jpi
                ztmp = zconvrad * gphit(ji,jj)
-               raa(ji,jj) = SIN( ztmp ) * SIN( zdecrad )
-               rbb(ji,jj) = COS( ztmp ) * COS( zdecrad )
+               raa(ji,jj) = SIN( ztmp ) * zsin
+               rbb(ji,jj) = COS( ztmp ) * zcos
             END DO  
          END DO  
@@ -140 +138 @@
          rdawn(:,:) = MOD((rdawn(:,:) + 1.), 1.)
          rdusk(:,:) = MOD((rdusk(:,:) + 1.), 1.)
-
 
          !     2.2 Compute the scalling function:
@@ -166 +163 @@
          END DO  
          !
-         ztmp = rday / rdt
+         ztmp = rday / ( rdttra(1) * REAL(nn_fsbc, wp) )
          rscal(:,:) = rscal(:,:) * ztmp
 
@@ -184 +181 @@
                   zupusd = MAX(zupusd, zlo)
                   ztmp = fintegral(zlousd, zupusd, raa(ji,jj), rbb(ji,jj), rcc(ji,jj)) 
-                  pqsr(ji,jj,1) = qsr_daily(ji,jj) * ztmp * rscal(ji,jj)
+                  zqsrout(ji,jj) = pqsrin(ji,jj) * ztmp * rscal(ji,jj)
                   !
                ELSE                                         ! day time in two parts
@@ -194 +191 @@
                   ztmp2 = fintegral(zlousd, zupusd, raa(ji,jj), rbb(ji,jj), rcc(ji,jj)) 
                   ztmp = ztmp1 + ztmp2
-                  pqsr(ji,jj,1) = qsr_daily(ji,jj) * ztmp * rscal(ji,jj)
+                  zqsrout(ji,jj) = pqsrin(ji,jj) * ztmp * rscal(ji,jj)
                ENDIF
             ELSE                                   ! 24h light or 24h night
                !
-               IF( raa(ji,jj) > rbb(ji,jj) ) THEN         ! 24h day
+               IF( raa(ji,jj) > rbb(ji,jj) ) THEN           ! 24h day
                   ztmp = fintegral(zlo, zup, raa(ji,jj), rbb(ji,jj), rcc(ji,jj)) 
-                  pqsr(ji,jj,1) = qsr_daily(ji,jj) * ztmp * rscal(ji,jj)
+                  zqsrout(ji,jj) = pqsrin(ji,jj) * ztmp * rscal(ji,jj)
                   !
                ELSE                                         ! No day
-                  pqsr(ji,jj,1) = 0.e0
+                  zqsrout(ji,jj) = 0.e0
                ENDIF
             ENDIF
@@ -209 +206 @@
       END DO  
       !
-   END SUBROUTINE sbc_dcy
+   END FUNCTION sbc_dcy
 
    !!======================================================================