Changeset 15027 for NEMO/branches/2021/dev_r14318_RK3_stage1/src/OCE/TRA/traqsr.F90

Timestamp:

2021-06-19T10:14:22+02:00 (3 years ago)

Author:

techene

Message:

2605: restartability OK reproducibility ~OK (when not using RGBc) for MLF, restartability OK reproducibility OK for RK3 (with appropriate CFGs)

File:

• NEMO/branches/2021/dev_r14318_RK3_stage1/src/OCE/TRA/traqsr.F90 (modified) (8 diffs)

NEMO/branches/2021/dev_r14318_RK3_stage1/src/OCE/TRA/traqsr.F90

@@ -143 +143 @@
          ELSE                                           ! No restart or Euler forward at 1st time step
             z1_2 = 1._wp
-            DO_3D( 0,0, 0,0, 1, jpk )
+            DO_3D( nn_hls, nn_hls, nn_hls, nn_hls, 1, jpk )
                qsr_hc_b(ji,jj,jk) = 0._wp
             END_3D
@@ -149 +149 @@
       ELSE                             !==  Swap of qsr heat content  ==!
          z1_2 = 0.5_wp
-         DO_3D( 0,0, 0,0, 1, jpk )
+         DO_3D( nn_hls, nn_hls, nn_hls, nn_hls, 1, jpk )
             qsr_hc_b(ji,jj,jk) = qsr_hc(ji,jj,jk)
          END_3D
@@ -212 +212 @@
          ENDIF
       END_2D
-!!st7-2 end
+      !                                       !===>>> CAUTION: lbc_lnk is required on fraqsr_lev since sea ice computes on the full domain
+      !                                       !                otherwise restartability and reproducibility are broken 
+      CALL lbc_lnk( 'tra_qsr', fraqsr_1lev(:,:), 'T', 1._wp )
+!!st      CALL lbc_lnk( 'tra_qsr', qsr_hc(:,:,:), 'T', 1._wp )
       !
       ! TEMP: [tiling] This change not necessary and working array can use A2D(nn_hls) if using XIOS (subdomain support)
@@ -384 +387 @@
       !
       DO jk = nk0+1, nkR                  !* down to Red extinction *!   (< ~71 meters : RGB , IR removed from calculation)
-         DO_2D( 0, 0, 0, 0 )
+          DO_2D( 0, 0, 0, 0 )
             !                                      !- inverse of RGB attenuation lengths
             zlogc     = zc(ji,jj,0)
@@ -429 +432 @@
       !
       DO jk = nkR+1, nkG                  !* down to Green extinction *!   (< ~350 m : GB , IR+R removed from calculation)
-         DO_2D( 0, 0, 0, 0 )
+          DO_2D( 0, 0, 0, 0 )
             !                                      !- inverse of RGB attenuation lengths
             zlogc     = zc(ji,jj,0)
@@ -470 +473 @@
       !
       DO jk = nkG+1, nkB                  !* down to Blue extinction *!   (< ~1300 m : B , IR+RG removed from calculation)
-         DO_2D( 0, 0, 0, 0 )
+          DO_2D( 0, 0, 0, 0 )
             !                                      !- inverse of RGB attenuation lengths
             zlogc     = zc(ji,jj,0)
@@ -558 +561 @@
       !
       DO jk = 1, nk0                      !* near surface layers *!   (< ~12 meters : IR + RGB )
-         DO_2D( 0, 0, 0, 0 )
+          DO_2D( 0, 0, 0, 0 )
             ze3t = e3t(ji,jj,jk,Kmm)
             zze0 = ze0(ji,jj) * EXP( - ze3t * r1_si0 )   ;   zzeR = zeR(ji,jj) * EXP( - ze3t * r1_LR )   ! IR    ; Red  at jk+1 w-level
@@ -591 +594 @@
       !
       DO jk = nk0+1, nkR                  !* down to Red extinction *!   (< ~71 meters : RGB , IR removed from calculation)
-         DO_2D( 0, 0, 0, 0 )
+          DO_2D( 0, 0, 0, 0 )
             ze3t = e3t(ji,jj,jk,Kmm)
             zzeR = zeR(ji,jj) * EXP( - ze3t * r1_LR )                                                 ! Red          at jk+1 w-level