Changeset 4206 for branches/2013/dev_LOCEAN_2013/NEMOGCM/NEMO/OPA_SRC/TRA/traqsr.F90

Timestamp:

2013-11-14T21:13:29+01:00 (10 years ago)

Author:

clem

Message:

 

File:

• branches/2013/dev_LOCEAN_2013/NEMOGCM/NEMO/OPA_SRC/TRA/traqsr.F90 (modified) (6 diffs)

branches/2013/dev_LOCEAN_2013/NEMOGCM/NEMO/OPA_SRC/TRA/traqsr.F90

@@ -163 +163 @@
          CALL iom_put( 'qsr3d', etot3 )   ! Shortwave Radiation 3D distribution
          ! clem: store attenuation coefficient of the first ocean level
-         IF ( ln_qsr_ice ) THEN
+         IF ( lk_lim3 .AND. ln_qsr_ice ) THEN
             DO jj = 1, jpj
                DO ji = 1, jpi
@@ -232 +232 @@
                END DO
                ! clem: store attenuation coefficient of the first ocean level
-               IF ( ln_qsr_ice ) THEN
+               IF ( lk_lim3 .AND. ln_qsr_ice ) THEN
                   DO jj = 1, jpj
                      DO ji = 1, jpi
@@ -256 +256 @@
                END DO
                ! clem: store attenuation coefficient of the first ocean level
-               IF ( ln_qsr_ice ) THEN
+               IF ( lk_lim3 .AND. ln_qsr_ice ) THEN
                   oatte(:,:) = etot3(:,:,1) / r1_rau0_rcp
                   iatte(:,:) = oatte(:,:)
@@ -280 +280 @@
                END DO
                ! clem: store attenuation coefficient of the first ocean level
-               IF ( ln_qsr_ice ) THEN
+               IF ( lk_lim3 .AND. ln_qsr_ice ) THEN
                   DO jj = 1, jpj
                      DO ji = 1, jpi
@@ -299 +299 @@
                END DO
                ! clem: store attenuation coefficient of the first ocean level
-               IF ( ln_qsr_ice ) THEN
+               IF ( lk_lim3 .AND. ln_qsr_ice ) THEN
                   oatte(:,:) = etot3(:,:,1) / r1_rau0_rcp
                   iatte(:,:) = oatte(:,:)
@@ -319 +319 @@
          !
       ENDIF
-      ! clem: store attenuation coefficient of the first ocean level
-      !IF (ln_traqsr) THEN
-      !   DO jj = 1, jpj
-      !      DO ji = 1, jpi
-      !         IF ( qsr(ji,jj) /= 0._wp ) THEN
-      !            oatte(ji,jj) = qsr_hc(ji,jj,1) / ( r1_rau0_rcp * qsr(ji,jj) )
-      !            iatte(ji,jj) = qsr_hc(ji,jj,1) / ( r1_rau0_rcp * qsr(ji,jj) )
-      !         ENDIF
-      !      END DO
-      !   END DO
-      !END IF
       !
       IF( lrst_oce ) THEN   !                  Write in the ocean restart file