Changeset 1462

Timestamp:

2009-06-08T17:40:59+02:00 (15 years ago)

Author:

rblod

Message:

Bugs in bdytides, see ticket #441

File:

• trunk/NEMO/OPA_SRC/BDY/bdytides.F90 (modified) (3 diffs)

trunk/NEMO/OPA_SRC/BDY/bdytides.F90

@@ -146 +146 @@
       CHARACTER(len=80) :: clfile         ! full file name for tidal input file 
       INTEGER ::   ipi, ipj, inum, idvar  ! temporary integers (netcdf read)
-      INTEGER, DIMENSION(3) :: lendta     ! length of data in the file (note may be different from nblendta!)
+      INTEGER, DIMENSION(3) :: lendta=0   ! length of data in the file (note may be different from nblendta!)
       REAL(wp) ::  z_arg, z_atde, z_btde, z1t, z2t           
       REAL(wp), DIMENSION(jpbdta,1) ::   zdta   ! temporary array for data fields
@@ -230 +230 @@
          ! This is the initialisation step, so nday, nmonth, nyear are the 
          ! initial date/time of the integration.
+           print *, nday,nmonth,nyear
+           nyear  = int(ndate0 / 10000  )                           ! initial year
+           nmonth = int((ndate0 - nyear * 10000 ) / 100 )          ! initial month
+           nday   = int(ndate0 - nyear * 10000 - nmonth * 100)
 
          CALL uvset( 0, nday, nmonth, nyear, z_ftc, z_vplu )
@@ -301 +305 @@
          z_arg = kt * rdt * rad / 3600.0
       ELSE                              ! we are in a barotropic subcycle (for timesplitting option)
-         z_arg = ( (kt-1) * rdt + jit * rdt / REAL(nn_baro,lwp) ) * rad / 3600.0
+!         z_arg = ( (kt-1) * rdt + jit * rdt / REAL(nn_baro,lwp) ) * rad / 3600.0
+         z_arg = ( (kt-1) * rdt + jit * rdt / REAL(nn_baro,wp) ) * rad / 3600.0
       ENDIF