Changeset 2187 for branches/dev_r2174_DCY/NEMO/OPA_SRC/SBC/sbcmod.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/sbcmod.F90 (modified) (3 diffs)

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

@@ -137 +137 @@
       IF( nn_ice == 2 .AND. .NOT.( ln_blk_clio .OR. ln_blk_core .OR. lk_cpl ) )   &
          &   CALL ctl_stop( 'sea-ice model requires a bulk formulation or coupled configuration' )
+      
+      IF( ln_dm2dc .AND. .NOT.( ln_flx .OR. ln_blk_core ) )   &
+         &   CALL ctl_stop( 'diurnal cycle into qsr field from daily values requires a flux or core-bulk formulation' )
+      
+      IF( ln_dm2dc .AND. ( ( NINT(rday) / ( nn_fsbc * NINT(rdt) ) )  < 8 ) )   &
+         &   CALL ctl_warn( 'diurnal cycle for qsr: the sampling of the diurnal cycle is too small...' )
       
       ! Choice of the Surface Boudary Condition (set nsbc)
@@ -216 +222 @@
       ! Misc. Options
       ! -------------
-
-!!gm  IF( ln_dm2dc       )   CALL sbc_dcy( kt )                 ! Daily mean qsr distributed over the Diurnal Cycle
       
       SELECT CASE( nn_ice )                                     ! Update heat and freshwater fluxes over sea-ice areas
@@ -239 +243 @@
          CALL iom_put( "emp"    , emp       )                   ! upward water flux
          CALL iom_put( "emps"   , emps      )                   ! c/d water flux
-         CALL iom_put( "qns+qsr", qns + qsr )                   ! total heat flux   (caution if ln_dm2dc=true, to be 
-         CALL iom_put( "qns"    , qns       )                   ! solar heat flux    moved after the call to iom_setkt)
-         CALL iom_put( "qsr"    ,       qsr )                   ! solar heat flux    moved after the call to iom_setkt)
+         CALL iom_put( "qns+qsr", qns + qsr )                   ! total heat flux
+         CALL iom_put( "qns"    , qns       )                   ! solar heat flux
+         CALL iom_put( "qsr"    ,       qsr )                   ! solar heat flux
          IF(  nn_ice > 0 )   CALL iom_put( "ice_cover", fr_i )  ! ice fraction 
       ENDIF