Changeset 5222 for branches/2015/dev_r5204_CNRS_PISCES_dcy/NEMOGCM/NEMO/LIM_SRC_3/limsbc.F90

Timestamp:

2015-04-17T12:24:12+02:00 (6 years ago)

Author:

cetlod

Message:

dev_r5204_CNRS_PISCES_dcy: some improvments

File:

• branches/2015/dev_r5204_CNRS_PISCES_dcy/NEMOGCM/NEMO/LIM_SRC_3/limsbc.F90 (modified) (1 diff)

branches/2015/dev_r5204_CNRS_PISCES_dcy/NEMOGCM/NEMO/LIM_SRC_3/limsbc.F90

@@ -219 +219 @@
 
 
+      IF( lk_cpl ) THEN 
+         !!! LIM2 version zqsr = qsr_tot(ji,jj) + ( fstric(ji,jj) - qsr_ice(ji,jj,1) ) * ( 1.0 - pfrld(ji,jj) )
+          zfcm1 = qsr_tot(ji,jj)
+          DO jl = 1, jpl
+             zfcm1 = zfcm1 + ( ftr_ice(ji,jj,jl) - qsr_ice(ji,jj,jl) ) * a_i_b(ji,jj,jl)
+          END DO
+      ELSE
+          !!! LIM2 version zqsr = pfrld(ji,jj) * qsr(ji,jj)  + ( 1.  - pfrld(ji,jj) ) * fstric(ji,jj)
+          zfcm1   = pfrld(ji,jj) * qsr(ji,jj)
+          DO jl = 1, jpl
+             zfcm1   = zfcm1 + a_i_b(ji,jj,jl) * ftr_ice(ji,jj,jl)
+          END DO
+      ENDIF
+
       !   daily mean qsr when diurnal cycle is applied on physics - for BGC models
       IF( l_trcdm2dc ) THEN
-         DO jj = 1, jpj
-            DO ji = 1, jpi
-               zfcm1 = pfrld(ji,jj) * qsr_mean(ji,jj)
-               DO jl = 1, jpl
-                  zfcm1 = zfcm1 + ftr_ice_mean(ji,jj,jl)
-               END DO
-               qsr_mean(ji,jj) = zfcm1                                      
+         IF( lk_cpl ) THEN 
+            DO jj = 1, jpj
+               DO ji = 1, jpi
+                  zfcm1 = qsr_mean(ji,jj)
+                  DO jl = 1, jpl
+                    zfcm1 = zfcm1 + ftr_ice_mean(ji,jj,jl) - qsr_ice_mean(ji,jj,jl) * a_i_b(ji,jj,jl)
+                 END DO
+                 qsr_mean(ji,jj) = zfcm1                                      
+               ENDDO
             ENDDO
-         ENDDO
+         ELSE
+            DO jj = 1, jpj
+               DO ji = 1, jpi
+                  zfcm1 = pfrld(ji,jj) * qsr_mean(ji,jj)
+                  DO jl = 1, jpl
+                     zfcm1 = zfcm1 + ftr_ice_mean(ji,jj,jl)
+                  END DO
+                  qsr_mean(ji,jj) = zfcm1                                      
+               ENDDO
+            ENDDO
+         ENDIF
       ENDIF