Changeset 2068 for branches/DEV_r1837_MLF/NEMO/OPA_SRC/TRA/traqsr.F90

Timestamp:

2010-09-06T17:56:51+02:00 (14 years ago)

Author:

mlelod

Message:

ticket: #663 ensuring restartability and conservation

File:

• branches/DEV_r1837_MLF/NEMO/OPA_SRC/TRA/traqsr.F90 (modified) (2 diffs)

branches/DEV_r1837_MLF/NEMO/OPA_SRC/TRA/traqsr.F90

@@ -112 +112 @@
          ztrdt(:,:,:) = ta(:,:,:) 
          ztrds(:,:,:) = 0.e0
+      ENDIF
+
+      !                                           ! ---------------------------------------- !
+      !                                           !          Swap of forcing field           !
+      !                                           ! ---------------------------------------- !
+      IF( kt /= nit000 ) qsr_trd_hc_b(:,:,:) = qsr_trd_hc_n(:,:,:)
+      !                                           ! ---------------------------------------- !
+      IF( kt == nit000 ) THEN                     !   set the forcing field at nit000 - 1    !
+         !                                        ! ---------------------------------------- !
+         IF( ln_rstart .AND.    &                               !* Restart: read in restart file
+            & iom_varid( numror, 'qsr_trd_hc_b', ldstop = .FALSE. ) > 0 ) THEN 
+            IF(lwp) WRITE(numout,*) '          nit000-1 surface forcing fields red in the restart file'
+            CALL iom_get( numror, jpdom_autoglo, 'qsr_trd_hc_b', qsr_trd_hc_b )   ! before heat content trend due to Qsr flux
+         ENDIF
       ENDIF
 
@@ -225 +239 @@
             END DO
          END DO
+      ENDIF
+
+      !                                            ! ---------------------------------------- !
+      IF( lrst_oce ) THEN                          !      Write in the ocean restart file     !
+         !                                         ! ---------------------------------------- !
+         IF(lwp) WRITE(numout,*)
+         IF(lwp) WRITE(numout,*) 'qsr : penetrative solar radiation forcing field written in ocean restart file ',   &
+            &                    'at it= ', kt,' date= ', ndastp
+         IF(lwp) WRITE(numout,*) '~~~~'
+         CALL iom_rstput( kt, nitrst, numrow, 'qsr_trd_hc_b', qsr_trd_hc_n )
+         !
       ENDIF