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 12724 for NEMO/branches/2020/dev_r12377_KERNEL-06_techene_e3/src/OCE/TRA/trabbc.F90 – NEMO

Ignore:
Timestamp:
2020-04-08T21:37:59+02:00 (4 years ago)
Author:
techene
Message:

branch KERNEL-06 : merge with trunk@12698 #2385 - in duplcated files : changes to comply to the new trunk variables and some loop bug fixes

Location:
NEMO/branches/2020/dev_r12377_KERNEL-06_techene_e3
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2020/dev_r12377_KERNEL-06_techene_e3

    • Property svn:externals
      •  

        old new  
        33^/utils/build/mk@HEAD         mk 
        44^/utils/tools@HEAD            tools 
        5 ^/vendors/AGRIF/dev_r11615_ENHANCE-04_namelists_as_internalfiles_agrif@HEAD      ext/AGRIF 
         5^/vendors/AGRIF/dev@HEAD      ext/AGRIF 
        66^/vendors/FCM@HEAD            ext/FCM 
        77^/vendors/IOIPSL@HEAD         ext/IOIPSL 
         8 
         9# SETTE 
         10^/utils/CI/sette@HEAD         sette 
  • NEMO/branches/2020/dev_r12377_KERNEL-06_techene_e3/src/OCE/TRA/trabbc.F90

    r12590 r12724  
    6767      !!       ocean bottom can be computed once and is added to the temperature 
    6868      !!       trend juste above the bottom at each time step: 
    69       !!            ta = ta + Qsf / (rau0 rcp e3T) for k= mbkt 
     69      !!            ta = ta + Qsf / (rho0 rcp e3T) for k= mbkt 
    7070      !!       Where Qsf is the geothermal heat flux. 
    7171      !! 
     
    104104      ENDIF 
    105105      ! 
    106       CALL iom_put ( "hfgeou" , rau0_rcp * qgh_trd0(:,:) ) 
     106      CALL iom_put ( "hfgeou" , rho0_rcp * qgh_trd0(:,:) ) 
    107107      IF(sn_cfctl%l_prtctl)   CALL prt_ctl( tab3d_1=pts(:,:,:,jp_tem,Krhs), clinfo1=' bbc  - Ta: ', mask1=tmask, clinfo3='tra-ta' ) 
    108108      ! 
     
    164164         CASE ( 1 )                          !* constant flux 
    165165            IF(lwp) WRITE(numout,*) '   ==>>>   constant heat flux  =   ', rn_geoflx_cst 
    166             qgh_trd0(:,:) = r1_rau0_rcp * rn_geoflx_cst 
     166            qgh_trd0(:,:) = r1_rho0_rcp * rn_geoflx_cst 
    167167            ! 
    168168         CASE ( 2 )                          !* variable geothermal heat flux : read the geothermal fluxes in mW/m2 
     
    181181 
    182182            CALL fld_read( nit000, 1, sf_qgh )                         ! Read qgh data 
    183             qgh_trd0(:,:) = r1_rau0_rcp * sf_qgh(1)%fnow(:,:,1) * 1.e-3 ! conversion in W/m2 
     183            qgh_trd0(:,:) = r1_rho0_rcp * sf_qgh(1)%fnow(:,:,1) * 1.e-3 ! conversion in W/m2 
    184184            ! 
    185185         CASE DEFAULT 
Note: See TracChangeset for help on using the changeset viewer.