Changeset 5207 for branches/2015/dev_r5204_CNRS_PISCES_dcy/NEMOGCM/NEMO/LIM_SRC_3

Timestamp:

2015-04-13T13:00:27+02:00 (9 years ago)

Author:

cetlod

Message:

dev_r5204_CNRS_PISCES_dcy:minor bugs corrections

Location:

branches/2015/dev_r5204_CNRS_PISCES_dcy/NEMOGCM/NEMO/LIM_SRC_3

Files:

• limsbc.F90 (modified) (1 diff)
• limthd.F90 (modified) (1 diff)
• limthd_dif.F90 (modified) (2 diffs)

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

@@ -220 +220 @@
 
       !   daily mean qsr when diurnal cycle is applied on physics - for BGC models
-      IF( l_trcdm2dc ) qsr_mean(:,:) =  pfrld(:,:) * qsr_mean(:,:) +  ftr_ice_mean(:,:)
+      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                                      
+            ENDDO
+         ENDDO
+      ENDIF
 
       !------------------------------------------------!

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

@@ -680 +680 @@
          !         
          IF( l_trcdm2dc )  THEN
-            CALL tab_1d_2d( nbpb, qns_ice_mean(:,:,jl), npb, qns_ice_mean_1d(1:nbpb) , jpi, jpj)
+            CALL tab_1d_2d( nbpb, qsr_ice_mean(:,:,jl), npb, qsr_ice_mean_1d(1:nbpb) , jpi, jpj)
             CALL tab_1d_2d( nbpb, ftr_ice_mean(:,:,jl), npb, ftr_ice_mean_1d(1:nbpb) , jpi, jpj )
          ENDIF

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

@@ -24 +24 @@
    USE wrk_nemo       ! work arrays
    USE lib_fortran    ! Fortran utilities (allows no signed zero when 'key_nosignedzero' defined)  
-   USE sbc_oce, ONLY : lk_cpl
+   USE sbc_oce, ONLY : lk_cpl, l_trcdm2dc
 
    IMPLICIT NONE
@@ -294 +294 @@
       IF( l_trcdm2dc ) THEN
          DO ji = kideb , kiut
-             ftr_ice_mean_1d(ji) =  ftr_ice_mean_1d(ji) & 
-                 &                + a_i_1d(ji) * zftrice_mean(ji) &
-                 &                             * EXP( - rn_kappa_i * ( MAX ( 0._wp , h_i_1d(ji) ) ) ) &
-                                               * EXP( - zraext_s   * ( MAX ( 0._wp , h_s_1d(ji) ) ) )
+             ftr_ice_mean_1d(ji) =  ftr_ice_mean_1d(ji)               & 
+                 &                + a_i_1d(ji) * zftrice_mean(ji)     &
+                 &                             * EXP( - rn_kappa_i * ( MAX ( 0._wp , ht_i_1d(ji) ) ) ) &
+                 &                             * EXP( - zraext_s   * ( MAX ( 0._wp , ht_s_1d(ji) ) ) )
          END DO
       ENDIF