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 9923 for NEMO/branches/2018/dev_r9838_ENHANCE04_MLF/src/OCE/SBC/sbccpl.F90 – NEMO

Ignore:
Timestamp:
2018-07-11T10:24:17+02:00 (6 years ago)
Author:
gm
Message:

#1911 (ENHANCE-04): step I.2: dev_r9838_ENHANCE04_MLF

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2018/dev_r9838_ENHANCE04_MLF/src/OCE/SBC/sbccpl.F90

    r9767 r9923  
    193193 
    194194   REAL(wp) ::   rpref = 101000._wp   ! reference atmospheric pressure[N/m2]  
    195    REAL(wp) ::   r1_grau              ! = 1.e0 / (grav * rau0)  
     195   REAL(wp) ::   r1_rhog              ! = 1 / (rho0*grav)  
    196196 
    197197   INTEGER , ALLOCATABLE, SAVE, DIMENSION(:) ::   nrcvinfo           ! OASIS info argument 
     
    11001100      LOGICAL  ::   llnewtx, llnewtau      ! update wind stress components and module?? 
    11011101      INTEGER  ::   ji, jj, jn             ! dummy loop indices 
    1102       INTEGER  ::   isec                   ! number of seconds since nit000 (assuming rdt did not change since nit000) 
     1102      INTEGER  ::   isec                   ! number of seconds since nit000 (assuming rn_Dt did not change since nit000) 
    11031103      REAL(wp) ::   zcumulneg, zcumulpos   ! temporary scalars      
    11041104      REAL(wp) ::   zcoef                  ! temporary scalar 
     
    11141114      !                                                      ! Receive all the atmos. fields (including ice information) 
    11151115      !                                                      ! ======================================================= ! 
    1116       isec = ( kt - nit000 ) * NINT( rdt )                      ! date of exchanges 
     1116      isec = ( kt - nit000 ) * NINT( rn_Dt )                   ! date of exchanges 
    11171117      DO jn = 1, jprcv                                          ! received fields sent by the atmosphere 
    11181118         IF( srcv(jn)%laction )   CALL cpl_rcv( jn, isec, frcv(jn)%z3, xcplmask(:,:,1:nn_cplmodel), nrcvinfo(jn) ) 
     
    12591259          IF( kt /= nit000 )   ssh_ibb(:,:) = ssh_ib(:,:)    !* Swap of ssh_ib fields  
    12601260 
    1261           r1_grau = 1.e0 / (grav * rau0)               !* constant for optimization  
    1262           ssh_ib(:,:) = - ( frcv(jpr_mslp)%z3(:,:,1) - rpref ) * r1_grau    ! equivalent ssh (inverse barometer)  
     1261          r1_rhog = 1.e0 / (grav * rho0)               !* constant for optimization  
     1262          ssh_ib(:,:) = - ( frcv(jpr_mslp)%z3(:,:,1) - rpref ) * r1_rhog    ! equivalent ssh (inverse barometer)  
    12631263          apr   (:,:) =     frcv(jpr_mslp)%z3(:,:,1)                         !atmospheric pressure  
    12641264     
     
    20472047      !!---------------------------------------------------------------------- 
    20482048      ! 
    2049       isec = ( kt - nit000 ) * NINT( rdt )        ! date of exchanges 
     2049      isec = ( kt - nit000 ) * NINT( rn_Dt )        ! date of exchanges 
    20502050 
    20512051      zfr_l(:,:) = 1.- fr_i(:,:) 
Note: See TracChangeset for help on using the changeset viewer.