Changeset 3858 for branches/2013/dev_r3856_MERCATOR3_QSRMEAN24H/NEMOGCM/NEMO/OPA_SRC/SBC/sbcblk_core.F90

Timestamp:

2013-04-03T11:59:02+02:00 (11 years ago)

Author:

cbricaud

Message:

commit first draft of MERCATOR3_QSRMEAN24H

File:

• branches/2013/dev_r3856_MERCATOR3_QSRMEAN24H/NEMOGCM/NEMO/OPA_SRC/SBC/sbcblk_core.F90 (modified) (2 diffs)

branches/2013/dev_r3856_MERCATOR3_QSRMEAN24H/NEMOGCM/NEMO/OPA_SRC/SBC/sbcblk_core.F90

@@ -186 +186 @@
       CALL fld_read( kt, nn_fsbc, sf )             ! input fields provided at the current time-step
 
+      ! If diurnal cycle is activated, compute a daily mean short waves flux for biogeochemistery 
+#if defined key_top
+      IF( ltrcdm2dc )CALL blk_bio_meanqsr
+#endif
+
       !                                            ! compute the surface ocean fluxes using CORE bulk formulea
       IF( MOD( kt - 1, nn_fsbc ) == 0 )   CALL blk_oce_core( kt, sf, sst_m, ssu_m, ssv_m )
@@ -423 +428 @@
       !
    END SUBROUTINE blk_oce_core
-   
+  
+   SUBROUTINE blk_bio_meanqsr
+      !!---------------------------------------------------------------------
+      !!                     ***  ROUTINE blk_bio_meanqsr
+      !!                     
+      !! ** Purpose :   provide daily qsr_mean for PISCES when
+      !!                analytic diurnal cycle is applied in physic
+      !!                
+      !! ** Method  :   add part where there is no ice
+      !! 
+      !!---------------------------------------------------------------------
+      INTEGER  :: ji,jj
+      REAL(wp) :: zfrld
+      !!---------------------------------------------------------------------
+      IF( nn_timing == 1 )  CALL timing_start('blk_bio_meanqsr')
+
+      qsr_mean(:,:) = (1. - albo ) *  sf(jp_qsr)%fnow(:,:,1)
+
+      IF( nn_timing == 1 )  CALL timing_stop('blk_bio_meanqsr')
+
+   END SUBROUTINE blk_bio_meanqsr
+ 
    
    SUBROUTINE blk_ice_core(  pst   , pui   , pvi   , palb ,   &