Changeset 194

Timestamp:

2004-12-22T16:43:16+01:00 (19 years ago)

Author:

opalod

Message:

CT : BUGFIX135 : correction of the zxday calculation using nday_year the number of days since the beginning of the run

Location:

trunk/NEMO/OPA_SRC/SBC

Files:

• bulk.F90 (modified) (1 diff)
• flxblk.F90 (modified) (4 diffs)

trunk/NEMO/OPA_SRC/SBC/bulk.F90

@@ -89 +89 @@
          zsst(:,:) = gsst(:,:) / REAL( nfbulk ) * tmask(:,:,1)
          zsss(:,:) = gsss(:,:) / REAL( nfbulk )
-         CALL flx_blk( kt, zsst )    
+         CALL flx_blk( zsst )    
      
          gsst(:,:) = 0.    

trunk/NEMO/OPA_SRC/SBC/flxblk.F90

@@ -77 +77 @@
 CONTAINS
 
-   SUBROUTINE flx_blk( kt, psst )
+   SUBROUTINE flx_blk( psst )
       !!---------------------------------------------------------------------------
       !!                     ***  ROUTINE flx_blk  ***
@@ -104 +104 @@
       !!----------------------------------------------------------------------
       !! * Arguments
-      INTEGER , INTENT( in  )  ::   kt        ! time step
       REAL(wp), INTENT( in ), DIMENSION(jpi,jpj) ::   &
          &                          psst      ! Sea Surface Temperature 
@@ -139 +138 @@
       REAL(wp)  ::            & 
          zxday             ,  &  ! day of year
-         ztsec             ,  &  ! time in seconds
          zdist             ,  &  ! distance between the sun and the earth during the year
          zdaycor           ,  &  ! corr. factor to take into account the variation of 
@@ -316 +314 @@
       !  compution of the day of the year at which the fluxes have to be calculate 
       !--The date corresponds to the middle of the time step.
-      ztsec = ( 2 * INT ( ( kt - 1 ) / nfbulk ) + 1 ) * rdtbs2
-      zxday = MOD( ztsec , raass ) / rday + 1.0  
+      zxday=nday_year + rdtbs2/rday
 
       iday   = INT( zxday )