- Timestamp:
- 2018-07-11T10:24:17+02:00 (6 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2018/dev_r9838_ENHANCE04_MLF/src/OCE/SBC/sbccpl.F90
r9767 r9923 193 193 194 194 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) 196 196 197 197 INTEGER , ALLOCATABLE, SAVE, DIMENSION(:) :: nrcvinfo ! OASIS info argument … … 1100 1100 LOGICAL :: llnewtx, llnewtau ! update wind stress components and module?? 1101 1101 INTEGER :: ji, jj, jn ! dummy loop indices 1102 INTEGER :: isec ! number of seconds since nit000 (assuming r dt did not change since nit000)1102 INTEGER :: isec ! number of seconds since nit000 (assuming rn_Dt did not change since nit000) 1103 1103 REAL(wp) :: zcumulneg, zcumulpos ! temporary scalars 1104 1104 REAL(wp) :: zcoef ! temporary scalar … … 1114 1114 ! ! Receive all the atmos. fields (including ice information) 1115 1115 ! ! ======================================================= ! 1116 isec = ( kt - nit000 ) * NINT( r dt )! date of exchanges1116 isec = ( kt - nit000 ) * NINT( rn_Dt ) ! date of exchanges 1117 1117 DO jn = 1, jprcv ! received fields sent by the atmosphere 1118 1118 IF( srcv(jn)%laction ) CALL cpl_rcv( jn, isec, frcv(jn)%z3, xcplmask(:,:,1:nn_cplmodel), nrcvinfo(jn) ) … … 1259 1259 IF( kt /= nit000 ) ssh_ibb(:,:) = ssh_ib(:,:) !* Swap of ssh_ib fields 1260 1260 1261 r1_ grau = 1.e0 / (grav * rau0) !* constant for optimization1262 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) 1263 1263 apr (:,:) = frcv(jpr_mslp)%z3(:,:,1) !atmospheric pressure 1264 1264 … … 2047 2047 !!---------------------------------------------------------------------- 2048 2048 ! 2049 isec = ( kt - nit000 ) * NINT( r dt ) ! date of exchanges2049 isec = ( kt - nit000 ) * NINT( rn_Dt ) ! date of exchanges 2050 2050 2051 2051 zfr_l(:,:) = 1.- fr_i(:,:)
Note: See TracChangeset
for help on using the changeset viewer.