Changeset 1274 for trunk/NEMO/OPA_SRC/SBC/sbcflx.F90

Timestamp:

2009-01-19T16:32:54+01:00 (15 years ago)

Author:

smasson

Message:

bugfix in sbcflx, see ticket:294

File:

• trunk/NEMO/OPA_SRC/SBC/sbcflx.F90 (modified) (1 diff)

trunk/NEMO/OPA_SRC/SBC/sbcflx.F90

@@ -143 +143 @@
       !                                          ! input fields at the current time-step
 
-      ! set the ocean fluxes from read fields
+      IF( MOD( kt-1, nn_fsbc ) == 0 ) THEN
+         !
+         ! set the ocean fluxes from read fields
 !CDIR COLLAPSE
-      DO jj = 1, jpj
-         DO ji = 1, jpi
-            utau(ji,jj) = sf(jp_utau)%fnow(ji,jj)
-            vtau(ji,jj) = sf(jp_vtau)%fnow(ji,jj)
-            qns (ji,jj) = sf(jp_qtot)%fnow(ji,jj) - sf(jp_qsr)%fnow(ji,jj)
-            qsr (ji,jj) = sf(jp_qsr )%fnow(ji,jj)
-            emp (ji,jj) = sf(jp_emp )%fnow(ji,jj)
+         DO jj = 1, jpj
+            DO ji = 1, jpi
+               utau(ji,jj) = sf(jp_utau)%fnow(ji,jj)
+               vtau(ji,jj) = sf(jp_vtau)%fnow(ji,jj)
+               qns (ji,jj) = sf(jp_qtot)%fnow(ji,jj) - sf(jp_qsr)%fnow(ji,jj)
+               qsr (ji,jj) = sf(jp_qsr )%fnow(ji,jj)
+               emp (ji,jj) = sf(jp_emp )%fnow(ji,jj)
+            END DO
          END DO
-      END DO
-
-      ! Initialization of emps (when no ice model)
-      emps(:,:) = emp (:,:) 
-
-      ! Estimation of wind speed as a function of wind stress ( |tau|=rhoa*Cd*|U|^2 )
-      CALL sbc_tau2wnd
-
-      ! control print (if less than 100 time-step asked)
-      IF( nitend-nit000 <= 100 .AND. lwp ) THEN
-         WRITE(numout,*) 
-         WRITE(numout,*) '        read daily momentum, heat and freshwater fluxes OK'
-         DO jf = 1, jpfld
-            IF( jf == jp_utau .OR. jf == jp_vtau )   zfact =     1.
-            IF( jf == jp_qtot .OR. jf == jp_qsr  )   zfact =     0.1
-            IF( jf == jp_emp                     )   zfact = 86400.
+         
+         ! Initialization of emps (when no ice model)
+         emps(:,:) = emp (:,:) 
+         
+         ! Estimation of wind speed as a function of wind stress ( |tau|=rhoa*Cd*|U|^2 )
+         CALL sbc_tau2wnd
+         
+         ! control print (if less than 100 time-step asked)
+         IF( nitend-nit000 <= 100 .AND. lwp ) THEN
             WRITE(numout,*) 
-            WRITE(numout,*) ' day: ', ndastp , TRIM(sf(jf)%clvar), ' * ', zfact
-            CALL prihre( sf(jf)%fnow, jpi, jpj, 1, jpi, 20, 1, jpj, 10, zfact, numout )
-         END DO
-         CALL FLUSH(numout)
+            WRITE(numout,*) '        read daily momentum, heat and freshwater fluxes OK'
+            DO jf = 1, jpfld
+               IF( jf == jp_utau .OR. jf == jp_vtau )   zfact =     1.
+               IF( jf == jp_qtot .OR. jf == jp_qsr  )   zfact =     0.1
+               IF( jf == jp_emp                     )   zfact = 86400.
+               WRITE(numout,*) 
+               WRITE(numout,*) ' day: ', ndastp , TRIM(sf(jf)%clvar), ' * ', zfact
+               CALL prihre( sf(jf)%fnow, jpi, jpj, 1, jpi, 20, 1, jpj, 10, zfact, numout )
+            END DO
+            CALL FLUSH(numout)
+         ENDIF
+         !
       ENDIF
       !