Ignore:
Timestamp:
2020-02-28T16:55:11+01:00 (9 months ago)
Author:
davestorkey
Message:

Preparation for new timestepping scheme #2390.
Main changes:

  1. Initial euler timestep now handled in stp and not in TRA/DYN routines.
  2. Renaming of all timestep parameters. In summary, the namelist parameter is now rn_Dt and the current timestep is rDt (and rDt_ice, rDt_trc etc).
  3. Renaming of a few miscellaneous parameters, eg. atfp → rn_atfp (namelist parameter used everywhere) and rau0 → rho0.

This version gives bit-comparable results to the previous version of the trunk.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/trunk/src/OCE/TRA/trabbc.F90

    r12377 r12489  
    6666      !!       ocean bottom can be computed once and is added to the temperature 
    6767      !!       trend juste above the bottom at each time step: 
    68       !!            ta = ta + Qsf / (rau0 rcp e3T) for k= mbkt 
     68      !!            ta = ta + Qsf / (rho0 rcp e3T) for k= mbkt 
    6969      !!       Where Qsf is the geothermal heat flux. 
    7070      !! 
     
    102102      ENDIF 
    103103      ! 
    104       CALL iom_put ( "hfgeou" , rau0_rcp * qgh_trd0(:,:) ) 
     104      CALL iom_put ( "hfgeou" , rho0_rcp * qgh_trd0(:,:) ) 
    105105      IF(sn_cfctl%l_prtctl)   CALL prt_ctl( tab3d_1=pts(:,:,:,jp_tem,Krhs), clinfo1=' bbc  - Ta: ', mask1=tmask, clinfo3='tra-ta' ) 
    106106      ! 
     
    162162         CASE ( 1 )                          !* constant flux 
    163163            IF(lwp) WRITE(numout,*) '   ==>>>   constant heat flux  =   ', rn_geoflx_cst 
    164             qgh_trd0(:,:) = r1_rau0_rcp * rn_geoflx_cst 
     164            qgh_trd0(:,:) = r1_rho0_rcp * rn_geoflx_cst 
    165165            ! 
    166166         CASE ( 2 )                          !* variable geothermal heat flux : read the geothermal fluxes in mW/m2 
     
    179179 
    180180            CALL fld_read( nit000, 1, sf_qgh )                         ! Read qgh data 
    181             qgh_trd0(:,:) = r1_rau0_rcp * sf_qgh(1)%fnow(:,:,1) * 1.e-3 ! conversion in W/m2 
     181            qgh_trd0(:,:) = r1_rho0_rcp * sf_qgh(1)%fnow(:,:,1) * 1.e-3 ! conversion in W/m2 
    182182            ! 
    183183         CASE DEFAULT 
Note: See TracChangeset for help on using the changeset viewer.