New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 2108 for branches/DEV_r1837_MLF/NEMO/OPA_SRC/TRA – NEMO

Ignore:
Timestamp:
2010-09-21T21:52:56+02:00 (14 years ago)
Author:
mlelod
Message:

ticket: #663 add heat/salt/volume budgets

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/DEV_r1837_MLF/NEMO/OPA_SRC/TRA/trasbc.F90

    r2091 r2108  
    127127         qns(:,:) = qns(:,:) + qsr(:,:)      ! total heat flux in qns 
    128128         qsr(:,:) = 0.e0                     ! qsr set to zero 
    129          IF( kt == nit000 ) THEN             ! idem on before field at nit000 
    130             qsr_b(:,:) = 0.e0                      
    131          ENDIF 
     129         ! IF( kt == nit000 ) THEN             ! idem on before field at nit000 
     130         !    qsr_b(:,:) = 0.e0                      
     131         ! ENDIF 
    132132      ENDIF 
    133133      !                                          Set before sbc tracer content fields 
     
    186186      !                                          Write in the ocean restart file 
    187187      !                                          ******************************* 
    188       IF(lwp) WRITE(numout,*) 
    189       IF(lwp) WRITE(numout,*) 'sbc : ocean surface tracer content forcing fields written in ocean restart file ',   & 
    190            &                    'at it= ', kt,' date= ', ndastp 
    191       IF(lwp) WRITE(numout,*) '~~~~' 
    192       CALL iom_rstput( kt, nitrst, numrow, 'sbc_hc_b', sbc_hc_n ) 
    193       CALL iom_rstput( kt, nitrst, numrow, 'sbc_sc_b', sbc_sc_n ) 
     188      IF( lrst_oce ) THEN 
     189         IF(lwp) WRITE(numout,*) 
     190         IF(lwp) WRITE(numout,*) 'sbc : ocean surface tracer content forcing fields written in ocean restart file ',   & 
     191            &                    'at it= ', kt,' date= ', ndastp 
     192         IF(lwp) WRITE(numout,*) '~~~~' 
     193         CALL iom_rstput( kt, nitrst, numrow, 'sbc_hc_b', sbc_hc_n ) 
     194         CALL iom_rstput( kt, nitrst, numrow, 'sbc_sc_b', sbc_sc_n ) 
     195      ENDIF 
    194196       
    195197      IF( l_trdtra ) THEN           ! save the sbc trends for diagnostic 
Note: See TracChangeset for help on using the changeset viewer.