Changes between Version 6 and Version 7 of 2009WP/2009Stream3/LeapFrogC

Timestamp:

2009-06-14T13:55:09+02:00 (15 years ago)

Author:

gm

Comment:

--

2009WP/2009Stream3/LeapFrogC

@@ -264 +264 @@
  applied the solar heat flux (qsr) as the mean of 2 consecutive time-steps (kt+1/2 and kt-1/2)
 
-'''to be explained'''
-
-{{{
-}}}
- becomes
-
-{{{
-}}}
-
-
-
-
-
-
-
+{{{
+            IF( nn_chldta ==1 ) THEN                             !*  Variable Chlorophyll
+...
+               zsi0r = 1.e0 / rn_si0
+               zcoef  = ( 1. - rn_abs ) / 3.e0                        ! equi-partition in R-G-B
+               ze0(:,:,1) = rn_abs  * qsr(:,:)
+               ze1(:,:,1) = zcoef * qsr(:,:)
+               ze2(:,:,1) = zcoef * qsr(:,:)
+               ze3(:,:,1) = zcoef * qsr(:,:)
+               zea(:,:,1) =         qsr(:,:)
+...
+                  ta(:,:,jk) = ta(:,:,jk) + etot3(:,:,jk) * qsr(:,:)
+
+         ENDIF
+         !                                                ! ------------------------- !
+         IF( ln_qsr_2bd ) THEN                            !  2 band light penetration !
+            !                                             ! ------------------------- !
+            !
+            DO jk = 1, nksr
+               DO jj = 2, jpjm1
+                  DO ji = fs_2, fs_jpim1   ! vector opt.
+                     ta(ji,jj,jk) = ta(ji,jj,jk) + etot3(ji,jj,jk) * qsr(ji,jj)
+                  END DO
+               END DO
+            END DO
+            !
+         ENDIF
+}}}
+ becomes
+
+{{{
+            IF( nn_chldta ==1 ) THEN                             !*  Variable Chlorophyll
+...
+               zsi0r = 1.e0 / rn_si0
+               zcoef  = ( 1. - rn_abs ) / 3.e0                        ! equi-partition in R-G-B
+               zea(:,:,1) =    0.5 * ( qsr_b(:,:) + qsr(:,:) )
+               ze0(:,:,1) = rn_abs * zea(:,:)
+               ze1(:,:,1) = zcoef  * zea(:,:)
+               ze2(:,:,1) = zcoef  * zea(:,:)
+               ze3(:,:,1) = zcoef  * zea(:,:)
+...
+            ELSE                                                 !*  Constant Chlorophyll
+               DO jk = 1, nksr
+                  ta(:,:,jk) = ta(:,:,jk) + etot3(:,:,jk) * 0.5 * ( qsr_b(:,:) + qsr(:,:) )
+               END DO
+            ENDIF
+
+         ENDIF
+         !                                                ! ------------------------- !
+         IF( ln_qsr_2bd ) THEN                            !  2 band light penetration !
+            !                                             ! ------------------------- !
+            !
+            DO jk = 1, nksr
+               ta(:,:,jk) = ta(:,:,jk) + etot3(:,:,jk) * 0.5 * ( qsr_b(:,:) + qsr(:,:) )
+            END DO
+            !
+         ENDIF
+}}}
+
+ '''CAUTION 1''' : the case where the light penetration is provided by the bio-model ( i.e. lk_qsr_bio = ln_qsr_bio =.true. is not yet treated. It probably requires changes in p4zopt.F90 and trcopt.F90. Not completely obvious since this depends on the TOP time-stepping strategy.