Changeset 2187 for branches/dev_r2174_DCY/NEMO/OPA_SRC/SBC/sbcblk_core.F90

Timestamp:

2010-10-07T17:44:48+02:00 (14 years ago)

Author:

smasson

Message:

add diurna cycle in dev_r2174_DCY, see ticket:730

File:

• branches/dev_r2174_DCY/NEMO/OPA_SRC/SBC/sbcblk_core.F90 (modified) (3 diffs)

branches/dev_r2174_DCY/NEMO/OPA_SRC/SBC/sbcblk_core.F90

@@ -26 +26 @@
    USE fldread         ! read input fields
    USE sbc_oce         ! Surface boundary condition: ocean fields
+   USE sbcdcy          ! surface forcing: diurnal cycle
    USE iom             ! I/O manager library
    USE in_out_manager  ! I/O manager
@@ -147 +148 @@
          READ  ( numnam, namsbc_core )
          !
+         ! do we plan to use ln_dm2dc with non-daily forcing?
+         IF( ln_dm2dc .AND. sn_qsr%nfreqh /= 24 )   &
+            &   CALL ctl_stop( 'sbc_blk_core: ln_dm2dc can be activated only with daily short-wave forcing' ) 
+         IF( ln_dm2dc .AND. sn_qsr%ln_tint ) THEN
+            CALL ctl_warn( 'sbc_blk_core: ln_dm2dc is taking care of the temporal interpolation of daily qsr',   &
+                 &         '              ==> We force time interpolation = .false. for qsr' )
+            sn_qsr%ln_tint = .false.
+         ENDIF
+         !
          ! store namelist information in an array
          slf_i(jp_wndi) = sn_wndi   ;   slf_i(jp_wndj) = sn_wndj
@@ -169 +179 @@
          !
          ! fill sf with slf_i and control print
-         CALL fld_fill( sf, slf_i, cn_dir, 'sbc_blk_core', 'flux formulattion for ocean surface boundary condition', 'namsbc_core' )
+         CALL fld_fill( sf, slf_i, cn_dir, 'sbc_blk_core', 'flux formulation for ocean surface boundary condition', 'namsbc_core' )
          !
       ENDIF
 
       CALL fld_read( kt, nn_fsbc, sf )                   ! input fields provided at the current time-step
+
+      IF( ln_dm2dc )   CALL sbc_dcy( kt , sf(jp_qsr)%fnow )   ! modify sf(jp_qsr)%fnow for diurnal cycle
 
 #if defined key_lim3