Changeset 2068 for branches/DEV_r1837_MLF/NEMO/OPA_SRC/TRA/trasbc.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/trasbc.F90 (modified) (5 diffs)

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

@@ -20 +20 @@
    USE trdmod          ! ocean trends 
    USE trdmod_oce      ! ocean variables trends
+   USE iom
    USE in_out_manager  ! I/O manager
+   USE restart         ! ocean restart
    USE prtctl          ! Print control
 
@@ -132 +134 @@
          ENDIF
       ENDIF
-
-      !                             ! ---------------------- !
-      IF( lk_vvl ) THEN             !  Variable Volume case  !
-         !                          ! ---------------------- !
+      !                                            ! ---------------------------------------- !
+      IF( kt /= nit000 ) THEN                      !          Swap of forcing fields          !
+         !                                         ! ---------------------------------------- !
+         sbc_trd_hc_b(:,:) = sbc_trd_hc_n(:,:)                         ! Swap the ocean forcing fields except at nit000
+         IF ( .NOT. lk_vvl ) sbc_trd_sc_b(:,:)   = sbc_trd_sc_n(:,:)
+      ENDIF
+      !                                            ! ---------------------------------------- !
+      IF( kt == nit000 ) THEN                      !   set the forcing field at nit000 - 1    !
+         !                                         ! ---------------------------------------- !
+         IF( ln_rstart .AND.    &                               !* Restart: read in restart file
+            & iom_varid( numror, 'sbc_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, 'sbc_trd_hc_b', sbc_trd_hc_b )   ! before heat content sbc trend
+            CALL iom_get( numror, jpdom_autoglo, 'qsr_trd_hc_b', qsr_trd_hc_b )   ! before heat content trend due to Qsr flux
+            IF ( .NOT. lk_vvl ) THEN
+               CALL iom_get( numror, jpdom_autoglo, 'sbc_trd_sc_b', sbc_trd_sc_b )   ! before salt content sbc trend
+            ENDIF
+         ENDIF
+      ENDIF
+      !                                            ! ---------------------- !
+      IF( lk_vvl ) THEN                            !  Variable Volume case  !
+         !                                         ! ---------------------- !
 !!gm BUG : in key_vvl emps must be modified to only include the salt flux due to sea-ice freezing/melting
 !!gm       otherwise this flux will be missing  ==> modification required in limsbc,  limsbc_2 and CICE interface.s
@@ -161 +181 @@
             END DO
          ENDIF
-         !                          ! ---------------------- !
-      ELSE                          !  Constant Volume case  !
-         !                          ! ---------------------- !
+         !                                         ! ---------------------- !
+      ELSE                                         !  Constant Volume case  !
+         !                                         ! ---------------------- !
          IF ( neuler == 0 .AND. kt == nit000 ) THEN
             DO jj = 2, jpj
@@ -197 +217 @@
       ENDIF
 
+      !                                            ! ---------------------------------------- !
+      IF( lrst_oce ) THEN                          !      Write in the ocean restart file     !
+         !                                         ! ---------------------------------------- !
+         IF(lwp) WRITE(numout,*)
+         IF(lwp) WRITE(numout,*) 'sbc : ocean surface forcing fields written in ocean restart file ',   &
+            &                    'at it= ', kt,' date= ', ndastp
+         IF(lwp) WRITE(numout,*) '~~~~'
+         CALL iom_rstput( kt, nitrst, numrow, 'sbc_trd_hc_b', sbc_trd_hc_n )
+         IF ( .NOT. lk_vvl ) CALL iom_rstput( kt, nitrst, numrow, 'sbc_trd_sc_b', sbc_trd_sc_n )
+         !
+      ENDIF
+
       IF( l_trdtra ) THEN           ! save the sbc trends for diagnostic
          ztrdt(:,:,:) = ta(:,:,:) - ztrdt(:,:,:)
@@ -205 +237 @@
       IF(ln_ctl)   CALL prt_ctl( tab3d_1=ta, clinfo1=' sbc  - Ta: ', mask1=tmask,   &
          &                       tab3d_2=sa, clinfo2=       ' Sa: ', mask2=tmask, clinfo3='tra' )
-      !
    END SUBROUTINE tra_sbc