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 3294 for trunk/NEMOGCM/NEMO/OPA_SRC/DIA/diahsb.F90 – NEMO

Ignore:
Timestamp:
2012-01-28T17:44:18+01:00 (12 years ago)
Author:
rblod
Message:

Merge of 3.4beta into the trunk

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMOGCM/NEMO/OPA_SRC/DIA/diahsb.F90

    r2528 r3294  
    1818   USE lib_mpp         ! distributed memory computing library 
    1919   USE trabbc          ! bottom boundary condition 
     20   USE obc_par         ! (for lk_obc) 
    2021   USE bdy_par         ! (for lk_bdy) 
    21    USE obc_par         ! (for lk_obc) 
     22   USE timing          ! preformance summary 
    2223 
    2324   IMPLICIT NONE 
     
    7273      REAL(dp)   ::   z_frc_trd_v                 !    -     - 
    7374      !!--------------------------------------------------------------------------- 
     75      IF( nn_timing == 1 )   CALL timing_start('dia_hsb') 
    7476 
    7577      ! ------------------------- ! 
     
    107109         ! heat content variation 
    108110         zdiff_hc = zdiff_hc + SUM( surf(:,:) * tmask(:,:,jk)          & 
    109             &                       * ( fse3t_n(:,:,jk) * tn(:,:,jk)   & 
     111            &                       * ( fse3t_n(:,:,jk) * tsn(:,:,jk,jp_tem)   & 
    110112            &                           - hc_loc_ini(:,:,jk) ) ) 
    111113         ! salt content variation 
    112114         zdiff_sc = zdiff_sc + SUM( surf(:,:) * tmask(:,:,jk)          & 
    113             &                       * ( fse3t_n(:,:,jk) * sn(:,:,jk)   & 
     115            &                       * ( fse3t_n(:,:,jk) * tsn(:,:,jk,jp_sal)   & 
    114116            &                           - sc_loc_ini(:,:,jk) ) ) 
    115117      ENDDO 
     
    139141      IF ( kt == nitend ) CLOSE( numhsb ) 
    140142 
     143      IF( nn_timing == 1 )   CALL timing_stop('dia_hsb') 
     144 
    1411459020  FORMAT(I5,11D15.7) 
    142146      ! 
     
    205209      WRITE(numout,*) "dia_hsb: heat salt volume budgets activated" 
    206210      WRITE(numout,*) "~~~~~~~  output written in the 'heat_salt_volume_budgets.txt' ASCII file" 
    207       IF( lk_obc .OR. lk_bdy) THEN 
     211      IF( lk_obc .or. lk_bdy ) THEN 
    208212         CALL ctl_warn( 'dia_hsb does not take open boundary fluxes into account' )          
    209213      ENDIF 
     
    248252      ! 4 - initial conservation variables ! 
    249253      ! ---------------------------------- ! 
    250       ssh_ini(:,:) = sshn(:,:)                               ! initial ssh 
     254      ssh_ini(:,:) = sshn(:,:)                                       ! initial ssh 
    251255      DO jk = 1, jpk 
    252          e3t_ini   (:,:,jk) = fse3t_n(:,:,jk)                ! initial vertical scale factors 
    253          hc_loc_ini(:,:,jk) = tn(:,:,jk) * fse3t_n(:,:,jk)   ! initial heat content 
    254          sc_loc_ini(:,:,jk) = sn(:,:,jk) * fse3t_n(:,:,jk)   ! initial salt content 
     256         e3t_ini   (:,:,jk) = fse3t_n(:,:,jk)                        ! initial vertical scale factors 
     257         hc_loc_ini(:,:,jk) = tsn(:,:,jk,jp_tem) * fse3t_n(:,:,jk)   ! initial heat content 
     258         sc_loc_ini(:,:,jk) = tsn(:,:,jk,jp_sal) * fse3t_n(:,:,jk)   ! initial salt content 
    255259      END DO 
    256260      frc_v = 0.d0                                           ! volume       trend due to forcing 
Note: See TracChangeset for help on using the changeset viewer.